Мы рады представить вам нашу очередную статью из цикла «Моделирование руками и глазами студента», посвящённую разработке контроллера для аппарата искусственной вентиляции лёгких, выполненной в программном обеспечении REPEAT. Автором этой работы является студентка РГУ имени Косыгина Ульяна, которую мы благодарим за её усердие и талант.
В статье описана математическая модель и результаты расчётов, которые показывают, как программное обеспечение REPEAT может быть использовано в области медицинской техники.
Ссылка на телеграм-канал REPEAT: https://t.me/repeatlab
Разработка контроллера для аппарата искусственной вентиляции легких
Гермогенова Ульяна Леонидовна, студент 4 курса, Института мехатроники и робототехники, кафедра Автоматики и промышленной электроники
РГУ им. А.Н. Косыгина, Москва
Назначение контроллера
Контроллер аппарата искусственной вентиляции позволяет настраивать и контролировать процесс вентиляции, чтобы обеспечить требуемую поддержку дыхательной функции пациента. Кроме того, контроллер позволяет наблюдать состояние пациента и адаптировать параметры вентиляции в соответствии с его потребностями [1].
Контроллер необходим для управления и регулирования следующими параметрами вентиляции:
1. Давление.
2. Расход.
3. Частота дыхания.
Для обеспечения работы контроллера требуется выполнить моделирование аппарата искусственной вентиляции с учётом математической модели лёгких пациента.
Описание математической модели
Аппарат искусственной вентиляции состоит из следующих частей:
- Центробежный насос (вентилятор)
- Воздушные соединители (шланги).
Модель центробежного насоса
Вентилятор нагнетает давление окружающего воздуха для вентиляции легких пациента и выдувает воздух с объёмным расходом Qout и давлением pout. Его источником является управляющий сигнал pcontrol.
Характеристики вентилятора назначены для установившегося режима. В REPEAT вентилятор моделируется как динамическая система с инерцией по следующей формуле:
где ωn = 2πn, n = 30 об/сек, а коэффициент затухания ζ = 1, соответствующий реальным экспериментальным данным, – комплексная переменная [2].
Для реализации модели центробежного насоса на платформе REPEAT используется передаточная функция с заранее вычисленными постоянными коэффициентами числителя и знаменателя.
Модель воздушных соединителей (шланги)
Воздух, выходящий из вентилятора, проходит через шланг, прежде чем попасть в лёгкие пациента. Шланг для утечки обеспечивает возможность выдыхаемому воздуху выйти из воздушных разъемов в условиях окружающей среды (pamb = 0). Также шланг для утечки допускает свежий воздух при вдохе. Давление на конце воздушных соединителей, т. е. давление, подаваемое в дыхательные пути легких, составляет paw. Объемный расход воздуха в легкие пациента составляет Qpat, а расход утечки Qleak.
Потери в шлангах моделируются линейными функциями сопротивления Rhose и Rleak для длинного шланга и шланга утечки, соответственно. Падение давления на этих участках вычисляется пропорционально объемному расходу, как показано ниже:
С учетом принципа сохранения потока выходной поток , поток пациента и поток утечки связаны следующим образом:
Модель лёгких пациента
Легкие реагируют на объемный приток/выход воздуха (Qpat) за счёт изменения легочного давления и за счёт объёмного расширения/сжатия.
Модель лёгких пациента последовательно состоит из следующих элементов:
- Rlung сопротивление для учёта линейных вязких потерь;
- Clung податливость для учёта эластичности легких.
Уравнения, описывающие динамику лёгких, представлены ниже:
где pext внешнее давление на лёгкие, которое принимается равным нулю. Clung рассчитывается по следующей формуле:
Таблица данных, по которым будут вычисляться и вводиться значения, представлена ниже (Таблица 1).
Стадия моделирования
Для ввода значения входного давления на аппарате искусственной вентиляции лёгких используем блоки “Константа”, “Время моделирования” и “Jython” (см. Рисунок 2). В константы соответственно задаются значения целевого давления плато, положительного давления в конце выдоха, целевого времени подъёма, времени дыхательного цикла, соотношения инспирации и экспирации.
Код блока “Jython” представлен ниже (см. Рисунок 3):
Для моделирования вентилятора используется блок “Передаточная функция” с заранее вычисленными коэффициентами из библиотеки “Автоматика” (см. Рисунок 4).
Для комбинированной модели для шланга и легких создан блок пространства состояний во внешнем проекте. Для его моделирования были использованы блоки “Усилитель”, “Сумматор”, “Интегрирующее звено” из библиотеки “Автоматика” (см. Рисунок 5).
Разработанный проект блока пространства состояний представлен ниже (см. Рисунок 6). Параметры блока задаются из значений A, B, C, D, которые вычисляются следующим образом:
Для ПИД-регулирования используется блок “ПИД-регулятор” и единичное усиление. В итоге получаем полную модель на общей схеме (см. Рисунок 7).
Результаты моделирования
Результаты моделирования были сняты с блоков “Усилитель” (номера 16 и 21) и представлены на графиках соответственно (Рисунок 8, Рисунок 9) для каждого параметра.
Регулировка коэффициента переменного усиления
Стратегия регулировки коэффициента переменного усиления может быть разработана на основе величины объемного расхода воздуха пациента Qpat. Основная идея заключается в следующем: когда величина расхода Qpat высока (скажем, выше порога δ – длины переключения), можно использовать контроллер с высоким коэффициентом усиления. Напротив, когда расход низкий, предпочтительнее использовать регулятор с низким коэффициентом усиления, который минимизирует колебания в расходе.
Нелинейный коэффициент усиления переключателя ф может быть описан следующим уравнением:
где α - высокий коэффициент усиления. Вышеуказанный коэффициент усиления переключателя может быть использован для измерения погрешности и передан в ранее описанный регулятор с низким коэффициентом усиления.
Проектируем нелинейный коэффициент усиления переключателя ф(Qpat) как новый подпроект. Схема показана ниже (см. Рисунок 10).
Вышеупомянутая подмодель теперь может быть включена в нашу систему для реализации схемы с коэффициентом усиления. Полная схема показана ниже (см. Рисунок 11).
Таким образом, программное обеспечение REPEAT обладает широким потенциалом для использования и имеет перспективы в области медицинской техники, и мы надеемся, что эта статья послужит стимулом для дальнейших исследований.
Список источников
B. Hunnekens, S. Kamps and N. Van De Wouw, "Variable-Gain Control for Respiratory Systems," in IEEE Transactions on Control Systems Technology, vol. 28, no. 1, pp. 163-171, Jan. 2020, doi: 10.1109/TCST.2018.2871002. keywords: {Lung;Hoses;Ventilation;Mathematical model;Hospitals;Mechanical ventilation;performance;respiratory systems;variable-gain control}.
Преобразование Лапаласа
Узнайте больше о преимуществах использования REPEAT попробовав демоверсию на нашем сайте!