CC BY
96
УДК 614.89
ИМИТАЦИОННОЕ МОДЕЛИРОВАНИЕ ВНЕШНЕГО ДЫХАНИЯ ЧЕЛОВЕКА
Г.Р. Камалетдинова1, П.М. Оневский2, В.А. Погонин2
ОАО «Научно-исследовательский и конструкторский институт химического машиностроения», г. Москва (1); кафедра «Информационные процессы и управление», ФГБОУ ВПО «ТГТУ» (2); onev1@mail.ru
Представлена членом редколлегии профессором Н.Ц. Гатаповой
Ключевые слова и фразы: газовая дыхательная смесь; искусственные легкие; концентрация диоксида углерода; потребление кислорода; система управления.
Аннотация: Рассмотрено моделирование функционирования системы управления стендом «Искусственные легкие» в пакете MatLab. Заданная концентрация диоксида углерода в легких при выходе на режим испытательного стенда достигается за счет управления частотой и глубиной дыхания. Потребление кислорода имитируется сбросом части газовой дыхательной смеси в атмосферу на фазе вдоха и подачей соответствующих количеств диоксида углерода и азота.
В настоящее время в России и за рубежом ведется интенсивная разработка технических средств для испытания и исследования средств индивидуальной защиты органов дыхания (СИЗОД) человека. В качестве СИЗОД широкое распространение получили индивидуальные дыхательные аппараты (ИДА) изолирующего типа с химически связанным кислородом, используемые в различных областях в экстремальных ситуациях: на земле и под землей, в космосе и на транспорте, на воде и под водой.
Поскольку испытание ИДА на людях в силу многих причин возможно лишь в ограниченном объеме и в ряде случаев небезопасно, такие испытания целесообразно проводить на специальных установках, называемых «Искусственные легкие» (ИЛ). В настоящее время установка ИЛ должна стать основным инструментом для определения характеристик ИДА и исследования его влияния на характер дыхания человека, что не требует привлечения людей-добровольцев.
Недостатком существующих зарубежных и отечественных установок ИЛ является невозможность изменения формы дыхательной кривой и реализации математическим и программным обеспечением автоматизированной системы управления установки дыхательного коэффициента меньше 1 (то есть при снижении производительности регенеративного патрона ИДА). В существующих установках не предусмотрено также изменение частоты и глубины дыхания для поддержания заданной фиксированной концентрации диоксида углерода в выдыхаемой воздушной смеси при изменении концентрации вдыхаемых газов из ИДА, что реально происходит при испытании ИДА на людях.
Целью данной работы является разработка и исследование Simulink-модели ИЛ, имитирующей внешнее дыхание человека.
Испытательный стенд ИЛ состоит из четырех основных блоков: подачи диоксида углерода и азота; имитации дыхания; имитации потребления кислорода (по массе и объему); управления [1].
Блок имитации дыхания создает пульсирующий поток газовой дыхательной смеси (ГДС), аналогичный потоку, формируемому легкими человека. Блок работает поочередно в режиме вдоха и выдоха.
Аналогично в двух режимах работает блок имитации потребления кислорода путем сброса части ГДС в атмосферу через соответствующие клапаны. Подача смеси диоксида углерода и азота в имитатор дыхания происходит на стадии вдоха.
Для проведения имитационных исследований функционирования стенда при различных психофизиологических состояниях человека необходимо использовать математическую модель потребления кислорода пользователем ИДА.
Основными входными параметрами модели являются: легочная вентиляция Жл, дм3/мин; глубина дыхания Уд, дм3; частота дыхания п, мин4.
Исходный режим для имитационного моделирования:
- глубина дыхания Уд = 1,75 дм3;
- частота дыхания п = 20 мин1;
- подача диоксида углерода ^Со2 (0) = 1,1 дм3/мин;
- коэффициент дыхания Кд = 1;
- объем системы ИЛ Уил = 10 дм3.
Задачей исследования является обеспечение заданной концентрации диоксида углерода на входе в ИДА (то есть на выдохе из ИЛ) и определение кривых дыхания, реализующих данную концентрацию, путем автоматизации процессов управления стендом ИЛ.
Основные соотношения математической модели:
3
- подача диоксида углерода за такт вдоха-выдоха, дм ,
Усо2 (0) = ^со2(0)/п, дм3; (1)
- потребление кислорода за такт вдоха-выдоха
Уо2(0) = ^со2(0)/Кд , дм3. (2)
Состав ГДС в испытательном стенде перед первым циклом вдоха-выдоха
ИЛ = 0,2 %; С™ = 21 %; С™
2 О2
принимался равным атмосферному: С^ИО = 0,2 %; С^1 = 21 %; СИ = 100 -
- СИЛ - сИЛ = 78,8 %. си2 и2
Объемы газов, поступающих в блок имитации дыхания ИЛ для схемы с имитацией потребления кислорода по массе и объему, дм3:
V™ = С™ ГД/100; (3)
СО 2 СО2
^2
ГИЛ = Сатм ГД /100 - Vq2; (4)
"О2 'Д/*"" И2 * С
n2 n2 д/
ГИЛ = Сатм ГД /100, (5)
где ССО , Са™ С^™ - объемные доли газов в атмосферном воздухе, %. На выдохе соответственно имеем:
Свыд (vCQ2(°) + СИЛ2 ГИл/100 + ГСИ0Л2 )'100 ; (g)
С02 VCQ2(0) + Гил + Гд ; ()
(СИЛ Гил/10° + У? ) 1°° ,
Свыд = 4 И2 ИЛ/_И2 ' ; (7)
02 Гсо2(0) + Гил + Гд '
С-Д = 10° - ССИД2 - СИ2Д, (8)
где СС0 , СО , С0ТЬ1Д - объемные доли газов на выдохе, %.
с02 02 ^2
Далее полученные значения С
выд
CÖ2'
C выд CO2 '
Свыд подставляются вместо
N2
Ссо2' СИЛ, С^^Л и так далее для других циклов.
На рисунке 1 представлена упрощенная модель стенда ИЛ, реализованная с помощью пакета моделирования динамических систем Simulink, входящий в состав пакета прикладных программ MatLab.
В блоках fcn и fcnl с помощью блоков с переменным временем задержки Variable Transport Delay реализована модель процесса дыхания, определяемая выражениями (1) - (8), с возможностью изменения формы дыхательной кривой и дыхательного коэффициента.
Система управления процессом дыхания построена по двухконтурной схеме. Первый контур обеспечивает управление заданной концентрацией диоксида углерода на выдохе путем изменения глубины дыхания, второй - изменением частоты дыхания. Влияние каждого контура на процесс дыхания можно изменять путем соответствующих настроек ПИД-регуляторов PID Controller. Ограничения по глубине и частоте дыхания задаются соответственно блоками Saturation и Saturationl.
Динамические свойства измерителя концентрации диоксида углерода (газоанализатора) на выдохе моделируются введением временной задержки т = 60 с в контуры обратных связей с помощью блоков Transport Delay и Transport Delay 1. Динамикой исполнительных механизмов пренебрегаем ввиду их слабого влияния на динамические свойства системы.
0.002
0.21
0.788
Vd Cco2vid
t Co2vid
n Cn2vid
Ceo 2at Vgds
Co2at fcn
Cn2at Wgds
Cn2il Wn2
Co2il Wco2
Cco2il tau
о
Function Function1
Cco2il Cco2vid Co2il fcnlCo2vid Cn2il Cn2vid
«-П
Vd
I
<^>Чрд>
о
0.02993
С C02
г*0.1761
С О^
Sbros GDS
С N2
^0.7939
Podacha N2
Podacha С02
os:
tibka
1.75 Vd nominal
■JH —
25 n nominal
Г?
0.03
С C02 vidox zadannaya
Рис. 1. Simulink-модель стенда ИЛ
На рисунках 2-4 представлены результаты имитационных исследований, которые характеризуют переходные режимы работы стенда ИЛ для различных заданных концентраций диоксида углерода на входе в ИДА.
СС02 0,05
0,045
0,04
0,035
0,03
0,025
0,02
0,015
0,01
0,005
/
\ч /
/ 1 х 4
2 3 /
4
0 100 200 300 400 500 600 Рис. 2. Концентрация СО2 в искусственных легких для различных заданных концентраций СО2
t, С
1 - 0,025; 2 - 0,03; 3 - 0,035; 4 - 0,04
Л
1,85 1,8 1,75 1,7 1,65 1,6 1,55 1,5 1,45 1,4
1
/7
/
/ 3
/ 4
0
100
200
300
400
500
600
t, c
Рис. 3. Глубина дыхания Уд в зависимости от заданной концентрации СО2 на входе в ИДА:
1 - 0,025; 2 - 0,03; 3 - 0,035; 4 - 0,04
п, мин
28
26
24
22
20
18
16
14
1
/
^—
{ /
1/ / 3 ....
IV ! /
V ■L-*
100
200
300
400
500
600
t, c
Рис. 4. Частота дыхания п при различных заданных концентрациях СОг:
1 - 0,025; 2 - 0,03; 3 - 0,035, 4 - 0,04
На рисунке 2 представлены кривые изменения концентрации диоксида углерода при выходе на режим испытательного стенда для различных заданных концентраций диоксида углерода.
На рисунках 3 и 4 представлены управляющие воздействия, реализующие заданные концентрации диоксида углерода.
Результаты моделирования показывают, что испытательный стенд выходит
на номинальный режим работы (С^О4 = 0,03) менее чем за 5 мин с точностью до
0,1 %, что подтверждается результатами испытаний реального стенда при «ручном» управлении. С увеличением концентрации диоксида углерода на выходе из блока имитации дыхания уменьшается глубина и частота дыхания. Минимальная частота дыхания ограничивалась величиной 15 мин-1.
Таким образом, автоматизация процессов управления стендом ИЛ позволяет повысить точность воспроизведения реальных дыхательных процессов, присущих человеку в различных психофизиологических состояниях. Полученные результаты могут быть использованы при принятии оптимальных проектных решений на всех этапах разработки и сопровождения ИЛ.
Работа выполнена в рамках соглашения №14.В37.21.2083 Федеральной целевой программы «Научные и научно-педагогические кадры инновационной России» на 2009-2013 годы.
0
Список литературы
1. Гудков, С.В. Совершенствование методики испытания изолирующих дыхательных аппаратов с химически связанным кислородом / С.В. Гудков, Д. С. Дворецкий, А.Ю. Хромов // Вестн. Тамб. гос. техн. ун-та. - 2009. - Т. 15, № 3. -С. 589-597.
Simulation of Human External Respiration G.R. Kamaletdinova1, P.M. Onevsky2, V.A. Pogonin2
Scientific Research and Design Institute of Chemical Engineering, Moscow (1);
Department "Information Processes and Control", TSTU (2); onev1@mail.ru
Key words and phrases: artificial lungs; carbon dioxide concentration; control system; oxygen consumption; respiratory gas.
Abstract: The paper examines the simulation of control over "Artificial lungs" facility in MatLab package. A given concentration of carbon dioxide in the lungs in testing mode is achieved by controlling the rate and depth of breathing. Oxygen consumption is simulated by respiratory gas discharge in the atmosphere in inspiratory phase and supply of the corresponding quantity of carbon dioxide and nitrogen.
Imitationsmodellierung der äußerlichen Atmung des Menschen
Zusammenfassung: Es ist die Modellierung des Funktionierens des Systems der Steuerung vom Stand «Künstliche Lungen» im Paket MatLab betrachtet. Die aufgegebene Konzentration des Kohlendioxides in den Lungen beim Ausgang auf das Regime des Prüfstandes wird auf Kosten von der Steuerung von der Frequenz und der Tiefe der Atmung erreicht. Der Verbrauch des Sauerstoffes wird von der Ableitung des Teiles der Gasatmungsmischung in die Atmosphäre auf der Phase des Einatmes und der Abgabe der entsprechenden Zahlen des Kohlendioxides und des Stickstoffes imitiert.
Modélage de simulation de la respiration extérieure de l'homme
Résumé: Est examiné le fonctionnement du système de la commande du stand «Poumons artificiels» dans le paquet de MatLab. La concentration donnée du dioxyde de l'hydrogène dans les poumons lors de la sortie sur le régime du stand d'essai est obtenue compte tenu de la commande de la fréquence et de la profondeur de la respiration. La consommation de l'oxygène est simulée par l'évacuation d'une partie du mélange gazeux de respiration dans l'atmosphère à la phase de l'aspiration et par le débit de la quantité correspodante du dyoxyde de l'hydrogène et de l'azote.
Авторы: Камалетдинова Гузель Ринатовна - научный сотрудник, ОАО «Научно-исследовательский и конструкторский институт химического машиностроения», г. Москва; Оневский Павел Михайлович - кандидат технических наук, доцент кафедры «Информационные процессы и управление»; Погонин Василий Александрович - доктор технических наук, профессор кафедры «Информационные процессы и управление», ФГБОУ ВПО «ТГТУ».
Рецензент: Матвейкин Валерий Григорьевич - доктор технических наук, профессор, заведующий кафедрой «Информационные процессы и управление» ФГБОУ ВПО «ТГТУ», заместитель генерального директора, ОАО «Корпорация «Росхимзащита», г. Тамбов.
