Классификация биотехнических систем. Биотехнические системы и технологии

1. Классификация биотехнических систем

КЛАССИФИКАЦИЯ
БИОТЕХНИЧЕСКИХ
СИСТЕМ
Биотехнические системы и
технологии

2. КЛАССИФИКАЦИЯ БТС

По характеру основной целевой функции БТС разделены на
3 группы:
- БТС медицинского назначения (БТС – М);
- БТС эргатического типа (БТС – Э) с человеком – оператором в качестве
управляющего звена;
- БТС целенаправленного управления поведением целостного организма (БТС
– У).
В настоящее время к этим группам дополнительно относятся:
БТС лабораторного анализа (БТС-ЛА);
Мониторные системы (МС);
БТС искусственного жизнеобеспечения.

3. БТС медицинского назначения (БТС – М)

БТС МЕДИЦИНСКОГО НАЗНАЧЕНИЯ (БТС – М)
БТС-М предназначены для использования в медицинских целях,
главным из которых являются:
Диагностика состояния живого организма (текущая),
проводимая
в
реальном
масштабе
времени,
и
дифференциальная, осуществляемая в процессе апостериорной
медико-биологической информации.
Управление состоянием организма для его нормализации
(методами дискретной или непрерывной коррекции).
Временная или длительная компенсация утраченных функций
органов или физиологических систем живого организма.
Протезирование и коррекция функций сенсорных систем или
двигательного аппарата.
Различные медико-биологические исследования и лечебные
процедуры, связанные с применением приборов активного
вмешательства, сочлененных с живым организмом в единую
БТС.

4. БТС медицинского назначения (БТС – М)

БТС МЕДИЦИНСКОГО НАЗНАЧЕНИЯ (БТС – М)
Вся медико-биологическая информация разделяется на медленно
изменяющиеся процессы (МИП) и быстроизменяющиеся процессы
(БИП). При этом к МИП относят процессы, частота изменения которых
ниже 1 Гц, а мгновенные значения могут быть выражены цифрой.
Пример:
частота сердечных сокращений;
частота дыхания;
температура тела и др.
К БИП относятся, главным образом, электрофизические процессы,
характеризуемые изменением электрических потенциалов на
поверхности кожи или на отдельных локальных участках организма под
ее покровами.
Пример:
электрокардиосигнал (ЭКС) (ЭКГ);
электроэнцефалограмма (ЭЭГ);
электромиограмма (ЭМГ) и др.
Частотные спектры БИП находятся в диапазоне 0,1…1000 Гц.

5. БТС медицинского назначения (БТС – М)

БТС МЕДИЦИНСКОГО НАЗНАЧЕНИЯ (БТС – М)
Схема получения медико-биологической информации
Биообъект
ДПИ
МИП
БИП
Блок предварительной
обработки
Передатчик ЛФ И

6. БТС медицинского назначения (БТС – М)

БТС МЕДИЦИНСКОГО НАЗНАЧЕНИЯ (БТС – М)
Структурная схема измерительно-информационной БТС – М
ДПИМИП
БУ
БП
ААС
П
ЛС
ДПИБИП
В
1
2
БУ
БСИ
3
СОИ
ДП
БР
ОП
БВИП

7. Биотехнические системы эргатического типа (БТС – Э)

БИОТЕХНИЧЕСКИЕ СИСТЕМЫ ЭРГАТИЧЕСКОГО ТИПА (БТС –
Э)
В подобных системах (БТС – Э) или (по общепринятой терминологии) в
системах “человек-машина” биологическое звено представлено
человеком-оператором, выполняющим различные функции в замкнутом
контуре управления технической системой.
В процессе анализа и синтеза БТС – Э практически приходится иметь
дело с моделированием на четырех уровнях:
тканевом (в случае моделирования процессов метаболизма при
определении энергозатрат организма человека в процессе
деятельности);
органном (при бионических исследованиях органов восприятия для
оптимального сопряжения их с техническими элементами
информационных каналов);
организменном (при исследовании поведения оператора в процессе
управления динамическими процессами);
популяционном (про разработке методов управления популяцией
живых организмов с целью организации их целенаправленного
поведения).

