MODELOWANIE KONSTYTUTYWNE TKANKI KOSTNEJ I ..., Ćwiczenia z Mechanika

Pobierz MODELOWANIE KONSTYTUTYWNE TKANKI KOSTNEJ I ... i więcej Ćwiczenia w PDF z Mechanika tylko na Docsity! Aktualne Problemy Biomechaniki, nr 4/2010 185 Cyprian SUCHOCKI, Zakład Konstrukcji Maszyn i Inżynierii Biomedycznej, Instytut Mechaniki i Poligrafii, Politechnika Warszawska, Warszawa MODELOWANIE KONSTYTUTYWNE TKANKI KOSTNEJ I METODY WYZNACZANIA STAŁYCH MATERIAŁOWYCH Streszczenie. Pracę rozpoczyna krótki rys historyczny dotyczący rozwoju mechaniki ośrodków ciągłych od początków XIX w., aż do dnia dzisiejszego. Omówiono kilka wariantów równań konstytutywnych, które mogą posłużyć do opisu własności mechanicznych tkanki kostnej. Uwzględniono przy tym najnowsze osiągnięcia w zakresie mechaniki ośrodków ciągłych. W dalszej części pracy opisano problem wyznaczania stałych materiałowych. Skupiono się przy tym na algorytmie aproksymacji Levenberg'a - Marquardt'a. 1. WSTĘP Teoria sprężystości rozwinęła się jako pierwszy z działów mechaniki ośrodków ciągłych (MOC) jeszcze na początku XIX w. Podstawowym modelem materiału wykorzystywanym w owej teorii jest model materiału o liniowej charakterystyce naprężenie - odkształcenie i jednakowych własnościach mechanicznych we wszystkich kierunkach. Materiał tego typu określa się mianem izotropowego materiału liniowo sprężystego, zaś jego model jest powszechnie znany pod nazwą „prawo Hooke'a". Modele materiałowe mechaniki ośrodków ciągłych określa się też często mianem równań konstytutywnych. Nazwa ta pochodzi od greckiego słowa constitue» - ustanawiać. Równanie konstytutywne jest zatem równaniem ustanawiającym materiał. Model materiału liniowo sprężystego i izotropowego znalazł bardzo szerokie zastosowanie i jest do dziś z powodzeniem wykorzystywany do opisu własności mechanicznych takich materiałów jak np. stal. W pierwszej połowie XX w. pojawiły się jednak materiały, których nie dało się opisać za pomocą standardowego modelu liniowego. Były to polimery, coraz szerzej stosowane jako materiały konstrukcyjne w rozmaitych gałęziach przemysłu. Doświadczenia wykazały, że reakcja materiałów polimerowych na wymuszenie (siłowe lub kinematyczne) jest zależna od czasu i co więcej często jest nieliniowa (elastomery). Pojawiła się konieczność zastosowania modeli materiałowych, które pozwoliłyby na dokładny opis własności nowych materiałów, a co za tym idzie także i na efektywne obliczenia inżynierskie w ramach prac projektowych. Pierwszy model tzw. materiału z pamięcią czyli materiału którego odpowiedź na obciążenie jest zależna od czasu, został opracowany jeszcze pod koniec XIX w. przez niemieckiego fizyka Ludwiga Boltzmanna. Stanowił on podwalinę nowego działu MOC, tj. teorii lepkosprężystości. Nazwa nowej teorii oddaje własności opisywanych przez nią materiałów. Materiały te wykazują zarówno własności sprężystego ciała stałego (magazynowanie energii) jak i lepkiej cieczy (rozpraszanie energii). Mówiąc obrazowo: są czymś pomiędzy stalą a wodą. Modele Boltzmanna opisywały materiały, których reakcja zależała liniowo od obciążenia, zaś nieliniowo od czasu. Teorię lepkosprężystości rozwinął w latach 20-tych i 30-tych XX w. włoski matematyk Vito Volterra, co umożliwiło opis materiałów o nieliniowej charakterystyce naprężenie - odkształcenie i równocześnie nieliniowej zmienności reakcji materiału w czasie. Teoria Volterry jest jednak bardzo 186 C. Suchocki skomplikowana i nawet dzisiaj w dobie powszechnego dostępu do komputerów osobistych znajduje ona ograniczone zastosowanie. W latach 40-tych XX w. amerykanie Melvin Mooney i Ronald Rivlin, opracowali niezależnie pierwsze modele nieliniowych materiałów sprężystych poddawanych dużym odkształceniom. W ten sposób powstał nowy, do dziś rozwijany dział MOC, tj. teoria hipersprężystości. Modele hipersprężystości znalazły w pierwszej kolejności zastosowanie do opisu własności mechanicznych gum i kauczuków. Druga połowa XX w. przyniosła gwałtowny rozwój inżynierii biomedycznej, w tym także biomechaniki. Badania własności mechanicznych materiałów biologicznych ujawniły, że prawie zawsze wykazują one mniejszą lub większą nieliniowość charakterystyki naprężenie - odkształcenie oraz zdolność do lepkiego rozpraszania energii (lepkosprężystość). Równolegle rozwinął się przemysł kompozytów polimerowych, coraz szerzej stosowanych jako materiały konstrukcyjne (np. w przemyśle samochodowym i lotniczym), które podobnie jak tkanki biologiczne wykazują nieliniowość materiałową i lepkosprężystość. W tej sytuacji na potrzeby obliczeń inżynierskich sięgnięto zarówno do teorii lepkosprężystości jak i hipersprężystości. Nowoczesna mechanika ośrodków ciągłych zajmuje się opracowywaniem równań konstytutywnych uwzględniających silnie nieliniowe charakterystyki naprężenie - odkształcenie, lepkosprężystość, zależność reakcji materiału od prędkości odkształcenia oraz anizotropię. Wykorzystuje się przy tym opracowaną pod koniec lat 70-tych metodę uwzględniania anizotropii materiału za pomocą tensorów struktury. W pracy zaprezentowano równania konstytutywne, które mogą posłużyć do opisu własności mechanicznych tkanki kostnej. Wskazano na rozwiązanie, które najwierniej oddałoby właściwości tkanki kostnej. Uwzględniono przy tym najnowsze trendy obowiązujące w mechanice ośrodków ciągłych. Omówiono również program komputerowy wykorzystujący metodę aproksymacji Levenberga - Marąuardta, który może służyć do wyznaczania stałych materiałowych dowolnego równania konstytutywnego. 2. MODELOWANIE KONSTYTUTYWNE I METODY WYZNACZANIA STAŁYCH MATERIAŁOWYCH 2.1. Równania konstytutywne tkanki kostnej stosowane do tej pory w symulacjach MES Jak wykazały badania tkanka kostna w zakresie sprężystym wykazuje zależność naprężenie - odkształcenie, którą w przybliżeniu można uznać za liniową (w wypadku kości gąbczastej widoczne jest jednak lekkie zakrzywienie). Tkanka kostna ma własności lepkosprężyste, a więc jest zdolna zarówno do magazynowania energii odkształcenia (jak sprężyna) jak i do lepkiego rozpraszania energii (zamianę jej na ciepło i oddanie do otoczenia). Oprócz tego stwierdzono, że tkanka kostna jest materiałem ortotropowy, tj. wykazującym odmienne własności mechaniczne w trzech prostopadłych kierunkach. Reakcja tkanki kostnej na obciążenie zależy też od prędkości odkształcenia. Wszystkie te właściwości, czynią modelowanie tkanki kostnej zadaniem trudnym. Nie powstało jeszcze jedno równanie konstytutywne modelujące tkankę kostną, które uwzględniałoby jednocześnie nieliniowość, lepkosprężystość, ortotropię i zależność od prędkości odkształcenia. A trzeba pamiętać jeszcze o zdolności tkanki kostnej do nieustannej przebudowy (remodelingu). Opracowane do tej pory równania konstytutywne tkanki kostnej cechują się mniejszymi lub większymi uproszczeniami. Nie znaczy to jednak, że równania te są bezużyteczne. Kluczowe pytanie odnoszące się do symulowania procesów deformacji tkanki kostnej brzmi: jakie cechy tkanki będą w danym wypadku odgrywać najważniejszą rolę? Niejednokrotnie może się okazać, że w danym wypadku uwzględnianie wszystkich własności mechanicznych Modelowanie konstytutywne tkanki kostnej i metody wyznaczania stałych.. 189 2.3. Metody wyznaczania stałych materiałowych Stałe materiałowe zaawansowanych modeli MOC wyznaczyć można posługując się metodami regresji nieliniowej. Algorytmem nieliniowej aproksymacji najodpowiedniejszym do wyznaczania stałych materiałowych jest algorytm Levenberg'a - Marquardt'a (rysunek 1). Na podstawie podanego algorytmu opracowano program napisany w języku MATLAB. Rys. 1. Algorytm aproksymacji metody Levenberg'a - Marquardt'a Optymalizacja modelu przebiegała w czterech etapach: 1) Optymalizacja współczynników ©i (wszystkich trzech w jednej pętli optymalizacji L - M); 2) Optymalizacja pary współczynników coi i t i; 3) Optymalizacja pary współczynników co2 i x i, 4) Optymalizacja pary współczynników 0)3 i T 3. Wyniki aproksymacji krzywej pełzania PVC przedstawiono na rysunku 2. 190 C. Suchocki e(t) [-]6 7.8 7.6 7.4 7.2 7 6.8 6.6 6 4 6 2 , r i 0 2000 4000 6000 8000 10000 12000 14000 f s Rys. 2. Wyniki aproksymacji krzywej pełzania PVC przy użyciu metody Levenberg'a - Marquardt'a W przyszłości planowane jest uzupełnienie programu o kryterium dodatniego znaku wyznaczanych stałych. LITERATURA [1] Dondurur D., Sari C.: A Fortran 77 computer code for damped least-squares inversion of Slingram electromagnetic anomalies over thin tabular. Computers & Geosciences 30 (2004), str. 591 - 5 9 9 . [2] Pioletti D. P.: Viscoelastic Constitutive Law Based on the Time Scale of the Mechanical Phenomena. Mechanics of Biological Tissue Part IV, Springer Berlin Heidelberg, 2006, 399-404. [3] Wilczyński A.: Mechanika polimerów w praktyce konstrukcyjnej. Wyd.l. Warszawa: WNT, 1984. [4] Fung Y. C.: Podstawy mechaniki ciała stałego. Wyd.l . Warszawa: PWN, 1969. CONSTITUTIVE MODELING OF BONE TISSUE AND METHODS FOR DETERMINING THE MATERIAL CONSTANTS Summary. This work begins with a short historical note on the development of continuum mechanics from the beginning of XlXth century to the present day. In the following part of the work different variants of constitutive equation that models bone tissue have been presented. At this point the newest achievements of the continuum mechanics have been taken into account. The rest of the work covers the problem of determining the material constants. Here special attention have been paid to the Levenberg - Marquardt nonlinear regression algorithm.