CC BY
48
УДК 531/534:57+612.7
МЕТОД ОПИСАНИЯ ТЕЧЕНИЯ И ОПРЕДЕЛЕНИЯ РЕОЛОГИЧЕСКИХ КОНСТАНТ ВЯЗКОПЛАСТИЧЕСКИХ БИОМАТЕРИАЛОВ. ЧАСТЬ 2
С.В. Шилько*, С.Л. Гавриленко*, В.Ф. Хиженок*, И.Н. Стакан**, С.П. Саливончик**
‘Институт механики металлополимерных систем им. В.А. Белого НАН Беларуси, 246050, Беларусь, г. Гомель, ул. Кирова, д. 32а, e-mail: [email protected]
‘‘Гомельский кардиологический диспансер, 246012, Беларусь, г. Гомель, Медгородок
Аннотация. Экспериментально изучено гемодинамическое сопротивление мелких вен и артерий с использованием имитатора кровеносного сосуда в виде полимерного капилляра. Показана возможность описания пуазейлевского течения крови при помощи модели вязкопластической среды типа Бингама. Определены реологические константы крови - пластическая вязкость, предельное напряжение сдвига, параметр скоростной чувствительности. Полученные результаты могут быть использованы для анализа динамических явлений в сердечно-сосудистой системе, диагностики заболеваний крови и при создании протезов сосудов.
Ключевые слова: кровь, вязкопластичность, течение Пуазейля, имитатор сосуда, гемореологические константы.
Введение
В первой части [1] данного исследования было показано, что для уточненного описания механических свойств ряда биологических материалов следует привлекать модель вязкопластической среды с нелинейной вязкостью. В частности, разработанная в [1] методика такого описания может быть использована для определения реологических констант крови при течении в артериях и капиллярах и пищевых субстратов при движении по кишечному тракту с образованием квазитвердой области, а также кинетики деформирования мягких тканей организма, демонстрирующих квазипластичность вследствие выраженного порогового эффекта по отношению к нагрузке (например, остеогенеза и резорбции зубных тканей при ортопедической коррекции зубного ряда).
Знание реологических констант крови необходимо для описания гемодинамических явлений [2], диагностики заболеваний [3] и разработки протезов сосудов. Целью настоящей работы было определение реологических констант крови при помощи эксперимента, имитирующего кровоток на масштабном уровне мелких вен и артерий, играющих важную роль в функционировании сердечно-сосудистой системы (табл. 1).
© С.В. Шилько, С.Л. Гавриленко, В.Ф. Хиженок, И.Н. Стакан, С.П. Саливончик, 2003
Таблица 1
Параметры изучаемых сосудов (средние значения)
Вид сосуда Диаметр, см Длина, см Количество
Артерии 0,1 10 600
Вены 0,06 1 1800
Таблица 2
Зависимость вязкости крови ц (с-П) от скорости сдвига при содержании гематокрита Н=44.8 %
Скорость сдвига, у , с-1
0,1 0,3 0,7 1,0 4,0 7,0 10,0 35,0 70,0 100,0
95 60 37 30 20 14 10 8 6 5
Таблица 3
Условия гемодинамических испытаний________________________________
Скорость поршня, мм/мин 60 120 180 240 360
Скорость течения в капилляре, м/с 0,7 1,4 2,1 2,8 4,2
Методика эксперимента
При проведении эксперимента учитывалась выраженная скоростная чувствительность механических характеристик крови как структурированной среды в виде суспензии эритроцитов в плазме [4]. Ранее было показано [2], что в зависимости от скорости сдвига наблюдается феномен упругого деформирования крови (при весьма малых значениях у), затем имеет место вязкопластическое течение и далее -ньютоновское течение с постоянной вязкостью. Трудности определения предела текучести крови при малых скоростях сдвига у ^ 0 обусловлены миграцией эритроцитов от стенок вискозиметра. При высокой скорости сдвига кровь ведет себя как ньютоновская жидкость с постоянной вязкостью, что описывается простейшей зависимостью
т = цу, ц = const. (1)
Эти особенности реологии крови были выявлены на ротационном вискозиметре, реализующем куэттовское течение гомогенной среды, когда ширина зазора между соосными цилиндрами прибора во много раз превышала размер эритроцитов. В соответствии с данными, приведенными в таблице 2, вязкость крови изменяется в широких пределах в зависимости от скорости сдвига.
Таким образом, в основном диапазоне скоростей сдвига уравнение состояния крови должно соответствовать вязкопластическому закону течения, причем этот диапазон расширяется при увеличении содержания гематокрита Н. Можно предположить, что вязкопластические свойства проявляются и в условиях течения Пуазейля через канал цилиндрической формы (капилляр), что более характерно для сосудистой системы.
Рис. 1. Способ реализации пуазейлевского течения крови
Тест на растяжение
Тест на растяжение
-10
-20
-30
и
<й
К
-40-
-50-
60 — 120180240 -360 -
Удлинение (мм) а
Тест на растяжение
20 40 60
Удлинение (мм) б
Тест на растяжение
80
20 40 60
Удлинение (мм) в
Удлинение (мм) г
Рис. 2. Зависимости гемодинамического сопротивления от перемещения: а - испытания при незаполненном цилиндре; б, в - испытания при постоянной скорости; г - испытания при
ступенчатом изменении скорости
В частности, в настоящей работе изучение реологических свойств крови проводилось по схеме, показанной на рис. 1.
В отличие от обычно используемых стандартных вискозиметров [5] в качестве нагружающего и силоизмерительного устройства применяли универсальный испытательный стенд ¡МБТЯОМ 5567 (ик, 2002), обеспечивающий прецизионное определение усилий и перемещений при стабильной скорости нагружения. Нагнетание крови в капилляр производилось при помощи медицинского шприца емкостью 60 мл и диаметром поршня 26,5 мм, установленного в зажимах стенда. Максимальный ход поршня устанавливался равным 90 мм.
Имитатор сосуда представлял собой прямолинейный отрезок тонкой трубки (капилляр) длиной 180 мм из бионейтрального полимера, используемого при изготовлении катетеров для ангиографии, и соединялся со шприцем при помощи стандартного адаптера. Внутренний диаметр имитатора составлял 1,0 мм.
Нужно отметить, что в силу степенной зависимости расхода от данного параметра важно обеспечить высокую точность измерения диаметра капилляра. Кроме того, материал имитатора имеет существенно более высокую жесткость по сравнению с биотканью сосудистой стенки и его дилатация (увеличение диаметра) под действием давления крови незначительна. Для более адекватного моделирования кровотока
необходимо корректировать диаметр, учитывая нелинейное деформационное поведение стенок эластичных кровеносных сосудов.
Для исследования использовалась К-масса (получаемая из цельной крови путем отделения плазмы методом центрифугирования), которая до измерения хранилась в холодильнике при низкой температуре. Измерение сопротивления кровотоку через капилляр для каждой скорости нагнетания производилось с использованием свежей порции крови. Испытания проводились при температуре 37° С для пяти значений скорости поршня и соответствующих значений скорости течения в капилляре, обратно пропорциональных отношению квадратов диаметров поршня и капилляра (табл. 3).
Результаты испытаний приведены на рис. 2 (отрицательный знак регистрируемого усилия соответствует сжатию шприца в стадии нагнетания). С целью устранения систематической погрешности, обусловленной трением поршня о стенки цилиндра, на каждой скорости предварительно проводилось нагружение в условиях смачивания стенок без заполнения цилиндра (рис. 2а). При этом имели место незначительные флуктуации усилия, вызванные, очевидно, изменениями условий трения поршня в цилиндре. Измеренное усилие трения вычиталось из регистрируемых усилий при заполненном цилиндре и тем самым выделялась компонента гемодинамического сопротивления. Для каждой скорости испытание повторялось 3 раза. Было также проведено испытание с изменением скорости нагнетания в середине рабочего хода, что позволяет выявить динамику переходного процесса (рис. 2г).
Зависимости гемодинамического сопротивления от времени для указанных скоростей приведены на рис. 2 б-г. Видно, что происходит достаточно быстрая стабилизации регистрируемого усилия, а с ростом скорости течения гемодинамическое сопротивление увеличивается.
Идентификация вязкопластической модели
Полученные экспериментальные данные были использованы для идентификации гемодинамической модели, основанной на решении задачи о течении вязкопластической среды через капилляр при стационарном перепаде давления для следующих определяющих соотношений [1]:
т = т0 + К
/ \ т
V £0 у
(2)
Н = 0, т < тп
где Зу - девиатор тензора напряжений, т - интенсивность напряжений сдвига, Н -интенсивность скоростей деформации; 4 - девиатор тензора скоростей деформации, К - пластическая вязкость, т0 - предельное напряжение сдвига, т - параметр скоростной чувствительности, ¿0 - характерная скорость деформации (параметр, определяющий размерность).
Расчет проводился путем решения системы нелинейных алгебраических уравнений
2л l
Kn\Api |(n + 1)(n + 2)(n + 3)
і Iap,
V n+1
2l
a -тп
n +1) (n + 2)a2 + 2r02 + 2ar0 (n +1)|= Qi, (3)
где r0 = - радиус квазитвердой области; Ap - перепад давления; /л, n.
Ap
реологические константы;
^0 -
dp
dz
- градиент давления, a, l - радиус и длина капилляра,
г
k = —; Qi - объемный расход среды в единицу времени, i = 1, 2, 3. a
В соответствии с полученными экспериментальными данными были найдены следующие значения реологических констант т0=0,58-10-3 Па, £=0,89-10-2 Па, и= 1,111. Таким образом, кровь является нелинейной вязкопластической средой, хотя и с весьма малым пределом текучести.
Оценка точности расчета была получена при помощи дополнительного эксперимента. В частности, проведение контрольного теста при скорости кровотока в капилляре 2,8 м/с показало, что относительная погрешность расчета т0, п, k не превышает 1%. Кроме того, результат расчета т0 согласуется с известным пределом текучести крови для высоких (более 40%) значений гематокрита.
Заключение
Проведенное расчетно-экспериментальное исследование подтверждает вязкопластическое поведение крови. Разработанная методика может быть использована для определения реологических констант крови в условиях артериального и венозного течения.
Список литературы
1. Gavrilenko S.L., Vasin R.A., Shilko S.V. A method for determining flow and rheological constants of viscoplastic biomaterials. Part 1 // Russian Journal of Biomechanics. 2002. V. 6. № 3. P. 92-98.
2. Регирер С.А. Биофизические основы движения крови по сосудам. // Руководство по кардиологии. Т. 1. Структура и функция сердечно-сосудистой системы в норме и при патологии / Под ред. Е.И. Чазова. М.: Медицина, 1982. С. 189-201.
3. Baev V.M. Analysis of dependence between the character of alcohol intake by adults and rheologic properties of the whole blood // Russian Journal of Biomechanics. 2001. V. 5. № 1. P. 95-100.
4. Шафранова Е.И., Снегирева Н.С. Кровь как структурированная среда. Особый структурный уровень пространственной организации крови - агрегаты эритроцитов в норме и патологии // Мех. композит. матер. и конструкций. 2001. Т. 6. № 3. С 359-371.
5. Малкин А.Я., ЧалыхА.Е. Диффузия и вязкость полимеров. Методы измерения. М.: Химия, 1979.
A METHOD FOR DEFINING FLOW AND RHEOLOGICAL CONSTANTS OF VISCOPLASTIC BIOMATERIALS. PART 2
S.V. Shilko, S.L. Gavrilenko, V.F. Khizhenok, I.N. Stakan, S.P. Salivonchik
(Gomel, Belarus)
Hemodynamic resistance of small vessels and arteries has been studied using a blood artery simulator in the form of a polymer capillary. It has been proved that Poiseuille’s flow can be described by a viscoplastic medium model of Bingham’s type. Rheological constants of blood, including plastic viscosity, limiting shear stress and velocity sensitivity have been estimated. The results derived are used for analyzing dynamic phenomena in the cardiovascular system, diagnostics of blood diseases and for designing vessel prostheses.
Keywords: blood, viscoplasticity, Poisseuille’s flow, vessel simulator, hemorheological constants.
Получено 24 апреля 2003