8. Биотехнические системы эргатического типа (БТС – Э)

БИОТЕХНИЧЕСКИЕ СИСТЕМЫ ЭРГАТИЧЕСКОГО ТИПА (БТС
– Э)
Техническая реализация метода поэтапного моделирования
воплощается в тренажерно – моделирующем комплексе (ТМК).
ТМК является многоцелевой универсальной системой,
предназначенной для:
учета человеческого фактора на ранних стадиях синтеза БТС – Э;
оптимизация режима работы оператора, выработки
количественных оценок эффективности их деятельности на
пультах управления различных конструкций и разработки
инженерных рекомендаций по усовершенствованию
существующих систем управления эргатического типа;
получение априорной информации с целью обоснования модели
идеального оператора применительно к системе данного типа;
разработки методик профотбора, обучения и тренировок
операторов, а так же для подготовки операторов методов
стадийной тренировки;
тренировки операторов комплексов и систем в период их
разработки задолго до их воплощения в действующие –
конструкции.

9. Биотехнические системы эргатического типа (БТС – Э)

БИОТЕХНИЧЕСКИЕ СИСТЕМЫ ЭРГАТИЧЕСКОГО ТИПА (БТС
– Э)
Структурная схема тренажерно – моделирующего комплекса
Моделирование внешней
обстановки
Среда
Моделирование
динамики объекта
управления
Анализ деятельности
оператора
Объект
воздействия
Анализ
психофизиологического
состояния оператора
Рабочее место оператора
Рабочее место инструктора
Управление
обстановкой
Управление
объектом
Управление
режимами
Блок
программ
Индикация
Регистрация
Оценка стоимости
результата

10. Биотехнические системы эргатического типа (БТС – Э)

БИОТЕХНИЧЕСКИЕ СИСТЕМЫ ЭРГАТИЧЕСКОГО ТИПА
(БТС – Э)
Структурная схема адаптивной БТС – Э
С
ОбУ
ПУ
ЛФП1
ИП
СУС
О
ОУ
СОДО
МА
МИП
ЛФП2
БИП
НСО
БВИП
МС
КС
ААС

11. Биотехнические системы управления (БТС-У)

БИОТЕХНИЧЕСКИЕ СИСТЕМЫ УПРАВЛЕНИЯ (БТС-У)
Экстракорпоральный искусственный орган
Fi
Фс
Fc i
i
Fdi
Фd
i
Мb- биологический организм, Fi-искусственный орган,
Fci-управляющий процессор, Fdi-исполнительный процессор,
Фci-управляющая среда,
Фdi-управляемая среда
Mb
Фb
i

12. Биотехнические системы управления (БТС-У)

БИОТЕХНИЧЕСКИЕ СИСТЕМЫ УПРАВЛЕНИЯ (БТС-У)
Интракорпоральный искусственный орган
Fi
Фс
Fc
i
Fd
i
Фd
i
i
Mb

13. Биотехнические системы управления (БТС-У)

БИОТЕХНИЧЕСКИЕ СИСТЕМЫ УПРАВЛЕНИЯ (БТС-У)
Искусственный
орган Mbi
Сердце
Управляющее воздействие
органа
Экстракорпоральный
Интракорпональный
Легкое
Экстракорпоральный
Интракорпональный
Почка
Экстракорпоральный
Интракорпональный
Поджелудочная
железа
Управляющий процессор Fci
Исполнительный процессор
Fdi
Не лимитированный источник
энергии
Неимплантированный
или
имплантированный источник
энергии
Неимплантированный насос
Имплантированный насос
Источник газов(баллоны
с О2,СО2 и воздухом)
Аппарат искусственного
дыхания
Оксигенатор
Источник
диализата
Источник
диализата
Диализатор
нестерильного
стерильного
Экстракорпоральный
Неимплантированный
дозатор инсулина
Интракорпональный
Неимплантированный
или
имплантированный
дозатор инсулина
Легкие
Брюшная
организма
полость
Магистраль
экстракорпоральной
кровопроводящей системы
Сосуд сердечно-сосудистой
системы организма

