Разработка экспериментальных методов исследования капиллярного движения жидкости применительно к задачам гемодинамики тема автореферата и диссертации по физике, 01.04.01 ВАК РФ

На правах рукописи

КОЗЛОВ Денис Юрьевич

РАЗРАБОТКА ЭКСПЕРИМЕНТАЛЬНЫХ МЕТОДОВ ИССЛЕДОВАНИЯ КАПИЛЛЯРНОГО ДВИЖЕНИЯ ЖИДКОСТИ ПРИМЕНИТЕЛЬНО К ЗАДАЧАМ ГЕМОДИНАМИКИ

01 04 01 - приборы и методы экспериментальной физики

Автореферат диссертации на соискание ученой степени кандидата физико-математических наук

Барнаул-2007

ооз

003174206

Работа выполнена в Алтайском государственном медицинском университете

Научный руководитель

доктор технических наук, профессор Волков Валерий Иванович

Официальные оппоненты

доктор физико-математических наук, профессор Шайдук Александр Михайлович

доктор физико-математических наук, старший научный сотрудник

Прокопьев Владимир Егорович

Ведущая организация

Алтайский государственный технический университет

им И И Ползунова

Защита состоится 2 ноября 2007г в 14— часов на заседании диссертационного совета Д 212 005 03 в Алтайском государственном университете по адресу 656049, г Барнаул, пр Ленина, 61

С диссертацией можно ознакомиться в библиотеке Алтайского государственного университета

Автореферат разослан «23» 2007 г

диссертационного совета С

Ученый секретарь

Рудер Д Д

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы

Капиллярно-гидродинамические течения жидкостей широко распространены в природе и технике, однако особую значимость капиллярные явления приобретают при исследовании задач гемодинамики живых организмов Характер кровотока определяет функциональность и жизнеспособность организма, и поэтому исследование капиллярно-гидродинамических течений имеет, помимо многих других областей знаний, важное прикладное значение для медицины

Несмотря на практическую значимость движения жидкостей в капиллярах, этот процесс недостаточно изучен Отчасти это связано с тем, что межфазные эффекты весьма чувствительны к примесям и физическому состоянию поверхности, различным неровностям и шероховатостям, что не всегда можно учесть в методике эксперимента Кроме того, течение биологических жидкостей имеет ряд существенных особенностей, и, если говорить о гемодинамике, то важен состав крови, диаметр и поверхность кровеносного сосуда Более того, по упругим стенкам крупных кровеносных сосудов от сердца распространяется пульсовая волна, что является еще одним требующим учета фактором Недостаточно полно развитые приборно-экспериментальная база и методы исследования подобных процессов позволяют говорить об актуальности предлагаемой работы

В настоящей работе описывается разработка математических моделей, а на их базе - экспериментальных методик изучения нестационарных течений жидкости в капиллярах С их помощью исследуется скачкообразный процесс заполнения пустого капилляра, изучаются пульсации давления в пульсовой волне, распространяющейся по стенке кровеносного сосуда На основе выработанной математической методики обработки данных сфигмографического обследования предлагается способ оценки упруго-эластических свойств артерий и ранней экспресс-диагностики атеросклероза, для чего создана измерительная приставка к компьютеру с программно-вычислительным комплексом обработки экспериментальных результатов

Приборы и методы исследования, применяемые в настоящей работе, актуальны в различных областях науки, особенно при изучении течений в закрытых и щелевых каналах, при исследовании фильтрации жидкостей в пористых структурах, в медицине, при исследовании состояния сердечно-сосудистой системы

Цель работы

Основной целью работы являлась разработка приборной и методической базы экспериментального исследования капиллярных течений различных жидкостей в нестационарном режиме При этом ставились следующие задачи

1 Создать экспериментальную базу и разработать методику исследования динамики заполнения капилляра, и провести экспериментальные исследования характера капиллярного течения

2 Разработать математическую модель процесса натекания жидкости в капилляр и проверить ее экспериментально

3 Применительно к задачам гемодинамики разработать физическую модель квазипериодического процесса распространения пульсовой волны в организме на базе параметрических дифференциальных уравнений Хилла и Матье

4 Разработать статистический критерий оценки степени хаотичности фазовых траекторий пульсовой волны применительно к диагностике атеросклероза, и на этой базе создать измерительную установку с программно-вычислительным комплексом обработки экспериментальных результатов

Практическая значимость

1 Разработана новая методика измерения вязкости жидкости и изготовлена измерительная установка на ее основе

2 Разработан статистический критерий оценки степени хаотичности фазовых траекторий пульсовой волны применительно к диагностике атеросклероза с математической обработкой данных сфигмографии, и созданы измерительная установка, сопряженная с компьютером, и программно-вычислительный комплекс обработки экспериментальных результатов

Научная новизна работы

1 На основе предложенной математической модели получено приближенное уравнение для процесса заполнения жидкостью капилляра, имеющее аналитическое решение, удовлетворительно согласующееся с численным решением точного нелинейного уравнения и экспериментальными данными

2 С использованием созданной экспериментальной методики обнаружен скачкообразный характер капиллярного движения жидко-

сти, впервые получены количественные результаты и предложено качественное описание этого явления

3 На основании применения разработанной экспериментальной методики сделан вывод о влиянии адгезионных свойств материала поверхностей растекания на капиллярное движение жидкостей, и предложена поправка для учета этого влияния

4 Для задач гемодинамики разработана физическая модель квазипериодических процессов организма на базе параметрических дифференциальных уравнений Хилла и Матье

5 Предложен новый способ оценки упруго-эластических свойств артерий и диагностики атеросклероза на основе анализа фазовых траекторий полученного сфигмографического сигнала

Апробация работы

Основные положения диссертационной работы докладывались и обсуждались на научных и научно-практических конференциях

1 Международной научной конференции «Устойчивость и турбулентность течений гомогенных и гетерогенных жидкостей» (г Новосибирск, 26-28 апреля 2004 г )

2 Международной научно-практической конференции «Измерительные информационные технологии и приборы в охране здоровья МЕТРОМЕД - 2007» (г Санкт-Петербург, 17-19 апреля 2007 г)

Публикации

По результатам диссертационной работы опубликовано 10 печатных работ

Структура и объем диссертации

Диссертационная работа состоит из введения, трех глав, заключения и списка литературы Она изложена на 108 страницах, включая 53 рисунка, 3 таблицы, список литературы из 102 названий

Основные положения, представляемые к защите

1 Экспериментальная методика для измерения вязкости жидкости, основанная на учете адгезионных свойств поверхности капилляров и коэффициентов поверхностного натяжения эталонной и исследуемой жидкостей

2 Новая физическая модель процесса натекания жидкости в капилляр, позволяющая адекватно описать динамику капиллярного течения

3 Экспериментальная методика обнаружения скачкообразного характера движения мениска жидкости в капилляре с количественными результатами описания этого явления

4 Физическая модель квазипериодической пульсовой волны, распространяющейся в организме, на базе параметрических дифференциальных уравнений Хилла и Матье

5 Программно-аппаратный комплекс обработки экспериментальных результатов, основанный на предложенном статистическом критерии оценки степени хаотичности фазовых траекторий пульсовой волны применительно к диагностике атеросклероза

КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ

Во введении обосновывается актуальность выбранной темы, обсуждаются цели и задачи работы, а также основные защищаемые положения

