ˇ ´ Vysoke ´ Uc ˇen´ı Technicke ´ Cesk e ´ Fakulta elektrotechnicka Katedra kybernetiky
´ pra ´ ce Diplomova
Modelov´an´ı pr˚ utoku krve a krevn´ı disperze ve vaskul´arn´ım ˇreˇciˇsti v z´avislosti na charakteru proudˇen´ı
Autor: Bc. Michal Pinc Vedouc´ı pr´ ace: MUDr. Ing. David Mack˚ u Praha, 2015
České vysoké učení technické v Praze Fakulta elektrotechnická Katedra kybernetiky
ZADÁNÍ DIPLOMOVÉ PRÁCE Student:
Bc. Michal P i n c
Studijní program:
Biomedicínské inženýrství a informatika (magisterský)
Obor:
Biomedicínské inženýrství
Název tématu:
Modelování průtoku krve a krevní disperze ve vaskulárním řečišti v závislosti na charakteru proudění
Pokyny pro vypracování: 1. Seznamte se s následujícími tématy: modelování, kardiovaskulární systém, disperze tekutin, pulsatilní vs. nepulsatilní krevní tok, perfuze orgánů. 2. Sestavte základní smyčku modelu vaskulárního řečiště s implementovanými sensory tlaku a průtoku, sestavte SW a HW rozhraní pro přenos dat ze senzorů tlaku a průtoku do PC, vytvořte algoritmy pro vizualizaci naměřených hodnot a ukládání dat. 3. Vytvořte model řečiště s větvením alespoň do IV st. a zařaďte do základní smyčky kardiovaskulárního systému. Zapojte senzory tlaku a průtoku do kmene systému a na různé pozice na posledním stupni větvení, propojte s vytvořeným SW a HW rozhraním. Vytvořte zařízení pro generování různého charakteru pulsatilního proudění (elektromagnetický ventil a inteligentní řízení) a zapojte do systému. 4. Sledujte průtoky na posledním stupni větvení v různých částech řečiště v závislosti na tvaru průtokové a tlakové křivky v kmeni řečiště. Detailně popište chování tekutin na posledním stupni větvení při nepulsatilním a při pulsatilním toku (tvar tlakové křivky v modelu podobný tlakové křivce v lidském těle na úrovni velkých cév). Výsledky měření vyhodnoťte a diskutujte. Seznam odborné literatury: [1] Nichols, W. W.; O'Rourke, M. F.; Hartley, C., & McDonald D. A.: McDonald's blood flow in arteries: Theoretic, experimental, and clinical principles. New York, 1998. [2] Batzel, J. J.; Kappel, F.; Schneditz, D.; Tran, H. T.: Cardiovascular and Respiratory Systems: Modeling, Analysis, and Control. Philadelphia, 2006. [3] Cobelli, C.; Carsin, E.: Introduction to modeling in physiology and medicine. Academic Press, Amsterdam, 2008. Vedoucí diplomové práce: MUDr. Ing. David Macků Platnost zadání: do konce zimního semestru 2015/2016 L.S.
doc. Dr. Ing. Jan Kybic vedoucí katedry
prof. Ing. Pavel Ripka, CSc. děkan V Praze dne 8. 7. 2014
Prohl´ aˇ sen´ı autora pr´ ace Prohlaˇsuji, ˇze jsem pˇredloˇzenou pr´aci vypracoval samostatnˇe a ˇze jsem uvedl veˇsker´e pouˇzit´e informaˇcn´ı zdroje v souladu s Metodick´ ym pokynem o dodrˇzov´an´ı etick´ ych princip˚ u pˇri pˇr´ıpravˇe vysokoˇskolsk´ ych z´avˇereˇcn´ ych prac´ı. V Praze dne 5.1.2015
Bc. Michal Pinc
Podˇ ekov´ an´ı Na tomto m´ıstˇe bych r´ad podˇekoval vedouc´ımu diplomov´e pr´ace, MUDr. Ing. Davidovi Mack˚ u, za jeho rady a pˇripom´ınky k vypracov´an´ı pr´ace a za vˇseobecnˇe pˇr´atelskou atmosf´eru naˇs´ı spolupr´ace. Tak´e bych r´ad podˇekoval cel´e katedˇre kybernetiky za financov´an´ı m´e pr´ace. D´ale bych r´ad podˇekoval Ing. Miroslavovi Loˇzkovi za cenn´e rady a zap˚ ujˇcen´ı nˇekolika d˚ uleˇzit´ ych komponent a Ing. Filipovi Jeˇzkovi za u ´vod do pr´ace s dan´ ym vybaven´ım. V neposledn´ı ˇradˇe patˇr´ı podˇekov´an´ı tak´e m´e rodinˇe, a to za umoˇznˇen´ı studia a podporu nejen bˇehem vypracov´av´an´ı t´eto pr´ace, ale bˇehem cel´eho m´eho studia.
Abstract We have created mechanical model of cardiovascular system. Liquid in the system is pumped by the peristaltic (or roller) pump, heart function is simulated by expansion vessel with solenoid valve, vascular tree is simulated by branching structure of silicon tubing with decreasing diameter. Mechanical model is able to work in pulsatile or continous mode, while maintaining same cumulative flow. On the created model we have meassured flow on the lowest level of vascular tree and on the second lowest level. We have found out, that in the pulsatile mode the flow is significantly higher on the lowest level of vascular tree (0,122 ± 0,005 ml/s, in comparison with 0,0763 ± 0,004 ml/s in continous mode), while flow through second lowest level is higher in continous flow (2,770 ± 0,125 ml/s, in comparison with 2,491 ± 0,066 ml/s in pulsatile mode). We conclude, that the pulsatility of flow is significant for redistribution of liquid in vascular tree and should be taken in consideration in future design of mechanical heart supports.
Abstrakt V t´eto pr´aci byl sestrojen mechanick´ y model kardiovaskul´arn´ıho syst´emu. Tekutina je poh´anˇena peristaltickou pumpou, funkce srdce je simulov´ana kompenzaˇcn´ı n´adobou s elektromagnetick´ ym ventilem, c´evn´ı ˇreˇciˇstˇe je nahrazeno vˇetvenou strukturou silikonov´ ych hadiˇcek s postupnˇe sniˇzuj´ıc´ımi se pr˚ umˇery. Model umoˇzn ˇuje vytv´aˇret kontinu´aln´ı a pulsatiln´ı proudˇen´ı pˇri zachov´an´ı stejn´eho minutov´eho pr˚ utoku. Na vytvoˇren´em modelu jsme zmˇeˇrili pr˚ utok hadiˇckou na nejniˇzˇs´ı a druh´e nejniˇzˇs´ı u ´rovni v pulsatiln´ım a kontinu´aln´ım reˇzimu, a tyto pr˚ utoky porovnali. Zjistili jsme, ˇze pro pulsatiln´ı proudˇen´ı je pr˚ utok nejniˇzˇs´ı u ´rovn´ı ˇreˇciˇstˇe v´ yraznˇe vyˇsˇs´ı (0,122 ± 0,005 ml/s, v porovn´an´ı s 0,0763 ± 0,004 ml/s pro kontinu´aln´ı reˇzim), zat´ımco pr˚ utok druhou nejniˇzˇs´ı u ´rovn´ı je vyˇsˇs´ı v kontinu´aln´ım reˇzimu (2,770 ± 0,125 ml/s, v porovn´an´ı s pr˚ utokem 2,491 ± 0,066 ml/s v pulsatiln´ım reˇzimu). To n´as vede k z´avˇeru, ˇze pulsatilita proudˇen´ı m´a velk´ y vliv na redistribuci tekutiny v c´evn´ım syst´emu a mˇela by b´ yt br´ana v potaz pˇri n´avrhu mechanick´ ych srdeˇcn´ıch podpor.
Kl´ıˇ cov´ a slova Kardiovaskul´arn´ı syst´em, modelov´an´ı, redistribuce tekutiny
Seznam pouˇ zit´ ych zkratek BSA CCP CI CO CPB CVP DP ECMO EEP FDA HR MAP MCL mmHG ODE PP SP SV TAH TPR VAD
Body surface area, celkov´ y povrch tˇela Critical closing pressure, kritick´ y tlak uzavˇren´ı Cardiac index, srdeˇcn´ı index Cardiac outout, srdeˇcn´ı v´ ydej Cardiopulmonary bypass, mimotˇeln´ı obˇeh Central venous pressure, centr´aln´ı ˇziln´ı tlak Diastolic pressure, diastolick´ y tlak Extracorporal membrane oxygenator, mimotˇeln´ı membr´anov´ y oxigen´ator Energy equivalent pressure, energeticky ekvivalentn´ı tlak ´ rad pro kontrolu potravin a Food and drugs administration, Uˇ l´eˇciv, vl´adn´ı agentura USA Heart rate, tepov´a frekvence Mean arterial pressure, stˇredn´ı arteri´aln´ı tlak Mock circulatory loop, umˇel´a obˇehov´a smyˇcka Milimetr rtut’ov´eho sloupce, 1 mmHg = 133,3 Pa Ordinary differential equation, obyˇcejn´a diferenci´aln´ı rovnice Pulse pressure, pulsn´ı tlak Systolic pressure, systolick´ y tlak Stroke volume, tepov´ y objem Total artifficial heart, tot´aln´ı umˇel´a srdeˇcn´ı n´ahrada Total peripherial resistance, celkov´ y odpor periferi´ı Ventricular assist device, mechanick´a srdeˇcn´ı podpora
Seznam pouˇ zit´ ych pojm˚ u Bridge-to-recovery Bridge-to-transplant Bridge-to-decision Destination therapy Ischemie in vitro in vivo Perfuze Transkut´ann´ı
l´ekaˇrsk´ y pˇr´ıstroj nebo postup pro zlepˇsen´ı ˇzivota do uzdraven´ı jedince l´ekaˇrsk´ y pˇr´ıstroj nebo postup pro zlepˇsen´ı ˇzivota do nalezen´ı vhodn´eho transplant´atu l´ekaˇrsk´ y pˇr´ıstroj nebo postup pro zlepˇsen´ı ˇzivota do rozhodnut´ı o dalˇs´ım postupu l´ekaˇrsk´ y pˇr´ıstroj nebo postup pro doˇzivotn´ı zlepˇsen´ı ˇzivota nedostateˇcn´e prokrven´ı tk´anˇe nebo org´anu v umˇel´ ych podm´ınk´ach v ˇziv´ ych organismech pr˚ utok krve org´anem pˇres/skrz k˚ uˇzi
Obsah ´ 1 Uvod
9
1.1
Struktura pr´ace . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
9
1.2
C´ıle
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2 Teoretick´ y rozbor problematiky 2.1
2.2
11
Struˇcn´a historie mimotˇeln´ıho obˇehu a mechanick´ ych srdeˇcn´ıch podpor . . . 11 2.1.1
Mimotˇeln´ı obˇeh a ECMO . . . . . . . . . . . . . . . . . . . . . . . . 11
2.1.2
Mechanick´e srdeˇcn´ı podpory . . . . . . . . . . . . . . . . . . . . . . 11
2.1.3
Tot´aln´ı srdeˇcn´ı n´ahrady . . . . . . . . . . . . . . . . . . . . . . . . 12
Pulsatiln´ı a kontinu´aln´ı proudˇen´ı . . . . . . . . . . . . . . . . . . . . . . . 12 2.2.1
Definice pulsatiln´ıho proudˇen´ı . . . . . . . . . . . . . . . . . . . . . 12
2.2.2
Studie vlivu pulsatiln´ıho a kontinu´aln´ıho proudˇen´ı krve v ˇziv´em organismu . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
3 Anatomie a fyziologie kardiovaskul´ arn´ıho syst´ emu 3.1
15
Srdce . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 3.1.1
Anatomie a fyziologie srdce . . . . . . . . . . . . . . . . . . . . . . 16
3.1.2
Srdeˇcn´ı revoluce . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
3.1.3
Mechanick´a pr´ace a v´ ykon srdce . . . . . . . . . . . . . . . . . . . . 18
3.2
Mal´ y plicn´ı obˇeh . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
3.3
Velk´ y syst´emov´ y obˇeh . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
3.4
3.5
3.6
3.3.1
Art´erie a arterioly . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
3.3.2
Kapil´ary a mikrocirkulace . . . . . . . . . . . . . . . . . . . . . . . 21
3.3.3
V´eny a venuly . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
Krev . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22 3.4.1
Proudˇen´ı krve . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
3.4.2
Fyzik´aln´ı vlastnosti krve . . . . . . . . . . . . . . . . . . . . . . . . 24
Regulaˇcn´ı mechanismy kardiovaskul´arn´ıho syst´emu . . . . . . . . . . . . . 25 3.5.1
M´ıstn´ı regulaˇcn´ı mechanismy
. . . . . . . . . . . . . . . . . . . . . 25
3.5.2
Syst´emov´e regulaˇcn´ı mechanismy . . . . . . . . . . . . . . . . . . . 26
Hemodynamick´e parametry . . . . . . . . . . . . . . . . . . . . . . . . . . 27 3.6.1
Tlak krve a jeho mˇeˇren´ı . . . . . . . . . . . . . . . . . . . . . . . . 27
3.6.2
Pr˚ utok krve a jeho mˇeˇren´ı . . . . . . . . . . . . . . . . . . . . . . . 28
3.6.3
Dalˇs´ı v´ yznamn´e mˇeˇren´e parametry . . . . . . . . . . . . . . . . . . 29
6
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
4 Teorie modelov´ an´ı kardiovaskul´ arn´ıho syst´ emu 30 ´ 4.1 Uvod do teorie modelov´an´ı . . . . . . . . . . . . . . . . . . . . . . . . . . . 30 4.1.1
Popis modelu . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
4.1.2
Obecn´ y postupy modelovan´ı . . . . . . . . . . . . . . . . . . . . . . 31
4.1.3
Kompartmentov´a anal´ yza . . . . . . . . . . . . . . . . . . . . . . . 31
4.2
Specifika biologick´ ych syst´em˚ u . . . . . . . . . . . . . . . . . . . . . . . . . 32
4.3
Modely kardiovaskul´arn´ı soustavy . . . . . . . . . . . . . . . . . . . . . . . 32
5 V´ ysledn´ y model kardiovaskul´ arn´ıho syst´ emu 34 ˇ 5.1 Cerpadlo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35 5.2
5.3
5.4
Kompenzaˇcn´ı n´adoba s elektromagnetick´ ym ventilem . . . . . . . . . . . . 36 5.2.1
Kompenzaˇcn´ı n´adoba . . . . . . . . . . . . . . . . . . . . . . . . . . 36
5.2.2
Elektromagnetick´ y ventil . . . . . . . . . . . . . . . . . . . . . . . . 36
C´evn´ı ˇreˇciˇstˇe . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38 5.3.1
Pouˇzit´e hadiˇcky . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
5.3.2
Pouˇzit´e spojky . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
5.3.3
Odporov´ y prvek . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
Senzory . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40 5.4.1
5.5
Senzory tlaku . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
5.4.2 Mˇeˇren´ı pr˚ utoku . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41 ˇ ızen´ı modelu a sbˇer dat . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42 R´ 5.5.1
Hardware pro ovl´ad´an´ı a sbˇer dat . . . . . . . . . . . . . . . . . . . 42
5.5.2
Software pro ovl´ad´an´ı modelu a sbˇer dat . . . . . . . . . . . . . . . 44
6 Experiment a namˇ eˇ ren´ e v´ ysledky 6.1
6.2
6.3
Pouˇzit´e statistick´e zpracov´an´ı dat . . . . . . . . . . . . . . . . . . . . . . . 47 6.1.1
Test rovnosti stˇredn´ıch hodnot - Student˚ uv t-test . . . . . . . . . . 47
6.1.2
Test rovnosti v´ ybˇer˚ u - Mann-Whitney˚ uv test . . . . . . . . . . . . . 49
Prerekvizity experimentu . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50 6.2.1
Pˇredpoklad zachov´an´ı minutov´eho pr˚ utoku pˇri zmˇenˇe charakteru proudˇen´ı . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
6.2.2
Pˇredpoklad spr´avn´e reprezentace pr˚ utoku pomoc´ı rychlosti kuliˇcky
51
Fin´aln´ı experiment . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54 6.3.1
V´ ysledky pro mˇeˇren´ı pr˚ utoku hadiˇckou o pr˚ umˇeru 4 mm . . . . . . 54
6.3.2
V´ ysledky pro mˇeˇren´ı pr˚ utoku hadiˇckou o pr˚ umˇeru 2 mm . . . . . . 55
6.3.3
Shrnut´ı v´ ysledk˚ u . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
7 Diskuze v´ ysledk˚ u
OBSAH
47
57
-7-
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
8 Z´ avˇ er
59
A Data pro mˇ eˇ ren´ı zachov´ an´ı minutov´ eho v´ ydeje ˇ cerpadla pˇ ri zmˇ enˇ e charakteru proudˇ en´ı 62 B Data pro mˇ eˇ ren´ı porovn´ an´ı pr˚ utoku mˇ eˇ ren´ eho pomoc´ı rychlosti kuliˇ cky a odmˇ ern´ e n´ adoby 63 C Data pro porovn´ an´ı redistribuce toku pˇ ri zmˇ enˇ e charakteru proudˇ en´ı
64
D Dodateˇ cn´ y obrazov´ y materi´ al
65
E Obsah pˇ riloˇ zen´ eho CD
68
OBSAH
-8-
Kapitola 1 ´ Uvod Kardiovaskul´arn´ı syst´em je org´anov´a soustava nezbytn´a pro ˇzivot ˇclovˇeka. Srdcem kardiovaskul´arn´ı soustavy je srdce. To ˇze tento idiom vznikl a je bˇeˇzn´ y v mnoha jazyc´ıch svˇeta naznaˇcuje, jak je tento org´an pro ˇclovˇeka d˚ uleˇzit´ y. Srdeˇcn´ı choroby dodnes patˇr´ı mezi nejvˇetˇs´ı pˇr´ıˇciny u ´mrt´ı, a to i pˇresto, ˇze prvn´ı u ´spˇeˇsn´a transplantace srdce byla provedena uˇz v roce 1967 v Jihoafrick´e republice doktorem Christiaanem Barnardem (za vyuˇzit´ı poznatk˚ u a technik vyvinut´ ych doktory Normanem Shumwayem a Richardem Lowerem). Pˇr´ıˇcinou je nedostatek vhodn´ ych d´arc˚ u org´an˚ u. Proto se jiˇz dlouh´a desetilet´ı usilovnˇe pracuje na v´ yvoji pˇr´ıstroj˚ u, kter´e by podpoˇrily ˇcinnost selh´avaj´ıc´ıho srdce, popˇr´ıpadˇe by cel´e srdce nahradily. Mezi tyto pˇr´ıstroje patˇr´ı mimotˇeln´ı obˇeh a ECMO, mechanick´e srdeˇcn´ı podpory (VAD) a tot´aln´ı srdeˇcn´ı n´ahrady (TAH). Tyto pˇr´ıstroje podporuj´ı, nebo zcela nahrazuj´ı, mechanickou ˇcerpac´ı funkci srdce. Probl´emem vˇsak je, ˇze fyziologick´ y tok krve je pulsatiln´ı. Generovat pulsatiln´ı tok je mnohem sloˇzitˇejˇs´ı neˇz generovat tok kontinu´aln´ı. Pulsatiln´ı pumpy se skl´adaj´ı z vˇetˇs´ıho mnoˇzstv´ı elektronick´ ych i mechanick´ ych prvk˚ u a tud´ıˇz jsou n´achylnˇejˇs´ı k poruˇse. Logick´ ym krokem tedy je pokusit se nahradit pulsatiln´ı tok tokem kontinu´aln´ım. Kontinu´aln´ı pumpy se bˇeˇznˇe a u ´spˇeˇsnˇe pouˇz´ıvaj´ı pro stimulaci toku krve v c´evn´ım ˇreˇciˇsti pˇri operac´ıch srdce, doba pouˇzit´ı se vˇsak omezuje na des´ıtky minut ˇci hodin. Ot´azkou z˚ ust´av´a, jak´ y bude m´ıt tento nefyziologick´ y tok vliv na zbytek kardiovaskul´arn´ı soustavy a cel´ y organismus ˇclovˇeka pˇri dlouhodob´em pouˇzit´ı. Jak ovlivn´ı perfuzi org´an˚ u a jak ovlivn´ı regulaˇcn´ı mechanismy? Zkoum´an´ı vlivu typu proudˇen´ı na organismus in vivo je velice problematick´e, znamen´a to velk´ y z´asah do organizmu a rozhodnˇe ho nelze prov´adˇet na lidsk´ ych bytostech. Proto je v tˇechto v´ yzkumech ˇclovˇek nejˇcastˇeji nahrazen modelem. Tyto modely se vˇsak ˇcasto soustˇred´ı na sledov´an´ı hemodynamick´ ych parametr˚ u na u ´rovn´ı velk´ ych tepen, zat´ımco mikrocirkulace b´ yv´a nahrazena jedn´ım funkˇcn´ım blokem (napˇr [10], [11], [12]). Tato pr´ace si klade za c´ıl vytvoˇrit mechanick´ y model kardiovaskul´arn´ı soustavy, v nˇemˇz by tok tekutiny bylo moˇzn´e mˇenit z pulsatiln´ıho na kontinu´aln´ı a naopak, a n´aslednˇe na vytvoˇren´em modelu zkoumat rozd´ıly v redistribuci toku na u ´rovni velmi mal´ ych c´ev. Nen´ı c´ılem vytvoˇrit model, kter´ y by pˇresnˇe simuloval tlakov´e pomˇery na u ´rovni velk´ ych c´ev.
1.1
Struktura pr´ ace
Prvn´ı tˇri kapitoly t´eto pr´ace slouˇz´ı jako u ´vod do problematiky. Najdete zde kr´atk´e shrnut´ı historick´eho v´ yvoje mimotˇeln´ıho obˇehu, mechanick´ ych srdeˇcn´ıch podpor a tot´aln´ıch srdeˇcn´ıch n´ahrad, shrnut´ı nˇekolika v´ yznamn´ ych studi´ı, struˇcn´ y popis anatomie a fyziologie kardiovaskul´arn´ı soustavy, u ´vod do teorie modelov´an´ı a zamyˇslen´ı nad existuj´ıc´ımi modely. Dalˇs´ı kapitola se vˇenuje popisu samotn´e konstrukce naˇseho modelu, jeho ˇr´ızen´ı a sbˇeru dat. N´asleduje kapitola statistick´eho vyhodnocen´ı namˇeˇren´ ych v´ ysledk˚ u a jejich 9
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
interpretace a na z´avˇer je pˇripraveno shrnut´ı pr´ace, v´ ysledk˚ u a doporuˇcen´ y postup do budoucna.
1.2
C´ıle
Mezi c´ıle t´eto pr´ace patˇr´ı: • Anal´ yza fyziologie kardiovaskul´arn´ıho syst´emu ˇclovˇeka a problematiky modelov´an´ı • Konstrukce fyzick´eho modelu pro mˇeˇren´ı vybran´ ych veliˇcin • Mˇeˇren´ı vybran´ ych veliˇcin na vytvoˇren´em modelu • Vyhodnocen´ı namˇeˇren´ ych v´ ysledk˚ u • Diskuze v´ ysledk˚ u
C´ıle
- 10 -
Kapitola 2 Teoretick´ y rozbor problematiky V t´eto kapitole jsou pokryty teoretick´a t´emata d˚ uleˇzit´a pro naˇs´ı pr´aci. Tato t´emata zahrnuj´ı historii mechanick´ ych srdeˇcn´ıch podpor a mimotˇeln´ıho obˇehu, je zde zm´ınˇena problematika definice pulsatility toku krve a jsou zde pops´any studie vlivu pulsatility (respektive nepulsatility) krve na ˇziv´ y organismus.
2.1 2.1.1
Struˇ cn´ a historie mimotˇ eln´ıho obˇ ehu a mechanick´ ych srdeˇ cn´ıch podpor Mimotˇ eln´ı obˇ eh a ECMO
Z´akladn´ımi syst´emy nahrazuj´ıc´ı mechanickou ˇcinnost srdce jsou syst´emy pro mimotˇeln´ı obˇeh a okysliˇcov´an´ı krve (v anglick´e literatuˇre oznaˇcov´ano jako Cardiopulmonary Bypass) a pˇr´ıstroje ECMO (extracorporeal membrane oxygenator). Rozd´ıl mezi mimotˇeln´ım obˇehem a ECMO spoˇc´ıv´a pˇredevˇs´ım v jejich indikaci. Mimotˇeln´ı obˇeh se pouˇz´ıv´a sp´ıˇse pˇri operac´ıch srdce, d˚ uleˇzitou roli hraje regulace teploty a jeho pouˇzit´ı vyˇzaduje celkovou anest´ezii. ECMO lze pouˇz´ıt i s lok´aln´ı anest´ezi´ı a doba jeho aplikace se pohybuje v ˇra´du dn˚ u (3-10). Funkce tˇechto pˇr´ıstroj˚ u je podobn´a, princip je zaloˇzen na vyveden´ı krve z tˇela, u ´pravy jej´ı teploty a jej´ıho okysliˇcen´ı a navr´acen´ı do tˇela. Mimotˇeln´ı obˇeh je nezbytn´ y pro operace, kter´e vyˇzaduj´ı zastaven´ı srdce. Pro tyto operace je typick´e, ˇze se cel´ y organismus v´ yraznˇe podchlad´ı - to je zajiˇst’ov´ano pr´avˇe postupn´ ym sniˇzov´an´ım teploty krve. Podchlazen´ı organismu vede ke sn´ıˇzen´ı metabolismu (pˇri teplotˇe tˇela 28◦ C se sn´ıˇz´ı metabolismus aˇz o 50%). Pohon krve je zajiˇst’ov´an peristaltickou ˇci centrifug´aln´ı pumpou. Peristaltick´a (nebo v´aleˇckov´a rotaˇcn´ı) pumpa m´a tu v´ yhodu, ˇze krev je uzavˇrena v pruˇzn´e trubici, tud´ıˇz nepˇrich´az´ı do kontaktu se samotnou pumpou a sterilizace je mnohem jednoduˇsˇs´ı a levnˇejˇs´ı. Krev je potom postupnou rotac´ı pˇr´ıtlaˇcn´ ych v´aleˇck˚ u poh´anˇena ve smˇeru proudˇen´ı. Nev´ yhodou je, ˇze v oblasti pˇrim´aˇcknut´ı trubice v´aleˇckem je krev stlaˇcena, coˇz m˚ uˇze zp˚ usobit poˇskozen´ı ˇcerven´ ych krvinek (hemol´ yzu). Okysliˇcen´ı krve se dnes prov´ad´ı nejˇcastˇeji pomoc´ı membr´anov´eho oxygen´atoru. Krev prot´ek´a velice tenk´ ymi kapil´arami se stˇenou prostupnou pro krevn´ı plyny, nikoli vˇsak pro erytrocyty. Na z´akladˇe rozd´ıln´ ych parci´aln´ıch tlak˚ u potom doch´az´ı k v´ ymˇenˇe plyn˚ u mezi krv´ı a plynn´ ym prostˇred´ım.
2.1.2
Mechanick´ e srdeˇ cn´ı podpory
Mechanick´e srdeˇcn´ı podpory (VAD) svoj´ı funkc´ı pom´ahaj´ı srdci pumpovat krev, nahrazuj´ı tak srdeˇcn´ı komoru. Z´akladem tˇechto pˇr´ıstroj˚ u je pulsn´ı nebo kontinu´aln´ı pumpa, 11
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
uloˇzen´a mimo tˇelo nebo implantovan´a. VAD m˚ uˇze nahrazovat levou (LVAD) nebo pravou (RVAD) komoru, popˇr´ıpadˇe obˇe souˇcasnˇe (BiVAD). Prvn´ı pˇr´ıpad uˇzit´ı VAD se pˇrisuzuje Americk´emu l´ekaˇri Michaelu DeBakeymu. Jednalo se o mimotˇeln´ı pumpu, kter´e po dobu 10 dn˚ u podporovala krevn´ı obˇeh ˇzeny po komplikovan´e operaci srdce[13]. N´asledoval nˇekolik des´ıtek let trvaj´ıc´ı v´ yvoj, kter´ y vedl k v´ yvoji nˇekolika nov´ ych pˇr´ıstroj˚ u, jako je napˇr´ıklad pˇr´ıstroj od firmy Novacor (pouˇzit v roce 1984 jako bridge-to-transplant) ˇci pˇr´ıstroj HeartMate od firmy Thoratec (implantovateln´ y, pouˇzit jako bridge-to-transplant v roce 1992) [13]. Od t´e doby se ale v´ yvoj a v´ yzkum soustˇred´ı sp´ıˇse na kontinu´aln´ı pumpy, a to z toho d˚ uvodu, ˇze syst´emy s kontinu´aln´ı pumpou jsou jednoduˇsˇs´ı, je tedy niˇzˇs´ı riziko poruchy, a tak´e jsou rozmˇerovˇe menˇs´ı. VAD se pouˇz´ıvaj´ı jak v kr´atkodob´em aˇz stˇrednˇedob´em horizontu, napˇr´ıklad jako bridge-to-recovery nebo bridge-to-transplant, tak i v dlouhodob´em horizontu jako destination therapy. VAD nenahrazuje cel´e srdce, ale podporuje jeho funkci. P˚ uvodn´ı srdce se st´ale do jist´e m´ıry u ´ˇcastn´ı ˇcerp´an´ı krve. T´ım m˚ uˇze v jinak kontinu´aln´ım proudˇen´ı dod´avan´em VAD stimulovat ˇc´asteˇcnou pulsaci [16]. Uˇzit´ı tˇechto pˇr´ıstroj˚ u je relativnˇe rozˇs´ıˇren´e, napˇr´ıklad pˇr´ıstroj Heart-Mate II od firmy Thoratec, spadaj´ıc´ı do kategorie LVAD, byl od roku 2005 implantov´an v´ıce neˇz 6000 pacient˚ u [16].
2.1.3
Tot´ aln´ı srdeˇ cn´ı n´ ahrady
Tot´aln´ı srdeˇcn´ı n´ahrady (TAH) maj´ı tak´e relativnˇe dlouhou historii, prvn´ı implantace probˇehla jiˇz v roce 1969 na dobu tˇr´ı dn˚ u, po kter´e bylo umˇel´e srdce nahrazeno transplantovan´ ym. V´ yzkum a v´ yvoj tot´aln´ım srdeˇcn´ıch n´ahrad posl´eze ustoupil do pozad´ı aˇz do roku 2001. TAH jsou u ´zce spjaty s VAD, princip jejich funkce je obdobn´ y. TAH vˇsak ˇcasto vyˇzaduj´ı resekci cel´eho srdce, m´a smysl je pouˇz´ıvat tedy tam, kdy nelze oˇcek´avat uzdraven´ı srdce nebo kdy je potˇreba srdce z d˚ uvodu onemocnˇen´ı odstranit. V souˇcasn´e dobˇe stoj´ı za zm´ınku pˇredevˇs´ım TAH CardioWest od firmy SynCardia (povolen´ y k uˇzit´ı v USA, Evropˇe a Kanadˇe). Tento pˇr´ıstroj je celkem rozmˇern´ y a proto je vhodn´ y jen pro vˇetˇs´ı pacienty. Celkem bylo tˇechto zaˇr´ızen´ı implantov´ano v´ıce neˇz 1350 [15]. V souˇcasn´e dobˇe se vˇsak pracuje tak´e na menˇs´ı verzi, urˇcen´e pro mladˇs´ı a menˇs´ı pacienty.
2.2
Pulsatiln´ı a kontinu´ aln´ı proudˇ en´ı
Pˇri n´avrhu mechanick´ ych srdeˇcn´ıch podpor je velk´a tendence nahradit fyziologick´ y pulsatiln´ı tok krve tokem kontinu´aln´ım, a to z toho d˚ uvodu, ˇze kontinu´aln´ı pumpy jsou jednoduˇsˇs´ı a maj´ı tedy vˇetˇs´ı ˇzivotnost. Velkou ot´azkou ale z˚ ust´av´a, jak kontinu´aln´ı tok ovlivn´ı cel´ y organismus, pˇredevˇs´ım na u ´rovni mikrocirkulace, kde m˚ uˇze doj´ıt k nespr´avn´e a nedostateˇcn´e redistribuci krevn´ıho toku, coˇz m˚ uˇze v´est k nedostateˇcn´e perfuzi org´an˚ u a dalˇs´ım neblah´ ym jev˚ um.
2.2.1
Definice pulsatiln´ıho proudˇ en´ı
Velk´ y probl´em nast´av´a uˇz u samotn´e definice pulsatiln´ıho proudˇen´ı krve v c´ev´ach, pro kterou jeˇstˇe nen´ı pˇrijata jednotn´a metodika. V nˇekter´ ych pˇr´ıpadech se pulsatilita toku
Tot´aln´ı srdeˇcn´ı n´ahrady
- 12 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
definuje jen podle pulsn´ıho tlaku. pulsn´ı tlak je rozd´ıl mezi maxim´aln´ı (systolickou) a minim´aln´ı (diastolickou) hodnotou tlakov´e kˇrivky. Proudˇen´ı s pulsn´ım tlakem vˇetˇs´ım neˇz 15 mmHg se povaˇzuje za pulsatiln´ı, s tlakem menˇs´ım neˇz 15 mmHg jako nepulsatiln´ı. Tento pˇr´ıstup zcela zanedb´av´a tvar pr˚ utokov´e kˇrivky krve a dalˇs´ı v´ yznamn´e hemodynamick´e parametry, jako je napˇr´ıklad tepov´ y objem (SV). Jedn´ım z postup˚ u jak l´epe definovat pulsatiln´ı proudˇen´ı je stanoven´ı energeticky ekvivalentn´ıho tlaku - EEP. Tuto metodu navrhl Richard B. Shepard a tuto hodnotu lze spoˇc´ıtat n´asleduj´ıc´ım zp˚ usobem: R f · p · dt (2.1) EEP = R f · dt kde EEP je energeticky ekvivalentn´ı tlak [mmHg], f je pr˚ utok dan´y pumpu [ml/s], p je arteri´aln´ı tlak [mmHg] EEP tedy d´av´a do vztahu tlakovou i tokovou kˇrivku. Ve v´ ysledku reprezentuje celkovou energii pulsn´ı vlny. Pˇri porovn´av´an´ı pulsatiln´ıho a kontinu´aln´ıho proudˇen´ı se ale nejˇcastˇeji operuje s veliˇcinou stˇredn´ı arteri´aln´ı tlak (MAP). Pro fyziologick´e pulsatiln´ı proudˇen´ı se d´a tento tlak vypoˇc´ıtat jako: SP + 2DP (2.2) M AP = 3 kde MAP je stˇredn´ı arteri´aln´ı tlak [mmHg], DP je diastolick´y arteri´aln´ı tlak [mmHg], SP je systolick´y arteri´aln´ı tlak [mmHg]
Toto je pˇribliˇzn´a a zjednoduˇsen´a metoda platn´a jen pro fyziologick´e tlakov´e kˇrivky. Pˇresnˇejˇs´ı metoda je tato: Z 1 T p · dt (2.3) M AP = T 0 kde MAP je stˇredn´ı arteri´aln´ı tlak [mmHg], p je okamˇzit´y aort´aln´ı tlak [mmHg], T je perioda pulsac´ı [s]
Pro kontinu´aln´ı proudˇen´ı plat´ı, ˇze EEP je stejn´ y jako MAP. Nicm´enˇe [9] uv´ad´ı, ˇze pro pulsatiln´ı tok jsou to jiˇz dost odliˇsn´a ˇc´ısla, konkr´etnˇe EEP, a tedy i energie nesen´a pulsatiln´ım tokem, je u koˇrene aorty 1.0 aˇz 2.3 kr´at vyˇsˇs´ı neˇz je MAP. Tato energie se m˚ uˇze pˇren´est aˇz do kapil´ar a zv´ yˇsit jejich pr˚ uchodnost. Je totiˇz potˇreba m´ıt na pamˇeti, ˇze c´evy nejsou ide´aln´ı rigidn´ı trubky, ale elastick´e tk´anˇe, na kter´e p˚ usob´ı tlak okoln´ıch tk´an´ı, coˇz ve v´ ysledku to znamen´a, ˇze arterioly a kapil´ary pˇri nulov´em tlaku krve nez˚ ust´avaj´ı otevˇreny. Tento parametr se oznaˇcuje jako kritick´ y tlak uzavˇren´ı (CCP). [2] uv´ad´ı ˇze hodnota CCP se pohybuje okolo 20 mmHg. Pokud tlak krve poklesne pod tuto hodnotu, c´evou nepoteˇce ˇza´dn´a krev.
2.2.2
Studie vlivu pulsatiln´ıho a kontinu´ aln´ıho proudˇ en´ı krve v ˇ ziv´ em organismu
Pˇrestoˇze sledov´an´ı vlivu kontinu´aln´ıho a pulsatiln´ıho proudˇen´ı in vivo nen´ı v˚ ubec jednoduch´e, je pˇresto provediteln´e, jak dokazuje Atsushi Baba a Jarom´ır Vaˇsk˚ u se sv´ ym t´ ymem Studie vlivu pulsatiln´ıho a kontinu´aln´ıho proudˇen´ı krve v ˇziv´em organismu
- 13 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
ve studii [5] a [6]. Tento t´ ym se zab´ yval sledov´an´ım fyziologick´ ych zmˇen vyvolan´ ych v tˇele kozi s implantovanou tot´aln´ı srdeˇcn´ı n´ahradou poh´anˇenou undulaˇcn´ı pumpou (umoˇzn ˇuj´ıc´ı vytv´aˇret libovoln´ y profil proudˇen´ı) v z´avislosti na charakteru proudˇen´ı krve. Prvn´ı ze studi´ı se vˇenuje akutn´ım zmˇen´am v tˇele pˇri zmˇenˇe profilu toku, konkr´etnˇe se zamˇeˇruje na rychlost pohybu erytrocyt˚ u v c´ev´ach a perfuzi oˇcn´ı spojivky. Tˇrem koz´am byla implantov´ana tot´aln´ı srdeˇcn´ı n´ahrada, zpoˇca´tku pracuj´ıc´ı v pulsatiln´ım reˇzimu s arteri´aln´ım tlakem dosahuj´ıc´ı hodnot 140/80 mmHg (stˇredn´ı arteri´aln´ı tlak stanoven jako 120 mmHg). Po zotaven´ı z operace (mezi 21. a 50. dnem po operaci) byla pomoc´ı digit´aln´ıho mikroskopu a kamery s vysok´ ym rozliˇsen´ım zmˇeˇrena rychlost proudˇen´ı erytrocyt˚ u v kapil´ar´ach a tak´e byl stanoven poˇcet perfundovan´ ych kapil´ar proch´azej´ıc´ı pomyslnou ˇca´rou o zn´am´e d´elce (oboj´ı mˇeˇreno na oˇcn´ı spojivce zv´ıˇrete, kter´emu byla pod´ana sedativa, ale nebylo v anest´ezii). N´aslednˇe byla pumpa pˇrepnuta do kontinu´aln´ıho reˇzimu s arteri´aln´ım tlakem 120 mmHg, ve kter´em pracovala 20 minut. Pot´e byla pumpa opˇet pˇrepnuta do pulsatiln´ıho reˇzimu. Bˇehem cel´e procedury byla pravidelnˇe mˇeˇrena rychlost erytrocyt˚ u a stanoven poˇcet perfundovan´ ych kapil´ar. P˚ uvodn´ı stˇredn´ı rychlost erytrocyt˚ u pˇri pulsatiln´ım proudˇen´ı byla stanovena jako 526 ± 83 µm/s, pˇri pˇrechodu na kontinu´aln´ı reˇzim klesla na 132 ± 41 µm/s a pˇri pˇrechodu zpˇet na pulsn´ı proudˇen´ı postupnˇe vystoupala aˇz na 433 ± 71 µm/s. Poˇcet perfundovan´ ych kapil´ar z 34,6 ± 6,3 kapil´ar/mm pˇri pulsatiln´ım proudˇen´ı klesl na 19,7 ± 4,1 kapil´ar/mm pˇri kontinu´aln´ım proudˇen´ı. Druh´a ze jmenovan´ ych studi´ı se potom zamˇeˇrila na chronick´e n´asledky kontinu´aln´ıho proudˇen´ı krve v tˇele zv´ıˇrete. Po ukonˇcen´ı pˇredchoz´ıho experimentu se koza nechala ˇz´ıt 33 dn´ı s pumpou pracuj´ıc´ı v kontinu´aln´ım reˇzimu (po 33 dnech doˇslo k mechanick´emu selh´an´ı pumpy). N´aslednˇe byly odebr´any histologick´e vzorky z ledvin, jater a plic. V odebran´ ych vzorc´ıch byly ve zv´ yˇsen´e m´ıˇre nalezeny stopy po chronick´e venost´aze (hromadˇen´ı krve chud´e na kysl´ık ve ven´ozn´ı ˇc´asti ˇreˇciˇstˇe), isch´emii a ischemick´e nekr´ozy, otoku a vnitˇrn´ıho krv´acen´ı. Atsushi Baba, Jarom´ır Vaˇsk˚ u a jejich t´ ym tvrd´ı, ˇze pulsatilita proudˇen´ı by mˇela b´ yt br´ana jako vit´aln´ı podm´ınka pro dlouhodob´e pˇreˇzit´ı organismu s tot´aln´ı srdeˇcn´ı n´ahradou. Zjiˇstˇen´e patologick´e zmˇeny jsou podle nich zp˚ usobeny hlavnˇe poˇskozen´ım funkce regulaˇcn´ıch mechanism˚ u kardiovaskul´arn´ı soustavy. Ty zajiˇst’uj´ı adekv´atn´ı perfuzi org´an˚ u podle aktu´aln´ıch potˇreb, ale z´aroveˇ n zachov´avaj´ı z´akladn´ı funkci kardiovaskul´arn´ı soustavy, tedy obˇeh krve v cel´em syst´emu (regulaˇcn´ı mechanismy m´ıstn´ı a syst´emov´e). Pokud bude tˇelem proudit krev kontinu´alnˇe, regulaˇcn´ı mechanismy nebudou dod´avat do centra ˇr´ızen´ı spr´avn´e informace, coˇz zp˚ usob´ı celkovou vazokonstrikci a t´ım i sn´ıˇzen´ı perfuze org´an˚ u. Toto je fyziologick´ y pohled na vˇec. Je ovˇsem moˇzn´e, ˇze samotn´a distribuce krevn´ıho toku z´aleˇz´ı na typu proudˇen´ı a vlastnostech c´evn´ıho ˇreˇciˇstˇe, i bez pˇriˇcinˇen´ı regulaˇcn´ıch mechanism˚ u. Ovˇeˇren´ı t´eto domnˇenky je c´ılem t´eto pr´ace.
Studie vlivu pulsatiln´ıho a kontinu´aln´ıho proudˇen´ı krve v ˇziv´em organismu
- 14 -
Kapitola 3 Anatomie a fyziologie kardiovaskul´ arn´ıho syst´ emu Kardiovaskul´arn´ı syst´em je vit´aln´ı org´anov´a soustava vyˇsˇs´ıch ˇzivoˇcich˚ u. Slouˇz´ı k distribuci krve, d´ıky n´ıˇz se po tˇele rozv´ad´ı d˚ uleˇzit´e plyny, ˇziviny, hormony, metabolick´e produkty a mnoh´e dalˇs´ı v´ yznamn´e l´atky. D´ale zajiˇst’uje termoregulaci organismu a je ned´ılnou souˇca´st´ı imunitn´ıho syst´emu. Lidsk´ y kardiovaskul´arn´ı syst´em se vyznaˇcuje t´ım, ˇze je uzavˇren´ y (krev nikdy neopouˇst´ı c´evy) a t´ım, ˇze se skl´ad´a ze dvou oddˇelen´ ych, s´eriovˇe ˇrazen´ ych syst´em˚ u. Mal´ y (plicn´ı) ’ obˇeh zajiˇst uje transport neokysliˇcen´e krve z prav´e komory do plicn´ıho ˇreˇciˇstˇe, kde je krev okysliˇcena, a d´ale pak do lev´e s´ınˇe. Velk´ y (syst´emov´ y obˇeh) potom pˇrepravuje krev z lev´e komory skrz syst´emov´e c´evy aˇz k jednotliv´ ym buˇ nk´am a pak zp´atky do prav´e s´ınˇe. Oba tyto syst´emy jsou spolu spojeny v srdci.
3.1
Srdce
Srdce je pumpa, kter´a m´a hlavn´ı pod´ıl na obˇehu krve v tˇele. Je to dut´ y, svalov´ y org´an pˇribliˇznˇe o velikosti pˇesti. Lidsk´e srdce se skl´ad´a ze dvou s´ın´ı a ze dvou komor. Prav´e srdce ˇcerp´a krev ze syst´emov´eho obˇehu do plicn´ıho, lev´e pak z plicn´ıho obˇehu dod´av´a krev do obˇehu syst´emov´eho.
15
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
3.1.1
Anatomie a fyziologie srdce
Anatomii srdce lze nejl´epe ilustrovat obr´azkem, d˚ uleˇzit´e pojmy jsou vysvˇetleny d´ale:
Obr´azek 3.1: Sch´ema srdce: 1 – endokard; 2 – myokard; 3 – epikard; 4 – perikard; 5 – srdeˇcn´ı dutiny; 6 – prav´a s´ıˇ n; 7 – lev´a s´ıˇ n; 8 – prav´a komora; 9 – lev´a komora; 10 – srdeˇcn´ı chlopnˇe; 11 – mitr´aln´ı chlopeˇ n; 12 – aort´aln´ı chlopeˇ n; 13 – pulmon´aln´ı chlopeˇ n; 14 – trikuspid´aln´ı chlopeˇ n; 15 – oblouk aorty; 16 – prav´a plicn´ı tepna; 17 – prav´e plicn´ı ˇz´ıly; 18 – horn´ı dut´a ˇz´ıla; 19 – doln´ı dut´a ˇz´ıla; 20 – s´ıˇ no-komorov´e septum. Obr´azek pˇrevzat z [14] Srdce se skl´ad´a ze ˇctyˇr vrstev [3], jedn´a se endokard, myokard, epikard a perikard. Endokard je nejvnitˇrnˇejˇs´ı vrstva srdce, je to hladk´a vrstva bunˇek, kter´a vyst´ yl´a vnitˇrn´ı prostory srdce (s´ınˇe a komory). Jsou j´ım tak´e tvoˇreny srdeˇcn´ı chlopnˇe. Myokard je prostˇredn´ı vrstva srdce a je tvoˇrena speci´aln´ım typem svalov´ ych bunˇek (srdeˇcn´ı svalovinou) a buˇ nkami pˇrevodn´ıho syst´emu srdeˇcn´ıho. Srdeˇcn´ı svalovina je jedineˇcn´a t´ım, ˇze dok´aˇze pracovat nepˇretrˇzitˇe cel´ y ˇzivot (jako hladk´a svalovina) a z´aroveˇ n je schopna po celou dobu poskytovat rychle a opakovanˇe relativnˇe velk´ y v´ ykon (jako ˇ ast myokardu, kter´e je tvoˇrena tˇemito buˇ pˇr´ıˇcnˇe pruhovan´e svalstvo). C´ nkami, se oznaˇcuje za pracovn´ı myokard a jej´ı hlavn´ı funkc´ı je kontrakce dutin. Buˇ nky pˇrevodn´ıho syst´emu srdeˇcn´ıho jsou potom zodpovˇedny za pravidelnou a automatickou srdeˇcn´ı ˇcinnost a rozvod vzruch˚ u v urˇcen´em poˇrad´ı po cel´em srdci. Pˇrevodn´ı syst´em srdeˇcn´ı charakterizuj´ı n´asleduj´ıc´ı vlastnosti: ˇ ızen´ı vzniku vzruch˚ • Autonomie. R´ u je na samotn´em pˇrevodn´ım syst´emu srdeˇcn´ım. Centr´aln´ı nervov´a soustava nen´ı schopna ovlivnit tvoˇren´ı puls˚ u, je schopna ovlivnit pouze jejich frekvenci. • Automacie. Vzruchy se v srdci vytv´aˇr´ı samoˇcinnˇe, bez vnˇejˇs´ıho ˇr´ızen´ı. • Rytmicita. Vzruchy se vytv´aˇr´ı pravidelnˇe s urˇcitou frekvenc´ı. Anatomie a fyziologie srdce
- 16 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
Pˇrevodn´ı syst´em srdeˇcn´ı se skl´ad´a z n´asleduj´ıc´ıch ˇc´ast´ı [3]: • Sinoatri´aln´ı uzel. Nach´az´ı se nad pravou s´ın´ı bl´ızko u ´st´ı horn´ı dut´e ˇz´ıly (Vena cava superior). Jedn´a se o prim´arn´ı pacemaker, tedy udavatele rytmu. • Atrioventrikul´arn´ı uzel. Nach´az´ı se ve stˇenˇe prav´e s´ınˇe pod sinoatri´aln´ım uzlem. Vede vzruch o pozn´an´ı pomaleji neˇz sinoatri´aln´ı uzel, d´ıky ˇcemuˇz maj´ı s´ınˇe dostatek ˇcasu na plnou kontrakci. V pˇr´ıpadˇe poˇskozen´ı sinoatri´aln´ıho uzlu pˇreb´ır´a jeho funkci atrioventrikul´arn´ı uzel, kter´ y se z tohoto d˚ uvodu oznaˇcuje za sekund´arn´ı pacemaker. • His˚ uv svazek. Zajiˇst’uje pˇrevod vzruchu ze s´ın´ı do komor, proch´az´ı skrz stˇenu oddˇeluj´ıc´ı s´ınˇe a komory. • Tawarova r´amenka. Uvnitˇr stˇeny komor je vzruch veden pomoc´ı prav´eho a lev´eho Tawarova ram´enka. • Purkyˇ nova vl´akna. Purkyˇ nova vl´akna se vˇetv´ı z Tawarov´ ych ram´enek a vedou vzruch k srdeˇcn´ı svalovinˇe prav´e a lev´e komory. Cel´e srdce je obaleno epikardem a uloˇzeno v perikardu, kter´ y tvoˇr´ı vakovit´ y obal srdce, jenˇz ho oddˇeluje od ostatn´ıch org´an˚ u. Epikard je pˇr´ımo pˇrirostl´ y k srdci. Mezi epikardem a perikardem je mal´e mnoˇzstv´ı tekutiny, kter´a zvyˇsuje kluzkost a t´ım usnadˇ nuje pohyb srdce.
3.1.2
Srdeˇ cn´ı revoluce
Jako srdeˇcn´ı revoluce se oznaˇcuje jeden cyklus srdce, ve kter´em dojde k systole (stahu myokardu) a diastole (uvolnˇen´ı myokardu). Srdeˇcn´ı revoluce sest´av´a ze ˇctyˇr f´az´ı [3]. Prvn´ı f´az´ı je pln´ıc´ı f´aze. Bˇehem t´eto f´aze doch´az´ı k tomu, ˇze se uvoln´ı s´ınˇe i komory. Nitrokomorov´ y tlak klesne k nule a krev vlivem setrvaˇcnosti proud´ı do komor. Druhou f´az´ı je izovolumetrick´a kontrakce. Dojde k systole komor, jejichˇz objem se nemˇen´ı, ale zvyˇsuje se vnitˇrn´ı tlak. Chlopnˇe mezi s´ınˇemi a komorami se uzavˇrou, coˇz se projev´ı vznikem prvn´ı srdeˇcn´ı ozvy (odposlechnuteln´e fonendoskopem). N´asleduje ejekˇcn´ı f´aze. Kdyˇz nitrokomorov´ y tlak stoupne nad hodnotu potˇrebnou pro otevˇren´ı aort´aln´ı (cca 80 mmHg) a pulmon´arn´ı (10 mmHg) chlopnˇe, tyto chlopnˇe se otevˇrou a dojde k vypuzen´ı krve z komor. Jakmile je krev vypuzena z komor, nast´av´a f´aze izovolumetrick´e relaxace. Aort´aln´ı a pulmon´arn´ı chlopnˇe se vlivem poklesu nitrokomorov´eho tlaku uzavˇrou, coˇz vyvol´a druhou srdeˇcn´ı ozvu. Jakmile tlak v komor´ach klesne pod tlak v s´ın´ıch, otevˇrou se chlopnˇe mezi s´ınˇemi a komorami a cel´ y cyklus se opakuje.
Srdeˇcn´ı revoluce
- 17 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
´ Obr´azek 3.2: PV diagram lev´e komory. Usek AB izovolumterick´a kontrakce, BC ejekˇcn´ı f´aze, CD izovolumetrick´a relaxace, DA pln´ıc´ı f´aze
3.1.3
Mechanick´ a pr´ ace a v´ ykon srdce
Srdce se vyznaˇcuje svoj´ı st´alou ˇcinnost´ı. Jedn´a se o ˇcerpadlo, kter´e z fyzik´aln´ıho pohledu kon´a pˇri kaˇzd´em stahu mechanickou pr´aci. Tato pr´ace se d´a rozdˇelit na kinetickou a tlakovˇe-objemovou sloˇzku [4]. Kinetick´a pr´ace je vynaloˇzena na udˇelen´ı urˇcit´eho zrychlen´ı tepov´emu objemu (stroke volume). Tlakovˇe-objemov´a pr´ace se vynakl´ad´a na posunut´ı dan´eho objemu proti dan´emu tlaku. Lze tedy ps´at, ˇze: W = Wk + Wp
(3.1)
kde W je celkov´a vykonan´a mechanick´a pr´ace [J], Wk je kinetick´a pr´ace [J] a Wp je tlakovˇe-objemov´a pr´ace [J] Pro kinetickou sloˇzku potom plat´ı: 1 Wk = ρv 2 V 2
(3.2)
kde Wk je kinetick´a pr´ace [J], ρ je hustota krve [kg/m3 ], v je rychlost proud´ıc´ı krve [m/s] a V je tepov´y objem [m3 ] Pro tlakovˇe-objemovou sloˇzku potom plat´ı: Wp = pV
(3.3)
kde Wp je tlakovˇe-objemov´a pr´ace [J], p je tlak, proti kter´emu se pumpuje [Pa] a V je tepov´y objem [m3 ] Pr´aci kon´a jak lev´a, tak prav´a komora. Hodnoty jednotliv´ ych promˇenn´ ych jsou uvedeny v n´asleduj´ıc´ı tabulce:
Mechanick´a pr´ace a v´ ykon srdce
- 18 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı Tepov´ y objem Hustota krve Rychlost proudˇen´ı krve Tlak pro pravou komoru Tlak pro levou komoru
70ml (7 ·10−5 m3 ) 1060 kg/m3 10,5 m/s 15 mmHg (1 999,5 Pa) 120 mmHg (15 996,5 Pa)
Tabulka 3.1: Hodnoty promˇenn´ ych uˇzit´e pro v´ ypoˇcet mechanick´e pr´ace srdce Po dosazen´ı hodnot promˇenn´ ych do uveden´ ych rovnic dost´av´ame n´asleduj´ıc´ı v´ ysledky: Wk Wk Wp Wp W
prav´e komory lev´e komory prav´e komory lev´e komory
0,009 J 0,009 J 0,140J 1,120 J 1,278 J
Tabulka 3.2: V´ ysledn´e d´ılˇc´ı mechanick´e pr´ace a jejich suma Dostali jsme tedy v´ yslednou mechanickou pr´aci srdce pˇri jedn´e systole. Pro v´ ypoˇcet v´ ykonu pouˇzijeme standardn´ı vzorec P =
W t
(3.4)
kde P je v´ykon [W], W je mechanick´a pr´ace [J] a t je ˇcas [s] Vezmeme-li v u ´vahu, ˇze klidov´a tepov´a frekvence ˇcin´ı 70 tep˚ u za minutu, dost´av´ame ˇze jeden tep trv´a 0,857 sekundy. Z toho dost´av´ame, ˇze v´ ykon srdce je pˇribliˇznˇe 1,5 wattu. To je vˇsak jen v´ ykon, kter´ y se spotˇrebuje pˇri stahu srdce. Celkov´ y v´ ykon srdce se uv´ad´ı jako 13 watt˚ u, z ˇcehoˇz 9/10 se spotˇrebuje na udrˇzen´ı svalov´eho tonu srdce.
3.2
Mal´ y plicn´ı obˇ eh
Mal´ y plicn´ı obˇeh slouˇz´ı k pˇrepravˇe krve do plic, kde dojde k jej´ımu okysliˇcen´ı a n´aslednˇe n´avratu do srdce. Tento syst´em se tak´e oznaˇcuje jako n´ızkotlak´ y, nebot’ hodnoty systolick´eho tlaku se v plicn´ı tepnˇe pohybuj´ı v rozmez´ı 15 aˇz 30 mmHg. Krev je puzena z prav´e komory do plicn´ı arterie, kter´a se postupnˇe ˇstˇep´ı a zmenˇsuje aˇz na u ´roveˇ n vl´aseˇcnic, kter´e obt´aˇc´ı jednotliv´e plicn´ı skl´ıpky. V tomto m´ıstˇe doch´az´ı k v´ ymˇenˇe d´ ychac´ıch plyn˚ u mezi krv´ı a okoln´ım prostˇred´ı. Zp´atky se krev vrac´ı ˇctyˇrmi plicn´ımi ˇzilami do lev´e s´ınˇe. Krev je tedy v plicn´ıch tepn´ach neokysliˇcena, v ˇzil´ach naopak okysliˇcena. Plicn´ı obˇeh je s´eriovˇe zaˇrazen za syst´emov´ y obˇeh, se kter´ ym je spojen v srdci. V naˇsem modelu plicn´ı obˇeh zanedb´av´ame a pˇredpokl´ad´ame jeho norm´aln´ı fyziologickou funkci.
3.3
Velk´ y syst´ emov´ y obˇ eh
Velk´ y syst´emov´ y obˇeh distribuuje krev ke vˇsem buˇ nk´am lidsk´eho tˇela. B´ yv´a tak´e oznaˇcov´an jako vysokotlak´ y, nebot’ hodnoty systolick´eho tlaku v arteri´ıch dosahuj´ı mnohem vyˇsˇs´ıch hodnot neˇz u arteri´ı plicn´ıho obˇehu, norm´aln´ı hodnota se pohybuje okolo 120 mmHg. Mal´ y plicn´ı obˇeh
- 19 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
Krev je z lev´e komory vypuzena do vzestupn´e aorty. Z koˇrene aorty vystupuj´ı koron´arn´ı tepny, kter´e z´asobuj´ı krv´ı samotn´e srdce. Vzestupn´a aorta se n´aslednˇe staˇc´ı doleva a dol˚ u a vytv´aˇr´ı typick´ y aort´aln´ı oblouk. Z aort´aln´ıho oblouku se vˇetv´ı tˇri tepny, truncus brachiocephalicus, lev´a tepna podkl´ıˇckov´a a lev´a krkavice. Truncus brachiocephalicus se d´ale vˇetv´ı na pravou podkl´ıˇckovou tepnu a pravou krkavici. Sestupn´a aorta, kter´a je pokraˇcov´an´ı aort´aln´ıho oblouku, pˇrech´az´ı v hrudn´ı aortu - aorta thoracica. Aorta thorarica z´asobuje krv´ı ˇca´st meziˇzebern´ıch sval˚ u, pˇredn´ı u ´sek bˇriˇsn´ıch sval˚ u, ˇca´st br´anice, p´ateˇrn´ı kan´al, m´ıchu a m´ıˇsn´ı obaly. D´ale potom pˇrech´az´ı v bˇriˇsn´ı aortu, kter´a rozv´ad´ı krev ke vˇsem bˇriˇsn´ım org´an˚ um a p´anevn´ım org´an˚ um, z´asobuje z´adov´e svaly, ˇca´st br´anice, bˇriˇsn´ı stˇeny a dalˇs´ı. N´asleduje aort´aln´ı bifurkace, kde se bˇriˇsn´ı aorta dˇel´ı na dvˇe vˇetvˇe, kaˇzd´a pokraˇcuje jako zevn´ı kyˇceln´ı tepna, a d´ale pak jako stehenn´ı tepna a pˇredn´ı a zadn´ı holenn´ı tepna, kter´e vyˇzivuj´ı doln´ı konˇcetiny. Vˇsechny tepny se postupnˇe d´ale dˇel´ı a zmenˇsuj´ı na tep´enky, neboli arterioly, kter´e se opˇet dˇel´ı a zmenˇsuj´ı aˇz na u ´roveˇ n vl´aseˇcnic - kapil´ar. Samotn´e vl´aseˇcnice jsou velice tenk´e c´evy, kter´e jsou tvoˇreny jen jedinou vrstvou bunˇek. Jejich pr˚ umˇer je tak´e minim´aln´ı, pohybuje se od 5 do 20 µm. Toto z´ uˇzen´ı vede k tomu, ˇze ˇcerven´e krvinky proud´ı v kapil´ar´ach v ˇradˇe za sebou, kaˇzd´a m´a maxim´aln´ı kontakt se stˇenami vl´aseˇcnice a doch´az´ı k nejefektivnˇejˇs´ı v´ ymˇenˇe l´atek mezi tk´anˇemi a krv´ı. Kapil´ary se na druh´e stranˇe opˇet spojuj´ı do venul (ˇzilek), kter´e se d´ale spojuj´ı do ˇ ıly z doln´ıch konˇcetin a trupu se nakonec spojuj´ı do doln´ı dut´e ˇz´ıly, ˇz´ıly z horn´ı ˇzil. Z´ ˇca´sti tˇela pak do horn´ı dut´e ˇz´ıly. Tyto ˇz´ıly pak vstupuj´ı do prav´e s´ınˇe srdeˇcn´ı, kde se cel´ y syst´em uzav´ıra. Za zm´ınku stoj´ı, ˇze ˇz´ıly doln´ıch konˇcetin obsahuj´ı chlopnˇe, kter´e zabraˇ nuj´ı zpˇetn´emu toku krve pˇri diastole v d˚ usledku vysok´eho hydrostatick´eho tlaku.
3.3.1
Art´ erie a arterioly
Anatomie art´eri´ı je zn´azornˇena na obr´azku 3.3a. Vnitˇrn´ı plochu arteri´ı vyst´ yl´a jednovrstevn´ y endotel. Ten je obalen vnitˇrn´ı elastickou vrstvou, tvoˇrenou pletivem elastinov´ ych a kolagenov´ ych bunˇek. N´asleduje relativnˇe tlust´a vrstva hladk´e svaloviny. Kontrakc´ı nebo relaxac´ı t´eto svaloviny lze ovlivˇ novat efektivn´ı pr˚ usvit art´erie, ˇcehoˇz se vyuˇz´ıv´a pro udrˇzen´ı st´al´eho stˇredn´ıho arteri´aln´ıho tlaku a poˇzadovan´eho toku. Vrstva svaloviny je obalena vnˇejˇs´ı elastickou vrstvou, opˇet tvoˇrenou elastinov´ ymi a kolagenov´ ymi vl´akny. Z vnˇejˇs´ı strany je potom art´erie chr´anˇena vnˇejˇs´ım obalem. Syst´emov´e art´erie lze podle [1] a [3] zjednoduˇsenˇe rozdˇelit do tˇr´ı skupin, z nichˇz kaˇzd´a m´a svoj´ı specifickou funkci, spoleˇcnou funkc´ı je potom distribuce okysliˇcen´e krve. Prvn´ı skupinou jsou velk´e art´erie v bl´ızkosti srdce (pruˇzn´ıkov´e tepny). Tyto art´erie maj´ı nejˇsirˇs´ı pr˚ umˇer (do 2 cm), jsou velice elastick´e, vrstva hladk´e svalovinu je u nich relativnˇe tenˇc´ı. Slouˇz´ı jako akumul´ator krve vypuzen´e pˇri systole. Energie krve uvolnˇen´e pˇri systole se nahromad´ı ve stˇen´ach art´eri´ı a ty se rozˇs´ıˇr´ı. Pˇri diastole se d´ıky svoj´ı elasticitˇe opˇet smrˇst´ı, ˇc´ımˇz doch´az´ı k toku krve i bˇehem diastoly. Druhou skupinou jsou pak dlouh´e svalnat´e art´erie a arterioly (rezistenˇcn´ı tepny a tep´enky). Tyto slouˇz´ı pˇredevˇs´ım k ˇr´ızen´ı redistribuci krve. D´ıky vyˇsˇs´ımu obsahu hladk´e svaloviny jsou tak´e schopny ovlivˇ novat profil toku a rozhodovat o distribuci minutov´eho v´ ydeje srdeˇcn´ıho k jednotliv´ ym org´an˚ um. Tˇret´ı skupinou jsou prekapil´arn´ı sfinktery, kter´e maj´ı minim´aln´ı pr˚ umˇer 20 µm. Maj´ı Art´erie a arterioly
- 20 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
tak´e v´ yraznou vrstvu hladk´eho svalstva, a svoj´ı konstrikc´ı ovlivˇ nuj´ı perfuzi jednotliv´ ych org´an˚ u, konkr´etnˇe ˇr´ıd´ı kter´e kapil´ary budou prokrveny a kter´e ne.
(a) Art´erie
(b) V´ena
Obr´azek 3.3: Anatomie art´erie a v´eny, pˇrevzato s u ´pravami z [14]
3.3.2
Kapil´ ary a mikrocirkulace
Vezmeme-li jako hlavn´ı funkci kardiovaskul´arn´ıho syst´emu dopravu l´atek a plyn˚ u z a do tk´an´ı, potom jsou kapil´ary nejd˚ uleˇzitˇejˇs´ı ˇca´st´ı kardiovaskul´arn´ıho syst´emu, nebot’ pr´avˇe v nich doch´az´ı k v´ ymˇenˇe l´atek. Kapil´ary jsou tvoˇreny jednovrstevn´ ym epitelem. Jejich pr˚ umˇer se pohybuje v rozmez´ı 5-20 µm, d´elka do 750 µm. Proudˇen´ı krve v kapil´ar´ach je konstantn´ı a relativnˇe pomal´e. Pro perfuzi kapil´ar maj´ı ale velk´ y v´ yznam i arterioly a venuly, dohromady se tedy ˇcasto mluv´ı o jednom funkˇcn´ım celku: mikrocirkulaci. Pro regulaci perfuze kapil´ar slouˇz´ı tzv. sfinktery, kter´e se nach´az´ı u u ´st´ı kapil´ar na arterioln´ı stranˇe. Jejich uzavˇren´ı nebo otevˇren´e pˇr´ımo rozhoduje o tom, zda krev danou kapil´arou poteˇce ˇci nikoli. D´ale se ˇcasto mezi arteriolami a venulami nach´az´ı arterioven´ozn´ı zkraty (anastom´ozy), kter´e pˇri otevˇren´ı umoˇzn ˇuj´ı obej´ıt celou kapil´arn´ı s´ıt’ nach´azej´ıc´ı se za anastom´ozou. Pro v´ ymˇenu l´atek mezi krv´ı a okoln´ı tk´an´ı je d˚ uleˇzit´a rychlost pr˚ utoku kapil´arou a plocha kapil´ary. Rychlost toku krve v kapil´ar´ach se pohybuje ˇra´dovˇe v mm/s. Velikost povrchu jedn´e kapil´ary (vezmeme-li v u ´vahu pr˚ umˇern´e hodnoty) potom dosahuje 14 000 µm2 ([3]). L´atkov´e v´ ymˇeny se u ´ˇcastn´ı i ˇca´st venuly, tud´ıˇz celkov´a plocha jednoho u ´seku 2 vzroste na 28 000 µm . Odhaduje se, ˇze celkov´ y poˇcet kapil´ar v lidsk´em tˇele se pohybuje okolo 40 miliard. Celkovou plochu pro v´ ymˇenu l´atek mez´ı krv´ı a tk´anˇemi tedy dost´av´ame 1000 m2 . Je vˇsak potˇreba si uvˇedomit, ˇze krev prot´ek´a souˇcasnˇe zhruba jednou tˇretinou vˇsech kapil´ar, obzvl´aˇstˇe je-li ˇclovˇek v klidu. Dalˇs´ı v´ yznamnou ˇca´st z mnoˇzstv´ı kapil´ar tvoˇr´ı kapil´ary k˚ uˇze, kter´e se vyuˇz´ıvaj´ı pˇredevˇs´ım pro termoregulaci organismu.
Kapil´ary a mikrocirkulace
- 21 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
3.3.3
V´ eny a venuly
V´eny a venuly maj´ı podobnou stavbu jako art´erie, ilustrace jejich anatomie je na obr´azku 3.3b. Narozd´ıl od art´eri´ı ale maj´ı vˇetˇs´ı pr˚ usvit, jsou m´enˇe elastick´e a nemaj´ı tak silnou vrstvu hladk´e svaloviny. V mnoha v´en´ach se tak´e nach´az´ı chlopnˇe, kter´e zabraˇ nuj´ı zpˇetn´emu toku krve, ke kter´emu by mohlo doch´azet v d˚ usledku hydrostatick´eho tlaku krve. Podle [3] lze i v´eny rozdˇelit do dvou kategori´ı. Prvn´ı jsou rezistenˇcn´ı venuly, kter´e maj´ı silnˇejˇs´ı vrstvu hladk´e svaloviny. Pomˇer mezi konstrikc´ı rezistenˇcn´ıch arteriol a venul ovlivˇ nuje rychlost pr˚ utoku krve kapil´arami a hydrostatick´ y tlak v kapil´ar´ach, ˇc´ımˇz se d´a ˇr´ıdit rychlost filtrace a resorpce l´atek v kapil´ar´ach, Druhou kategori´ı jsou potom kapacitn´ı v´eny. Jsou m´enˇe elastick´e neˇz art´erie, jejich pr˚ usvit je vˇetˇs´ı (maj´ı slabˇs´ı vrstvu hladk´e svaloviny) a dok´aˇzou pojmout velk´e mnoˇzstv´ı krve a slouˇz´ı jako jej´ı rezervo´ar. Tlak ve v´en´ach je mnohem niˇzˇs´ı (kolem 20 mmHg) neˇz v art´eri´ıch. Proudˇen´ı je lehce pulsatiln´ı, coˇz je zp˚ usobeno ˇcinnost´ı prav´eho srdce.
3.4
Krev
Krev je suspenze skl´adaj´ıc´ı se z tekut´e plasmy (54.3%) a pevn´e smˇesi erytrocyt˚ u (ˇcerven´e krvinky), leukocyt˚ u (b´ıl´e krvinky) a krevn´ıch destiˇcek (45,7%). Tato smˇes se oznaˇcuje jako pln´a krev, pˇriˇcemˇz erytrocyty tvoˇr´ı pˇribliˇznˇe 45% pln´e krve. Plasma je z 92% tvoˇrena voˇ dou, zbyl´ ych 8% tvoˇr´ı proteiny a rozpuˇstˇen´e miner´aly. Cerven´ e krvinky slouˇz´ı pˇredevˇs´ım k pˇrepravˇe kysl´ıku ke tk´an´ım. Maj´ı diskovit´ y bikonk´avn´ı tvar, jejich pr˚ umˇer se pohybuje v rozmez´ı 6-8 µm, jejich tlouˇst’ka v nejtlustˇs´ım m´ıstˇe (na okraji) dosahuje 2 µm. Nemaj´ı bunˇeˇcn´e j´adro, jsou relativnˇe elastick´e a tak se mohou v pˇr´ıpadˇe potˇreby prot´ahnout (pro lepˇs´ı prostupnost kapil´arami). Leukocyty jsou o nˇeco vˇetˇs´ı neˇz erytrocyty, zpravidla kulovit´eho tvaru. Jsou souˇc´ast´ı imunitn´ıho syst´emu. Rozliˇsujeme pˇet typ˚ u leukocyt˚ u, pˇriˇcemˇz kaˇzd´ y kon´a specifickou funkci. Krevn´ı destiˇcky slouˇz´ı pˇredevˇs´ım k zastaven´ı krv´acen´ı.
3.4.1
Proudˇ en´ı krve
Charakter proudˇen´ı krve v c´evn´ım ˇreˇciˇsti se postupnˇe mˇen´ı z pulsn´ıho na kontinu´aln´ı. To je nejl´epe ilustrov´ano n´asleduj´ıc´ım obr´azkem:
V´eny a venuly
- 22 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
Obr´azek 3.4: Charakter proudˇen´ı na jednotliv´ ych u ´rovn´ıch c´evn´ıho ˇreˇciˇstˇe. Pˇrevzato s u ´pravami z [17] pulsn´ı charakter ve velk´ ych c´ev´ach je d´an ˇcinnost´ı srdce. Art´erie a arterioly jsou elastick´e, pˇri systole akumuluj´ı ve sv´ ych stˇen´ach energii pulsu, kterou uvolˇ nuj´ı v diastolick´e f´azi. S klesaj´ıc´ım pr˚ umˇerem c´evy se zvyˇsuje odpor kladen´ y tekutinˇe. Nav´ıc celkov´a plocha, kterou krev prot´ek´a, se na niˇzˇs´ıch u ´rovn´ıch zvyˇsuje (pr˚ umˇer vˇsech kapil´ar je 700kr´at vˇetˇs´ı neˇz pr˚ umˇer vˇsech art´eri´ı), ˇc´ımˇz se zpomaluje tok a sniˇzuje se pulsn´ı tlak aˇz do t´e m´ıry, ˇze v kapil´ar´ach je proudˇen´ı jiˇz kontinu´aln´ı. Kontinu´aln´ı tok krve se potom spojuje ve venul´ach a v´en´ach a vrac´ı do srdce. Proudˇen´ı ve v´en´ach je lehce pulsatiln´ı. Tato pulsatilita je zp˚ usobena ˇcinnost´ı chlopnˇe prav´e s´ınˇe (tok krve je v pravideln´ ych intervalech zastavov´an). pulssatilita ve ven´ozn´ım ˇreˇciˇsti je ale v porovn´an´ı s pulsatilitou v arteri´aln´ım ˇreˇciˇsti zanedbateln´a. Mimo celkov´ y charakter proudˇen´ı krve je i d˚ uleˇzit´e, zda krev proud´ı lamin´arnˇe ˇci turbulentnˇe. Rozd´ıl mezi lamin´arn´ım a turbulentn´ım proudˇen´ım ilustruje n´asleduj´ıc´ı obr´azek:
Obr´azek 3.5: Ilustrace lamin´arn´ıho a turbulentn´ıho proudˇen´ı. Pˇrevzato z [14] Za norm´aln´ıch podm´ınek proud´ı krev v c´ev´ach lamin´arnˇe, coˇz znamen´a, ˇze vˇsechny Proudˇen´ı krve
- 23 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
molekuly se pohybuj´ı rovnobˇeˇznˇe s pod´elnou osou c´evy (aˇckoli se z d˚ uvodu viskozity pohybuj´ı jinou rychlost´ı). Tento typ proudˇen´ı je v´ yhodn´ y z nˇekolika d˚ uvod˚ u. Pˇrednˇe pr˚ utok je pˇr´ımo u ´mˇern´ y tlakov´emu gradientu v kapalinˇe. Dalˇs´ı v´ yhodou je, ˇze jednotliv´e sloˇzky krve spolu tolik nekoliduj´ı, stejnˇe tak nedoch´az´ı k u ´der˚ um v´ıru na stˇeny c´evy, tud´ıˇz nedoch´az´ı k mechanick´emu poˇskozen´ı krve a c´ev. Za urˇcit´ ych podm´ınek pˇrech´az´ı lamin´arn´ı proudˇen´ı na turbulentn´ı. Tento typ proudˇen´ı je charakteristick´ y t´ım, ˇze molekuly se nepohybuj´ı po rovnobˇeˇzn´ ych drah´ach, naopak, mohou se pohybovat vˇsemi smˇery. Pr˚ utok krve jiˇz nen´ı pˇr´ımo u ´mˇern´ y tlakov´emu gradientu, ale sp´ıˇse jeho druh´e odmocninˇe [3]. Tud´ıˇz na dvojn´asobn´e zv´ yˇsen´ı pr˚ utoku je potˇreba tlak zv´ yˇsit ˇctyˇrikr´at. To pˇredstavuje zbyteˇcnou z´atˇeˇz pro cel´ y organismus. Zda bude proudˇen´ı sp´ıˇse lamin´arn´ı nebo sp´ıˇse turbulentn´ı se d´a zjistit pomoc´ı Reynoldsova ˇc´ısla. Jedn´a se o bezrozmˇernou veliˇcinu, kter´a se vypoˇc´ıt´a n´asleduj´ıc´ım zp˚ usobem [1]: vs · D · ρ (3.5) Re = µ kde vs je pr˚ umˇern´a rychlost proudˇen´ı [m/s], D je pr˚ umˇer c´evy [m], ρ je hustota krve 3 [kg/m ], aµ je dynamick´a viskozita krve [Pa/s] Z parametr˚ u, kter´e ovlivˇ nuj´ı Reynoldsovo ˇc´ıslo m˚ uˇze organismus ovlivˇ novat jen pr˚ umˇer c´evy a rychlost proudˇen´ı. Jako kritick´a hodnota Reynoldsova ˇc´ısla se nejˇcastˇeji ud´av´a Re=2320. Do t´eto hranice lze povaˇzovat proudˇen´ı za lamin´arn´ı. Pro hodnoty Reynoldosva ˇc´ısla vˇetˇs´ı neˇz 2320 a menˇs´ı neˇz 4000 se hovoˇr´ı o tzv. pˇrechodov´e oblasti, kdy proudˇen´ı nen´ı ˇcistˇe lamin´arn´ı, tvoˇr´ı se v nˇem v´ıry, ale jeˇstˇe se ned´a povaˇzovat za turbulentn´ı proudˇen´ı. Pokud Reynoldsovo ˇc´ıslo pˇres´ahne 4000, d´a se jiˇz mluvit o turbulentn´ım proudˇen´ı. V lidsk´em tˇele doch´az´ı k turbulentn´ımu proudˇen´ı ve velk´ ych c´ev´ach, kdy krev teˇce vysokou rychlost´ı, popˇr´ıpadˇe se m˚ uˇze vyskytnout v c´ev´ach postiˇzen´ ych ateroskl´er´ozou, kdy se do stˇeny c´evy ukl´adaj´ı tukov´e l´atky, coˇz vyvol´a z´ uˇzen´ı c´evy, zv´ yˇsen´ı rychlosti pr˚ utoku a vznik v´ır˚ u v toku.
3.4.2
Fyzik´ aln´ı vlastnosti krve
Na plnou krev je moˇzno z hlediska dynamiky kapalin nahl´ıˇzet jako na ne-newtonovskou kapalinu, a to zejm´ena kv˚ uli dynamick´e viskozitˇe krve. Viskozita pˇredstavuje odpor, kter´ ym kapalina p˚ usob´ı proti s´ıle snaˇz´ıc´ı se ji uv´est do pohybu. Jak bylo vysvˇetleno v pˇredchoz´ı podkapitole, lze pˇredpokl´adat, ˇze proudˇen´ı krve v c´ev´ach m´a lamin´arn´ı charakter, to znamen´a ˇze krev se pohybuje v soustˇredn´ ych kruhov´ ych vrstv´ach. Jednotliv´e vrstvy se vˇsak pohybuj´ı r˚ uznou rychlost´ı, vznik´a mezi nimi tˇren´ı. Pokud toto tˇren´ı vzt´ahneme na jednotku plochy, dostaneme smykov´e napˇet´ı. Rozd´ıl rychlosti sousedn´ıch vrstev potom pˇredstavuje smykov´a rychlost. Viskozitu kapalin definuje vztah: viskozita = smykov´e napˇet´ı/smykov´a rychlost (3.6) Pro newtonovsk´e kapaliny plat´ı, ˇze viskozita je konstantn´ı. Pro ne-newtonovsk´e kapaliny tento pˇredpoklad nelze pˇrijmout a mus´ıme uvaˇzovat viskozitu v urˇcit´em rozpˇet´ı. Z hlediska stavby modelu n´as u krve zaj´ım´a pˇredevˇs´ım jej´ı hustota (kter´a ovlivˇ nuje setrvaˇcnost toku) a viskozita. V n´asleduj´ıc´ı tabulce jsou porovn´any tyto vlastnosti krve s vodou:
Fyzik´aln´ı vlastnosti krve
- 24 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
3
Hustota [kg/m ] Dynamick´a viskozita [Pa/s]
Krev 1060 (3-4)·10−3
Voda pˇri 20◦ C 998,2 1·10−3
Tabulka 3.3: Vlastnosti krve Jak je vidˇet, hustota krve je podobn´a hustotˇe vody, nicm´enˇe dynamick´a viskozit je tˇrikr´at aˇz ˇctyˇrikr´at vˇetˇs´ı. Pˇresto v naˇsem modelu budeme jako medium pouˇz´ıvat pˇrev´aˇznˇe vodu.
3.5
Regulaˇ cn´ı mechanismy kardiovaskul´ arn´ıho syst´ emu
N´aroky jednotliv´ ych org´an˚ u a tk´an´ı na z´asobov´an´ı krv´ı se velice liˇs´ı a v´ yraznˇe se mˇen´ı v z´avislosti na aktu´aln´ım vyt´ıˇzen´ı dan´eho org´anu (tk´anˇe). Kardiovaskul´arn´ı syst´em mus´ı zaˇr´ıdit dostateˇcnou perfuzi org´an˚ u, ale z´aroveˇ n mus´ı udrˇzet svoj´ı z´akladn´ı funkci, kterou je pr˚ utok krve c´evami. Pr˚ utok krve c´evou (nebo c´evn´ım syst´emem) lze popsat analogi´ı Ohmova z´akona takto: ∆P (3.7) Q= R kde Q je je tok krve c´evou [m3 /s], ∆P je rozd´ıl tlak˚ u na zaˇc´atku a konci c´evy [Pa], R je periferi´aln´ı odpor, kter´y c´eva klade toku krve [P a · s/m3 ] Je tedy jasn´e, ˇze tlak a tok krve a odpor c´ev na sobˇe z´avis´ı a zmˇena jednoho parametru vyvol´a zmˇenu ostatn´ıch. Z tˇechto tˇr´ı parametr˚ u mohou ˇr´ıd´ıc´ı mechanismy kardiovaskul´arn´ıho syst´emu pˇr´ımo ovlivnit tok, a to zv´ yˇsen´ım (minutov´eho) srdeˇcn´ıho v´ ydeje, a odpor c´ev. Srdeˇcn´ı v´ ydej lze zmˇenit srdeˇcn´ım tepem ˇci zmˇenou tepov´eho objemu (nebo kombinac´ı ob´eho), odpor c´ev potom jejich konstrikc´ı ˇci dilatac´ı. Regulaˇcn´ı mechanismy m˚ uˇzeme z hlediska jejich p˚ usoben´ı rozdˇelit na m´ıstn´ı a syst´emov´e [3].
3.5.1
M´ıstn´ı regulaˇ cn´ı mechanismy
M´ıstn´ı regulaˇcn´ı mechanismy (ˇcasto oznaˇcov´any jako autoregulaˇcn´ı mechanismy) maj´ı na starosti pˇredevˇs´ım zajiˇstˇen´ı dostateˇcn´e perfuze tk´anˇe podle jej´ıch aktu´aln´ıch metabolick´ ych poˇzadavk˚ u, a d´ale udrˇzen´ı lamin´arn´ıho profilu a konstantn´ıho objemu toku krve pˇri zmˇen´ach tlaku krve. Mezi m´ıstn´ı mechanismy patˇr´ı myogenn´ı autoregulace, endotelov´a regulace a metabolick´a regulace. Myogenn´ı autoregulace spoˇc´ıv´a v kontrakci hladk´e svaloviny c´evy pˇri zv´ yˇsen´ı krevn´ıho tlaku (naopak, relaxaci pˇri sn´ıˇzen´ı tlaku). Kontrakce hladk´e svaloviny vede ke zmenˇsen´ı pr˚ usvitu c´evy a t´ım zv´ yˇsen´ı jej´ıho odporu (relaxace naopak odpor sniˇzuje). Z rovnice 3.7 potom vypl´ yv´a, ˇze pˇri odpov´ıdaj´ıc´ı zmˇenˇe odporu c´evy pˇri zmˇenˇe tlakov´eho gradientu z˚ ustane pr˚ utok c´evou stejn´ y. Endotelov´a regulace je zaloˇzena na sekreci oxidu dusnat´eho (NO) endotelem pˇri zv´ yˇsen´ı rychlosti proudˇen´ı krve v c´evˇe. Oxid dusnat´ y m´a vazodilataˇcn´ı u ´ˇcinek. Zvˇetˇsen´ım pr˚ uˇrezu plochy, kterou krev prot´ek´a, se sn´ıˇz´ı jej´ı rychlost. Tato reakce m´a za c´ıl zabr´anit vzniku turbulentn´ıho proudˇen´ı, kter´e by mohlo mechanicky poˇskodit endotel. Tento jev popisuje
Regulaˇcn´ı mechanismy kardiovaskul´arn´ıho syst´emu
- 25 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
rovnice kontinuality: v1 S2 = v2 S1 kde v1 , respektive v2 , je rychlost toku v trubici o ploˇse pr˚ uˇrezu S1 , respektive S2
(3.8)
Metabolick´e regulaˇcn´ı mechanismy jsou zaloˇzeny na zmˇenˇe koncentrace r˚ uzn´ ych l´atek z d˚ uvodu zv´ yˇsen´eho metabolismu org´anu. Pokud v c´ev´ach dan´eho org´anu klesne koncentrace kysl´ıku nebo stoupne koncentrace oxidu uhliˇcit´eho ˇci jin´ ych metabolick´ ych produkt˚ u, zp˚ usob´ı to vazodilataci c´ev a t´ım i zv´ yˇsen´ı perfuze org´anu.
3.5.2
Syst´ emov´ e regulaˇ cn´ı mechanismy
Syst´emov´e regulaˇcn´ı mechanismy ˇr´ıd´ı celkov´ y perif´ern´ı odpor a minutov´ y srdeˇcn´ı v´ ydej srdce tak, aby byl zachov´an tlakov´ y gradient nutn´ y pro tok krve v c´ev´ach. Z hlediska rychlosti a doby trv´an´ı u ´ˇcinku je dˇel´ıme na rychl´e (kr´atkodob´e) a pomal´e (dlouhodob´e). Mezi rychl´e regulaˇcn´ı mechanismy patˇr´ı pˇredevˇs´ım baroreflex a hormon´aln´ı syst´emy adrenalinnoradrenalin a renin-angiotenzin. Baroreflex (jak n´azev napov´ıd´a) je vyvol´an zmˇenou tlaku krve. Nejv´ıce baroroceptor˚ u je uloˇzeno ve velk´ ych tepln´ach, pˇredevˇs´ım aort´aln´ım oblouku. Tyto receptory vys´ılaj´ı do center v mozkov´em kmeni sign´al, jehoˇz frekvence z´avis´ı na tlaku krve. Zv´ yˇs´ı-li se krevn´ı tlak, zv´ yˇs´ı se i frekvence sign´alu. To vyvol´a odpovˇed’ tˇechto center ve formˇe sn´ıˇzen´ı tepov´e frekvence, sn´ıˇzen´ı tepov´eho objemu a vazodilatace, ˇc´ımˇz krevn´ı tlak opˇet klesne. Adrenalin a noradrenalin patˇr´ı mezi hormony vyluˇcovan´e dˇren´ı nadledvinek (katecholaminy), kter´e ovlivˇ nuj´ı jak tonus c´ev, tak ˇcinnost srdce. Reakce c´ev a srdce na tyto hormony je vˇsak pomˇernˇe sloˇzit´a, nebot’ receptory na nˇe reaguj´ıc´ı jsou dva (alfa a beta), pˇriˇcemˇz kaˇzd´ y m´a jeˇsta dva podtypy (alfa1 , alfa2 , beta1 , beta2 ). Aktivace alfa receptor˚ u vyvol´av´a vazokonstrikci, aktivace beta receptor˚ u potom vazodilataci. Zat´ımco na noradrenalin reaguj´ı jen alfa receptory, na adrenalin reaguj´ı oba receptory. Z´aleˇz´ı tedy na koncentraci receptor˚ u v dan´e c´evˇe, jak´ ym smˇerem se zmˇen´ı jej´ı tonus. Zat´ımco v nˇekter´ ych c´ev´ach se vyvol´a vazokonstrikce (c´evy k˚ uˇze, tr´avic´ı soustavy) v jin´ ych dojde naopak k vazodilataci (kostern´ı svalstvo). Jednou z funkc´ı tohoto syst´emu je totiˇz pˇripravit tˇelo na svalov´ y v´ ykon. Adrenalin tak´e ovlivˇ nuje ˇcinnost srdce, a to tak, aby doˇslo ke zv´ yˇsen´ı minutov´eho v´ ydeje srdeˇcn´ıho. Pomalejˇs´ım regulaˇcn´ım mechanismem je syst´em renin-angiotenzin. Hlavn´ı u ´lohou tohoto syst´emu je zajistit dostateˇcn´ y tlak krve v ledvin´ach (tlak krve zde pˇr´ımo souvis´ı s efektivitou ledvin pˇri filtraci krve). Poklesne-li tlak krve v ledvin´ach, uvoln´ı se do krve hormon renin. Ten ˇstˇep´ı angiotenzinogen, uvolˇ novan´ y j´atry, na angiotenzin I, kter´ y je d´ale pˇremˇenˇen na angiotenzin II (k tomu doch´az´ı zejm´ena v plic´ıch). Angiotenzin II m´a vazokonstrikˇcn´ı u ´ˇcinek, coˇz vede ke zv´ yˇsen´ı perifern´ıho odporu a t´ım i zv´ yˇsen´ı tlaku krve v ledvin´ach. Tohoto efektu je dosaˇzeno zhruba po dvaceti minut´ach. Nejvˇetˇs´ı v´ yznam m´a tento syst´em pˇri celkov´em poklesu syst´emov´eho tlaku, napˇr´ıklad v d˚ usledku krv´acen´ı. Nejv´ yznamnˇejˇs´ım pomal´ ym regulaˇcn´ım mechanismem je ˇr´ızen´ı celkov´eho objemu krve v tˇele, souvis´ı tedy s ˇr´ızen´ım objemu vody v tˇele. Pˇri n´ızk´em tlaku krve dojde ke zv´ yˇsen´ı produkce antidiuretick´eho hormonu, ˇc´ımˇz vzroste resorpce vody v glomerulech ledvin. T´ım se zvˇetˇsuje objem krve (respektive krevn´ı plasmy). Antidiuretick´ y hormon m´a nav´ıc vazokonstrikˇcn´ı u ´ˇcinky. Kombinac´ı tˇechto faktor˚ u doch´az´ı k zv´ yˇsen´ı tlaku krve. Naopak, pˇri vyˇsˇs´ım tlaku se zvˇetˇsuje i filtraˇcn´ı tlak ledvin, ˇc´ımˇz doch´az´ı k vˇetˇs´ımu vyluˇcov´an´ı Syst´emov´e regulaˇcn´ı mechanismy
- 26 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
vody, omezen´ı sekrece antidiuretick´eho hormonu, sn´ıˇzen´ı objemu krve a poklesu krevn´ıho tlaku.
3.6
Hemodynamick´ e parametry
Aktivitu kardiovaskul´arn´ıho syst´emu prov´azej´ı ˇcetn´e projevy, kter´e lze mˇeˇrit a pozorovat. Tyto projevy lze mˇeˇrit bud’ neinvazivnˇe, nebo invazivnˇe. Nejjednoduˇsˇs´ım neinvazivn´ım mˇeˇren´ım je mˇeˇren´ı tepov´e frekvence. Dalˇs´ımi neinvazivn´ımi metodami lze mˇeˇrit elektrickou aktivitu srdce (elektrokardiografie), systolick´ y a diastolick´ y krevn´ı tlak (auskultaˇcn´ı ˇci oscilometrick´a metoda), pomoc´ı sledov´an´ı srdeˇcn´ıch ozev lze pozorovat funkci srdeˇcn´ıch chlopn´ı (fonokardiografie), pr˚ utok krve danou c´evou (sonografie) ˇci prokrven´ı dan´e oblasti (pletysmografie), a mnoho dalˇs´ıch. Invazivn´ımi metodami lze zpravidla dos´ahnout pˇresnˇejˇs´ıch v´ ysledk˚ u, ovˇsem za cenu z´asahu do organismu. Invazivn´ı metody zahrnuj´ı mˇeˇren´ı krevn´ıho tlaku v jinak nedostupn´ ych m´ıstech (komorov´ y tlak), diluˇcn´ımi metodami lze mˇeˇrit srdeˇcn´ı v´ ydej, tepov´ y objem, rychlost proudˇen´ı, d´ale lze mˇeˇrit pr˚ utok krve (indukˇcn´ı pr˚ utokomˇery) a mnoho dalˇs´ıch parametr˚ u.
3.6.1
Tlak krve a jeho mˇ eˇ ren´ı
Jedn´ım ze z´akladn´ıch a relativnˇe snadno mˇeˇriteln´ ych parametr˚ u kardiovaskul´arn´ıho syst´emu je krevn´ı tlak na u ´rovni velk´ ych tepen. Nejrozˇs´ıˇrenˇejˇs´ı metodou neinvazivn´ıho mˇeˇren´ı krevn´ıho tlaku je auskultaˇcn´ı metoda. Toto mˇeˇren´ı se prov´ad´ı pomoc´ı rtut’ov´eho manˇzetov´eho tonometru a fonendoskopu. Rtut’ se v tonometru pouˇz´ıv´a proto, ˇze m´a v´ yraznˇe vyˇsˇs´ı hustotu neˇz krev (pˇribliˇznˇe 13,6 kr´at vyˇsˇs´ı, pokud by se tedy pouˇz´ıvala tekutina s hustou bl´ızkou vody, mˇel by sloupec tekutiny odpov´ıdaj´ıc´ı systolick´emu tlaku v´ yˇsku 155 cm). Z tohoto d˚ uvodu se zmˇeˇren´e v´ ysledky uv´adˇej´ı ve speci´aln´ı jednotce mmHg (milimetr rtut’ov´eho sloupce), kter´a odpov´ıd´a pˇribliˇznˇe 133,3 Pa. Mˇeˇren´ı se nejˇcastˇeji prov´ad´ı na lev´e brachi´aln´ı tepnˇe (paˇzn´ı tepna) pˇribliˇznˇe v u ´rovni srdce. Manˇzeta tonometru se obtoˇc´ı okolo paˇze, n´aslednˇe se zv´ yˇs´ı tlak v manˇzetˇe v´ yraznˇe nad oˇcek´avan´ y systolick´ y tlak (200 mmHg). T´ım dojde k zaˇskrcen´ı tepny a zastaven´ı pr˚ utoku krve. N´aslednˇe se tlak manˇzety sniˇzuje aˇz k m´ıstu, kdy opˇet zaˇcnˇe krev tepnou proudit. Protoˇze m´a ale tepna v´ yraznˇe sn´ıˇzen´ y pr˚ umˇer, proudˇen´ı je turbulentn´ı. V´ıry v proudu se projev´ı tzv. Korotkovov´ ymi zvuky, kter´e se daj´ı zachytit fonendoskopem. Tlak manˇzety v dobˇe, kdy se objev´ı Korotkovovy zvuky odpov´ıd´a systolick´emu tlaku. Tlak manˇzety se n´aslednˇe st´ale sniˇzuje aˇz do chv´ıle, kdy proudˇen´ı nabyde sv´e fyziologick´e lamin´arn´ı podoby. T´ım odezn´ı Korotkovovy zvuky, tlak v manˇzetˇe tonometru odpov´ıd´a diastolick´emu tlaku.
Hemodynamick´e parametry
- 27 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
Obr´azek 3.6: Typick´ y pr˚ ubˇeh tlakov´e kˇrivky mˇeˇren´ y na brachi´aln´ı art´erii. Obr´azek pˇrevzat su ´pravami z [14] Na obr´azku v´ yˇse je zn´azornˇena typick´a kˇrivka arteri´aln´ıho tlaku. Dobˇre viditeln´ y je tzv. dikrotick´ y z´aˇrez (tlakov´ y pokles mezi prvn´ı a druhou vlnou tlakov´e kˇrivky). Druh´a vlna n´asleduj´ıc´ı za dikrotick´ ym z´aˇrezem je zp˚ usobena smrˇstˇen´ım tepen, kter´e d´ıky sv´e elasticitˇe akumuluj´ı energii pulsu, kterou n´aslednˇe uvolˇ nuj´ı a t´ım pom´ahaj´ı k zachov´an´ı tlakov´eho gradientu v c´evˇe. T´ım doch´az´ı k toku krve i mimo ejekˇcn´ı f´azi srdce.
3.6.2
Pr˚ utok krve a jeho mˇ eˇ ren´ı
Mˇeˇren´ı pr˚ utoku krve c´evou jiˇz nen´ı tak jednoduch´e. Neinvazivnˇe lze mˇeˇrit pr˚ utok pomoc´ı ultrazvukov´e sondy za vyuˇzit´ı Dopplerova efektu. Mechanick´e vlnˇen´ı o vysok´e frekvenci (ultrazvuk) je vys´ıl´ano vys´ılaˇcem sondy, vlnˇen´ı se odraz´ı od mˇeˇren´eho objektu a vrac´ı se zp´atky k pˇrij´ımaˇci sondy. Frekvence pˇrijat´eho vlnˇen´ı je vˇsak v d˚ usledku pohybu mˇeˇren´eho objektu odliˇsn´a od vyslan´eho zvuku. Z tohoto rozd´ılu se d´a vypoˇc´ıtat rychlost a smˇer proud´ıc´ıch ˇc´astic v mˇeˇren´e oblasti. Invazivn´ı zp˚ usob mˇeˇren´ı pr˚ utoku krve je moˇzn´e prov´est pomoc´ı indukˇcn´ıch pr˚ utokomˇer˚ u. Ty jsou zaloˇzeny na principu Faradayova z´akona elektromagnetick´e indukci, podle nˇehoˇz (zjednoduˇsenˇe) se v d˚ usledku pohybuj´ıc´ıho se vodiˇce (krev) v homogenn´ım magnetick´em poli (vytv´aˇren´e sondou) bude indukovat elektrick´e napˇet´ı na mˇeˇr´ıc´ıch c´ıvk´ach sondy. Tyto pr˚ utokomˇery maj´ı srpkovit´ y tvar, kter´ y jim umoˇzn ˇuje nasazen´ı na danou c´evu, ve kter´e si pˇrejeme mˇeˇrit pr˚ utok.
Obr´azek 3.7: Typick´ y rychlostn´ı profil toku krve aortou Vzhledem k tomu, ˇze ultrazvukov´e i indukˇcn´ı senzory pr˚ utoku jsou pro naˇse u ´ˇcely velmi drah´e, nebudeme v naˇsem modelu mˇeˇrit okamˇzit´ y pr˚ utok, ale jen pr˚ utok pr˚ umˇern´ y.
Pr˚ utok krve a jeho mˇeˇren´ı
- 28 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
3.6.3
Dalˇ s´ı v´ yznamn´ e mˇ eˇ ren´ e parametry
R˚ uzn´ ymi metodami lze mˇeˇrit mnoho dalˇs´ıch parametr˚ u kardiovaskul´arn´ıho syst´emu. V t´eto podkapitole jsou uvedeny jen nˇekter´e z´akladn´ı. Jak jiˇz bylo uvedeno, jednou ze z´akladn´ıch mˇeˇren´ ych veliˇcin je krevn´ı tlak. Je moˇzno mˇeˇrit systolick´ y (SP) a diastolick´ y (DP) tlak, ze kter´ ych je moˇzn´e odvodit dalˇs´ı parametry jako pulsn´ı tlak (PP) a stˇredn´ı arteri´aln´ı tlak (MAP). Pulsn´ı tlak pˇredstavuje rozd´ıl mezi tlakem systolick´ ym a diastolick´ ym: P P = SP − DP (3.9) Stˇredn´ı arteri´aln´ı tlak lze vypoˇc´ıtat nˇekolika zp˚ usoby, napˇr´ıklad jako: M AP =
SP + 2DP 3
M AP = (CO · SV R) + CV P
(3.10)
(3.11)
kde CO je srdeˇcn´ı v´ydej (cardiac output), SV R je odpor syst´emov´ych c´ev (systemic vascular resistance), CV P je centr´aln´ı ˇziln´ı tlak (central venous pressure) Srdeˇcn´ı v´ ydej (CO) pˇredstavuje celkov´e mnoˇzstv´ı krve, kter´e srdce vypumpuje za jednu minutu. Z´avis´ı na tepov´em objemu (SV) a tepov´e frekvenci (HR): CO = SV · HR
(3.12)
kde CO je srdeˇcn´ı v´ydej [l/min], SV je tepov´y objem [l] a HR je tepov´a frekvence [min−1 ] V klidov´em stavu se jedn´a pˇribliˇznˇe o 5 litr˚ u za minutu, tedy cel´ y objem krve v lidsk´em tˇele. V pˇr´ıpadˇe potˇreby ale m˚ uˇze v´ yraznˇe vzr˚ ust (aˇz pˇetkr´at), a to pˇredevˇs´ım kv˚ uli zv´ yˇsen´ı tepov´e frekvence, ale i tepov´eho objemu. Srdeˇcn´ı v´ ydej se bude samozˇrejmˇe liˇsit v z´avislost na rozmˇerech lidsk´eho tˇela. Normalizovan´ y parametr se naz´ yv´a srdeˇcn´ı index (CI), a jedn´a se o srdeˇcn´ı v´ ydej vztaˇzen´ yk rozmˇer˚ um ˇclovˇeka (udan´em pomoc´ı celkov´e plochy tˇela BSA): CI =
CO BSA
(3.13)
kde CI je srdeˇcn´ı index [l/min−1 ·m2 ], CO je srdeˇcn´ı v´ydej [l/min] a BSA je celkov´ a 2 plocha tˇela[m ] Celkovou plochu tˇela lze vypoˇc´ıtat empiricky napˇr´ıklad pomoc´ı Mostellerovy rovnice: r h·m BSA = ( ) (3.14) 3600 kde BSA je celkov´a plocha tˇela [m2 ], h je v´yˇska ˇclovˇeka [m], m je tˇelesn´a hmotnost [kg] Dalˇs´ım zaj´ımav´ ym parametrem je celkov´ y odpor periferi´ı (TPR). Ten lze zmˇeˇrit pomoc´ı modifikovan´eho Ohmova z´akona: TPR =
M AP − CV P CO
(3.15)
kde T P R je celkov´y odpor periferi´ı [mmHg · min / l], M AP − CV P je rozd´ıl stˇredn´ıho arteri´aln´ıho tlaku a centraln´ıho ˇziln´ıho tlaku [mmHg], CO je srdeˇcn´ı v´ydej [l/min] Dalˇs´ı v´ yznamn´e mˇeˇren´e parametry
- 29 -
Kapitola 4 Teorie modelov´ an´ı kardiovaskul´ arn´ıho syst´ emu Tato kapitola slouˇz´ı jako struˇcn´ yu ´vod do problematiky modelov´an´ı, pˇredevˇs´ım pak kardiovaskul´arn´ıho syst´emu. Naleznete zde struˇcn´ yu ´vod do teorie modelov´an´ı, popis probl´em˚ u spojen´ y s modelov´an´ım biologick´ ych soustav a kr´atk´e shrnut´ı existuj´ıc´ıch model˚ u.
4.1
´ Uvod do teorie modelov´ an´ı
Na zaˇca´tek je vhodn´e definovat nˇekolik pojm˚ u a obecn´ y postup modelov´an´ı. Z´akladem pro modelov´an´ı je re´aln´ y objekt (syst´em). Tento re´aln´ y objekt (at’ uˇz je to org´anov´a soustava, pˇr´ıstroj, poˇcas´ı) b´ yv´a zpravidla sloˇzit´ y, sest´av´a se z mnoha prvk˚ u, kter´e se vz´ajemnˇe ovlivˇ nuj´ı, chov´an´ı objektu b´ yv´a do jist´e m´ıry z´avisl´e na okol´ı a na ˇcase. Vˇsechny vztahy, kter´e ovlivˇ nuj´ı chov´an´ı tohoto syst´emu, nemus´ı b´ yt zn´amy, nebo mohou b´ yt velice komplexn´ı. Modelov´an´ı si klade za c´ıl vytvoˇrit jednoduˇsˇs´ı syst´em, na kter´em lze studovat poˇzadovan´e vlastnosti a z jeho chov´an´ı vyvozovat z´avˇery o p˚ uvodn´ım syst´emu - model1 . Model je syst´em, skl´adaj´ıc´ı se ze zn´am´ ych prvk˚ u a vazeb mezi nimi, kter´ y se sv´ ym chov´an´ım bl´ıˇz´ı dan´emu re´aln´emu syst´emu, respektive pro n´as zaj´ımav´e ˇc´asti chov´an´ı re´aln´eho syst´emu. Z principu vˇeci vypl´ yv´a, ˇze model nebude nikdy pˇresnou reprezentac´ı p˚ uvodn´ıho syst´emu. Kvalita modelu se tedy neposuzuje podle toho, jak vˇernˇe kop´ıruje p˚ uvodn´ı syst´em, ale podle toho jak je pro n´as uˇziteˇcn´ y a jak´e n´am pˇrin´aˇs´ı nov´e poznatky o p˚ uvodn´ım syst´emu. M´ıru toho, jak je pro n´as model uˇziteˇcn´ y, naz´ yv´ame platnost´ı modelu. Podle [8] Rozliˇsujeme tˇri stupnˇe platnosti modelu: • Replikaˇcn´ı platnost. Model generuje stejn´a data jako re´aln´ y syst´em dosud. • Predikˇcn´ı platnost. Model generuje spr´avn´a data dˇr´ıve, neˇz re´aln´ y syst´em. • Struktur´aln´ı platnost. Model generuje stejn´a data jako re´aln´ y syst´em, nav´ıc ale tato data generuje stejn´ ym zp˚ usobem, jako re´aln´ y syst´em. Modelovat lze i syst´em, u kter´eho nezn´ame vnitˇrn´ı strukturu. V tomto pˇr´ıpadˇe bereme modelovan´ y syst´em jako blackbox, kter´ y pˇrij´ım´a vstupn´ı data a vrac´ı data v´ ystupn´ı. Pakliˇze data generovan´a modelem se bl´ıˇz´ı dat˚ um generovan´ ym re´aln´ ym syst´emem, lze vytvoˇrit hypot´ezu o vnitˇrn´ı struktuˇre p˚ uvodn´ıho syst´emu. 1
Slovo model m˚ uˇze v ˇceˇstinˇe oznaˇcovat i repliku re´aln´eho objektu, kter´ y je se sv´ ym vzorem takˇrka identick´ y a slouˇz´ı napˇr´ıklad pro destruktivn´ı zkouˇsky. Tento v´ yznam vˇsak v tomto textu nen´ı uvaˇzov´ an
30
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
4.1.1
Popis modelu
Model lze popsat dvˇema zp˚ usoby, a to form´alnˇe ˇci neform´alnˇe. Neform´aln´ı popis m˚ uˇze b´ yt v pˇrirozen´em jazyce, m˚ uˇze to b´ yt blokov´e sch´ema, ˇci kombinace oboj´ıho. Neform´aln´ı popis obsahuje v´ yˇcet prvk˚ u, ze kter´ ych se skl´ad´a model, z popisu promˇenn´ ych veliˇcin, kter´e se k jednotliv´ ym prvk˚ um v´aˇzou, z definice parametr˚ u (obvykle konstantn´ı charakteristiky prvk˚ u), z definice vazeb mezi prvky, tedy pravidel, podle kter´ ych se prvky navz´ajem ovlivˇ nuj´ı, a z pˇredpokl´adan´ ych poˇc´ateˇcn´ıch podm´ınek. Pro v´ ybˇer prvk˚ u, promˇenn´ ych, parametr˚ u a poˇca´teˇcn´ıch podm´ınek neexistuj´ı ˇza´dn´a jednotn´a pravidla. Obecnˇe se d´a pˇr´ıstup k jejich v´ ybˇeru rozdˇelit na deduktivn´ı a induktivn´ı. Deduktivn´ı pˇr´ıstup je zaloˇzen na znalosti obecnˇe platn´ ych z´akonitost´ı dan´ ych jev˚ u, uplatn´ı se tedy u modelov´an´ı jev˚ u, jejichˇz podstata je zn´am´a a u kter´ ych v´ıme, ˇze se chovaj´ı podle dan´ ych pravidel. Induktivn´ı pˇr´ıstup se potom uplatˇ nuje u jev˚ u, kter´e nejsou pops´any zn´am´ ymi z´akony (tento pˇr´ıpad je u biologick´ ych syst´emu nejˇcastˇejˇs´ı). Form´aln´ı popis vyuˇz´ıv´a r˚ uzn´ ych matematick´ ych prostˇredk˚ u, je vˇzdy jednoznaˇcn´ y, konzistentn´ı a u ´pln´ y.
4.1.2
Obecn´ y postupy modelovan´ı
Obecnˇe postup modelov´an´ı podle [8] z´avis´ı na tom, zda modelujeme jiˇz teoreticky popsan´ y syst´em ˇci vych´az´ıme z namˇeˇren´ ych dat, na tom, zda struktura modelu odpov´ıd´a struktuˇre re´aln´eho syst´emu a na rozsahu modelu. Podle zvolen´eho postupu potom dˇel´ıme modely na: • Apriorn´ı a aposteriorn´ı. Apriorn´ı (teoretick´ y) model vych´az´ı ze znalosti teoretick´eho fungov´an´ı syst´emu a je sestrojen na z´akladˇe pˇredpokl´adan´eho zp˚ usobu generov´an´ı v´ ysledn´ ych dat. Aposteriorn´ı (empirick´ y) model vych´az´ı z namˇeˇren´ ych dat pˇri zn´am´em vstupu, tedy bez znalosti vnitˇrn´ı struktury syst´emu. Aposteriorn´ı model mus´ı generovat data co nejvˇernˇejˇs´ı dat˚ um namˇeˇren´ ym. • Popisn´e a mechanistick´e. Popisn´ y model je sloˇzen z prvk˚ u, jejichˇz parametry a promˇenn´e nemaj´ı vztah k re´aln´emu syst´emu, nicm´enˇe model generuje spr´avn´a data. Mechanistick´ y model je sloˇzen z prvk˚ u, jejichˇz parametry a promˇenn´e odpov´ıdaj´ı re´aln´emu syst´emu, tud´ıˇz v´ ystupn´ı data modelu jsou generov´ana obdobn´ ym zp˚ usobem, jako v re´aln´em syst´emu. • Redukovan´ y a rozs´ahl´ y. Redukovan´e modely popisuj´ı co nej´ uspornˇeji re´aln´ y syst´em. Rozs´ahl´e modely zahrnuj´ı i dalˇs´ı charakteristiky re´aln´eho syst´emu. Rozs´ahl´e modely mohou pˇrispˇet k vytv´aˇren´ı nov´ ych hypot´ez o chov´an´ı re´aln´eho syst´emu. Nicm´enˇe model mus´ı b´ yt maxim´alnˇe tak sloˇzit´ y, aby bylo moˇzn´e interpretovat a pochopit jeho v´ ysledky, popsat jeho vztah k re´aln´emu syst´emu, a aby bylo v˚ ubec moˇzn´e (nebo ekonomicky v´ yhodn´e) dan´ y model sestrojit.
4.1.3
Kompartmentov´ a anal´ yza
Aˇckoli neexistuje jednotn´ y n´avod jak vytvoˇrit (at’ uˇz form´aln´ı nebo neform´aln´ı) model, jedna metoda pˇresto dosti zautomatizovala postup tvoˇren´ı modelu. Jedn´a se o kompartmentovou anal´ yzu. Tato metoda spoˇc´ıv´a v rozkladu modelovan´eho syst´emu na kompartPopis modelu
- 31 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
menty, urˇcen´ı vztah˚ u mezi nimi (ˇcasto oznaˇcovan´e jako kan´aly) a definov´an´ı poˇca´teˇcn´ıch podm´ınek. Kompartment je urˇcit´a (koncepˇcn´ı) oblast modelovan´eho syst´emu, kterou je moˇzno popsat urˇcit´ ymi fyzik´aln´ımi veliˇcinami (napˇr´ıklad objemem, koncentrac´ı l´atek, tlakem....). Pˇredpokl´ad´a se, ˇze veˇsker´e fyzik´aln´ı vlastnosti jsou uvnitˇr kompartmentu uniformn´ı. Napˇr´ıklad pokud se v kompartmentu nach´az´ı smˇes l´atek, je koncentrace tˇechto l´atek vˇsude stejn´a. Pokud se koncentrace l´atek v kompartmentu mˇen´ı, mˇen´ı se stejnˇe v cel´em objemu ve stejn´em ˇcase. Vstup kompartmentu je pomysln´ y bod, skrz kter´ y m˚ uˇze sledovan´a veliˇcina do kompartmentu vstupovat (at’ uˇz z okol´ı, jin´ ych kompartment˚ u, nebo synt´ezou uvnitˇr kompartmentu). Obdobnˇe v´ ystup pˇredstavuje m´ısto, kudy sledovan´a veliˇcina kompartment opouˇst´ı (do okol´ı, jin´eho kompartmentu, nebo pˇremˇenou v jin´e l´atky). I jedn´ım kompartmentem se vstupy a v´ ystupy do okol´ı se daj´ı modelovat urˇcit´e syst´emy, ale ˇcastˇeji se konstruuj´ı multikompartmentov´e modely. Tyto modely jsou propojeny kan´aly, kter´ ymi se kompartmenty navz´ajem ovlivˇ nuj´ı. Kan´aly slouˇz´ı pouze pro spojen´ı kompartment˚ u, nemaj´ı ˇz´adn´ y objem ani jin´e vlastnosti. V´ ysledkem kompartmentov´e anal´ yzy je kompartmentov´ y model, jenˇze se d´a zcela popsat obyˇcejn´ ymi diferenci´aln´ımi rovnicemi (ODE). Pro ˇreˇsen´ı tˇechto rovnic je jeˇstˇe nutn´e definovat poˇc´ateˇcn´ı podm´ınky. Pokud pˇrijmeme pˇredpoklad, ˇze model nen´ı spojit´ y, ale diskr´etn´ı (tzn. ˇze existuje jen v pˇresnˇe definovan´ ych ˇcasov´ ych bodech), pˇrech´az´ı diferenci´aln´ı rovnice v rovnice diferenˇcn´ı.
4.2
Specifika biologick´ ych syst´ em˚ u
Biologick´e syst´emy se vyznaˇcuj´ı nˇekolika specifick´ ymi vlastnostmi, kter´e je potˇreba pˇri modelov´an´ı m´ıt na pamˇeti. Pˇredevˇs´ım se ˇcasto jedn´a o rozs´ahl´e syst´emy se sloˇzitou strukturou, kter´e obsahuj´ı mnoho promˇenn´ ych, pˇriˇcemˇz vztahy mezi jednotliv´ ym promˇenn´ ymi ani nemus´ı b´ yt zn´amy. Z d˚ uvodu mnoha zpˇetn´ ych vazeb a regulaˇcn´ıch mechanism˚ u nen´ı ˇcasto z´avislost mezi jednotliv´ ymi promˇenn´ ymi line´arn´ı. Biologick´e syst´emy se tak´e vyznaˇcuj´ı velkou interindividu´aln´ı a intrainduvidu´aln´ı variabilitou. Interindividualn´ı variabilita znamen´a, ˇze stejn´ y syst´em se mezi jednotliv´ ymi jedinci populace bude chovat jinak (napˇr´ıklad krevn´ı tlak dvou r˚ uzn´ ych lid´ı bude jin´ y). Intraindividualn´ı variabilita pak znamen´a, ˇze stejn´ y syst´em se bude u stejn´eho jedince chovat jinak v z´avislosti na ˇcase (napˇr´ıklad krevn´ı tlak jedince bude jin´ y r´ano a veˇcer).
4.3
Modely kardiovaskul´ arn´ı soustavy
Jelikoˇz je kardiovaskul´arn´ı syst´em jedn´ım z prim´arn´ıch syst´em˚ u lidsk´eho tˇela, je mu vˇseobecnˇe vˇenov´ana velk´a pozornost, a to i na poli modelov´an´ı. Za zm´ınku stoj´ı pˇredevˇs´ım tzv. Mock Circulatory Loop - MCL. Tyto modely (at’ uˇz fyzick´e nebo poˇc´ıtaˇcov´e, napˇr [10] nebo [11]) se vesmˇes soustˇred´ı na modelov´an´ı hlavn´ıch hemodynamick´ ych parametr˚ u jako je (minutov´ y) srdeˇcn´ı v´ ydej, systolick´ y a diastolick´ y arteri´aln´ı tlak, centr´aln´ı ven´ozn´ı tlak a fyziologick´ y pr˚ utok krve hlavn´ımi c´evami. Jejich hlavn´ım c´ılem b´ yv´a demonstrovat vz´ajemn´e z´avislosti hemodynamick´ ych parametr˚ u nebo testovat a kalibrovat pˇr´ıstroje ovlivˇ nuj´ıc´ı nebo nahrazuj´ıc´ı pr´aci srdce, jako jsou umˇel´e chlopnˇe, komorov´e asistenˇcn´ı pˇr´ıstroje (VAD) ˇci tot´aln´ı srdeˇcn´ı n´ahrady (TAH). Specifika biologick´ ych syst´em˚ u
- 32 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
Pro tyto modely je typick´e, ˇze sdruˇzuj´ı jednotliv´e c´evy do funkˇcn´ıch blok˚ u, ˇcasto je pr´avˇe cel´a mikrocirkulace nahrazena jedn´ım blokem s danou resistenc´ı a complienc´ı. Tento pˇr´ıstup vˇsak nen´ı pro n´aˇs model pˇrijateln´ y, nebot’ n´as zaj´ım´a distribuce a disperze krve ve velmi mal´ ych c´ev´ach, tud´ıˇz v naˇsem pˇr´ıpadˇe bude c´evn´ı ˇreˇciˇstˇe modelov´ano rozvˇetven´ ym stromem sloˇzen´ ym z hadiˇcek s postupnˇe se sniˇzuj´ıc´ımi pr˚ umˇery.
Modely kardiovaskul´arn´ı soustavy
- 33 -
Kapitola 5 V´ ysledn´ y model kardiovaskul´ arn´ıho syst´ emu Sestaven´ y model kardiovaskul´arn´ı soustavy lze rozdˇelit na nˇekolik ˇc´ast´ı, konkr´etnˇe na ˇca´st mechanickou, ˇc´ast elektronickou a ˇc´ast softwarovou. Mechanick´a ˇc´ast se skl´ad´a z peristaltick´eho ˇcerpadla, poh´anˇej´ıc´ıho tekutinu v ˇreˇciˇsti, z kompenzaˇcn´ı n´adoby s elektromagnetick´ ym ventilem, vyrovn´avaj´ıc´ı tok a simuluj´ıc´ı ˇcinnost srdce, z c´evn´ıho ˇreˇciˇstˇe a z rezervo´aru pro koluj´ıc´ı tekutinu. Elektronick´a ˇca´st zahrnuje pouˇzit´e tlakov´e senzory, zdroj nap´ajen´ı a ˇr´ıd´ıc´ı a mˇeˇr´ıc´ı obvody. Software potom obstar´av´a samotn´e ˇr´ızen´ı modelu a sbˇer dat. Celkov´e sch´ema je zobrazeno na n´asleduj´ıc´ım obr´azku:
Obr´azek 5.1: Sch´ema vytvoˇren´eho modelu: 1 - Z´asobn´ık tekutiny; 2 - Pumpa; 3 - Kompenzaˇcn´ı n´adoba; 4 - Ventil; 5 - C´evn´ı ˇreˇciˇstˇe; P - tlakomˇer
34
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
5.1
ˇ Cerpadlo
Pro poh´anˇen´ı kapaliny je v modelu pouˇzito peristaltick´e ˇcerpadlo. Tento typ ˇcerpadla se skl´ad´a z rotoru osazen´eho pˇr´ıtlaˇcn´ ymi v´alci, elastick´e trubice a kruhovit´eho pevn´eho bloku. Elastick´a trubice je vsazena mezi v´aleˇcky rotoru a pevnou stˇenu. Ot´aˇcen´ım rotoru doch´az´ı k okluzi trubice v m´ıstˇe pˇr´ıtlaku, t´ım je kapalina pˇred v´alcem hn´ana ve smˇeru rotace a d´ale z ˇcerpadla ven. Po uvolnˇen´ı pˇr´ıtlaku se trubice vr´at´ı do p˚ uvodn´ıho stavu, ˇc´ımˇz se v trubici vyvol´a podtlak a dojde k nas´at´ı kapaliny z rezervo´aru. Peristaltick´e ˇcerpadlo je samonas´avac´ı a poskytuje v´ ybornou kontrolu ˇcerpan´eho ob’ jemu, nebot poˇcet ot´aˇcek pˇresnˇe odpov´ıd´a pˇreˇcerpan´emu objemu (viz obr´azek 5.2). Nav´ıc nedoch´az´ı ke styku ˇcerpan´e kapaliny a sloˇzitˇejˇs´ıch ˇca´st´ı ˇcerpadla, protoˇze kapalina je zcela uzavˇrena ve snadno vymˇeniteln´e a relativnˇe levn´e trubici. To sniˇzuje n´aklady spojen´e s u ´drˇzbou ˇcerpadla. Pro budouc´ı u ´ˇcely naˇseho modelu to pak znamen´a, ˇze jako analogie krve je moˇzn´e pouˇz´ıt i jin´e tekutiny, neˇz jen ˇcist´e vody. V naˇsem modelu je pouˇzito ˇcerpadlo Verderflex Smart S40 se ˇsesti pˇr´ıtlaˇcn´ ymi v´alci. Toto ˇcerpadlo je vybaveno rozhran´ım pro ovl´ad´an´ı extern´ımi sign´aly. Pro naˇse u ´ˇcely je vyuˇzit logick´ y sign´al pro okamˇzit´e spuˇstˇen´ı/zastaven´ı ˇcerpadla a napˇet’ov´ y sign´al 0-10V, ˇ kter´ y ˇr´ıd´ı poˇcet ot´aˇcek za minutu (tud´ıˇz i pˇresnˇe stanovuje ˇcerpan´ y objem). Cerpadlo je nap´ajeno z rozvodn´e s´ıtˇe. Trubice ˇcerpadla je z materi´alu Verderprene, kter´ y poskytuje vysokou mechanickou a chemickou odolnost, s tlouˇst’kou stˇeny 4,8 mm a celkov´ ym pr˚ umˇerem 19,0 mm. Na n´ıˇze uveden´em obr´azku je uvedena z´avislost dod´avan´eho pr˚ utoku ˇcerpadla na poˇctu ot´aˇcek ˇcerpadla:
Obr´azek 5.2: Z´avislost pr˚ utoku ˇcerpadla na poˇctu jeho ot´aˇcek. Pro n´aˇs pˇr´ıpad plat´ı kˇrivka s rozmˇery 19,0x4,8. Obr´azek pˇrevzat z manu´alu ˇcerpadla. ˇ Cerpadlo
- 35 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
Vytvoˇren´ y model m´a m´a slouˇzit ke studiu kontinu´aln´ıho a pulsatiln´ıho proudˇen´ı tekutiny v c´evn´ım ˇreˇciˇsti. K tomu vˇsak nelze pouˇz´ıt samotn´e ˇcerpadlo, nebot’ neposkytuje kontinu´aln´ı tok ani konstantn´ı tlak. Dalˇs´ım probl´emem, kter´ y by hrozil pˇri napojen´ı v´ ystupu ˇcerpadla pˇr´ımo na elektromagnetick´ y ventil, simuluj´ıc´ı srdeˇcn´ı chlopeˇ n, je efekt vodn´ıho kladiva. Tento jev nast´av´a v situaci, kdy se proud´ıc´ı tekutinˇe rychle do cesty postav´ı pˇrek´aˇzka (zavˇren´ y ventil). Proud´ıc´ı voda m´a totiˇz velkou inertanci, snaˇz´ı se z˚ ustat v pohybu. Tato kinetick´a energie se mus´ı nˇekde vyb´ıt, a pokud by se vyb´ıjela o ventil, mohlo by doj´ıt k jeho poˇskozen´ı. Proto je mezi v´ ystup ˇcerpadla a ventil vˇrazena kompenzaˇcn´ı n´adoba.
5.2
Kompenzaˇ cn´ı n´ adoba s elektromagnetick´ ym ventilem
Kompenzaˇcn´ı n´adoba s elektromagnetick´ ym ventilem pln´ı troj´ı funkci v naˇsem syst´emu. Prvn´ı je zabr´anˇen´ı vzniku efektu vodn´ıho kladiva. Druhou je vyhlazen´ı kmitaj´ıc´ıho v´ ystupu ˇcerpadla a pˇremˇena tohoto proudˇen´ı na proudˇen´ı pulsatiln´ı nebo kontinu´aln´ı. Tˇret´ı je potom role akumul´atoru energie dod´avan´e ˇcerpadlem, kter´a se uvoln´ı pˇri otevˇren´ı ventilu (systola) a hromad´ı pˇri zavˇren´ı ventilu (diastola). Vˇsechny tyto role jsou zp˚ usoben´e jevem zn´am´ ym jak windkessel efekt (z nˇemˇciny, doslova efekt vˇetrn´eho kotle). Tento jev nast´av´a v situaci, kdy je do uzavˇren´e n´adoby naplnˇen´e stlaˇciteln´ ym plynem pulsatilnˇe vh´anˇena nestlaˇciteln´a tekutina. Plyn v sobˇe akumuluje energii pulsac´ı, kterou v opaˇcn´e f´azi uvolˇ nuje. Z n´adoby potom tekutina proud´ı kontinu´alnˇe. Podobn´ y efekt (i kdyˇz v menˇs´ı m´ıˇre) nast´av´a i v tepn´ach kardiovaskul´arn´ıho syst´emu, konkr´etnˇe v pruˇzn´ıkov´ ych tepn´ach, kde energie puls˚ u nen´ı akumulov´ana plynem, ale elastickou stˇenou tepny, ˇc´ımˇz doch´az´ı k postupn´e pˇremˇenˇe pulsatiln´ıho proudˇen´ı v arteri´ıch a arteriol´ach na kontinu´aln´ı proudˇen´ı v kapil´ar´ach.
5.2.1
Kompenzaˇ cn´ı n´ adoba
Jako samotn´a kompenzaˇcn´ı n´adoba slouˇz´ı uzavˇren´ y plastov´ y kanystr se dvˇema kohouty pro vtok a v´ ytok kapaliny. Tento kanystr je z ˇca´sti naplnˇen vodou, z ˇca´sti vzduchem. Kanystr je oznaˇcen certifik´atem UN 3H1/Y/120 , kter´ y zaruˇcuje (mimo jin´e) to, ˇze kanystr byl testov´an na vnitˇrn´ı pˇretlak o hodnˇe 120 kPa (pˇribliˇznˇe 900 torr). Vzhledem k tomu, ˇze tlak tekutin v n´adobˇe dosahuje maxim´alnˇe 33 kPa (pˇribliˇznˇe 250 torr), lze pˇredpokl´adat, ˇze nehroz´ı exploze n´adoby. Pro monitorov´an´ı tlaku v n´adobˇe je za jej´ı v´ ystup zaˇrazen tlakov´ y senzor.
5.2.2
Elektromagnetick´ y ventil
Elektromagnetick´ y ventil je z´asadn´ı pro vytvoˇren´ı pulsatiln´ıho toku v naˇsem modelu. Pouˇzit je ventil s n´asleduj´ıc´ımi parametry:
Kompenzaˇcn´ı n´adoba s elektromagnetick´ ym ventilem
- 36 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı Nap´ajen´ı Pr˚ umˇer u ´st´ı Pracovn´ı tlak Dalˇs´ı vlastnosti
12V DC 1/2” 0-7500 mmHg bez napˇet´ı zavˇren´ y, bez minim´aln´ıho tlaku
Tabulka 5.1: Parametry elektromagnetick´eho ventilu Pro naˇse u ´ˇcely je d˚ uleˇzit´e, ˇze ventil nepotˇrebuje pro otevˇren´ı ˇz´adn´ y minim´aln´ı tlak. Doba otevˇren´ı/zavˇren´ı ventilu se pohybuje v ˇr´adu des´ıtek milisekund, coˇz je tak´e dostaˇcuj´ıc´ı. Dalˇs´ı v´ yhodou je, ˇze bez napˇet´ı je ventil zavˇren´ y. V pˇr´ıpadˇe poruchy se tedy nem˚ uˇze st´at, ˇze ventil nep˚ ujde zavˇr´ıt. ˇ Cerpadlo, kompenzaˇcn´ı n´adoba a elektromagnetick´ y ventil simuluj´ı dohromady funkci srdce. V´ ystupem t´eto ˇca´sti je kontinu´aln´ı nebo pulsatiln´ı tok, jehoˇz tlakov´a sloˇzka vypad´a takto:
(a) Arteri´ aln´ı tlak pro kontinu´ aln´ı tok
(b) Ven´ozn´ı tlak pro kontinu´aln´ı tok
(c) Arteri´ aln´ı tlak pro pulsatiln´ı tok
(d) Ven´ozn´ı tlak pro pulsatiln´ı tok
Z poˇca´tku v´ yvoje modelu se st´ale vyskytoval probl´em, kter´ y neumoˇzn ˇoval generovat zcela vyhlazen´ y tok. Tento probl´em byl zp˚ usoben samotn´ ym ˇcerpadlem, kter´e se pˇri pr´aci (pˇredevˇs´ım pˇri vysok´ ych ot´aˇck´ach) dosti chvˇeje. Toto chvˇen´ı se pˇren´aˇselo na kompenzaˇcn´ı n´adobu a n´aslednˇe na cel´ y model. Tento jev byl odstranˇen v´ ymˇenou trubice spojuj´ıc´ı ˇcerpadlo a kompenzaˇcn´ı n´adobu, respektive v´ yrazn´ ym prodlouˇzen´ım jej´ı d´elky, a to z 50 cm na 250 cm.
Elektromagnetick´ y ventil
- 37 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
5.3
C´ evn´ı ˇ reˇ ciˇ stˇ e
K modelov´an´ı c´evn´ıho ˇreˇciˇstˇe jsou pouˇzity silikonov´e hadiˇcky r˚ uzn´ ych pr˚ umˇer˚ u, plastov´e spojky, spojky vyroben´e ze speci´aln´ıho teplem tvrditeln´eho polymeru (hmota FIMO) a odporov´e prvky vlastn´ı v´ yroby. Silikonov´e hadiˇcky jsou nˇekolikan´asobnˇe draˇzˇs´ı neˇz hadiˇcky z PVC obdobn´ ych rozmˇer˚ u, ale maj´ı v´ yraznˇe vˇetˇs´ı elasticitu ˇc´ımˇz se v´ıce bl´ıˇz´ı fyziologick´ ym ˇ c´ev´am. Reˇciˇstˇe je vˇetveno do p´at´e u ´rovnˇe, ˇsestou u ´roveˇ n pak pˇredstavuj´ı odporov´e prvky. Pr˚ umˇer hadiˇcky na nejniˇzˇs´ı u ´rovni je 2 mm, nedost´av´ame se tedy aˇz na u ´roveˇ n kapil´ar. Nicm´enˇe proudˇen´ı v kapil´ar´ach je pˇrev´aˇznˇe kontinu´aln´ı, tud´ıˇz pro n´as nem´a aˇz takov´ y smysl. D˚ uleˇzitˇejˇs´ı je pro n´as odpor kapil´ar, kter´ y je ˇc´asteˇcnˇe simulov´an speci´aln´ım odporov´ ym prvkem. Sch´ema vytvoˇren´eho c´evn´ıho ˇreˇciˇstˇe vypad´a takto:
Obr´azek 5.4: Sch´ema vˇetven´ı ˇreˇciˇstˇe, neodpov´ıd´a pˇresn´e geometrii v´ ysledn´eho modelu. ˇ Cervenˇ e - arteri´aln´ı ˇc´ast; fialovˇe - kapil´arn´ı odpor; modˇre - ven´ozn´ı ˇc´ast; b´eˇzovˇe - mˇeˇr´ıc´ı okruhy ˇ Celkov´ y objem c´evn´ıho ˇreˇciˇstˇe je pˇribliˇznˇe 300 ml. Cerpadlo je nastaveno tak, aby tento objem pˇreˇcerpalo za jednu minutu. Vych´az´ıme pˇritom z re´aln´eho kardiovaskul´arn´ıho syst´emu, kdy srdce v klidu pˇreˇcerp´a (pˇribliˇznˇe) 5 l/min, coˇz je (pˇribliˇznˇe) objem krve v tˇele.
5.3.1
Pouˇ zit´ e hadiˇ cky
Pro arteri´aln´ı a ven´ozn´ı ˇc´ast ˇreˇciˇstˇe jsou pouˇzity stejn´e hadiˇcky stejn´ ych vlastnost´ı. Pro jednotliv´e u ´rovnˇe vˇetven´ı (viz obr´azek 5.4) jsou zvoleny n´asleduj´ıc´ı d´elky a pr˚ umˇery hadiˇcek:
C´evn´ı ˇreˇciˇstˇe
- 38 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
´ Uroveˇ n L1 L2 L3 L4 L5
Vnitˇrn´ı pr˚ umˇer/ vnˇejˇs´ı [mm] 10/15 8/12 6/9 4/6 2/4
Pr˚ umˇern´a d´elka [cm] 90 50 20 10 5
Poˇcet hadiˇcek v u ´rovni 1 2 4 7 19
Tabulka 5.2: Parametry pouˇzit´ ych hadiˇcek pro jednotliv´e u ´rovnˇe V modelu jsou tak´e zaˇrazeny dva mˇeˇr´ıc´ı okruhy (pro mˇeˇren´ı rychlosti pohybu vloˇzen´e kuliˇcky). Tyto okruhy jsou na obr´azku 5.4 zobrazeny b´eˇzovou barvou. Parametry tˇechto hadiˇcek jsou: ´ Uroveˇ n L4 L5
Vnitˇrn´ı pr˚ umˇer/ vnˇejˇs´ı pr˚ umˇer [mm] 4/6 2/4
D´elka [cm] 150 60
Tabulka 5.3: Parametry hadiˇcek pouˇzit´ ych v mˇeˇr´ıc´ıch okruz´ıch
5.3.2
Pouˇ zit´ e spojky
Ke spojov´an´ı hadiˇcek r˚ uzn´ ych pr˚ umˇer˚ u jsou pouˇzity standardn´ı pˇr´ım´e kuˇzelovit´e polypropylenov´e redukˇcn´ı spojky. Pro vˇetven´ı syst´emu jsou pouˇzity standardn´ı polypropylenov´e spojky typu y“ a tak´e spojky vlastn´ı v´ yroby typu dvojit´e ypsilon“. Tyto vlastn´ı spojky ” ” jsou vyrobeny z polymerov´e teplem tvrditeln´e hmoty, a byly navrˇzeny tak, aby pˇr´ıtok z jedn´e hadiˇcky o vnitˇrn´ım pr˚ umˇeru 4 mm byl distribuov´an do ˇctyˇr vˇetv´ı o pr˚ umˇeru 2 mm.
Obr´azek 5.5: Pˇr´ıklad spojky typu ”dvojit´e ypsilon”
Pouˇzit´e spojky
- 39 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
5.3.3
Odporov´ y prvek
Odporov´ y prvek do jist´e m´ıry simuluje odpor kapil´ar. Je vyroben ze silikonov´e trubiˇcky, ve kter´e je um´ıstˇen mal´ y kousek obyˇcejn´e houbiˇcky na myt´ı n´adob´ı. Se zbytkem ˇreˇciˇstˇe je tento prvek spojen vlastn´ımi spojkami. Zaˇrazen´ım tˇechto odpor˚ u do modelu se povedlo zv´ yˇsit diastolick´ y tlak v pulsatiln´ım m´odu o 15 mmHg a tak´e se zabr´anilo zpˇetn´emu proudˇen´ı. Samotn´ y odporov´ y prvek je zobrazen na n´asleduj´ıc´ı fotografii:
Obr´azek 5.6: Odporov´ y prvek
5.4
Senzory
V naˇsem modelu jsou vyuˇzity jen senzory tlaku. Zcela jistˇe by bylo vhodn´e pouˇz´ıt i senzory pr˚ utoku, probl´em je vˇsak v jejich cenˇe. Turb´ınkov´e a vrtulkov´e pr˚ utokomˇery jsou relativnˇe levn´e (v ˇra´du tis´ıc˚ u korun), ale pro naˇse mˇeˇren´ı skoro nepouˇziteln´e, nebot’ turb´ınka m´a relativnˇe velkou setrvaˇcnost. Ide´aln´ı by bylo pouˇz´ıt ultrazvukov´ y sn´ımaˇc pr˚ utoku, ale u nˇej je probl´em vysok´a poˇrizovac´ı cena (v ˇr´adu desetitis´ıc˚ u). Proto nakonec od pˇr´ım´eho mˇeˇren´ı pr˚ utoku bylo upuˇstˇeno a pˇreˇslo se jen k nepˇr´ım´ ym metod´am mˇeˇren´ı pr˚ utoku.
5.4.1
Senzory tlaku
K mˇeˇren´ı tlaku se vyuˇz´ıvaj´ı dva typy senzor˚ u. Prvn´ım je sn´ımaˇc tlaku s keramickou mˇeˇr´ıc´ı membr´anou od firmy Cressto (k´od v´ yrobku TM G617A3F), kter´ y mˇeˇr´ı relativn´ı pˇretlak vzhledem k atmosf´eˇre. Tento senzor je nap´ajen stejnosmˇern´ ym napˇet´ım o hodnotˇe 12 V. V´ ystupem senzoru je proudov´ y sign´al v rozsahu 4 - 20 mA. Pˇrevodn´ı funkce senzoru je: PP a =
(2500 · I − 10000) 133, 32
(5.1)
kde PmmHg je tlak [mmHg], I je proud [mA] Senzor se k ˇreˇciˇsti pˇripojuje z´avitem o pr˚ umˇeru 1/2”, je jej tedy moˇzn´e pˇripojit jen k patˇriˇcnˇe velk´ ym hadic´ım. Proto jsou tyto senzory vyuˇzity k mˇeˇren´ı tlaku uvnitˇr kompenzaˇcn´ı n´adoby (senzor um´ıstˇen mezi kompenzaˇcn´ı n´adobu a elektromagnetick´ y ventil) a pro mˇeˇren´ı arteri´aln´ıho tlaku (senzor um´ıstˇen mezi elektromagnetick´ y ventil a art´erii). Nap´ajen´ı i pˇrenos dat prob´ıh´a pomoc´ı jedn´e dvojvodiˇcov´e linky.
Odporov´ y prvek
- 40 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
Obr´azek 5.7: Tlakov´ y senzor od firmy Cressto Druh´ ym typem je potom sn´ımaˇc tlaku od firmy Argon Medical Devices (oznaˇcen´ı v´ yrobku DTX+, TNF-R, ˇrada 682021). Nap´ajen je tak´e stejnosmˇern´ ym napˇet´ım o hodnotˇe 12 V a vrac´ı rovnˇeˇz proudov´ y sign´al. Tento sn´ımaˇc se bˇeˇznˇe pouˇz´ıv´a v medic´ınsk´e praxi. Hlavn´ı v´ yhodou pro n´aˇs model je moˇznost sn´ımaˇc pohodlnˇe pˇripojit i na niˇzˇs´ı vˇetve c´evn´ıho stromu (pomoc´ı hadiˇcky s pr˚ umˇerem 4 mm). V naˇsem syst´emu je vyuˇzit pro monitorov´an´ı tlaku ve ven´ozn´ı ˇc´asti ˇreˇciˇstˇe. Druh´ y senzor je potom moˇzno pˇripojit na libovoln´e m´ısto v ˇreˇciˇsti.
Obr´azek 5.8: Tlakov´ y senzor od firmy Argon Medical Devices
5.4.2
Mˇ eˇ ren´ı pr˚ utoku
Jak jiˇz bylo zm´ınˇeno, nen´ı v modelu implementov´ano pˇr´ım´e mˇeˇren´ı pr˚ utoku. K dispozici byl jeden turb´ınkov´ y pr˚ utokomˇer s pr˚ umˇerem mˇeˇr´ıc´ıho otvoru 1/2”. Pˇri implementaci tohoto senzoru v pˇredchoz´ıch prac´ıch vˇsak bylo zjiˇstˇeno, ˇze j´ım nelze mˇeˇrit okamˇzit´ y pr˚ utok, nebot’ setrvaˇcnost turb´ınky je pro tak mal´e pr˚ utoky pˇr´ıliˇs velk´a. Z toho d˚ uvodu se pˇreˇslo k nepˇr´ım´e metodˇe mˇeˇren´ı pr˚ utoku. Mˇeˇren´ı pr˚ utoku
- 41 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
Celkov´ y pr˚ utok syst´emem lze mˇeˇrit pomoc´ı odmˇern´e n´adoby um´ıstˇen´e na konec c´evn´ıho stromu. Na u ´rovni posledn´ı a pˇredposledn´ı vˇetve c´evn´ıho stromu (´ uroveˇ n L5, respektive L4 podle obr´azku 5.4) se nemˇeˇr´ı pˇr´ımo pr˚ utok, ale pr˚ umˇern´a rychlost proudˇen´ı. Toto mˇeˇren´ı je realizov´ano pomoc´ı sledov´an´ı pohybu kuliˇcky v mˇeˇr´ıc´ı hadiˇcce o zn´am´em objemu. Kuliˇcky jsou vyrobeny ze stejn´eho materi´alu jako vyroben´e spojky. Pr˚ umˇer kuliˇcky je takov´ y, aby kuliˇcka vyplnila cel´ y pr˚ uˇrez hadiˇcky (a nedoch´azelo tak k ovlivnˇen´ı mˇeˇren´ı pˇr´ıpadn´ ymi v´ıry v proudˇen´ı a tak´e aby bylo moˇzn´e zanedbat hustotu kuliˇcky), ale z´aroveˇ n takov´ y, aby se kuliˇcka mohla v hadiˇcce volnˇe pohybovat (a nezvyˇsovala tak odpor kladen´ y proudˇen´ı tekutiny nad u ´nosnou mez). Kuliˇcky jsou zn´azornˇeny na n´asleduj´ıc´ım obr´azku:
Obr´azek 5.9: Kuliˇcky a vzorky hadiˇcek, pouˇzit´e k mˇeˇren´ı pr˚ umˇern´e rychlosti proud´ıc´ı kapaliny
5.5
ˇ ızen´ı modelu a sbˇ R´ er dat
Pro ˇr´ızen´ı modelu a sbˇer mˇeˇren´ ych dat je vyuˇzito pr˚ umyslov´e vstupnˇe v´ ystupn´ı zaˇr´ızen´ı DataLab IO4 s USB rozhran´ım od spoleˇcnosti Moravsk´e pˇr´ıstroje a.s. Cel´ y model je ˇr´ızen´ y pomoc´ı aplikace vyvinut´e v prostˇred´ı Control Web od t´eˇze spoleˇcnosti.
5.5.1
Hardware pro ovl´ ad´ an´ı a sbˇ er dat
DataLab IO4 je pr˚ umyslov´e vstupnˇe v´ ystupn´ı zaˇr´ızen´ y, skl´adaj´ıc´ı se z nˇekolika (voliteln´ ych) modul˚ u. V´ yhodou tohoto syst´emu je mnoˇzstv´ı vstupn´ıch a v´ ystupn´ıch obvod˚ u, velice snadn´e ovl´ad´an´ı, zabudovan´a komunikace s poˇc´ıtaˇcem pomoc´ı bˇeˇzn´eho rozhran´ı USB, moˇznost nap´ajen´ı zaˇr´ızen´ı pˇres USB a snadn´e ovl´ad´an´ı pomoc´ı aplikac´ı vyvinut´ ych v prostˇred´ı Control Web.
ˇ ızen´ı modelu a sbˇer dat R´
- 42 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
Obr´azek 5.10: Zaˇr´ızen´ı DataLab IO4 Pouˇzit´ y syst´em DataLab IO4 se skl´ad´a z centr´aln´ıho modulu CPU/USB, kter´ y obstar´av´a komunikaci s poˇc´ıtaˇcem pomoc´ı USB a nap´ajen´ı cel´e desky (voliteln´e mezi nap´ajen´ım z USB nebo z extern´ıho zdroje). K tomuto centr´aln´ımu prvku jsou pˇripojeny n´asleduj´ıc´ı moduly: • DO1–modul rel´eov´ ych v´ ystup˚ u • AO1–modul analogov´ ych v´ ystup˚ u • AIO1–modul analogov´ ych vstup˚ u/v´ ystup˚ u - v naˇsem modelu nevyuˇzit • AI3–modul analogov´ ych vstup˚ u DO1 - modul rel´eov´ ych v´ ystup˚ u obsahuje 8 rel´eov´ ych sp´ınac´ıch kontakt˚ u. V naˇsem modelu se vyuˇz´ıvaj´ı dva kan´aly, jeden pro sp´ın´an´ı ventilu, druh´ y pro sp´ın´an´ı ˇcerpadla. Konstrukce rel´eov´ ych sp´ınaˇc˚ u umoˇzn ˇuje sp´ınat stejnosmˇern´e napˇet´ı aˇz 30 V pˇri 3 A, coˇz je pro naˇse u ´ˇcely v´ıce neˇz dostaˇcuj´ıc´ı (pouˇz´ıv´ame 12 V stejnosmˇern´e napˇet´ı). Mezn´ı kmitoˇcet cyklu (sepnut´ı/rozepnut´ı) je aˇz 50 Hz, coˇz je opˇet zcela dostaˇcuj´ıc´ı, nebot’ maxim´aln´ı frekvence sp´ın´an´ı elektromagnetick´eho ventilu (tepov´a frekvence) je 2 Hz. AO1 - obsahuje 8 analogov´ ych v´ ystup˚ u s 12 bitov´ ym D/A pˇrevodn´ıkem. Tyto v´ ystupy mohou pracovat jak v napˇet’ov´em (0-10 V), tak v proudov´em reˇzimu (4-20 mA). V naˇsem pˇr´ıpadˇe je vyuˇzit jeden napˇet’ov´ y kan´al pro ovl´ad´an´ı rychlosti ˇcerpadla. Maxim´aln´ı poˇcet ot´aˇcek ˇcerpadla za minutu je 240. Zmˇena o jednu ot´aˇcku za minutu vyˇzaduje zmˇenu vstupn´ıho napˇet’ov´eho sign´alu o 41,67 mV. Vzhledem k tomu, ˇze jeden krok D/A pˇrevodn´ıku odpov´ıd´a 2,5 mV, je tento sign´al pro pˇresn´e ovl´ad´an´ı ˇcerpadla zcela dostaˇcuj´ıc´ı. AI3 je modul 8 analogov´ ych vstup˚ u. Tyto vstupy mohou b´ yt jak napˇet’ov´e, tak proudov´e (nastaven´ı lze prov´est pro kaˇzd´ y kan´al samostatnˇe zaˇrazen´ım propojky na desce modulu. Propojka zaˇrad´ı do vstupu pˇresn´ y sn´ımac´ı odpor 10 Ω). Vzorkovac´ı frekvence modulu je 50 Hz, kter´a se d´al dˇel´ı na jednotliv´e kan´aly. V naˇsem pˇr´ıpadˇe se vyuˇz´ıvaj´ı 4 kan´aly pro mˇeˇren´ı tlaku, ostatn´ı jsou vypnut´e, to znamen´a ˇze vzorkovac´ı frekvence je 12,5 Hz. V modelu se poˇc´ıt´a s typickou tepovou frekvenc´ı 1 Hz, nav´ıc nepotˇrebujeme pˇresn´e vyhlazen´e tlakov´e kˇrivky, proto je tato vzorkovac´ı frekvence pro n´as dostaˇcuj´ıc´ı. Pˇrevod Hardware pro ovl´ad´an´ı a sbˇer dat
- 43 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
analogov´eho sign´alu na ˇc´ıslicov´ y je realizov´an pomoc´ı 16 bitov´eho sigma-delta pˇrevodn´ıku. Kaˇzd´emu kan´alu lze nav´ıc stanovit vlastn´ı napˇet’ov´ y nebo proudov´ y rozsah. To zajiˇst’uje celkovˇe vysokou pˇresnost mˇeˇren´ı.
5.5.2
Software pro ovl´ ad´ an´ı modelu a sbˇ er dat
Aplikace pro ovl´ad´an´ı modelu a sbˇer dat byla vytvoˇrena v prostˇred´ı Control Web 6.1. Programov´ y syst´em Control Web je v´ıce neˇz 15 let vyv´ıjen´ y n´astroj, umoˇzn ˇuj´ıc´ı v´ yvoj aplikac´ı pro kontrolu a ˇr´ızen´ı pr˚ umyslov´ ych syst´em˚ u, vzd´alen´ y sbˇer dat, modelov´an´ı a simulace a mnoho dalˇs´ıch funkc´ı. Programov´an´ı aplikace se d´a z velk´e ˇc´asti oznaˇcit za grafick´e programov´an´ı (podobnˇe jako napˇr´ıklad v´ yvojov´e prostˇred´ı LabVIEW). Uˇzivatel navrhuje v´ yslednou aplikaci vkl´ad´an´ım virtu´aln´ıch pˇr´ıstroj˚ u, definov´an´ım komunikaˇcn´ıch kan´al˚ u, konstant a promˇenn´ ych a vytv´aˇren´ım procedur pro jednotliv´e pˇr´ıstroje. V naˇsem pˇr´ıpadˇe byla vyvinuta aplikace pro ˇr´ızen´ı a sbˇer dat za vyuˇzit´ı zaˇr´ızen´ı DataLab IO. Pro u ´spˇeˇsnou komunikaci s t´ımto zaˇr´ızen´ım je potˇreba do programu zav´est patˇriˇcn´ y mapovac´ı (DMF) a parametrick´ y soubor (PAR). Oba soubory jsou souˇc´ast´ı syst´emu Control Web. Mapovac´ı soubor definuje vlastnosti ovl´adan´eho zaˇr´ızen´ı a zp˚ usob komunikace. Parametrick´ y soubor potom definuje jednotliv´e pˇripojen´e moduly, ˇc´ısluje jejich kan´aly, umoˇzn ˇuje vybrat, v jak´ ych jednotk´ach se bude dan´ y kan´al mˇeˇrit, popˇr´ıpadˇe umoˇzn ˇuje nˇekter´e kan´aly vypnout. N´ami pouˇzit´ y parametrick´ y soubor vypad´a n´asledovnˇe:
Obr´azek 5.11: Nastaven´ı parametrick´eho souboru V´ yznam d˚ uleˇzit´ ych kl´ıˇcov´ ych slov je uveden na pˇr´ıkladu modulu analogov´ ych vstup˚ u AI3. Kl´ıˇcov´e slovo module type definuje o jak´ y typ modulu se jedn´a, first input channel urˇcuje adresu, kde zaˇc´ın´a mapov´an´ı jednotliv´ ych kan´al˚ u (v tomto pˇr´ıpadˇe m´a kan´al jedna adresu 110, kan´al osm adresu 117), mode definuje mˇeˇr´ıc´ı rozsah kan´al˚ u (jedn´a se o hexadecim´aln´ı ˇc´ıslo, v naˇsem pˇr´ıpadˇe jsou prvn´ı ˇctyˇri kan´aly vypnut´e, zbyl´e ˇctyˇri jsou nastaveny na proudov´ y vstup ± 20 mA), a kl´ıˇcov´e slovo unit definuje, v jak´ ych jednotk´ach dan´ y kan´al mˇeˇr´ı (v tomto pˇr´ıpadˇe pro vˇsechny kan´aly nastaveny mA). O dalˇs´ıch moˇznostech Software pro ovl´ad´an´ı modelu a sbˇer dat
- 44 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
nastaven´ı parametrick´eho souboru se nem´a cenu zde rozepisovat, nebot’ syst´em Control Web je vybaven pˇrehlednou a jasnou n´apovˇedou. Po spr´avn´em zaveden´ı ovladaˇce zaˇr´ızen´ı DataLab IO je jiˇz moˇzn´e definovat pouˇzit´e kan´aly, vytvoˇrit si vlastn´ı promˇenn´e a konstanty, a n´aslednˇe je jiˇz moˇzno pˇristoupit k n´avrhu grafick´eho uˇzivatelsk´eho prostˇred´ı a k implementaci vlastn´ı funkcionality. Vytvoˇren´ y program vypad´a n´asledovnˇe:
Obr´azek 5.12: Grafick´e uˇzivatelsk´e prostˇred´ı aplikace ovl´adaj´ıc´ı model V horn´ı ˇca´sti GUI jsou ovl´adac´ı prvky modelu. Posuvn´ıkem Pump voltage[mV] se nastavuje napˇet´ı pˇriveden´e na ˇr´ıd´ıc´ı vstup pumpy, a t´ım i poˇcet jejich ot´aˇcek. Pˇrep´ınaˇc Pump ON/OFF slouˇz´ı k okamˇzit´emu spuˇstˇen´ı nebo zastaven´ı pumpy. Pˇrep´ınaˇc Pulsation ON/OFF rozhoduje, zda bude tok kontinu´aln´ı (OFF) nebo pulsatiln´ı (ON). V pˇr´ıpadˇe, ˇze model pracuje v kontinu´aln´ım reˇzimu, lze pˇrep´ınaˇcem Valve OPEN/CLOSE otev´ırat a zav´ırat ventil. V pˇr´ıpadˇe, ˇze je zapnut pulsatiln´ı reˇzim, je ventil automaticky otev´ır´an a zav´ır´an. Pomˇer doby, po kterou je ventil otevˇren/zavˇren, ˇr´ıd´ı posuvn´ık Valve opening ratio (kdy 0 znamen´a ventil otevˇren´ y po celou dobu periody, 10 znamen´a ventil zavˇren´ y po celou dobu periody). Posuvn´ık Heart rate [BPM] pak ud´av´a poˇcet pulsac´ı za minutu (tedy srdeˇcn´ı tep). V druh´e ˇc´asti GUI jsou pak informaˇcn´ı prvky. Graf Pressure after valve vykresluje tlakovou kˇrivku mˇeˇrenou za ventilem, tud´ıˇz arteri´aln´ı tlak. Graf Pressure in compensator ud´av´a hodnotu tlaku uvnitˇr kompenzaˇcn´ı n´adoby. Tuto hodnotu je tˇreba sledovat a zajistit, aby nevzrostla nad u ´nosnou mez. V programu je jako mez nastavena hodnota 150 Software pro ovl´ad´an´ı modelu a sbˇer dat
- 45 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
mmHg, coˇz je hluboko pod maxim´aln´ım povolen´ ym tlakem v kompenzaˇcn´ı n´adobˇe (750 mmHg). Graf Venous pressure monitoruje tlak na konci c´evn´ıho stromu, tedy odpov´ıd´a ven´ozn´ımu tlaku. Graf Auxillary Pressure Sensor vykresluje tlakovou kˇrivku mˇeˇrenou ˇctvrt´ ym tlakov´ ym senzorem, kter´ y je moˇzno um´ıstit na libovoln´e m´ısto ˇreˇciˇstˇe. Vˇsechny v´ yˇse uveden´e grafy vykresluj´ı hodnoty v jednotk´ach mmHg. Pod grafem vykresluj´ıc´ım arteri´aln´ı tlak se nach´az´ı jeˇstˇe jeden pomocn´ y graf, kter´ y zobrazuje otevˇren´ı/zavˇren´ı ventilu. V lev´em doln´ım rohu jsou zobrazeny hodnoty systolick´eho tlaku (SP), diastolick´eho tlaku (DP) a stˇredn´ıho arteri´aln´ıho tlaku (MAP). Systolick´ y tlak pˇredstavuje maxim´aln´ı arteri´aln´ı tlak zmˇeˇren´ y v dan´e periodˇe, diastolick´ y tlak potom pˇredstavuje minim´aln´ı arteri´aln´ı tlak zmˇeˇren´ y v dan´e periodˇe. Stˇredn´ı arteri´aln´ı tlak se potom poˇc´ıt´a z pr˚ ubˇehu arteri´aln´ıho tlaku podle pomoc´ı aritmetick´eho pr˚ umˇeru, konkr´etnˇe podle n´asleduj´ıc´ıho vzorce: Σi Part (5.2) M AP = 0 i kde M AP je stˇredn´ı arteri´aln´ı tlak [mmHg], Part je arteri´aln´ı tlak vzorku i [mmHg], i je poˇcet vzork˚ u v periodˇe Periodou se mysl´ı ˇcasov´ yu ´sek 5 sekund, po kter´em jsou tyto hodnoty vynulov´any a zmˇeˇreny znovu. Bˇehem periody jsou tyto hodnoty kaˇzdou sekundu aktualizov´any. Ukl´ad´an´ı dat nen´ı nijak softwarovˇe ˇreˇseno, a to z toho d˚ uvodu, ˇze Control Web umoˇzn ˇuje nastaven´ı periodick´eho z´alohov´an´ı dat do vlastn´ı datab´aze. V pˇr´ıpadˇe potˇreby tedy lze tuto funkci v aplikaci nastavit a zapnout.
Software pro ovl´ad´an´ı modelu a sbˇer dat
- 46 -
Kapitola 6 Experiment a namˇ eˇ ren´ e v´ ysledky C´ılem t´eto pr´ace je porovnat redistribuci vodn´ıho toku v z´avislosti na charakteru proudˇen´ı. Samozˇrejmˇe nejpˇresvˇedˇcivˇejˇs´ı by bylo mˇeˇrit pˇr´ımo pr˚ utok na nejniˇzˇs´ı u ´rovni kaˇzd´e vˇetve, ale vzhledem k tomu, ˇze nen´ı k dispozici vhodn´a metoda pˇr´ım´eho mˇeˇren´ı pr˚ utoku, museli jsme pˇristoupit k nepˇr´ım´e metodˇe mˇeˇren´ı pr˚ utoku.
6.1
Pouˇ zit´ e statistick´ e zpracov´ an´ı dat
Po proveden´ı pokus˚ u z´ısk´ame nˇekolik skupin dat, kter´e budeme cht´ıt mezi sebou nˇejak porovnat, k ˇcemuˇz vyuˇzijeme statistiku. Pˇrijmeme-li pˇredpoklad, ˇze z´ıskan´a data odpov´ıdaj´ı realizaci n´ahodn´eho v´ ybˇeru z norm´aln´ıho rozdˇelen´ı, je moˇzn´e prov´est Student˚ uv t-test rovnosti stˇredn´ıch hodnot tˇechto v´ ybˇer˚ u. Pro jistotu ale provedeme i Mann-Whitney˚ uv test rovnosti rozdˇelen´ı, kter´ y nevyˇzaduje pˇredpoklad norm´alnosti rozdˇelen´ı realizace n´ahodn´eho v´ ybˇeru. Pro pˇrehlednost si statistick´e zpracov´an´ı dat definujeme jiˇz na zaˇca´tku, a pˇri pozdˇejˇs´ım zpracov´an´ı se jiˇz budeme jen odkazovat na toto m´ısto.
6.1.1
Test rovnosti stˇ redn´ıch hodnot - Student˚ uv t-test
Student˚ uv nep´arov´ y t-test rovnosti stˇredn´ıch hodnot patˇr´ı do skupiny parametrick´ ych test˚ u a lze ho aplikovat na norm´aln´ı rozdˇelen´ı. Pˇredpoklad norm´alnosti rozdˇelen´ı, ze kter´ ych poch´az´ı n´ami mˇeˇren´a data (realizace v´ ybˇeru z rozdˇelen´ı) pˇrij´ım´ame na z´akladˇe typu mˇeˇren´ ych dat (ˇcasy do urˇcit´e ud´alosti) 1 . Nulov´a hypot´eza v tomto testu n´am ˇr´ık´a, ˇze obˇe rozdˇelen´ı ze kter´ ych poch´azej´ı realizace n´ahodn´ ych v´ ybˇer˚ u maj´ı stejnou stˇredn´ı hodnotu, neboli: H0 : µ1 = µ2
(6.1)
kde µ1 je stˇredn´ı hodnota prvn´ıho v´ybˇeru (realizovan´a aritmetick´ym pr˚ umˇerem), µ2 je stˇredn´ı hodnota druh´eho v´ybˇeru (realizovan´a aritmetick´ym pr˚ umˇerem) Nejprve si pro kaˇzd´ y v´ ybˇer urˇc´ıme v´ yznamn´e statistick´e parametry, a to stˇredn´ı hodnotu a rozptyl dat. Jelikoˇz m´ame k dispozici jen namˇeˇren´a data, nahrad´ıme tyto ukazatele aritmetick´ ym pr˚ umˇerem, respektive v´ ybˇerov´ ym rozptylem. Aritmetick´ y pr˚ umˇer spoˇc´ıt´ame jako: Σ n xi (6.2) x¯ = i=1 n kde x¯ je aritmetick´y pr˚ umˇer, xi je konkr´etn´ı jedno mˇeˇren´ı, n je poˇcet mˇeˇren´ı 1
Tento pˇredpoklad byl tak´e pro namˇeˇren´e datov´e ˇrady ovˇeˇren proveden´ım Shapiro-Wilkova testu normality v programu R-studio. Nicm´enˇe z d˚ uvodu relativnˇe n´ızk´eho poˇctu vzork˚ u provedeme pozdˇeji i Mann-Whitneyho test rovnosti v´ ybˇer˚ u.
47
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
D´ale si vypoˇcteme v´ ybˇerov´ y rozptyl dat: s2 =
Σni=1 |xi − x¯|2 n−1
(6.3)
kde s2 je v´ybˇerov´y rozptyl, x¯ je aritmetick´y pr˚ umˇer, xi je konkr´etn´ı jedno mˇeˇren´ı, n je poˇcet mˇeˇren´ı Pro kaˇzd´e mˇeˇren´ı tedy z´ısk´ame n´asleduj´ıc´ı statistick´e parametry: Aritmetick´ y pr˚ umˇer V´ ybˇerov´ y rozptyl
Set dat 1 x¯1 s21
Set dat 2 x¯1 s22
Tabulka 6.1: Vybran´e statistick´e parametry Dalˇs´ı postup z´avis´ı na tom, zda rozdˇelen´ı, ze kter´ ych poch´azej´ı naˇse realizovan´e v´ ybˇery, maj´ı stejn´ y rozptyl nebo ne. Otestujeme tedy rovnost rozptyl˚ u pomoc´ı Fisherova F-testu. Vypoˇcteme si testovac´ı statistiku F : F =
max(s21 , s22 ) min(s21 , s22 )
(6.4)
kde s21 , s22 jsou patˇriˇcn´e v´ybˇerov´e rozptyly Toto ˇc´ıslo porovn´ame s tabulkovou hodnotou F1−α/2(n1 −1,n2 −1) , tedy kvantilem Fisherova rozdˇelen´ı na hladinˇe v´ yznamnosti α se stupni volnosti (n1 − 1,n2 − 1). • Pakliˇze dost´av´ame, ˇze F < F1−α/2(n1 −1,n2 −1) , znamen´a to, ˇze mezi v´ ybˇerov´ ymi rozptyly nen´ı statisticky v´ yznamn´ y rozd´ıl. Vypoˇcteme si tedy testovac´ı statistiku t pro ovˇeˇren´ı stejnosti aritmetick´ ych pr˚ umˇer˚ u obou v´ ybˇer˚ u jako: t= q
|¯ xpuls − x¯kont | (n1 −1)·(s21 )+(n2 −1)·(s22 ) n1 +n2 −2
(6.5) ·
n1 +n2 n1 ·n2
kde x¯1 a x¯kont jsou aritmetick´e pr˚ umˇery, s21 a s22 jsou v´ybˇerov´e rozptyly, n1 a n2 jsou poˇcty mˇeˇren´ı Pro porovn´an´ı s tabulkovou hodnotou mus´ıme zn´at stupnˇe volnosti v, kter´e se v tomto pˇr´ıpadˇe vypoˇctou jako: v = n1 + n2 − 2
(6.6)
kde n1 a n2 jsou poˇcty mˇeˇren´ı ybˇerov´ ymi rozptyly statisticky • Pakliˇze dost´av´ame, ˇze F > F1−α/2(n1 −1,n2 −1) , je mezi v´ v´ yznamn´ y rozd´ıl. Vypoˇcteme si tedy testovac´ı statistiku t pro ovˇeˇren´ı stejnosti aritmetick´ ych pr˚ umˇer˚ u obou v´ ybˇer˚ u jako: t=
|¯ xpuls − x¯kont | q 2 s1 s2 + n22 n1
(6.7)
kde x¯1 a x¯kont jsou aritmetick´e pr˚ umˇery, s21 a s22 jsou v´ybˇerov´e rozptyly, n1 a n2 jsou poˇcty mˇeˇren´ı Test rovnosti stˇredn´ıch hodnot - Student˚ uv t-test
- 48 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
Pro porovn´an´ı s tabulkovou hodnotou mus´ıme zn´at stupnˇe volnosti v, kter´e se v tomto pˇr´ıpadˇe vypoˇctou jako: s2
v=
( n11 + s2
( n1 )2 1
n1 −1
s22 2 ) n2 s2
+
(6.8)
( n2 )2 2
n2 −1
kde s21 a s22 jsou v´ybˇerov´e rozptyly, n1 a n2 jsou poˇcty mˇeˇren´ı Vypoˇctenou testovac´ı statistiku porovn´ame s tabulkovou hodnotou t1−α/2(v) , tedy s hodnotou kvantilu Studentova rozdˇelen´ı na hladinˇe v´ yznamnosti α a pro stupeˇ n volnosti v. • Je-li t < t1−α/2(v) , nen´ı mezi stˇredn´ımi hodnotami (aritmetick´ ymi pr˚ umˇery) statisticky v´ yznamn´ y rozd´ıl a nem˚ uˇzeme tedy zam´ıtnou hypot´ezu H0 , tedy ˇze stˇredn´ı hodnoty se rovnaj´ı. • Je-li t > t1−α/2(v) , je mezi stˇredn´ımi hodnotami (aritmetick´ ymi pr˚ umˇery) statisticky v´ yznamn´ y rozd´ıl a m˚ uˇzeme tedy zam´ıtnou hypot´ezu H0 , tud´ıˇz stˇredn´ı hodnoty se statisticky v´ yznamnˇe liˇs´ı.
6.1.2
Test rovnosti v´ ybˇ er˚ u - Mann-Whitney˚ uv test
Mann-Whitney˚ uv test umoˇzn ˇuje porovn´avat skupiny dat poch´azej´ıc´ı z libovoln´eho rozdˇelen´ı. Pro nenorm´aln´ı rozdˇelen´ı dosahuje lepˇs´ıch v´ ysledk˚ u neˇz Student˚ uv t-test, pro norm´aln´ı rozdˇelen´ı potom obdobn´ ych. Je to velice elegantn´ı a lehce pochopiteln´ y test. Nulov´a hypot´eza H0 je v tomto testu tvrzen´ı, ˇze obˇe mˇeˇren´ı poch´az´ı ze shodn´ ych rozdˇelen´ı. V prvn´ım kroce testu prvn´ı skupinu namˇeˇren´ ych dat (napˇr´ıklad data namˇeˇren´a v kontinu´aln´ım reˇzimu) oznaˇc´ıme jako skupinu A. Druhou skupinu dat oznaˇc´ıme jako B. N´aslednˇe namˇeˇren´a data slouˇc´ıme do jedn´e skupiny (pˇri zachov´an´ı informace z kter´e skupiny kter´ y vzorek poch´az´ı) a seˇrad´ıme ve vzestupn´em poˇrad´ı. V dalˇs´ım kroku pˇriˇrad´ıme ke kaˇzd´emu vzorku jeho poˇradov´e ˇc´ıslo ve sm´ıˇsen´e skupinˇe. Pakliˇze se ve sm´ıˇsen´e skupinˇe vyskytnout stejn´e vzorky, o jejichˇz poˇrad´ı nelze jednoznaˇcnˇe rozhodnout, pˇriˇrad´ıme tˇemto stejn´ ym vzork˚ um pr˚ umˇern´e poˇrad´ı. N´aslednˇe urˇc´ıme sumu RA poˇradov´ ych ˇc´ısel pro skupinu A a sumu RB poˇradov´ ych ˇc´ısel pro skupinu B. Pˇritom plat´ı: RA + RB =
n · (n + 1) 2
(6.9)
kde n je celkov´y poˇcet vzork˚ u v obou skupin´ach N´aslednˇe vypoˇcteme testovac´ı statistiky UA a UB jako: U A = n1 · n2 +
n1 · (n1 + 1) − RA 2
(6.10)
n2 · (n2 + 1) − RB (6.11) 2 kde n1 je poˇcet vzork˚ u ve skupinˇe A, n2 je poˇcet vzork˚ u ve skupinˇe B, RA je suma poˇradov´ych ˇc´ısel vˇsech vzork˚ u ve skupinˇe A, RB je suma poˇradov´ych ˇc´ısel vˇsech vzork˚ u ve skupinˇe B UB = n1 · n2 +
Test rovnosti v´ ybˇer˚ u - Mann-Whitney˚ uv test
- 49 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
Menˇs´ı z hodnot UA a UB pouˇzijeme jako testovac´ı krit´erium. Tuto hodnotu porovn´ame s tabulkovou kritickou hodnotou Mann-Whitneyova testu pro pˇr´ısluˇsn´e n1 , n2 a na hladinˇe v´ yznamnosti α. ych • Je-li U < U(α,n1 ,n2 ) , zam´ıt´ame nulovou hypot´ezu, ˇze obˇe mˇeˇren´ı poch´az´ı ze shodn´ rozdˇelen´ı, tud´ıˇz mezi nimi je statisticky v´ yznamn´ y rozd´ıl. • Je-li U > U(α,n1 ,n2 ) , nem˚ uˇzeme zam´ıtnout nulovou hypot´ezu, ˇze obˇe mˇeˇren´ı poch´az´ı ze shodn´ ych rozdˇelen´ı, tud´ıˇz mezi nimi nen´ı statisticky v´ yznamn´ y rozd´ıl.
6.2
Prerekvizity experimentu
Proveden´ y experiment je postaven na dvou z´akladn´ıch pˇredpokladech, a sice ˇze minutov´ y v´ ydej srdeˇcn´ı (celkov´e mnoˇzstv´ı tekutiny, kter´e proteˇce syst´emem za jednu minutu) je pro kontinu´aln´ı a pulsatiln´ı proudˇen´ı stejn´e, a ˇze stˇredn´ı rychlost pohybu kuliˇcky v uzavˇren´e trubici je pˇr´ımo u ´mˇern´a pr˚ utoku touto trubic´ı. Tyto pˇredpoklady jsou dok´az´any samostatn´ ymi pokusy.
6.2.1
Pˇ redpoklad zachov´ an´ı minutov´ eho pr˚ utoku pˇ ri zmˇ enˇ e charakteru proudˇ en´ı
V koneˇcn´em experimentu nepˇr´ımo mˇeˇr´ıme pr˚ utok jednotliv´ ymi ˇc´astmi ˇreˇciˇstˇe a porovn´av´ame tyto pr˚ utoky pro pulsatiln´ı a kontinu´aln´ı proudˇen´ı. Aby toto mˇeˇren´ı mˇelo smysl, mus´ıme m´ıt zaruˇceno, ˇze pr˚ utok cel´ ym syst´emem je pro pulsatiln´ı a kontinu´aln´ı proudˇen´ı stejn´ y. Proto byl nejprve provedeno mˇeˇren´ı pr˚ utoku v z´avislosti na charakteru proudˇen´ı. Postup ovˇeˇren´ı tohoto pˇredpokladu byl n´asleduj´ıc´ı. Na zaˇca´tku mˇeˇren´ı se model pˇrepnul do kontinu´aln´ıho m´odu s v´ ydejem pˇribliˇznˇe 300 ml za minutu. Po deseti minut´ach (doba, za kterou je zaruˇceno ˇze odezn´ı vˇsechny pˇrechodov´e jevy, kompenzaˇcn´ı n´adoba ust´al´ı sv˚ uj vnitˇrn´ı pˇretlak a stˇredn´ı arteri´aln´ı tlak bude konstantn´ı) se pomoc´ı odmˇern´e n´adoby a stopek zmˇeˇril pr˚ utok syst´emem. Mezi jednotliv´a mˇeˇren´ı byla pro jistotu zaˇrazena kr´atk´a minutov´a pauza. N´aslednˇe se model pˇrepnul do pulsatiln´ıho reˇzimu pˇri zachov´an´ı stejn´eho v´ ydeje ˇcerpadla (opˇet tedy pˇribliˇznˇe 300 ml) a po deseti minut´ach probˇehla druh´a s´erie mˇeˇren´ı. Kompletn´ı namˇeˇren´a data jsou k dispozici v pˇr´ıloze A. Namˇeˇren´a data nyn´ı zpracujeme podle statistick´eho postupu, uveden´eho v ˇc´asti 6.1. Jako prvn´ı provedeme Student˚ uv t-test rovnosti stˇredn´ıch hodnot. Stanov´ıme tedy nulovou hypot´ezu H0 jako: H0 : µpuls = µkont (6.12) kde µpuls je stˇredn´ı hodnota pr˚ utoku pro pulsatiln´ı reˇzim (realizovan´a aritmetick´ym pr˚ umˇerem), µkont je stˇredn´ı hodnota pr˚ utoku pro kontinu´aln´ı reˇzim (realizovan´a aritmetick´ym pr˚ umˇerem) Z namˇeˇren´ ych dat si vypoˇcteme d˚ uleˇzit´e statistick´e parametry:
Prerekvizity experimentu
- 50 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
Aritmetick´ y pr˚ umˇer [ml/s] V´ ybˇerov´ y rozptyl[(ml/s)2 ]
Kontinu´aln´ı reˇzim
pulsatiln´ı reˇzim
302,78
300,69
21,67
18,77
Tabulka 6.2: Vybran´e statistick´e parametry mˇeˇren´ı rovnosti pr˚ utok˚ u Nejprve provedeme Fisher˚ uv F-testu rovnosti rozptyl˚ u. Dost´av´ame, ˇze: F = 1, 155 F0,975(20,20) = 2, 465 Jelikoˇz F < F0,975(20,20) , poch´az´ı obˇe mˇeˇren´ı z dat se stejn´ ym v´ ybˇerov´ ym rozptylem. M˚ uˇzeme tedy pˇrej´ıt k samotn´emu v´ ypoˇctu testovac´ı statistiky t.
t = 1, 506 t0,975(40) = 2, 021 Jelikoˇz dost´av´ame, ˇze t < t0,975(40) , nelze zam´ıtnout na hladinˇe v´ yznamnosti α = 0, 05 hypot´ezu H0 a nelze tedy vylouˇcit, ˇze jak pulsatiln´ı, tak kontinu´aln´ı proudˇen´ı m´a stejnou stˇredn´ı hodnotu pr˚ utoku. Jako druh´ y provedeme Mann-Whitney˚ uv test rovnosti v´ ybˇer˚ u. Skupina A reprezentuje mˇeˇren´ı prov´adˇen´e v kontinu´aln´ım reˇzimu, skupina B reprezentuje mˇeˇren´ı v pulsatiln´ım reˇzimu. Dost´av´ame n´asleduj´ıc´ı v´ ysledky testovac´ıch statistik:
UA = 164, 5 UB = 276, 5 U(0,95,21,21) = 113 Jelikoˇz UA > U(0,95,21,21) , nem˚ uˇzeme zam´ıtnout nulovou hypot´ezu, tud´ıˇz nem˚ uˇzeme vylouˇcit, ˇze mˇeˇren´ı poch´az´ı ze stejn´ ych v´ ybˇer˚ u. Proveden´ ymi statistick´ ymi testy se n´am nepodaˇrilo vylouˇcit moˇznost, ˇze pro kontinu´aln´ı a pulsatiln´ı reˇzim je pr˚ utok syst´emem stejn´ y, proto budeme d´ale uvaˇzovat, ˇze stejn´ y je.
6.2.2
Pˇ redpoklad spr´ avn´ e reprezentace pr˚ utoku pomoc´ı rychlosti kuliˇ cky
Druh´ y pˇredpoklad, nutn´ y pro spr´avnost v´ ysledn´eho mˇeˇren´ı, je myˇslenka, ˇze pohyb kuliˇcky, respektive pr˚ umˇern´a rychlost jej´ıho pohybu hadiˇckou, spr´avnˇe reprezentuje pr˚ utok danou hadiˇckou. Pro ovˇeˇren´ı tohoto pˇredpokladu byl model rozpojen a za elektromagnetick´ y ventil a senzor tlaku byla um´ıstˇena hadiˇcka dan´eho rozmˇeru. Do hadiˇcky byla vloˇzena
Pˇredpoklad spr´avn´e reprezentace pr˚ utoku pomoc´ı rychlosti kuliˇcky
- 51 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
kuliˇcka patˇriˇcn´e velikosti a byla zmˇeˇrena doba t, za kterou kuliˇcka projde hadiˇckou, budeli hn´ana vodou proud´ıc´ı pod zn´am´ ym tlakem. Jelikoˇz zn´ame d´elku hadiˇcky s a jej´ı plochu pr˚ uˇrezu S, vypoˇcteme si podle zn´am´ ych vzorc˚ u pr˚ utok Q hadiˇckou: v=
s t
(6.13)
Q=S·v
(6.14)
kde t je doba pohybu kuliˇcky [s], s je d´elka hadiˇcky (dr´aha pohybu) [cm], v je pr˚ umˇern´ a 2 rychlost pohybu kuliˇcky [cm/s], S je plocha pr˚ uˇrezu hadiˇcky [cm ], Q je pr˚ utok hadiˇckou [ml/s] N´aslednˇe byl pomoc´ı odmˇern´e n´adoby zmˇeˇren re´aln´ y pr˚ utok touto hadiˇckou [ml/s]. ’ Mˇeˇren´ı probˇehlo jen pro hadiˇcku o pr˚ umˇeru 4 mm, nebot nebyla k dispozici dostateˇcnˇe dlouh´a hadiˇcka o pr˚ umˇeru 2 mm a nebylo tedy moˇzn´e zmˇeˇrit rychlost pohybu kuliˇcky v t´eto u ´rovni. Mˇeˇren´ı bylo provedeno pro ˇctyˇri tlakov´e u ´rovnˇe v kontinu´aln´ım reˇzimu (30 mmHg, 50 mmHg, 60 mmHg a 70 mmHg). Namˇeˇren´e v´ ysledky jsou k dispozici v pˇr´ıloze B, zde jsou uvedeny jen nejd˚ uleˇzitˇejˇs´ı v´ ysledky:
Tlak [mmHg] 70 60 50 30
Pr˚ umˇern´ y pr˚ utok mˇeˇren´ y pomoc´ı pohybu kuliˇcky [ml/s] 10,728 9,163 7,551 5,197
Pr˚ umˇern´ y pr˚ utok mˇeˇren´ y pomoc´ı odmˇern´e n´adoby [ml/s] 9,753 8,423 6,996 4,699
Tabulka 6.3: Pr˚ utoky hadiˇckou mˇeˇren´e pomoc´ı rychlosti pohybu kuliˇcky a odmˇern´e n´adoby pro r˚ uzn´e tlakov´e u ´rovnˇe a kontinu´aln´ı tok Z namˇeˇren´ ych dat vid´ıme, ˇze mˇeˇren´ı pr˚ umˇern´eho pr˚ utoku pomoc´ı rychlosti kuliˇcky nen´ı zcela shodn´e s mˇeˇren´ım pr˚ utoku pomoc´ı odmˇern´e n´adoby. Co je ale d˚ uleˇzitˇejˇs´ı, tento rozd´ıl je line´ arnˇ e z´ avisl´ y na tlaku, pod kter´ ym tekutina proud´ı (alespoˇ n pro tlaky, ve kter´ ych se bude pohybovat naˇse mˇeˇren´ı). To ilustruje n´asleduj´ıc´ı graf:
Pˇredpoklad spr´avn´e reprezentace pr˚ utoku pomoc´ı rychlosti kuliˇcky
- 52 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
Obr´azek 6.1: Graf z´avislosti zmˇeˇren´eho pr˚ utoku pomoc´ı odmˇern´e n´adoby a vypoˇcten´eho pr˚ utoku pomoc´ı mˇeˇren´ı rychlosti pro kontinu´aln´ı tok D´ıky t´eto line´arn´ı z´avislosti lze udˇelat z´avˇer, ˇze mˇeˇren´ı pr˚ utoku pomoc´ı rychlosti pohybu kuliˇcky v hadiˇcce vˇernˇe reprezentuje skuteˇcn´ y pr˚ utok touto hadiˇckou a tud´ıˇz je tento nepˇr´ım´ y zp˚ usob mˇeˇren´ı pr˚ utoku pro naˇse u ´ˇcely pouˇziteln´ y.
Pˇredpoklad spr´avn´e reprezentace pr˚ utoku pomoc´ı rychlosti kuliˇcky
- 53 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
6.3
Fin´ aln´ı experiment
Jelikoˇz m´ame dok´az´any dva hlavn´ı pˇredpoklady, m˚ uˇzeme smˇele pˇrikroˇcit k vlastn´ımu experimentu. Jak je vidˇet na obr´azku 5.4, v modelu jsou zaˇrazeny dva mˇeˇr´ıc´ı okruhy, jeden tvoˇren´ y hadiˇckou o pr˚ umˇeru 4 mm, druh´ y hadiˇckou o pr˚ umˇeru 2 mm. Naˇse teorie prav´ı, ˇze redistribuce toku z´avis´ı na charakteru proudˇen´ı. Aby se tato teorie potvrdila, mus´ıme pˇri zmˇenˇe proudˇen´ı z kontinu´aln´ıho na pulsn´ı namˇeˇrit r˚ uzn´ y pr˚ utok tˇemito hadiˇckami. J´adrem experimentu je mˇeˇren´ı rychlosti pohybu kuliˇcek v mˇeˇr´ıc´ıch okruz´ıch. Na zaˇca´tku experimentu se nastav´ı konstantn´ı v´ ydej ˇcerpadla (pˇribliˇznˇe 300 ml). Pot´e se nech´a model pracovat cca 10 minut v kontinu´aln´ım reˇzimu (bˇehem t´eto doby se vyrovn´a tlak v kompenzaˇcn´ı n´adobˇe, odstran´ı se vˇsechny pˇr´ıpadn´e vzduchov´e bublinky v c´evn´ım syst´emu a dojde k vˇseobecn´emu ust´alen´ı). N´aslednˇe se vloˇz´ı do hadiˇcky o pr˚ umˇeru 4 mm kuliˇcka o stejn´em rozmˇeru a zmˇeˇr´ı se za jak dlouho uraz´ı zn´amou trasu. Pot´e se ˇcerpadlo pˇrepne do pulsatiln´ıho m´odu a cel´ y model se nech´a opˇet ust´alit. N´asleduje dalˇs´ı s´erie mˇeˇren´ı na hadiˇcce o pr˚ umˇeru 4 mm. Nakonec se cel´ y postup zopakuje pro hadiˇcku o pr˚ umˇeru 2 mm. Mezi jednotliv´ ymi mˇeˇren´ımi n´ahodnˇe zmˇeˇr´ıme celkov´ y pr˚ utok syst´emem (pomoc´ı odmˇern´e n´adoby), pro kontrolu, ˇze celkov´ y pr˚ utok syst´emem je st´ale stejn´ y, Z´ısk´ame tedy sadu 4 mˇeˇren´ı ˇcasu kter´ y kuliˇcka potˇrebuje k proplut´ı hadiˇckou. Veˇsker´e namˇeˇren´e v´ ysledky jsou k dispozici v pˇr´ıloze C.
6.3.1
V´ ysledky pro mˇ eˇ ren´ı pr˚ utoku hadiˇ ckou o pr˚ umˇ eru 4 mm
Namˇeˇren´a data zpracujeme podle statistick´eho postupu, uveden´eho v ˇc´asti 6.1. Jako prvn´ı provedeme Student˚ uv t-test rovnosti stˇredn´ıch hodnot. Stanov´ıme tedy nulovou hypot´ezu H0 jako: H0 : µpuls = µkont (6.15) kde µpuls je stˇredn´ı hodnota pr˚ utoku pro pulsatiln´ı reˇzim (realizovan´a aritmetick´ym pr˚ umˇerem), µkont je stˇredn´ı hodnota pr˚ utoku pro kontinu´aln´ı reˇzim (realizovan´a aritmetick´ym pr˚ umˇerem) Z namˇeˇren´ ych dat si vypoˇcteme d˚ uleˇzit´e statistick´e parametry: Pr˚ umˇer [s] V´ ybˇerov´ y rozptyl [s2 ]
Kontinu´aln´ı reˇzim 6,815
pulsatiln´ı reˇzim 7,570
0,094
0,042
Tabulka 6.4: Vybran´e statistick´e parametry mˇeˇren´ı pr˚ utoku hadiˇckou o pr˚ umˇeru 4 mm Nejprve provedeme Fisher˚ uv F-testu rovnosti rozptyl˚ u. Dost´av´ame, ˇze: F = 2, 238 F0,975(20,20) = 2, 465 Jelikoˇz F < F0,975(20,20) , poch´az´ı obˇe mˇeˇren´ı z dat se stejn´ ym v´ ybˇerov´ ym rozptylem. M˚ uˇzeme tedy pˇrej´ıt k samotn´emu v´ ypoˇctu testovac´ı statistiky t. t = 9, 382 t0,975(40) = 2, 021 Fin´aln´ı experiment
- 54 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
Jelikoˇz dost´av´ame, ˇze t > t0,975(40) , zam´ıt´ame na hladinˇe v´ yznamnosti α = 0, 05 hypot´ezu H0 , ˇze pr˚ utok touto u ´rovn´ı m´a stejnou stˇredn´ı hodnotu pro arteri´aln´ı a pulsatiln´ı proudˇen´ı. Jako druh´ y provedeme Mann-Whitney˚ uv test rovnosti v´ ybˇer˚ u. Skupina A reprezentuje mˇeˇren´ı prov´adˇen´e v kontinu´aln´ım reˇzimu, skupina B reprezentuje mˇeˇren´ı v pulsatiln´ım reˇzimu. Dost´av´ame n´asleduj´ıc´ı v´ ysledky testovac´ıch statistik:
UA = 1 UB = 440 U(0,95,21,21) = 113 Jelikoˇz UA < U(0,95,21,21) , zam´ıt´ame nulovou hypot´ezu, ˇze mˇeˇren´ı z´ıskan´a v kontinu´aln´ım a pulsatiln´ım reˇzimu jsou shodn´a. Obˇema testy se podaˇrilo zam´ıtnout nulovou hypot´ezu, jin´ ymi slovy lze ˇr´ıct, ˇze na hladinˇe v´ yznamnosti α=0,05 existuje rozd´ıl v dobˇe, kterou kuliˇcka potˇrebuje k proplut´ı danou hadiˇckou, v z´avislosti na pulsatilitˇe proudˇen´ı.
6.3.2
V´ ysledky pro mˇ eˇ ren´ı pr˚ utoku hadiˇ ckou o pr˚ umˇ eru 2 mm
Namˇeˇren´a data zpracujeme podle statistick´eho postupu, uveden´eho v ˇc´asti 6.1. Jako prvn´ı provedeme Student˚ uv t-test rovnosti stˇredn´ıch hodnot. Stanov´ıme tedy nulovou hypot´ezu H0 jako: H0 : µpuls = µkont (6.16) kde µpuls je stˇredn´ı hodnota pr˚ utoku pro pulsatiln´ı reˇzim (realizovan´a aritmetick´ym pr˚ umˇerem), µkont je stˇredn´ı hodnota pr˚ utoku pro kontinu´aln´ı reˇzim (realizovan´a aritmetick´ym pr˚ umˇerem) Z namˇeˇren´ ych dat si vypoˇcteme d˚ uleˇzit´e statistick´e parametry: Pr˚ umˇer [s] V´ ybˇerov´ y rozptyl [s2 ]
Kontinu´aln´ı reˇzim 24,782
pulsatiln´ı reˇzim 15,467
1,764
0,354
Tabulka 6.5: Vybran´e statistick´e parametry mˇeˇren´ı pr˚ utoku hadiˇckou o pr˚ umˇeru 2 mm Nejprve provedeme Fisher˚ uv F-testu rovnosti rozptyl˚ u. Dost´av´ame, ˇze: F = 4, 983 F0,975(20,20) = 2, 465 Jelikoˇz F > F0,975(20,20) , nepoch´azej´ı obˇe mˇeˇren´ı z dat se stejn´ ym v´ ybˇerov´ ym rozptylem. Testovac´ı statistiku t tedy budeme poˇc´ıtat podle postupu pro data poch´azej´ıc´ı z rozdˇelen´ı s r˚ uzn´ ym rozptylem.
t = 29, 331 t0,975(40) = 2, 021 V´ ysledky pro mˇeˇren´ı pr˚ utoku hadiˇckou o pr˚ umˇeru 2 mm
- 55 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
Jelikoˇz dost´av´ame, ˇze t > t0,975(40) , zam´ıt´ame na hladinˇe v´ yznamnosti α = 0, 05 hypot´ezu H0 , ˇze pr˚ utok hadiˇckou o pr˚ umˇeru 2 mm m´a stejnou stˇredn´ı hodnotu pro pulsatiln´ı a kontinu´aln´ı proudˇen´ı. Jako druh´ y provedeme Mann-Whitney˚ uv test rovnosti v´ ybˇer˚ u. Skupina A reprezentuje mˇeˇren´ı prov´adˇen´e v kontinu´aln´ım reˇzimu, skupina B reprezentuje mˇeˇren´ı v pulsatiln´ım reˇzimu. Dost´av´ame n´asleduj´ıc´ı v´ ysledky testovac´ıch statistik: UA = 0 UB = 441 U(0,95,21,21) = 113 Jelikoˇz UA < U(0,95,21,21) , zam´ıt´ame nulovou hypot´ezu, ˇze mˇeˇren´ı z´ıskan´a v kontinu´aln´ım a pulsatiln´ım reˇzimu jsou shodn´a. Obˇema testy se podaˇrilo zam´ıtnout nulovou hypot´ezu, jin´ ymi slovy lze ˇr´ıct, ˇze na hladinˇe v´ yznamnosti α=0,05 existuje rozd´ıl v dobˇe, kterou kuliˇcka potˇrebuje k proplut´ı danou hadiˇckou, v z´avislosti na pulsatilitˇe proudˇen´ı.
6.3.3
Shrnut´ı v´ ysledk˚ u
Mˇeˇren´ım jsme dok´azali, ˇze pr˚ utok jednotliv´ ymi u ´rovnˇemi c´evn´ıho stromu nen´ı pro kontinu´aln´ı a pulsatiln´ı proudˇen´ı stejn´ y. Pˇrepoˇcteme-li si dobu pohybu kuliˇcky v hadiˇck´ach na jejich rychlost (d´ıky zn´am´e d´elce hadiˇcky) a d´ale na pr˚ utok hadiˇckou (d´ıky zn´am´emu pr˚ umˇeru), m˚ uˇzeme v´ ysledky shrnout do n´asleduj´ıc´ıch dvou graf˚ u:
Obr´azek 6.2: Stˇredn´ı hodnoty pr˚ utoku se smˇerodatnou odchylkou pro pulsatiln´ı a kontinu´aln´ı proudˇen´ı, vlevo pro hadiˇcku 4 mm, vpravo pro hadiˇcku 2 mm Nejenˇze se n´am tedy podaˇrilo pozorovat zv´ yˇsen´ı pr˚ utoku nejniˇzˇs´ı vˇetv´ı c´evn´ıho ˇreˇciˇstˇe pˇri pulsatiln´ım toku v porovn´an´ı s kontinu´aln´ım (z 0,0763 ± 0,004 ml/s na 0,122 ± 0,005 ml/s), ale tak´e se n´am podaˇrilo pozorovat sn´ıˇzen´ı pr˚ utoku druhou nejniˇzˇs´ı u ´rovn´ı pˇri pulsatiln´ım toku v porovn´an´ı s kontinu´aln´ım (z 2,770 ± 0,125 ml/s na 2,491 ± 0,066 ml/s). To n´as vede k z´avˇeru, ˇze charakter proudˇen´ı skuteˇcnˇe ovlivˇ nuje redistribuci tekutiny v c´evn´ım ˇreˇciˇsti. Shrnut´ı v´ ysledk˚ u
- 56 -
Kapitola 7 Diskuze v´ ysledk˚ u Proveden´ ym experimentem jsme z´ıskali siln´e d˚ ukazy pro domnˇenku, ˇze pulsatilita proudˇen´ı tekutiny a vlastnosti c´evn´ıho ˇreˇciˇstˇe maj´ı vliv na redistribuci toku v jednotliv´ ych u ´rovn´ıch ˇreˇciˇstˇe. Ke stejn´emu z´avˇeru doˇsel i Atsushi Baba se sv´ ym t´ ymem ve studii [5], kteˇr´ı zkoumali vliv pulsatility na distribuci krve in vivo. Ti vˇsak vysvˇetlovali pokles rychlosti erytrocyt˚ u v kapil´ar´ach a sn´ıˇzen´ı poˇctu perfundovan´ ych kapil´ar vlivem regulaˇcn´ıch mechanism˚ u kardiovaskul´arn´ı soustavy. Toto tvrzen´ı se ned´a vylouˇcit, nebot’ regulaˇcn´ı mechanismy kardiovaskul´arn´ı soustavy jsou sloˇzit´e syst´emy, kter´e jsou nastaveny na fyziologick´ y pulsatiln´ı tok. V pˇr´ıpadˇe, ˇze je tento pulsatiln´ı tok nahrazen tokem kontinu´aln´ım, mohou b´ yt regulaˇcn´ı mechanismy zmateny, mohou zp˚ usobit vazokonstrikci c´evn´ıho syst´emu periferi´ı nebo dokonce sn´ıˇzit perfuzi vit´aln´ıch org´an˚ u (j´atra, ledviny, pl´ıce). V t´eto souvislosti tak´e stoj´ı za zm´ınku (nejzn´amˇejˇs´ı) pˇr´ıpad dvou osob, kter´ ym byla implantov´ana srdeˇcn´ı n´ahrada s ˇcistˇe kontinu´aln´ım tokem. Jedn´a se o Ameriˇcana Craiga ˇ Lewise a Cecha Jakuba Hal´ıka. Obˇema muˇz˚ um muselo b´ yt kv˚ uli v´aˇzn´e nemoci odstranˇeno srdce, kter´e bylo nahrazeno pumpou s ˇcistˇe kontinu´aln´ım tokem. Craig Lewis zemˇrel po 6 t´ ydnech na selh´an´ı ledvin a jater (podle vˇsech zdroj˚ u zp˚ usoben´e amyloidosou, tedy chorobou, kter´a vznik´a v d˚ usledku ukl´ad´an´ı ˇspatn´ ych protein˚ u v org´anech, a kter´a p˚ uvodnˇe napadla i jeho srdce). Jakub Hal´ık trpˇel vz´acn´ ym n´adorov´ ym onemocnˇen´ım, kv˚ uli kter´emu mu bylo vyjmutu srdce, kter´e bylo nahrazeno dvˇema LVAD pˇr´ıstroji. Po operaci ˇzil Jakub Hal´ık 6 mˇes´ıc˚ u, n´aslednˇe vˇsak zemˇrel na selh´an´ı jater, zp˚ usoben´e houbovou infekc´ı. N´am se podaˇrilo naj´ıt d˚ ukazy pro tvrzen´ı, ˇze stejn´ y jev, jako vypozoroval t´ ym A. Baby (tedy pokles rychlosti proudˇen´ı a pr˚ utoku v menˇs´ıch c´ev´ach v porovn´an´ı se ˇsirˇs´ımi) m˚ uˇze b´ yt zp˚ usoben i ˇcistˇe mechanick´ ymi vlastnostmi proudˇen´ı a vlastnostmi c´evn´ıho ˇreˇciˇstˇe. V t´eto pr´aci jsme vych´azeli z bˇeˇzn´e praxe, kdy je pulsatiln´ı a kontinu´aln´ı proudˇen´ı porovn´av´ano z hlediska stˇredn´ıho arteri´aln´ıho tlaku. D˚ uleˇzit´ ym faktem je, ˇze pˇri stejn´em stˇredn´ım arteri´aln´ım tlaku je zachov´an minutov´ y pr˚ utok syst´emem, nez´avisle na charakteru proudˇen´ı. Nicm´enˇe pˇri kontinu´aln´ım proudˇen´ı stoupne pr˚ utok v ˇsirˇs´ıch c´ev´ach a klesne v c´ev´ach nejmenˇs´ıch. ¨ Jak vˇsak p´ıˇse A. Undar ve sv´e studii [9], porovn´avat pulsatiln´ı a kontinu´aln´ı tok z hlediska stˇredn´ıho arteri´aln´ıho tlaku zˇrejmˇe nen´ı zcela vhodn´e. Optim´alnˇe by se totiˇz mˇely typy proudˇen´ı mezi sebou porovn´avat pomoc´ı energeticky ekvivalentn´ıho tlaku, kter´ y je pro pulsatiln´ı proudˇen´ı vyˇsˇs´ı neˇz pro kontinu´aln´ı proudˇen´ı stejn´eho stˇredn´ıho arteri´aln´ıho tlaku. Dalˇs´ım d˚ uleˇzit´ ym jevem, kter´ y byl v syst´emu vypozorov´an, je vliv pulsatiln´ıho proudˇen´ı na uvolˇ nov´an´ı pˇrek´aˇzek v ˇreˇciˇsti. Nejv´ıce se tento efekt projevoval na mal´ ych vzduchov´ ych bublink´ach. Zat´ımco pˇri pulsatiln´ım proudˇen´ı se vzduchov´e bublinky ˇcasem uvolnily samy, v kontinu´aln´ım reˇzimu se pˇrilepily na stˇenu hadiˇcky a bylo nutn´e je manu´alnˇe (poklepem na hadiˇcku, tedy vytvoˇren´ım pulsu) uvolnit. Samozˇrejmˇe ve fyziologick´em ˇreˇciˇsti se takto velk´e bublinky nevyskytuj´ı, na druhou stranu krev je suspenze rozliˇcn´ ych prvk˚ u r˚ uzn´ ych hustot, a nab´ız´ı se tedy ot´azka, zda by kontinu´aln´ı proudˇen´ı nevedlo k usazov´an´ı v´ yraznˇe 57
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
tˇeˇzˇs´ıch nebo lehˇc´ıch ˇc´ast´ı u stˇeny c´evy. N´ami vytvoˇren´ y model byl sestrojen na z´akladˇe mnoha zjednoduˇsen´ı a pˇredpoklad˚ u. Krev byla nahrazena vodou, kter´a m´a v´ yraznˇe niˇzˇs´ı viskozitu. C´evy byly nahrazeny silikonov´ ymi hadiˇckami, kter´e jsou vˇseobecnˇe rigidnˇejˇs´ı a z˚ ust´avaj´ı otevˇreny i pˇri nulov´em tlaku (to rozhodnˇe pro lidsk´e c´evy neplat´ı, obzvl´aˇstˇe pro c´evy menˇs´ıch pr˚ umˇer˚ u a pro v´eny, kter´e by se pˇri nulov´em vnitˇrn´ım tlaku uzavˇrely v d˚ usledku p˚ usoben´ı tlaku okoln´ıch tk´an´ı - tzv. critical closing pressure). Pˇresto si dovol´ıme tvrdit, ˇze vytvoˇren´ y model dos´ahl minim´alnˇe replikaˇcn´ı platnosti pro studium vlivu redistribuce tekutiny.
- 58 -
Kapitola 8 Z´ avˇ er V t´eto pr´aci byl sestrojen mechanick´ y model kardiovaskul´arn´ıho syst´emu. Tekutina je poh´anˇena peristaltickou pumpou, funkce srdce je simulov´ana kompenzaˇcn´ı n´adobou s elektromagnetick´ ym ventilem, c´evn´ı ˇreˇciˇstˇe je nahrazeno vˇetvenou strukturou silikonov´ ych hadiˇcek s postupnˇe sniˇzuj´ıc´ımi se pr˚ umˇery. Na nˇekolik m´ıstech c´evn´ıho ˇreˇciˇstˇe jsou inˇ ızen´ı a sbˇer dat na hardwarov´e u tegrov´any tlakov´e senzory. R´ ´rovni je realizov´ano pomoc´ı zaˇr´ızen´ı DataLab IO4, na softwarov´e u ´rovni potom pomoc´ı programu vyvinut´eho v prostˇred´ı ControlWeb. Model umoˇzn ˇuje generovat pulsatiln´ı a kontinu´aln´ı proudˇen´ı. Tlakov´a sloˇzka generovan´eho toku v pulsatiln´ım modu nekop´ıruje pˇresnˇe fyziologick´e hodnoty mˇeˇren´e ve velk´ ych tepn´ach lidsk´eho tˇela, ale jej´ı tvar je podobn´ y. Rozd´ıl (jedn´a se pˇredevˇs´ım o v´ yˇsku systolick´eho a diastolick´eho tlaku) si vysvˇetlujeme n´ızk´ ym perif´ern´ım odporem modelu. Pˇr´ım´e mˇeˇren´ı pr˚ utoku nebylo v modelu z d˚ uvodu vysok´ ych finanˇcn´ıch n´aklad˚ u implementov´ano. Nicm´enˇe byla navrˇzena metoda nepˇr´ım´eho mˇeˇren´ı pr˚ utoku pomoc´ı sledov´an´ı stˇredn´ı rychlosti pohybu kuliˇcky v hadiˇcce. Pr˚ utok vypoˇc´ıtan´ y pomoc´ı stˇredn´ı rychlosti pohybu kuliˇcky je line´arnˇe z´avisl´ y na skuteˇcn´em pr˚ utoku (minim´alnˇe pro naˇse tlakov´e rozsahy), a tud´ıˇz je tato metoda platn´a. Hlavn´ım experimentem se povedlo potvrdit teorii, ˇze redistribuce tekutiny v jednotliv´ ych u ´rovn´ıch c´evn´ıho stromu je z´avisl´a na charakteru proudˇen´ı. Pro kontinu´aln´ı proudˇen´ı byl namˇeˇren vyˇsˇs´ı pr˚ utok ve vyˇsˇs´ıch u ´rovn´ıch syst´emu (2,770 ± 0,125 ml/s, v porovn´an´ı s pr˚ utokem 2,491 ± 0,066 ml/s v pulsatiln´ım reˇzimu). Pro pulsatiln´ı proudˇen´ı byl zjiˇstˇen n´ar˚ ust pr˚ utoku na nejniˇzˇs´ı u ´rovni c´evn´ıho stromu (0,122 ± 0,005 ml/s, v porovn´an´ı s 0,0763 ± 0,004 ml/s pro kontinu´aln´ı reˇzim). Pro pˇresvˇedˇcivˇejˇs´ı v´ ysledky by bylo vhodn´e do modelu zaˇradit syst´em pro automatick´e mˇeˇren´ı pr˚ utoku (napˇr. pomoc´ı teplotn´ı nebo barvivov´e diluce), a nahradit vodu tekutinou, kter´e by se sv´ ymi vlastnostmi v´ıce bl´ıˇzila krvi (napˇr. zˇredˇen´ y propylenglykol) a samozˇrejmˇe prov´est v´ıce mˇeˇren´ı.
59
Literatura [1] Nichols, W., O’Rourke, M., (2011): McDonal’s blood flow in arteries: theoretic, experimental, and clinical principles. (6. vyd´an´ı) London: Hodder Arnold. ISBN 978-034-0985-014. [2] Sembulingam, K., (2006): Essentials of medical physiology. (4. vyd´an´ı) New Delhi: Jaypee Brothers Medical Publishers. ISBN 978-034-0985-014. [3] Trojan, S. et al., (2003): L´ekaˇrsk´a fyziologie. (4. vyd´an´ı) Praha: Gradda Publishing a.s. ISBN 978-818-0618-260 [4] Kittnar, O. et al., (2003): L´ekaˇrsk´a fyziologie. Praha: Gradda Publishing a.s. ISBN 978-80-247-3068-4 [5] Baba, A. et al., (2004): Microcirculation of the bulbar conjunctiva in the goat implanted with a total artificial heart: effects of pulsatile and nonpulsatile Flow. ASAIO journal. [6] Vaˇsk˚ u, J. et al., (2007): Acute and chronic consequences of non-pulsatile blood flow pattern in longterm total artificial heart experiment. Pathophysiology, vol. 14, issue 2, s. 327-330. DOI 10.1201/b14731-34. ˇ e Vysok´e Uˇcen´ı Tech[7] Eck, V., Raz´ım, M., (1996): Biokybernetika. Praha: Cesk´ nick´e. [8] Holˇc´ık, J., Fojt, O., (2001): Modelovn´ı biologick´ych syst´em˚ u (vybran´e kapitoly). Brno: Vysok´e Uˇcen´ı Technick´e. ISBN 80-214-2023-5 ¨ [9] Undar, A. et al., (1999): Defining pulsatile perfusion: quantification in terms of energy equivalent pressure. International Society for Artificial Organs. [10] Shaun, D., (2009): Simulation and development of mock circulation loop with variable resistance. Brisbane: Quennsland University of Technology. [11] Loˇzek, M., Havl´ık, J., Horˇc´ık, Z., (2013): N´avrh mechanick´eho modelu kardioˇ e Vysok´e Uˇcen´ı Technick´e. vaskul´arn´ıho syst´emu. Praha: Cesk´ [12] Mushi, S., Yu, Y., (2008): Control of a mock circulatory system to simulate the short-term baroreflex. Washington: American Control Conference. [13] Kirklin, J., Naftel, D., (2008): Mechanical circulatory support: registering a therapy in evolution. Uveˇrejnˇeno v: Circulation: Heart Failure [online] DOI: 10.1161/CIRCHEARTFAILURE.108.782599
60
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
[14] commons.wikimedia.org (listopad 2014) http://commons.wikimedia.org/ [15] Web spoleˇcnosti SynCardia (listopad 2014) http://www.syncardia.com/ [16] Web spoleˇcnosti Thoratec (listopad 2014) http://www.thoratec.com/ [17] Vietnam open educational resources (listopad 2014, licence Creative Commons Attribution 3.0) http://voer.edu.vn/ [18] Greenspon, A. et al., (2012) Trends in permanent pacemaker implantation in the United States from 1993 to 2009. Journal of the American college of cardiology, online (prosinec 2014) http://content.onlinejacc.org/article.aspx?articleid=1358197
LITERATURA
- 61 -
Pˇ r´ıloha A Data pro mˇ eˇ ren´ı zachov´ an´ı minutov´ eho v´ ydeje ˇ cerpadla pˇ ri zmˇ enˇ e charakteru proudˇ en´ı
Mˇeˇren´ı 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Pr˚ umˇer [ml/s] V´ ybˇerov´ y rozptyl [(ml/s)2 ]
Pr˚ utok pro kontinu´aln´ı reˇzim [ml/min] 302,00 305,14 303,90 302,57 311,10 293,97 299,40 303,44 305,81 297,77 306,75 299,60 302,98 309,70 299,40 303,34 303,29 293,26 303,44 299,95 309,65 302,78
Pr˚ utok pro pulsatiln´ı reˇzim [ml/min] 308,48 300,25 305,14 305,03 293,59 297,37 298,11 300,70 302,78 305,34 303,13 295,91 304,47 301,15 298,61 302,78 295,03 297,91 301,15 305,29 292,40 30,69
21,67
18,77
Tabulka A.1: Minutov´ y pr˚ utok syst´emem pro pulsatiln´ı a kontinu´aln´ı mˇeˇren´ı
62
Pˇ r´ıloha B Data pro mˇ eˇ ren´ı porovn´ an´ı pr˚ utoku mˇ eˇ ren´ eho pomoc´ı rychlosti kuliˇ cky a odmˇ ern´ e n´ adoby
Tlak [mmHg]
70
60
50
30
Pr˚ utok mˇeˇren´ y pomoc´ı pohybu kuliˇcky [ml/s] 10,87 10,44 10,82 10,82 10,74 9,38 8,88 9,18 9,28 9,11 7,32 7,61 7,63 7,56 7,65 5,07 5,31 5,17 5,20 5,23
Pr˚ utok mˇeˇren´ y pomoc´ı odmˇern´e n´adoby [ml/s] 9,74 9,74 9,76 9,76 7,76 8,36 8,53 8,36 8,42 8,35 6,96 7,04 7,03 6,93 7,02 4,52 4,71 4,77 4,77 4,75
Tabulka B.1: Pr˚ utoku hadiˇckou o pr˚ umˇeru 4 mm v kontinu´aln´ım reˇzimu
63
Pˇ r´ıloha C Data pro porovn´ an´ı redistribuce toku pˇ ri zmˇ enˇ e charakteru proudˇ en´ı ˇ ıslo mˇeˇren´ı C´ 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Pr˚ umˇer [s] V´ ybˇerov´ y rozptyl [s2 ]
Hadiˇcka o pr˚ umˇeru 4 mm Kontinu´aln´ı pulsatiln´ı reˇzim [s] reˇzim [s] 7,28 7,49 7,31 7,72 6,75 7,93 6,45 7,48 6,72 7,84 6,98 7,81 6,35 7,65 7,01 7,69 6,35 7,38 6,52 7,36 6,85 7,29 6,96 7,46 7,25 7,43 6,35 7,36 6,68 7,51 6,66 7,79 6,81 7,35 7,23 7,42 7,02 7,48 6,75 7,56 6,84 7,96 6,82 7,57 0,094
0,042
Hadiˇcka o pr˚ umˇeru 2 mm Kontinu´aln´ı pulsatiln´ı reˇzim [s] reˇzim [s] 26,06 16,11 27,65 15,09 24,58 14,69 21,62 15,04 26,18 14,56 23,26 16,30 24,56 15,69 26,53 14,72 23,69 15,98 25,53 16,12 24,69 15,42 26,53 16,39 24,36 14,96 25,01 15,15 24,69 15,78 23,99 16,42 24,13 14,99 24,89 15,36 23,59 15,85 24,10 15,21 24,79 14,96 24,78 15,47 1,764
0,355
Tabulka C.1: Doba pohybu kuliˇcky v hadiˇcce o dan´em pr˚ umˇeru s dan´ ym charakterem prouden´ı
64
Pˇ r´ıloha D Dodateˇ cn´ y obrazov´ y materi´ al
Obr´azek D.1: Fotografie sestaven´eho c´evn´ıho stromu
65
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
Obr´azek D.2: Fotografie sestaven´eho c´evn´ıho stromu
Obr´azek D.3: Detail kompenzaˇcn´ı n´adoby
- 66 -
Model kardiovaskul´arn´ıho syst´emu a jeho chov´an´ı v z´avislosti na typu proudˇen´ı
Obr´azek D.4: Detail kompenzaˇcn´ı n´adoby
Obr´azek D.5: Fotografie pouˇzit´e pumpy
- 67 -
Pˇ r´ıloha E Obsah pˇ riloˇ zen´ eho CD Na pˇriloˇzen´em CD najdete elektronickou verzi t´eto pr´ace a archiv obsahuj´ıc´ı program pro ovl´ad´an´ı modelu. Struktura je n´asleduj´ıc´ı: • DP2015 Michal Pinc • Ovladani modelu KVS.zip
68