14. Биотехнические системы управления (БТС-У)

БИОТЕХНИЧЕСКИЕ СИСТЕМЫ УПРАВЛЕНИЯ (БТС-У)
Биотехнический организм
Ma
Fc i
Fd i
Kd
Fc m
Fdm
Кровь
Mb
Fc
m+1
Fc n
Fd
m+1
Fdn
Ма - биотехнический организм, Kd - исполнительный канал связи

15. ПРИНЦИПЫ СИНТЕЗА и моделирования БИОТЕХНИЧЕСКИХ СИСТЕМ

ПРИНЦИПЫ СИНТЕЗА И
МОДЕЛИРОВАНИЯ
БИОТЕХНИЧЕСКИХ
СИСТЕМ
Биотехнические системы
и технологии

16. ПРИНЦИПЫ СИНТЕЗА БИОТЕХНИЧЕСКИХ СИСТЕМ

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

17. ПРИНЦИПЫ СИНТЕЗА БИОТЕХНИЧЕСКИХ СИСТЕМ

Основные принципы сопряжения технических и
биологических элементов
1.
Принцип адекватности, требующий согласования основных
конструктивных
параметров
и
"управленческих
характеристик" биологических и технических элементов
БТС.
2.
Принцип единства информационной среды, требующий
согласования
свойств
информационных
потоков,
циркулирующих между техническими и биологическими
элементами как в афферентных, так и в эффекторных цепях
БТС.

18. ПРИНЦИПЫ СИНТЕЗА БИОТЕХНИЧЕСКИХ СИСТЕМ

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

19. ОБЩИЕ ВОПРОСЫ МОДЕЛИРОВАНИЯ В БТС

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

20. ОБЩИЕ ВОПРОСЫ МОДЕЛИРОВАНИЯ В БТС

Классификация математических моделей
Одномерные вероятностные модели : параметрические; непараметрические.
При построении параметрической модели распределение случайных величин
аппроксимируется известным вероятными законами распределения (чаще всего
нормальным законом).
Пример: когда в справочнике указано, что у здорового человека в состоянии покоя объем
вдыхаемого воздуха составляет 400-800 мл, это значит, что распределение величины
объема моделируются равномерным законом.
Непараметрическими моделями называются те, у которых закон распределения случайных
величин неизвестен.
Многомерные вероятностные модели
Представляют собой многомерные описания состояния биотехнических систем и процессов
в них протекающих, к признакам этих систем относят существенные параметры,
характеризующие определенные единицы наблюдения.
Двухмерное пространство признаков, образованное двумя переменными, содержит
скопление точек на плоскости.
Трехмерное пространство, образованное тремя переменными содержит облака точек наблюдений.
Пример: геометрические размеры и биохимические характеристики клетки, интервалы
между зубцами ЭКГ и т. д.

21. ОБЩИЕ ВОПРОСЫ МОДЕЛИРОВАНИЯ В БТС

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

22. ОБЩИЕ ВОПРОСЫ МОДЕЛИРОВАНИЯ В БТС

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

23. ОБЩИЕ ВОПРОСЫ МОДЕЛИРОВАНИЯ В БТС

Основные этапы моделирования
1. Постановка цели моделирования
Для одной и той же системы можно составить множество моделей. Они будут
отличаться степенью детализации и учета тех или иных особенностей системы и
режимов ее функционирования.
2.Создание концептуальной модели
Концептуальная (или содержательная) модель – это словесное описание системы –
ее состав и структура, свойства входящие в нее элементов и причинно–
следственные связи.
3. Формирование информационного обеспечения модели
При создании концептуальной модели определяются качественные и
количественные параметры БТС и внешних воздействий на нее. К качественным
параметрам относят те, которые нельзя оценить в каких- то единицах измерений.
Для количественных параметров необходимо определить их конкретные значения,
которые будут использованы в виде исходных данных при моделировании.
4. Разработка математической модели
5. Выбор метода моделирования
Разработанная математическая модель может быть исследована различными
методами: аналитическими; имитационными.