В первой главе содержится литературный обзор методик экспериментальных и теоретических исследований поверхностных свойств жидкостей на границах раздела гетерогенных сред, капиллярных течений и тонких пленок Также рассматривается литература, посвященная методам и приборам изучения квазипериодических процессов организма, а также математическим моделям, лежащим в основе таких изысканий Течение биологических жидкостей имеет ряд существенных особенностей, и если говорить о гемодинамике, то важен состав крови, диаметр и поверхность кровеносного сосуда Более того, по упругим стенкам крупных кровеносных сосудов от сердца распространяется пульсовая волна, что является еще одним требующим учета фактором Поэтому для мелких кровеносных сосудов определяющими исследование параметрами являются вязкость и поверхностное натяжение крови, а для крупных артерий - пульсации стенки сосуда, имеющие квазипериодический характер

Вышеуказанным темам посвящено большое количество работ, в которых описываются контактные и бесконтактные методы экспериментального исследования поверхностного натяжения, напряжения начального сдвига и вязкости Из приведенных работ можно заметить, что все эти методы основаны на экспериментальном исследовании жидкостей при их движении, что приводит к необходимости учета сил трения и поведения жидкости под действием этих сил Обзор измерительных устройств вязкости, поверхностного натяжения и напряжения начального сдвига показал, что до настоящего времени существовало сравнительно мало разработок, которые бы позволили просто, а глав-

ное оперативно определять факторы, влияющие на капиллярные течения Расчет коэффициента вязкости в ряде методов опирается на табличные значения физических характеристик жидкости (коэффициент поверхностного натяжения, плотность), которые меняются в зависимости от влажности и температуры Кроме того, свойства биологических жидкостей являются индивидуальными для каждого организма, поэтому для них известные методы измерения вязкости могут приводить к существенным погрешностям

Математическое описание различных квазипериодических процессов организма, например, распространения пульсовой волны по сосудам кровеносной системы, в настоящее время не достаточно разработано Существует ряд работ, в которых предлагаются попытки моделирования таких процессов при помощи различных математических методов, в частности, с использованием параметрических дифференциальных уравнений Существующие методики диагностики состояния организма связаны с субъективностью врачебной оценки Поэтому актуальна задача разработки способов диагностики, в которых критерием дифференциации норма/патология выступали бы объективные числовые показатели Главным образом, в разрабатываемых на этом направлении методах используются статистические оценки и численный Фурье-анализ сигнала, с последующим исследованием его спектральных характеристик Однако на этом пути из-за нестрогой периодичности сердечного ритма есть ряд технических и вычислительных трудностей Вторая глава посвящена разработке экспериментальной базы и соответствующих методик для изучения течения жидкостей в тонких капиллярах Эти тонкие трубки выбраны в качестве физической модели мелких кровеносных сосудов, а для них важнейшими параметрами, определяющими движение, являются вязкость и поверхностное натяжение крови В рамках разработки математической модели для процесса заполнения жидкостью пустого капилляра, расположенного под углом Р к горизонтали (рис 1) записано нелинейное уравнение

здесь а — коэффициент поверхностного натяжения, К - радиус капилляра, г - длина растекания жидкости в капилляре, - длина погруженного в жидкость участка, V - коэффициент кинематической вязкости жидкости, g — ускорение свободного падения, 0 — параметр, определяемый из усреднения градиента скорости по сечению канала (в частности, для установившегося пуазейлевского течения 0 = 8/3)

Ж

