Biomechanical musculoskeletal model

Předkládaná práce je zaměřena na svalově-kosterní modelování, především pak na výpočet svalových sil a ramen momentů při libovolném pohybu s využitím nové metody určení průběhů svalů. Jejím hlavním přínosem je vývoj unikátní metody založené na svalovém obepínání anuloidů, která výrazné snižuje nedostatky již existujících metod pro určení svalových trajektorií. Metoda je vyvinuta pro výpočet korektního tvaru svalu při jakékoliv konfiguraci kloubů. Je založena na obecné známé metodě svalového obepínání série překážek tvořených tuhými geometrickými tvary a nahrazujících okolní tkáně, obecné známá jako metoda obstacle-set. Z důvodu vylepšení původní metody byly překážky tvaru koule či válce nahrazeny anuloidy. Nové vzniklá metoda dále umožňuje automatický výpočet umístění svalových úponů; pozic, natočení a poloměrů jednotlivých anuloidů; nárůst aktuálního fyziologického průřezu svalu během kontrakce či změny tvaru svalu s ohledem na sousedící svalové skupiny. Veškerá geometrie metody je založena na MRI a počtu uvažovaných svalových vláken. Dílčím cílem studie je vytvořit jednoduchý model ramene v programu MATLAB, který obsahuje pouze dvojhlavý sval pažní a je založen na nové vyvinuté metodě obepínání anuloidu. Touto cestou je prezentována implementace, použití, výhody a nevýhody této metody. Kosti modelu jsou nahrazeny tuhými tělesy spojenými reálnými klouby; skutečné chování svalů je simulováno modelem Hillova typu. Pro potřeby této práce jsou pohyby lopatky a klíční kosti zanedbány. Svalový komplex je prezentován elastickými svalovými vlákny zanedbatelného tření generující stejnou sílu po celé své délce a obepínající sousední struktury nahrazené anuloidy. Pro validaci modelu a metody obepínání anuloidu je simulován pohyb pažní kosti - abdukce a přední flexe do úhlu 90°. Trajektorie svalových vláken, síly ve svalech, aktuální délka a momentová ramena svalů jsou poté porovnána s výsledky obdobných modelů prezentovaných v literatuře, s elektromyografickým měřením a se dvěma modely ramene sestavených v programu AnyBody Modeling System. Výsledky prokazují úspěšnou validaci hlavních akčních členů abdukce a přední flexe ramene. Nová metoda svalového obepínání anuloidů je vhodnou metodou pro simulaci všech kloubů lidského těla - především pro komplikované klouby jako je např. ramenní komplex, či pro všechny typy svalů - silný, slabý, plochý, dlouhý, krátký, aj. Prezentovaná studie také stručné představuje anatomii a fyziologii ramenního komplexu, nabízí rešerši existujících ramenních modelů a metod pro výpočet svalové trajektorie a do větších detailů popisuje dynamiku vázaných mechanických systémů v prostoru. Závěrem lze říci, že metoda svalového obepínání anuloidů je užitečným nástrojem při svalové-kosterním modelování.NeobhájenoPresented thesis work is focused on musculoskeletal modeling, especially on muscle forces and moment arms calculation using the new method for muscle path determination. This method is based on obstacle-set method. However, the new torus obstacles are implemented instead of standard obstacles such as spheres and cylinders to improve the original process of muscle wrapping. This method also enables the automatic calculation of muscle lines attachments; positions, rotation and radius of torus obstacles originated from MRI and respecting the input number of muscle lines set by the user. The torus-obstacle method also considers the muscle bulging up as well as changes of muscle shapes influenced by surrounding muscles. The case of this study is to create the simple shoulder model in MATLAB including the deltoid muscle and using developed torus-obstacle method. Thanks that, the implementation, usage, advantages and disadvantages of presented method are shown. The bones are modeled by rigid bodies connected by real joints; the real muscle behavior is simulated by Hill-type model. For purpose of this work, the scapula and the clavicle are fixed. The muscle complex is replaced by elastic frictionless muscle lines of action generating the same force along the whole band and wrapping around the neighboring structures replaced by torus obstacles. The humeral abduction and forward flexion till 90° are simulated to validate the model and also the wrapping method. The paths of muscle lines, muscle forces, actual lengths and the muscle moment arms are compared to the similar models published in literature, to the electromyography measurement and to two shoulder models built in AnyBody Modeling System. The results show the successful validation of major actuators of abduction and forward flexion. In addition, the method is absolutely not time-consuming. The new torus-obstacle method is suitable for all human body joints - especially for complicated joints as shoulder complex, for all muscles - thick, thin, shallow, long, short etc. Presented study also introduces briefly the anatomy and physiology of the shoulder complex, offers the research of existing shoulder models and methods for muscle path definition and describes the multibody spatial dynamics in more details. In conclusion, developed torus-obstacle method designed for muscle trajectory computation in musculoskeletal modeling seems to be useful tool