24. ОБЩИЕ ВОПРОСЫ МОДЕЛИРОВАНИЯ В БТС

Основные этапы моделирования
6. Технические средства моделирования
7. Корректировка модели
При выявлении недопустимого рассогласования модели и системы возникает
необходимость в корректировке модели. При этом могут быть выделены следующие
типы изменений:
глобальные (требуют для устранения разработки новой модели);
локальные (требуют частичного изменения математической модели);
параметрические (изменения некоторых параметров (калиброванных).
8. Проверка адекватности модели
9. Планирование экспериментов с моделью
а) Стратегическое планирование.
б) Тактическое планирование.
10. Обработка измерений имитационного эксперимента
11. Использование результатов моделирования

25. ОБЩИЕ ВОПРОСЫ МОДЕЛИРОВАНИЯ В БТС

Метод поэтапного моделирования
Подготовительный (этап 1)
Разрабатывается
структурно
–
функциональная
схема
БТС,
конкретизируются её целевая функция (ЦФ) и возможные режимы работы.
Определяется биологический объект и предварительный алгоритм его
функционирования в БТС.
На основании априорных данных создается модель БТС с математической
моделью биологического элемента.
Управленческое согласование характеристик элементов БТС (этап 2)
Осуществляются итерационные процедуры согласования характеристик
элементов БТС в едином контуре управления.
Практическое значение результатов этого этапа заключается в том, что
можно отобрать операторов для БТС эргатического типа или подобрать
прямой биологический аналог человеческому организму при экспериментах.

26. ОБЩИЕ ВОПРОСЫ МОДЕЛИРОВАНИЯ В БТС

Метод поэтапного моделирования
Информационное согласование (этап 3)
Исследуются информационные процессы, обеспечивающие соблюдение
принципов адекватности и идентификации информационной среды.
Корректируются решающие правила, заложенные в виде программ в
системы обработки информации о состоянии биологического объекта;
Разрабатываются
специальные
требования
к
техническим
устройствам, согласующим информационные и управленческие
характеристики технической и биологической частей БТС, получившие
названия логических фильтров – преобразователей (ЛФП).
Заключительный (этап 4)
Проводится исследование БТС в полунатурных (модельных) и натурных
условиях.
Идет обработка данных эксперимента и окончательная корректировка
математической модели.
Подготавливаются задания на инженерную разработку БТС.

27.

Примеры моделирования БТС

28.

Основные направления математического моделирования гемодинамики
Моделирование течения крови в отдельном сосуде
2D и 3D модели течения в крупных сосудах (уравнения Навье-Стокса)
(моделирование упругости стенки сосуда, турбулентность течения, многокомпонентность
крови (не ньютоновская жидкость), взаимодействие со стенками сосуда, области
бифуркации сосудов, моделирование тромбообразования, стенозов, аневризмов и т.д.)
2D и 3D модели течения крови в мелких сосудах (капиллярах)с учетом реологии.
Моделирование течения крови в сердце (2D и 3D модели)
Моделирование течения крови в сети сосудов (дерево сосудов, замкнутая
система) для исследования общих закономерностей течения крови.
- На основе балансных соотношений
- На аналогиях с «электрической цепью»
- Квазиодномерное приближение:
- соответствует типу сосудистой сети
- дает возможность описать систему кровообращения в целом
- позволяет отслеживать параметры течения крови вдоль сосуда
- позволяет учесть особенности каждого сосуда
- является основой для построения разномасштабных моделей
-предоставляет возможность расчета переноса веществ кровью
-предоставляет возможность использовать различные модели сосудов и органов
-использует доступные физиологические данные
-обладает хорошей точностью
-предъявляет разумные требования к вычислительным мощностям

29.