2у6 (к

Для горизонтально расположенного капилляра уравнение упрощается, его удается решить аналитически. Таким образом, получена линейная зависимость квадрата смещения мениска жидкости от времени.

Для записи капиллярного движения жидкости была разработана следующая экспериментальная методика; производилась цифровая видеосъемка с разрешением 1024x768 пикселей и частотой 15 кадров в секунду. Затем происходила обработка видеофайла на компьютере для извлечения экспериментальных данных. Сумма случайной и систематической экспериментальной погрешности не превосходила 10%. В экспериментах со спиртом и трансформаторным маслом подтверждена ожидавшаяся закономерность (рис. 2), кроме того, получено, что при заполнении капилляров профиль скорости течения при малых числах Вебера может существенно отличаться от пуазейлевского, и величина этого отличия зависит от адгезионных свойств материала. Так, для двух разных стекол (в первом содержится 14.5% №20 и 1% К20, а во втором сорте стекла - 9.4% 1Ча20 и 7.1% К20) получены различия в динамике растекания, что иллюстрируется на рис. 2. Скорость течения определяется тангенсом угла наклона экспериментальной прямой

. Дополнительные эксперименты по вертикальному капиллярному поднятию показали, что высота поднятия в капиллярах из стекла большим содержанием К20 оказывается в 1.2 раза выше, чем в капиллярах из стекла первого сорта. Поэтому, сделан вывод о влиянии адгезионных свойств материала поверхностей растекания на движение

8

"П ---"I— т--1

жидкостей и оценена поправка, позволяющая учитывать это обстоятельство.

Рис. 2. Зависимость квадрата координаты мениска спирта от времени. Сплошная линия соответствует капиллярам из стекла первого сорта с диаметрами 1 - с!=9.4- 10"4 м, 3 - сНЗ.4- 10"4м; пунктир с точкой - капиллярам из стекла второго сорта с диаметром 2 - <1=5.9-10~4 м, 4 —

<1=з. ноЛ*.

На основании полученных закономерностей и отработанных экспериментальных методик была создана новая установка для измерения вязкости биологической жидкости, использующая два горизонтальных капилляра. Вязкость исследуемой жидкости устанавливается по характеристикам ее движения относительно эталонной. Методика измерений состоит в следующем. Предварительно определяется высота капиллярного поднятия эталонной жидкости в поле тяжести для обоих капилляров, затем высота поднятия исследуемой жидкости в поле тяжести, таким образом, отпадает необходимость знать коэффициент поверхностного натяжения обеих жидкостей. Далее измеряются длины одновременного растекания эталонной и исследуемой жидкостей по горизонтальным капиллярам и определяют кинематическую вязкость из соотношения:

V = (XV..

Л Я2

¿1 I2

где уэ, V - кинематические вязкости эталонной и исследуемои жидкостей, Ьэ, И - высоты капиллярного поднятия эталонной и исследуемой жидкостей в поле тяжести, Иэ, Я - радиусы капилляров по которым текут соответственно эталонная и исследуемая жидкость, 1.,, Ь - длины растекания эталонной и исследуемой жидкостей, а - коэффициент, учитывающий возможную разницу адгезионных свойств стекол, и определяющийся с помощью эталонной жидкости. Т.о., по сравнению с вискозиметром Гесса, увеличивается диапазон измерения относительной кинематической вязкости, а учет сил адгезии поверхности капилляра позволяет определить это значение более достоверно.

Рис. 3. Безразмерная длина растекания в зависимости от безразмерного времени в стеклянном капилляре диаметром 9.4-10 4 м и экспериментальные точки с учетом погрешностей.

Затем рассматривался процесс затекання жидкости в капилляр, расположенный под некоторым углом. Аналитическое решение нелинейного уравнения (I) найти не удается, поэтому была сделана попытка изменения исходного уравнения, с тем, чтобы получить приемлемое аналитическое решение, которое было бы близко как к численному решению, так и соответствовало экспериментальным данным. Для чис-

ленного решения дифференциального уравнения (1) использовался метод Булирша-Штёра. Модификация исходного уравнения приводит к безразмерному уравнению

(12м>

. 2 п 2v9 х с1и>

■ = 1 - XV БШ р-----

(2)

Л* Я'

где и' имеет смысл безразмерной разницы между квадратами длин заполненной и погруженной частей капилляра (рис. 1).

Решение (2) известно и описывает апериодические затухающие колебания. При некотором выборе параметра 6 модифицированного уравнения удается получить в качестве его решения кривую, которая очень близка к численному решению уравнения (1) и не противоречит экспериментальным данным (рис. 3). В ходе работы изучалась зависимость динамики движения этилового спирта в капилляре от глубины его погружения в широкий сосуд. Результаты представлены на рис. 4. Очевидно, что скорость заполнения капилляра падает с ростом глубины погружения.

Рис. 4. Зависимость длины растекания спирта от времени; сверху - по-

гружение на 9.1 10° м, в центре Капилляр с диаметром 3.4-10'* м.

на 4.5-10"2 м, внизу - на 0.2-КГ

1 1

Т

Из сравнения точного решения и экспериментальных данных следует, что профиль скорости при затекании жидкости в капилляр существенно отличается от пуазейлевского Так, значение максимальной скорости превышает расходную скорость более чем в 6 раз (точнее в 68 раз) С другой стороны, единственной силой толкающей жидкость вдоль капилляра является сила адгезии между жидкостью и стенкой капилляра, вследствие короткодействия эта сила распределена вблизи поверхности Отсюда можно предположить, что скорость жидкости вблизи стенки цилиндрического капилляра имеет максимальное значение, существенно превышающее и среднее значение скорости, и скорость в центре капилляра

Таким образом, разработана физическая модель капиллярного на-текания, в рамках которой предложено приближенное уравнение, имеющее аналитическое решение, позволяющее удовлетворительно аппроксимировать точное численное решение нелинейного уравнения процесса заполнения жидкостью капилляра и получить согласие с экспериментальными данными

Однако при детальном исследовании динамики затекания жидкости в вертикальный стеклянный капилляр было обнаружено, что заполнение капилляра происходит далеко не монотонным образом Так, во время заполнения капилляра скорость движения мениска жидкости несколько раз то увеличивалась, то уменьшалась Амплитуда и количество пульсаций по мере подъема мениска жидкости менялись, и максимальная высота поднятия жидкости также могла отличаться в отдельных экспериментах Вероятно, причина этого обусловлена гистерезисом смачивания Отметим, что в экспериментах с капиллярами, расположенными под различными углами к горизонту, движение также имело скачкообразный характер Особенно отчетливо неравномерность движения выявилась с помощью следующей методики Производилась видеосъемка теневого изображения капилляра с жидкостью В ходе экспериментов использовались стеклянные цилиндрические капилляры диаметрами от 3 10'4 до 10~3 м и капиллярная щель, образованная двумя близко расположенными параллельными стеклянными пластинами, расстояние между которыми можно было варьировать при помощи микровинта Для визуализации процесса заполнения капилляра была собрана оптическая система, проецирующая на экран теневое изображение капиллярного движения, что увеличивало капилляр в 35 раз Затем происходила обработка видеофайла на компьютере для извлечения экспериментальных данных Фиксировалось положение мениска жидкости через каждые три кадра (0 2 с) и, по изменению координаты, вы-

числялась скорость движения мениска. Рис. 5 иллюстрирует случайные скачкообразные изменения скорости.

Одно из возможных объяснений обнаруженного эффекта неравномерности капиллярного движения может заключаться в следующем. Как только жидкость соприкасается с твердой вертикальной стенкой, смачиваемой жидкостью, на ней образуется тонкий адгезионный слой жидкости, который устремляется вдоль свободной поверхности капилляра. При этом поверхность жидкости натягивается подобно тетиве лука, и, когда сила, обусловленная поверхностным натяжением, сравнивается с силой начального напряжения сдвига в жидкости, происходит срыв основной массы относительно жидкостной пленки вдоль стенки, что фиксируется в виде скачка столбика жидкости.

1 1 ■

■ 1 I,

0123456789 10

Ъс

Рис. 5. Скорость движения мениска воды по капилляру диаметром 0.34 мм.

Другими словами, по мерс капиллярного заполнения вдоль стенки капилляра скользит вверх, значительно опережая мениск, тонкая пленка жидкости - так называемая прекурсиопная пленка (в обзоре де Жена есть указания на ряд экспериментов по обнаружению этой пленки). А так как пленка тонкая, то она достаточно быстро испаряется, поэтому влияние ее на процесс заполнения капилляра будет тем меньше, чем меньше скорость движения жидкости. Таким образом, эффект нерав-

номерности капиллярного заполнения будет менее выражен в конце капиллярного движения при малой скорости жидкости, особенно для вязких жидкостей, скорость капиллярного заполнения для которых невелика

Для обнаружения неразличимой визуально прекурсионной пленки жидкости - «предвестника» первого скачка столбика жидкости, была разработана следующая методика измерялось электрическое сопротивление между жидкостью в кювете и стенкой капилляра, на которую крепился электрод из тонкой медной проволоки Было обнаружено, что после добавления жидкости в кювету сопротивление в цепи менялось в течение нескольких секунд, а затем восстанавливалось первоначальное значение При этом на теневой картинке с большим увеличением удалось обнаружить, что изменение сопротивления возникало на значительном удалении мениска от электрического контакта на стенке капилляра Это можно интерпретировать как появление прекурсионной пленки, а затем, поскольку капиллярное поднятие не продолжалось, ее испарение Так, в экспериментах с 96%м этиловым спиртом мениск в среднем отстоял от электрического контакта на 0,01 ±0,002 м При капиллярном поднятии дистиллированной воды также была обнаружена прекурсионная пленка, и ее длина составила 0,026±0,002 м

Третья глава посвящена разработке методик исследования гемодинамики крупных артерий, в частности, колебательных процессов жизнедеятельности организма, имеющих квазипериодический характер, с помощью дифференциальных параметрических уравнений Представляется важным создание физических моделей таких явлений, и построение на их основе математических моделей функционирования организма Тогда появляется универсальная возможность разработки различных методов компьютерной диагностики состояния организма Для крупных кровеносных сосудов требующими учета факторами являются характеристики пульсовой волны, распространяющейся от сердца по упругим стенкам артерий Для описания этих факторов в работе предлагается использовать дифференциальное уравнение Ма-тье, описывающее параметрические колебания

ф + (а + Ь 81П( ЮГ))ф = 0; (5)

в качестве способа моделирования квазипериодических колебательных процессов в организме параметрического движения крови и соответствующих ему пульсаций давления, пульсовой волны Метод диагностики состояния сердечно-сосудистой системы по параметрам пульсовой волны называется сфигмографией

Сначала разрабатывается физическая модель параметрических колебаний пульсовой волны, затем строится математическая модель, в

которой предлагается алгоритм построения уравнения Матье по более сложному, имеющему квазипериодическое решение уравнению Хилла, в предположении, что уравнение Хилла ближе к описанию сердечной деятельности. Показано, что в этом случае удается аппроксимировать решение уравнения Хилла решением уравнения Матье с удовлетворительной точностью. Обнаружено, что близость фазовой траектории, получаемой при решении какого-либо уравнения, к фазовой траектории уравнения Матье, может служить основанием для попытки аппроксимации решения такого уравнения предложенным методом.

^ : ^ Рис. 6. Фрагмент сфигмограммы и наложенное решение уравнения Матье (темная линия). Временная и амплитудная разница не превосходят 10%.

Далее, рассматриваются колебания пульсовой волны, идущей по кровеносной системе человека в результате сердечных сокращений. На рис. 6 представлены фрагмент сфигмограммы, которая представляют собой зависимости давления от времени в пульсовой волне, полученной на артерии стопы пациента. Получено качественное сходство сфигмограмм с численными решениями уравнения Матье (5) при специально заданной форме возбуждающего сигнала.

-1.2 -Об -0 4 О С 4 О в

•эо -го -ю

(вив

Рис. 7. Фазовые траектории уравнения Матье. а —устойчивое решение, Ь-биения вблизи предела устойчивости по параметру а, с расходящееся решение.

Далее исследуется устойчивость относительно параметров о и Л полученных решений уравнения Матье (5), при помощи их фазовых траекторий (рис. 7). С увеличением какого-либо из коэффициентов а или Ь происходит переход к неустойчивости через биения. Математи-

чески это проявляется как возрастание площади, охватываемой фазовой траекторией Такое возрастание можно трактовать как переход от нормального состояния сердечно-сосудистой системы к критическому Известно, что увеличение скорости пульсовой волны при таком переходе составляет 1 5-2 5 раза Пределы устойчивости решения уравнения Матье (5), полученные для рассматриваемой задачи по параметру а - 1 93 раза, а по Ъ - 1 42 раза Состояние сосудов в предлагаемой модели ассоциировалось именно с параметром а, относительное изменение которого как раз находится в указанном выше интервале Отметим возможность использования фазовой траектории для диагностики, как это было продемонстрировано при рассмотрении устойчивости решения уравнения Матье Вероятно, что предложенная схема может быть распространена на различные методы исследования периодических процессов в организме, которые, в большинстве своем, относятся к классу параметрических колебаний

Атеросклероз является одной из ведущих причин смерти и ин-валидизации взрослого населения Нарушение упруго-эластических свойств артериальной стенки - одно из важнейших проявлений атеросклероза Цель этой части работы состояла в выработке критерия экспресс-анализа атеросклеротического поражения сонных артерий по данным сфигмографии В предлагаемом исследовании для снятия сфигмограмм использовался пьезоэлемент РМЬ-12Т-9 2А1-50 с резонансной частотой 1500 Гц, сигнал с которого после схемной обработки и амплитудной модуляции подавался на линейный вход звуковой платы компьютера Частота несущего сигнала составляла в экспериментах 500 Гц, а частота дискретизации записи информации составляла 10000 Гц Таким образом, на входе звуковой платы компьютера имелся сигнал, пропорциональный давлению на пьезоэлектрическом датчике Значение производной от давления вычислялось компьютерной программой Затем на компьютере выполнялись демодуляция и ограничение полосы спектра сигнала сверху частотой среза 70 Гц с помощью фильтра Баттерворта Поскольку пульсовые колебания более полно описываются как значением давления, так и скоростью его изменения, после компьютерной демодуляции строились и анализировались фазовые траектории процесса в координатах (Р, Р ) Здесь и далее, обозначение Р соответствует безразмерному значению приращения давления, отнесенному к максимальному значению этого приращения Обозначение Р имеет аналогичный смысл отношение производной приращения давления к ее максимальному значению

С помощью сканера MEDISON 8000ЕХ (диагностический центр АГМУ, В П Куликов, С В Засорин) в двухмерном ультразвуковом ре-

жиме оценивали состояние стенки сонных артерий в области наложения датчика давления. В качестве критериев атеросклероза использовали толщину комплекса интима-медиа более 1 мм и/или наличие ате-росклеротических бляшек. Выявленные таким образом больные атеросклерозом подвергались сфигмографическому обследованию, данные которого обрабатывались при помощи описываемой ниже методики. На рис. 8 и рис. 9 слева приведены характерные фазовые траектории в зависимости от времени сонных артерий, соответственно здорового человека и пациента с атеросклерозом. Ультразвуковое сканирование выявило значительное (1.3 мм при верхней границе нормы в 1 мм) утолщение стенки артерии и наличие множественных мелких атером.

Рис. 9. Пациент с атеросклерозом. Слева фазовая траектория сигнала от правой общей сонной артерии в зависимости от времени. Справа сечение трехмерной поверхности количества фазовых состояний при уровне числа состояний я = 15

Рис. 8. Здоровый. Слева фазовая траектория сигнала от левой общей сонной артерии в зависимости от времени. Справа сечение трехмерной поверхности количества фазовых состояний при уровне числа состояний п = 15 •

Сравнивая рис. 8 и рис. 9 можно заметить, что фазовая траектория сигнала, полученного от пациента с атеросклеротически измененными сосудами, более хаотична. В ней практически отсутствовала цикличность, которая характерна для здоровых. Как было установлено

при исследовании стенок сосудов в двухмерном ультразвуковом режиме, этот пациент имел множество мелких атером на правой общей сонной артерии и протяженные атеросклеротические бляшки на левой общей сонной артерии Вследствие атеросклеротических отложений стенки сосудов отчасти утрачивали свои упруго-эластические свойства, становясь более жесткими, что повышало скорость распространения пульсовой волны Кроме того, отложения на границе сосуда турбу-лизировали кровоток, что проявлялось хаотическими искажениями формы сфигмографического сигнала (рис 9 слева)

Для выработки численного критерия степени хаотичности фазовых траекторий была выполнена следующая процедура Вся фазовая площадь (Р, Р') разбивалась на 400 (20x20) одинаковых прямоугольников, и подсчитывалось количество попавших в каждый прямоугольник состояний п, образующих фазовые траектории Далее строилась трехмерная поверхность, у которой две координаты М|(Р) и М2(Р') соответствовали фазовым координатам (Р, Р'), а значение третьей координаты п равнялось количеству фазовых состояний с определенными величинами давления и производной от давления

Большая хаотичность фазовой траектории у больных (рис 9 слева), проявлялась на трехмерной поверхности приближением распределения фазовых состояний к нормальному распределению Максимальное значение и величину дисперсии этого распределения можно использовать для оценки состояния пациента Для здоровых людей это распределение фазовых состояний характеризовалось сложными фигурами, в которых встречались как спирали, так и несколько мало отличающихся по величине вершин Следовательно, распределения для пациентов с атеросклерозом и здоровых людей отличались степенью заполнения прямоугольника, описанного вокруг фигуры, полученной в определенном сечении этой трехмерной поверхности Например, на рис 8-9 справа представлены сечения таких поверхностей при п — 15 Степень заполнения прямоугольника, полученного в сечении трехмерной поверхности при заданном значении уровня числа состояний п, определится следующим образом С{п) = N/где д<- количество

занятых ячеек, /V - количество ячеек, необходимых для того, чтобы

заполнить вышеозначенный прямоугольник целиком, т е все возможные состояния в рамках области, ограниченной прямоугольником Для пациента с атеросклерозом в большинстве сечений была гораздо большая степень заполнения ограничивающего прямоугольника, чем для здорового человека На основании этого предложен численный критерий оценки состояния сосудов

J(0,5 -G{n))dn + 7(0,5 - G(n)\dn S = -2-.

2

где ит1ч - точка, в которой G(«) достигает 1

Для правой общей сонной артерии пациента с утолщенной стенкой и множественными мелкими атеромами S — 0 В тоже время, для здорового - S — 24 Дальнейшие исследования показали, что для пациентов с атеросклеротически измененными сосудами указанное число либо равнялось нулю, либо оказывалось гораздо меньше, чем у здоровых Так, у здоровых самое меньшее полученное значение S составило 22, а иногда выражалось и трехзначным числом Для больных самое высокое значение не превышало 13, а в 80 % случаев - принимало нулевое значение

Расчетный метод тестировался для решений уравнения Матье с разной устойчивостью Так, для устойчивого решения (рис 7а) S = 125, далее с уменьшением устойчивости и, соответственно, с повышением степени хаотичности, значение числового критерия ожидаемо быстро снижалось Это позволяет говорить как о корректности введенного критерия для соответствующей степени устойчивости уравнения Матье, так и о том, что решение уравнения Матье приближается к описанию сфигмографического сигнала

Итак, отметим, что определен сфигмографический критерий оценки степени атеросклеротического поражения сонных артерий, основанный на функции распределения фазовых состояний сигнала, устойчивый к небольшим ошибкам позиционирования датчика давления Предложенная методика исследования фазовой траектории сигнала может быть использована и для некоторых других методов изучения периодических процессов организма Так, можно с определенностью утверждать, что введенный критерий будет работать для реографии (регистрации изменения импеданса биологической ткани вследствие сердечных сокращений), поскольку морфология реографического сигнала качественно очень близка к морфологии сфигмограммы

Основные выводы и результаты работы

1 Разработана новая методика определения вязкости жидкости, основанная на учете адгезионных свойств поверхности капилляров и коэффициентов поверхностного натяжения эталонной и исследуемой жидкостей, и создана установка для измерения вязкости

2 Предложена новая физическая модель, в рамках которой записано нелинейное уравнение для процесса натекания жидкости в капилляр и получено приближенное уравнение, имеющее аналитическое решение, удовлетворительно согласующееся с численным решением и экспериментальными данными

3 Разработана экспериментальная методика и создан экспериментальный стенд для обнаружения скачкообразного характера капиллярного движения жидкости, впервые получены количественные характеристики этого движения и предложено его качественное описание

4 Разработана физическая модель квазипериодического процесса распространения пульсовой волны в организме на базе параметрических дифференциальных уравнений Хилла и Матье применительно к задачам гемодинамики

5 Разработан статистический критерий оценки степени хаотичности фазовых траекторий пульсовой волны применительно к диагностике атеросклероза с математической обработкой данных сфигмографии, и созданы измерительная установка, сопряженная с компьютером, и программно-вычислительный комплекс обработки экспериментальных результатов

Основное содержание диссертации опубликовано в следующих работах

1 Волков В И , Козлов Д Ю , Шеломанов А И Близость колебаний параметрических систем // Устойчивость и турбулентность течений гомогенных и гетерогенных жидкостей Тезисы доклада Международной конференции Вып 9 / Под ред В В Козлова Новосибирск, 2004, с 48-49

2 Козлов Д Ю Моделирование периодических процессов с помощью уравнений Матье и Хилла // Физика, радиофизика - Новое поколение в науке, Межвузовский сборник научных статей молодых ученых, аспирантов, студентов Вып 4/ под ред В В Полякова, Барнаул изд-во АлтГУ, 2004 - с 67-72

3 Волков В И , Козлов Д Ю Математические модели сфигмографии // Сборник научных работ «Актуальные проблемы фармации», Вып IV, Барнаул Изд-во АГМУ, 2007, с 55-62

4 Волков В И , Козлов Д Ю , Кирколуп Е Р Исследование пульсаций при капиллярном натекании // Сборник научных работ «Актуальные проблемы фармации», Вып IV, Барнаул Изд-во АГМУ, 2007, с 63-68

5 Волков В И , Козлов Д Ю , Кирколуп Е Р Исследование динамики движения жидкости по капилляру // Известия АГУ, № 1 (53), - Барнаул изд-во АлтГУ, 2007 - с 100-104

6 Волков В И , Козлов Д Ю , Оценка капиллярного поднятия // Известия АГУ, № 1 (53), - Барнаул изд-во АлтГУ, 2007 -с 105 - 107

7 Волков В И, Козлов Д Ю и Кирколуп Е Р Динамика капиллярного натекания // Известия вузов Физика, №5, Томск - 2007 - с 52-57

8 Волков В И , Козлов Д Ю . Кирколуп Е Р Исследование затекания жидкости в тонкие капилляры при малых числах Вебера // Журнал технической физики, 2007, т 77, Вып 7, с 24-27

9 Волков В И , Козлов Д Ю , Кирколуп Е Р К вопросу о затекании жидкости в тонкие капилляры // Труды международной научной конференции "Измерительные и информационные технологии в охране здоровья МЕТРОМЕД-2007" СПБ Изд-во Политехи ун-та, 2007, с 56-57

10 Волков В И , Козлов Д Ю Использование параметрических уравнений для пульсовой диагностики // Труды международной научной конференции "Измерительные и информационные технологии в охране здоровья МЕТРОМЕД-2007" СПБ Изд-во Политехи ун-та, 2007, с 58-59

Подписано в печать 26 09 2007

Объем 1 пл

Бесплатно

Печать ризографическая Бумага офсетная Тираж 100 экз

Отпечатано РА «ПАРАГРАФ» г Барнаул, пр Ленина 40 каб 334, тел (385-2)366-143 Лицензия на полиграфическую деятельность ПД 12-061 от 04 01 2002 г

ВВЕДЕНИЕ.-3

ГЛАВА 1. МАТЕМАТИЧЕСКИЕ МОДЕЛИ И МЕТОДЫ ИССЛЕДОВАНИЯ КАПИЛЛЯРНЫХ ТЕЧЕНИЙ ЖИДКОСТЕЙ.- 7

1Л. Оценка физических свойств жидкостей.- 7

1.2. Анализ и обзор технических решений по измерению поверхностного натяжения жидкостей и краевого угла.- 9

1.3. Обзор устройств по измерению напряжения начального сдвига.- 14

1.4. Устройства для измерения вязкости жидкости.-181.5. Обзор методов и приборов изучения квазипериодических процессов организма.- 20

ГЛАВА 2. СРЕДСТВА ИССЛЕДОВАНИЯ ТЕЧЕНИЯ ЖИДКОСТИ В ТОНКИХ КАПИЛЛЯРАХ.-222 Л. Капиллярные течения в гемодинамике.- 22

2.2. Уравнение заполнения пустого капилляра.- 23

2.2.1. Случай горизонтально расположенного капилляра.- 25

2.2.2. Методика эксперимента.- 27

2.2.3. Разработка устройства для измерения вязкости.- 37

2.2.4. Расположение капилляра под углом к горизонту.-412.3. Оценка высоты капиллярного поднятия.-512.4. Неравномерность капиллярного течения.- 56

ГЛАВА 3. МЕТОДИКА ДИАГНОСТИКИ АТЕРОСКЛЕРОЗА АРТЕРИЙ. - 66 -3.1. Средства и модели для изучения колебательных процессов в организме

663.2. Уравнения Матье и Хилла.- 66

3.2.1. Вывод уравнения Матье.- 68

3.2.2. Построение дифференциального уравнения Матье по известному дифференциальному уравнению Хилла.- 68

3.2.3. Проверка предложенного алгоритма.- 71

3.3. Прикладная задача с использованием предложенного алгоритма. - 75

3.4. Критерий оценки состояния сердечно-сосудистой системы по фазовой траектории сигнала.- 82

Актуальность:

Капиллярно-гидродинамические течения жидкостей широко распространены в природе и технике, однако особую значимость капиллярные явления приобретают при исследовании задач гемодинамики живых организмов. Характер кровотока определяет функциональность и жизнеспособность организма, и поэтому исследование капиллярно-гидродинамических течений имеет, помимо многих других областей знаний, важное прикладное значение для медицины.

Несмотря на практическую значимость движения жидкостей в капиллярах, этот процесс недостаточно изучен. Отчасти это связано с тем, что межфазные эффекты весьма чувствительны к примесям и физическому состоянию поверхности, различным неровностям и шероховатостям, что не всегда можно учесть в методике эксперимента. Течение биологических жидкостей имеет ряд существенных особенностей, и, если говорить о гемодинамике, то важен как состав крови, так и диаметр и поверхность кровеносного сосуда. Более того, по упругим стенкам крупных кровеносных сосудов от сердца распространяется пульсовая волна, что является еще одним требующим учета фактором. Недостаточно полно развитые приборно-экспериментальная база и методы исследования подобных процессов позволяют говорить об актуальности предлагаемой работы.

В настоящей работе описывается разработка математических моделей, а на их базе - экспериментальных методик изучения нестационарных течений жидкости в капиллярах. С их помощью исследуется скачкообразный процесс заполнения пустого капилляра, изучаются пульсации давления в пульсовой волне, распространяющейся по стенке кровеносного сосуда. На основе выработанной математической методики обработки данных сфигмографического обследования предлагается способ оценки упруго-эластических свойств артерий и ранней экспресс-диагностики атеросклероза, для чего создана измерительная приставка к компьютеру с программно-вычислительным комплексом обработки экспериментальных результатов.

Приборы и методики исследования, применяемые в настоящей работе, актуальны в различных областях науки, особенно при изучении течений в закрытых и щелевых каналах, при исследовании фильтрации жидкостей в пористых структурах, в медицине, при исследовании состояния сердечнососудистой системы.

Цель работы.

Основной целью работы являлась разработка приборной и методической базы исследования капиллярных течений различных жидкостей в нестационарном режиме. При этом ставились следующие задачи:

1. Создать экспериментальную базу и разработать методику исследования динамики заполнения капилляра, и провести экспериментальные исследования характера капиллярного течения.

2. Разработать математическую модель процесса натекания жидкости в капилляр и проверить ее экспериментально.

3. Применительно к задачам гемодинамики разработать физическую модель квазипериодического процесса распространения пульсовой волны в организме на базе параметрических дифференциальных уравнений Хилла и Матье.

4. Разработать статистический критерий оценки степени хаотичности фазовых траекторий пульсовой волны применительно к диагностике атеросклероза, и на этой базе создать измерительную установку с программно-вычислительным комплексом обработки экспериментальных результатов.

Практическая значимость

1. Разработана новая методика измерения вязкости жидкости и изготовлена измерительная установка на ее основе.

2. Разработан статистический критерий оценки степени хаотичности фазовых траекторий пульсовой волны применительно к диагностике атеросклероза с математической обработкой данных сфигмографии, и созданы измерительная установка, сопряженная с компьютером, и программно-вычислительный комплекс обработки экспериментальных результатов.

Научная новизна работы

1. На основе предложенной математической модели получено приближенное уравнение для процесса заполнения жидкостью капилляра, имеющее аналитическое решение, удовлетворительно согласующееся с численным решением точного нелинейного уравнения и экспериментальными данными.

2. С использованием созданной экспериментальной методики обнаружен скачкообразный характер капиллярного движения жидкости, впервые получены количественные результаты и предложено качественное описание этого явления.

3. На основании применения разработанной экспериментальной методики сделан вывод о влиянии адгезионных свойств материала поверхностей растекания на капиллярное движение жидкостей, и предложена поправка для учета этого влияния.

4. Для задач гемодинамики разработана физическая модель квазипериодических процессов организма на базе параметрических дифференциальных уравнений Хилла и Матье.

-65. Предложен новый способ оценки упруго-эластических свойств артерий и диагностики атеросклероза на основе анализа фазовых траекторий полученного сфигмографического сигнала.

Основные положения, представляемые к защите

1. Экспериментальная методика для измерения вязкости жидкости, основанная на учете адгезионных свойств поверхности капилляров и коэффициентов поверхностного натяжения эталонной и исследуемой жидкостей.

2. Новая физическая модель процесса натекания жидкости в капилляр, позволяющая адекватно описать динамику капиллярного течения.

3. Экспериментальная методика обнаружения скачкообразного характера движения мениска жидкости в капилляре с количественными результатами описания этого явления.

4. Физическая модель квазипериодической пульсовой волны, распространяющейся в организме, на базе параметрических дифференциальных уравнений Хилла и Матье.

5. Программно-аппаратный комплекс обработки экспериментальных результатов, основанный на предложенном статистическом критерии оценки степени хаотичности фазовых траекторий пульсовой волны применительно к диагностике атеросклероза.

Основные выводы и результаты работы:

Разработана новая методика определения вязкости жидкости, основанная на учете адгезионных свойств поверхности капилляров и коэффициентов поверхностного натяжения эталонной и исследуемой жидкостей, и создана установка для измерения вязкости.

Предложена математическая модель, в рамках которой записано нелинейное уравнение для процесса натекания жидкости в капилляр и получено приближенное уравнение, имеющее аналитическое решение, удовлетворительно согласующееся с численным решением и экспериментальными данными.

Разработана экспериментальная методика и создан экспериментальный стенд, с помощью которого удалось обнаружить скачкообразный характер капиллярного движения жидкости, впервые получены количественные характеристики этого движения и предложено его качественное описание.

Разработана физическая модель квазипериодического процесса распространения пульсовой волны в организме на базе параметрических дифференциальных уравнений Хилла и Матье применительно к задачам гемодинамики.

Разработан статистический критерий оценки степени хаотичности фазовых траекторий пульсовой волны применительно к диагностике атеросклероза с математической обработкой данных сфигмографии, и созданы измерительная установка, сопряженная с компьютером, и программно-вычислительный комплекс обработки экспериментальных результатов.

ЗАКЛЮЧЕНИЕ

1. Адамсон А. Физическая химия поверхностей. М.: Мир. 1979. 568 с.

2. Моррисон С. Химическая физика поверхности твердого тела. М.: Мир. 1980. 488 с.

3. Физика поверхности полупроводников // Сб. статей. Перевод под ред. Г.Е.Пикуса. М.: Иностранная литература. 1959. 423 с.

4. Электронные процессы на поверхности и в монокристаллических слоях полупроводников // Труды симпозиума. Под ред. A.B. Ржанова. Новосибирск: Наука. 1967. 240 с.

5. Физико-химические процессы в полупроводниках и на их поверхности. Воронеж: Изд-во ВГУ. 1981. 200 с.

6. Полтавцев Ю.Г., Князев A.C. Технология обработки поверхностей в микроэлектронике. Киев: Техника. 1990. 206 с.

7. Власов В.М. Работоспособность упрочненных трущихся поверхностей. М.: Машиностроение. 1987. 304 с.

8. Hauksbee F. Physico-mechanical experiments, London, 1709, pp. 139-169.

9. Maxwell J.C. Capillary action. / The encyclopaedia britannica, 11th edition, cambrige: at the university press, 1910, vol. 5, p. 256

10. Гиббс Дж. В. Термодинамические работы. / Пер. с англ. Под ред. В.К. Семенченко. M. JL, Гостехиздат, 1950. - 492 с.

11. Русанов А.И. 100 лет теории капиллярности Гиббса. // В сборнике: Современная теория капиллярности. JL: Химия, 1980.

12. Ван-дер-Ваальс И.Д., Констам Ф. Курс термостатики. 4.I.II. М.:ОНТИ, 1936.

13. Bakker С. Kapillarität und Oberflächenspannung. / Handb. der exper. phys. bd. vi. Leipzig, Wien-harms, 1928.

14. Гугенгейм Э. Современная термодинамика. / Пер. с англ. М.: Госхимиздат, 1941.

15. Русанов А. И. Термодинамика поверхностных явлений. JL, Изд. ЛГУ, 1960.

16. Русанов А. И. Фазовые равновесия и поверхностные явления. JL, «Химия», Ленингр. отд-ние 1967. 388 с.

17. Эткинс П. Физическая химия. Т.2. М.: Мир. 1980. 584 с.

18. Френкель Я.И. Кинетическая теория жидкостей. М. 1983.

19. Стасенко Л. Молекулы, сосиски и алмазы // жур. Квант. 2003. №1.

20. Дедков Г.В. Нанотрибология: экспериментальные факты и теоретические модели // Успехи физических наук. 2000. Т. 170. №6. С.585-618.

21. Харламов В.Ф., Ануфриев K.M., Крутовский Е.П., Мосин Ю.В., Злоткин Е.А., Емельянов И.В. Кинетика адсорбции и рекомбинации атомов водорода на поверхности твердых тел // Письма в ЖТФ. 1998. Т. 24. №5. С.23-27.

22. Харламов В.Ф., Васильев Н.Ф., Иващук O.A., Крутовский Е.П., Мосин Ю.В., Злоткин Е.А. Начальный пик на зависимости от времени скорости гетерогенной рекомбинации атомов водорода на поверхности кристаллофосфоров // Письма в ЖТФ. 1998. Т. 24. №3. С.54-59.

23. Вайнштейн Э.Ф. Зависимость поверхностного натяжения от размера системы и диаметра волокна // Институт химии растворов РАН, 1998.

24. Гегузин Я.Е. Капля./ Я.Е. Гегузин; АН СССР. 2-ое доп. изд., М., "Наука", 1977. - 176 е.: ил. (науч. попул. серия)

25. Сульман Э.М., Семагина Н.В., Матвеева В.Г. Влияние поверхностных характеристик участников реакции на скорость селективного каталитического гидрирования// Вестник ТГТУ. №1 (1). 2002, с. 28-31.

26. Оно С., Кондо С. Молекулярная теория поверхностного натяжения в жидкостях./ Под ред. И.З. Фишера. М., Изд. иностр. лит., 1978, 291 с.

27. Захарченко В.Н. Коллоидная химия: Учеб. для медико-биолог. спец. вузов. 2-е изд., перераб. и доп. - М.: Высш. Шк., 1989. - 238 е.: ил.

28. Лабораторные работы и задачи по коллоидной химии / Под ред. Ю.Г. Фролова и A.C. Гродского. М.: Химия, 1986. 216 с.

29. Синяченко О.В., Губергриц Н.Б., Череватская Е.Ю., Иванова Е.В. перспективы применения динамической межфазной тензиометрии в клинической гематологии. // Украиский медицинский журнал. № 1 (15), 2000, с. 16-20.

30. Голубев М.В., Халикова А.И. Оптимизация процесса подачи деэмульгатора. // Нефтегазовое дело, 2004, http://www.ogbus.ru

31. Физические основы измерений: Метод, указ. / Сост. Д.М. Мордасов, М.М. Мордасов. Тамбов: Изд-во Тамб. гос. ун-та., 2002, 32 с.

32. Русанов А. И., Прохоров В. А. Межфазная тензиометрия.— СПб: Химия.— 1994.— 400 с.

33. Соболева O.A. Влияние pH, добавок карбамида и бромида натрия на гистерезис смачивания в системе водный раствор цетилтриметиламмоний бромида - стекло. // Жур. Вестн. Моск. Ун-та. сер.2. Химия. Т. 40, № 4, 1999, с. 267-269.

34. Сумм Б.Д. Гистерезис смачивания // Соросовский образовательный журнал, № 7, 1999, с. 98 102.-10238. Де Жен П.Ж., Смачивание: статика и динамика. // Успехи физических наук, т. 151, вып. 4, 1987-с. 619-681.

35. Дерягин Б.В., Карасев В.В. Изучение граничной вязкости органических веществ по кинетике утоныпения их смачивающих слоев в процессе сдувания. // Доклады АН СССР, т. 101, №2, 1955, с. 289-292.

36. Блинов В.И., Худяков Г.Н., Петров И.И., Реутт В.Ч. О движении жидкости в резервуаре при перемещении ее струей воздуха // Инж.-физ. журнал, т. 1, №11, 1958, с. 6 13.

37. Карасев В.В., Дерягин Б.В. Измерения граничной вязкости по кинетике утоныпения смачивающих пленок жидкостей в процессе сдувания // ЖФХ, т. 33, вып. 1, 1959, с. 100- 106.

38. Голик А.З., Соломко В.П. Исследование физических свойств системы вода ацетон - спирты. // Укр. хим. журнал, т. 24, вып. 6, 1958, с. 594 - 600.

39. Дерягин Б.В., Захаваева H.H., Карасев В.В., Лазарев В.И. Механизм граничной смазки и свойства граничного смазочного слоя // ЖТФ, т. 27, вып. 5, 1957, с. 1076- 1086.

40. Дерягин Б.В., Зорин З.М. Оптическое исследование адсорбции и поверхностной конденсации паров вблизи насыщения. II Доклады АН СССР, т. 98, № 1,1954, с. 93 -96.

41. Базарон У.Б., Булгадаев A.B., Дерягин Б.В. Измерение сдвиговой упругости жидкостей и их граничных слоев резонансным методом. // ЖТФ, Т. 51, Вып. 4(10), 1966.

42. Дерягин Б.В., Чураев Н.В. Новые свойства жидкостей. Сверхплотная вода. М.: Наука, 1971.- 176 с.

43. Федякин H.H. О температурном расширении воды в макрокапиллярах. // Доклады АН СССР, т. 138, № 6,1961, с. 1389- 1391.

44. Дерягин Б.В., Нерпин С.В., Арутюнян М.А. О механокалорическом эффекте при обычных температурах. // Доклады АН СССР, т. 160, № 2, 1965, с. 387-389.

45. Дерягин Б.В. Определение структурной характеристики и удельной поверхности пористого тела по изотерме десорбции // ЖФХ, т. 31, вып. 2, 1957, с. 516 518.

46. Мезон У. Пьезоэлектрические кристаллы и их применение в ультраакустике, ИИЛ, 1952.

47. Мандельштам Л.И. Собр.соч., 4, Л., АН СССР, 1955, с. 217-218.

48. Хайкин С.Э., Лисовский Л.П., Саломонович А.Е. О силах «сухого трения». // Доклады АН СССР, т. 24, № 2,1939, с. 134 138.

49. Сабсай О.Ю., Дубинский М.Б., Барштейн Г.Р. Условия перехода к нерегулярному режиму течения расплавов термопластов // Пластические массы, 1991, № 9, с. 40 41.

50. Сухарев Ю.И., Лужнова О.В., Крупнова Т.Г., Васильева A.B. Эффект сильного влияния малых воздействий на свойства неравновесной гелевой системы оксигидрата ниобия // Известия Челябинского научного центра. Вып. 3 (16), 2002, с. 98-101

51. Нуген Х.Л., Сагитов Р.Ф., Полищук В.Ю. Моделирование течения материала высокой вязкости для исследования его свойств при больших скоростях сдвига. // Вестник ОГУ, № 4, 2004, с. 142 145.

52. Ремизов А.Н. Медицинская и биологическая физика: Учеб. для мед. спец. вузов. 3-е изд. испр. - М.: Высш. школа, 1999. - 616 с.

53. Гатчек Э. Вязкость жидкостей / Пер. с англ. М.П. Воларовича, Д.Н. Толстого. М., 1932.-215 с.

54. Гуляев Ю.В., Обухов Ю.В., Сударев A.M. Методы, аппаратура и информационные технологии телемедицинской функциональной диагностики. // Фундаментальная наука медицине. Материалы конференции, г. Москва, 2526.11.2002, -М.: Фирма "Слово", 2002. стр. 90-91.

55. Бурцев Е.В., Райков А.Н., Сасорова Е.В., Мультипараметрическая оценка изменения кровообращения при стрессе и коррекционных процедурах с помощью тахоосциллографии. // Вестник новых медицинских технологий. 1998, т. V, 1,ч. 3

56. Гаврилушкин А.П., Киселев С.В., Медведев А.П., Шелепнев A.B., Маслюк А.П. Геометрический анализ нелинейных хаотических колебаний в оценке вариабельности ритма сердца. Учеб.-метод. пособие для студентов и врачей. Нижний Новгород; 2001.

57. Баевский P.M., Кириллов О.И., Клецкин С.З. Математический анализ изменений сердечного ритма при стрессе. М: Наука; 1984; 221 с.

58. Борисов В.И. Анализ вариабельности сердечного ритма в оценке состояний больных инфарктом миокарда. Нижний Новгород: Изд-во НГМА; 1997.

59. Матусова А.П. Клиническое значение статистического анализа сердечного ритма у больных с острым инфарктом миокарда. Кардиология 1989; 29(1): с. 29—32.

60. Степура О.Б., Остроумова О.Д., Курильченко И.Т., Мартынов А.И. Клиническая значимость сигналусредненной электрокардиографии // Российские Медицинские Вести №1, 1997.

61. Латфуллин И.А., Тептин Г.М. Возможности электрокардиографии высокого разрешения. //Казанский медицинский журнал, 1998, 79, 2.

62. Латфуллин И.А., Тептин Г.М., Контуров С.В. Сравнительный анализ погрешностей в аналоговой и компьютерной электрокардиографии. // Вестник аритмологии. 1999, № 13, с. 54-58.

63. Писарук A.B. Анализ механизмов возрастных изменений системы барорефлекторной регуляции с помощью математической модели // Пробл. старения и долголетия. 1999, № 2.

64. В.А. Лукшин, С.И. Мухин, Т.В. Соколова, Н.В. Соснин, А.П. Фаворский, Математическое моделирование церебральной гемодинамики в квазипериодическом режиме // Препринт. Москва, МАКС Пресс, 2003, 20 с.

65. Хидирова М.Б. Об одной замкнутой модели сердечно-сосудистой системы: Основные уравнения и результаты вычислительных экспериментов // Проблемы механики, №5, 1998, 26-30.

66. Полубаринова-Кочина П.Я. Теория движения грунтовых вод. Изд. 2-е. М.: Наука, 1977.-664 с.

67. Сумм Б.Д., Горюнов Ю.В. Физико-химические основы смачивания и растекания. М.: Химия, 1976. 324 с.

68. Дерягин Б.В., Чураев Н.В., Муллер В.М. Поверхностные силы. М.: Наука, 1985.-398 с.

69. Штеренлихт Д.В. Гидравлика, М.: Энергоатомиздат, 1984, 640 с.

70. Левич В.Г. Физико-химическая гидродинамика. М.: Издательство академии наук, 1952, 538 с.

71. Лескова С.С. Диагностика свойств жидкости на границах раздела гетерогенных сред: Автореф. дис. к-та ф.-м. наук. Барнаул, 2006. - 20 с.

72. Волков В.И., Козлов Д.Ю., Кирколуп Е.Р. Исследование затекания жидкости в тонкие капилляры при малых числах Вебера // Журнал технической физики, 2007, т. 77, Вып. 7, с. 24-27.

73. В. 3. Канчукоев, Б. С. Карамурзов, В. А. Созаев, В. В. Чернов. Определение начальной скорости и продолжительности движения жидкости в капиллярах. ИФЖ, Том 76, Номер 1, 2003.

74. Волков В.И., Козлов Д.Ю. и Кирколуп Е.Р. Динамика капиллярного натекания // Известия вузов. Физика, №5, Томск 2007 - с. 52-57.

75. Волков В.И., Козлов Д.Ю., Кирколуп Е.Р. Исследование динамики движения жидкости по капилляру // Известия АТУ, № 1 (53), Барнаул: изд-во АлтГУ, 2007.-с. 100- 104.

76. Волков В.И., Лескова С.С., Кирколуп Е.Р. Экспериментальное исследование гистерезиса смачивания с помощью плоского капилляра // Вестник уральской медицинской академической науки, 2006, №2, с. 48-50.

77. Волков В.И., Козлов Д.Ю., Оценка капиллярного поднятия. // Известия АГУ, № 1 (53), Барнаул: изд-во АлтГУ, 2007. - с. 105 - 107.

78. Волков В.И., Кадышева С.С. Экспериментальное исследование напряжения сдвига у некоторых жидкостей. // Вестник Уральской Медицинской Академической Науки, Екатеринбург, 2005, № 1, с. 46 48.

79. Ghiradella Н., Radigan W., Frish Н. L. // 1975, J. Colloid and Interface Sci. V. 51, P. 522.

80. Мак-Лахлан Н.В. Теория и приложения функций Матье, М., ИЛ, 1951

81. Григорьев А.И., Голованов A.C. Деформации и перекрытия зон неустойчивости уравнения Матье-Хилла, Письма в ЖТФ, 1999, т. 25, вып. 20, с. 13-18.

82. Ландау Л.Д., Лифшиц Е.М. Механика. М.: Наука, 1965, 204 с.

83. Андронов A.A., Витт A.A., Хайкин С.Э., Теория колебаний, изд. 2, М.: Физматгиз, 1959.

84. Safar H., Mourad J.-J., Safar M., Blacher J. Aortic Pulse Wave Velocity, an Independent Marker of Cardiovascular Risk // Arch Mal Coeur 2002; 95:1215-1218.

85. Волков В.И., Козлов Д.Ю. Математические модели сфигмографии // Сборник научных работ «Актуальные проблемы фармации», Вып. IV, Барнаул: Изд-во АГМУ, 2007, с. 55-62.

86. Волков В.И., Куликов В.П., Останин С.А., Колосов Д.А. Пьезоэлектрический метод исследования пульсаций кровяного давления // Сборник научных работ «Актуальные проблемы фармации», Вып. IV, Барнаул: Изд-во АГМУ, 2007, с. 84-90.