Комплексная нелокальная математическая модель
сердечно-сосудистой системы
1. Создание математической
модели течения крови в
замкнутой системе сосудов
(графе сосудов) произвольной
топологии
2. Разработка эффективных
моделей различных органов,
сопряженных с работой
сердечно-сосудистой системы (в
том числе – точечной модели
сердца )
Базовая модель
описание течения крови в сосуде
квазиодномерными уравнениями гемодинамики +
нелокальные граничные условия +
условия сопряжения в точках бифуркации
Модели сосудов
Модели почки
S
q
kid
q
stat
exp( p p
stat
)
Ki
dn
Модели сердца Renal
p
ey
outflow
qkid
Q
t
3. Создание эффективных
однородных методов описания
графа сосудов и численного
решения глобальной
математической модели
Однородная консервативная неявная
разностная схема для системы уравнений
на графе

30.

Комплексная нелокальная математическая модель сердечно-сосудистой системы
4. Создание интерактивного
программного комплекса со
средствами подготовки и
обработки данных.
5. Построение и анализ точных
решений системы уравнений
гемодинамики на графе
6. Проведение вычислительных
экспериментов в интересах
фундаментальной и
практической физиологии и
медицины
Моделирование
церебральной гемодинамики
Моделирование влияния
гравитационных нагрузок на
сердечно-сосудистую
систему
g
Моделирование влияния
физических нагрузок
Моделирование переноса веществ
кровью по графу сосудов с учетом
процессов сорбции-десорбции
Моделирование почечной
регуляции давления
...

31.

Адекватность свойств модели законам функционирования сердечнососудистой системы.
Формализация системы
кровообращения в граф
+
+
Использование
различных моделей
элементов системы
кровообращения
Q
t
q
kid
q
stat
exp( p p
Renal outflow qkid
Kidne
y
stat
)

32.

Обозначения
S(t,x) –площадь сечения сосуда
u(t, x) –скорость течения крови
p(t,x) -давление
Q(t,x) – поток крови (Q=Su)
t - время
x – локальная пространственная координата
- плотность крови ( = const).
Локальная
координата
u(t,x)
x
Сосуд
L
D(t,x), S=D2 /4

33.

Модель гемодинамики в одном сосуде
Одиночный сосуд рассматриваем как трубку кругового сечения, протяженную по
сравнению со своими поперечными размерами. Под эластичностью стенок
понимается возможность изменения сечения сосуда под действием давления.
x- координата вдоль оси сосуда, t - время,
S(x,t) - площадь сечения сосуда, u(x,t) скорость движения крови (направлена
вдоль оси сосуда), p(x,t) - давление в
крови, =const - плотность крови, L- длина
сосуда, F(x,t) объемная плотность силы.
S(x2)
S(x1)
X
X2
X1
В основу описания движения крови в сосуде положены законы
сохранения массы и импульса (количества движения).
Интегральный вид уравнений
x2
t2
2
2
S
(
x
,
t
)
dx
S
(
x
,
t
)
u
(
x
,
t
)
dt 0
t
x
t
x
1
x1
x2
S=S(p) – Эмпирическая
1
t2
t1
x2
t2
t2 x2
t2 x2
t1 x1
t1 x1
2
Su
dx
u
S
dt
pS
t
x dt pdS FSdxdt
t2
2
1
x1
t1
S S ( p)
x1
x
1
t1
0 x L, t 0
зависимость площади
сечения сосуда от
давления.

34.

Последовательность математических
моделей
Модели кровеносных сосудов
Сосуд
S=const
Жесткая трубка
S const
Эластичная трубка
S=S(p)
S=S(p,u,Q, … )
Q(S,p,u)=const
Диаметр сосуда может быть постоянным или переменным и может
зависеть от любого числа физических и физиологических
параметров, таких как давление, коэффициент эластичности,
гравитация и т.д. Эту зависимость называют
“уравнением
состояния” . Кроме того, считаются стенки сосуда тонкими.

35.

Уравнение состояния
Характерный экспериментальный вид S(p)
S
Smax
Простейшее приближение
S
Smax
Smin
pmin
dS
0
dp
Smin
pmax
p
S ( p) Smax
p
S ( p) Smin
p
Нелинейное приближение
s
s m ax
P
Pmax
Pmin
=8
=2
=1
= 0 .5
smax smin
p p
s s
arctg
pmax pmin
s
s m in
p
p

36.

Типичное «уравнение состояния» сосуда
S
В области нормальных
значений давления
зависимость площади
поперечного сечения
сосудов от давления
близка к линейной.
Smax
Smin
pmin
pmax
Простейшая математическая модель
“уравнения состояния”
S
S
S ( p) Smin max min ( p pmin )
pmax pmin
p
S
Smax
Smin
pmin
pmax
p
Зависимость площади поперечного сечения от давления для большинства
сосудов может моделироваться этим уравнением как в нормальном состоянии,
так даже и в случае некоторых патологий. Заметим, что параметры Smin , Smax, Pmin,
Pmax могут быть функциями времени t и пространственной переменной x.
Например, стеноз, вызванный атеросклерозом, может быть промоделирован
таким образом.

37.

Пример специального «уравнения состояния»
Q
Эффект Остроумова-Балисса
в церебральных артериях
Q const
pmin
p0
p1
pmax
Необходимо иметь в виду,
что различные формы
уравнений состояния могут
порождать специфические
математические проблемы.
p
Этот эффект можно
моделировать уравнением
состояния следующего вида
S S ( p),
Q0
S (t , x) S ,
u
S S ( p),
pmin p p0
p0 p p1
p1 p pmax

38.

Элементы модели
1. Система сосудов (вся сердечнососудистая система или ее часть) граф сосудов.
2. Вершины графа :
- области бифуркаций
сосудов;
- мышечные ткани;
- отдельные органы.
Сосуды - ребра графа
вершины
Ткань, мышцы
Почка, печень,
кишечник, селезенка,
легкие, …
•Области бифуркации сосудов моделируются законом сохранения потока
вещества и условием непрерывности давления или интеграла Бернулли.
Области фильтрации крови через ткани моделируются законом
сохранения потока вещества и законом фильтрации Дарси.
•Каждый отдельный орган описывается специальной моделью, которая в
простейшем случае является точечной.
•Модели вершин могут быть любого уровня сложности.

39.

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

40.

Элементы модели
QA
QV
3. Сердце описывается
двух или четырех
камерной моделью.
предсердие
желудочек
«Двухкамерная» модель сердца состоит из
двух элементов- предсердия и желудочка и
работает как насос. В течение систолы
кровь из желудочка поступает в аорту. Этот
процесс регулируется набором
параметров: ударный объем, текущий
объем предсердия и желудочка, давление в
аорте и т.д. В течение диастолы желудочек
наполняется.
Большой круг
кровообращения
легкие
«Четырех камерная» модель объединяет
две «двухкамерные» модели со своими
параметрами.
p
Пример «двухкамерной» модели сердца
1
(qtop qbot )(t t sys ) 2 , 0 t t sys
q
top
2
t sys
Q(t )
q
, t sys t tdias
bot
t
Vk VKD QA ( )d ,0 t tsys
0
t
pbot
Vk VKS QV ( )d , tsys t tdias
t sys
mm Hg
ptop
tsys
tdias
t, сек

41.

Неконсервативная модель
Модели сердца
График функции сердечного выброса приближает экспериментальную кривую потока или
давления
Q
t S – продолжительность систолы
t D – продолжительность диастолы
Vsurg – ударный выброс
VKD , VKS – объемы желудочка в конце диастолы и систолы
Консервативная модель
QV
предсердие
0
tS
(VK – текущий объем крови в желудочке)
QA
Желудочек
VK
3t
Vsurg (4t
QA (tS , t D )
tS
0, t t t
S
S
2
)
1
, 0 t tS
tS2
tD
t
VKD QA ( )d , 0 t tS
0
Vk
t
V Q ( )d , t t t t
KS
S
D
S
V
0
tD
t
заданная
параметрическая
функция выброса
вычисляемый
текущий объем
желудочка
Поток крови QV, поступающий в предсердие, определяется сердечным
выбросом и потоком крови во всей системе. Такая модель позволяет
сохранять объем циркулирующей крови, исследовать механизмы
регуляции.
Регуляция по величине сердечного
Регуляция по времени систолы и
выброса по конечнодиастолическому
диастолы : Vsurg=const
объему VKD : Vsurg= Kf VKD

42.

Математическая модель на графе
1. Каждому ребру графа сопоставлена система уравнений гемодинамики
S uS
0
v
t x
u u 2 1 p
( )
FT FTP
t x 2
x
2. Каждой вершине графа,
соответствующей области бифуркации
z
сосудов, сопоставлено уравнение
неразрывности и условие непрерывности
интеграла Бернулли.
3. Каждой вершине графа,
соответствующей тканям, сопоставлено
уравнение сохранения вещества и
уравнение фильтрации Дарси.
z i si u i z j s j u j 0 , z i si u i k d ( pi p j ) ,
Kd- коэффициент фильтрации
zi si ui 0,
i
s s( p)
.
dS
0.
dp
2
u i2 pi u j p j
, i j
2 2
i,j –номера всех ребер, соединенных с
каждой вершиной бифуркации, zi=±1
ОБОЗНАЧЕНИЯ
S(t,x) –площадь поперечного сечения
u(t, x) -скорость потока крови
p(t,x) -давление
t - время
x - локальная пространственная
координата
- плотность крови ( = const).
FT – сила трения
FT – сила тяжести

43.

Свойства уравнений гемодинамики ( ГД )
Система уравнений гемодинамики ( ГД )
имеет гиперболический тип при условии,
что для уравнения состояния выполнено
неравенство dS/dp>0. В этом случае
имеется
две
характеристики,
два
соотношения на характеристиках, скорость
распространения малых возмущений.
Так как давление в
кровеносной системе
мало отклоняется от
своего среднего
значения, в ряде
случаев поведение
системы
удовлетворительно
описывается
линеаризованными
уравнениями
гемодинамики ( ЛГД ).
Характеристики
Скорость
малых
возмущений
Соотношения на
характеристиках
s us
0 ,
t x
u u 2 p
( ) FT FTP ,
t x 2
s=s(p)
ds
0
dp
dx
u с,
dt
s
с
,
u
c
M 1
2
d( s)
du ( FT FTP )dt

44.

4. В граничных вершинах графа должны быть заданы краевые условия. Это,
например, самосогласованная модель сердца или некоторые законы изменения
потоков или давления.
Вход ( Q(t) )
Граничные вершины
Выход( p(t) )
Пример графа сосудов головного мозга
5. При моделировании процессов переноса растворенных в крови
веществ в вершинах графа формулируются дополнительные уравнения,
описывающие эти процессы.

45.

Перенос веществ по графу сосудов
Перенос вещества с массовой концентрацией C с учетом диффузии (D=constкоэффициент диффузии) по сосуду с переменным сечением описывается
дифференциальным уравнением
C
C D
C
u
(S
)
t
x S x x
число Маха М<<1 ,
сечение сосудов
ограничено снизу
C
C
2C
u
D
t
x
x 2
Уравнения переноса вещества для графа
Cl , k
t
uk
Cl , k
x
D
2Cl , k
x 2
Cl ,i
z
S
C
u
D
i i l ,i i
x
i Ш ( m )
fl , k ( x, t ) k (C1, k ,..., ClC , k )
f m 0
Сl,k – концентрация l-ого вещества в k –ом
сосуде, f – источник , - химические
реакции.
В вершинах ветвления - закон Киргхофа с
учетом диффузии и равенство
соответствующих концентраций.
C l , k , i C l , k , j , i j , i , j Ш ( m)
Уравнения для концентраций в сердце
Рассматривается случай полного перемешивания веществ Cl в желудочке. Позволяют
рассматривать замкнутую систему кровообращения с сохранением количества веществ
во всей системе.
Cl (t t )
t t
t
( S a Cl ,a u a DS a
Cl ,a
)dt Cl (t )V (t ) S a u a dt V (t )
x
t
t t
1
Здесь t – время начала диастолы, а t изменяется от нуля до величины равной
продолжительности диастолы. Индекс a обозначает аорту. Аналогичное выражение
выписывается для концентраций во время систолы.

46.

Моделирование барорецепторной нейрогенной регуляции
Принципиальная схема нейрогенной регуляции
Повышение
артериального
давления
Увеличение
импульсации от
барорецепторов
Снижение тонуса сосудов.
Увеличение заполненности
тканей кровью.
Уменьшение частоты
сокращений сердца.
Реакция
ЦНС
Понижение давления: обратные процессы
Механизм нейрорегуляции настроен на поддержание характерного давления
pbar
в сосуде с барорецепторами, которые реагируют на отклонение
pbar (t ) pсред (t )
1
pbar , где p сред (t )
T
Модель изменения тонуса сосудов
Площадь сечения сосуда
S
S max, 1
S max, 0
S max, 2
t
p( )d
t T
Повышение давления приводит к
увеличению площади сечения и к
уменьшению жесткости
p max
pmax,0 1,5 pr ( pmax,0 pmin,0 )
p max,0
pmax,0 1,5 pr ( pmax,0 pm in,0 )
S min, 2
Нормальный тонус
p max, 1
Пониженный тонус
p max, 0
p max, 2 p
Повышенный тонус
,
– среднее давление в сосуде.
S min,1
S min, 0
p min, 2 p min, 0 p min,1
pbar (t )
0
pbar

47.

Моделирование барорецепторной нейрогенной регуляции
Модель изменения наполненности капилляров (тканей)
Повышение (уменьшение) давления приводит к увеличению (уменьшению)
количества капилляров в тканях, заполненных кровью. В рамках точечных
моделей это можно интерпретировать как увеличение (уменьшение)
коэффициента фильтрации в законе Дарси uS K D ( p1 p 2 )
pbar
(1 tis )k D,0 ,
p
k D (1 tis bar )k D,0 , pbar
(1 tis )k D,0 ,
pbar
kD
(1 tis )k D,0
k D,0
(1 tis ) k D, 0
0
pbar
Модель изменения частоты сокращения сердца
Повышение (уменьшение) давления приводит к увеличению (уменьшению)
продолжительности сердечного цикла tprd .
t prd
pbar
(1 prd )t prd ,0 ,
pbar
(1 prd
)t prd ,0 , pbar
pbar
(1 prd )t prd ,0 ,
t p rd
(1 prd )t prd , 0
t prd , 0
(1 prd )t prd , 0
0
pbar

48.

Моделирование барорецепторной нейрогенной регуляции
P, мм рт.ст.
Давление в аорте
После кратковременного
повышения артериального
давления нейрогенная
регуляция приводит к
возврату давления в норму.
Снижается среднее
давление в аорте и в
артериях руки.
250
210
170
130
90
50
10.45
13.65
16.85
20.05
Расчет А
Расчет С
23.25
время, с
Расчет А – течение без регуляции, B – течение с частичной регуляцией, С – течение с
полной регуляцией
Среднее давление в аорте
P, мм рт.ст.
150
140
130
120
110
100
90
80
10,45
13,65
16,85
Расчет А'
20,05
Расчет В'
23,25
26,45
Расчет С'
Среднее давление на ткани руки
29,65
32,85
время, с
P, мм рт.ст.
55
50
45
40
35
30
10,45
13,65
16,85
Расчет А'
20,05
23,25
Расчет В'
26,45
29,65
Расчет С'
32,85
время, с

49.

Моделирование гравитационного воздействия
Исследуется влияние гравитации на
гемодинамику человека.
g
Моделирование на полном графе (большой
круг кровообращения + церебральное
кровообращение) позволяет установить
изменения в гемодинамике, например, в
случае быстро меняющейся гравитации.
Объем кровотока в мозге сильно
падает под действием нарастающей
гравитации.
Меняется распределение потоков
крови по отделам головного мозга.

50.

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