litceysel.ru 1

Труды МФТИ, т. 1, № 1, 2009г., с. 5-16.


Петров И.Б.

Математическое моделирование в медицине и биологии на основе моделей механики сплошных сред

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

Разработка этого аппарата предполагает:

-построение замкнутой механико-математической модели процесса, описывающей поведение биологической среды на основе системы уравнений в частных производных механики сплошных сред (МСС);

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

-корректную математическую постановку задачи, т.е. представление замкнутой системы МСС, постановка необходимых для ее решения начальных и граничных условий, условий на контактных границах (если они есть);

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

-разработка алгоритма численного решения задачи и его программная реализация;

-численное решение задачи и визуализация полученных результатов.


Разумеется, при исследовании биомедицинских проблем встречаются процессы, для математического описания которых используются аппараты обыкновенных дифференцированных уравнений (ОДУ), систем алгебраических нелинейных уравнений, разностные отображения, теории бифуркций, хаоса и порядка. Примеры успешного использования таких математических аппаратов представлены в [1] для прогнозирования развития болезни, в [2-5] – для решения задач нелинейной динамики в биологии, химической кинетике и др.

Кроме того, при изучении решений некоторых медицинских процессов необходимо численно решать жесткие системы ОДУ, например, при моделировании протекания химических реакций, что представляет собой самостоятельную проблему, которой посвящена обширная литература (см., например, [6-8]).

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

Постановки биологических и медицинских задач, которые приводят к необходимости численного решения систем уравнений в частных производных (МСС, уравнения параболического и эллиптического типа) появились относительно недавно. Они представлены в работах [9-11]. Для численного решения таких задачи использовались методы, ранее применявшиеся для решения задач газогидродинамики [12-16]. Реологические соотношения для биологических сплошных сред разрабатывались в работах [17-19].

Круг рассматриваемых задач их этой области достаточно широк.

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


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

Исследование распространения импульсов Пуркинье (структуры кабельного типа, проводящие импульсы от предсердий до желудочков в сердце) проводится в работах [23,24]. Основываясь на численном эксперименте с использованием математической модели МакАллистера-Нобла-Тасена, описывающей динамику распространения импульсов в волокнах Пуркинье, авторы этих работ обнаружили различные режимы распространения этих импульсов. В частности оказалось, что при различных режимах возбуждения возможно появление таких явлений, как отражения, автоколебания, аннигиляция, солитоноподобные режимы. Интересно, что подобные явления имеют место в нервных волокнах, возбуждение импульсов в которых описывается в рамках известной математической модели Ходжкина–Хаксли [25,26].

Описание простейших математических моделей работы систем кровообращения и сердца можно найти, например, в [27,28]. Функции кровеносной системы человека, которая состоит из малого и большого кругов кровообращения, очень важны и разнообразны, поэтому их моделирование, как в нормальных, так и в патологических условиях, представляет одну из важнейших задач медицины.

По-видимому, на сегодняшней день наиболее адекватными представляются динамические модели пульсирующих течений несжимаемой жидкости в системе растяжимых трубок. В работах [29-31] использовалась квазиодномерная гидравлическая модель несжимаемой жидкости в деформируемом кровеносном сосуде переменного сечения, обобщенная на случай иерархической ветвящейся системы кровеносных сосудов факториальной структуры. Подобная же иерархическая (или сетевая) квазиодномерная модель использовалась и для математического описания работы дыхательной системы на участке трахея – бронхи [31,33].


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

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

По инициативе нейрохирургов Главного военного клинического госпиталя им. Бурденко и Института скорой помощи им. Склифасовского была поставлена задача о расчете последствий черепно-мозговых травм. Надо сказать, что в рамках модели сплошной среды реальным представляется постановка задач о расчете травм костей, грудной клетки, суставов, появлении гематом в мягких тканях тела, что также рассматривалось в работе. Экспериментальным данным, описывающим последствия черепно-мозговых травм, посвящены хорошо известные в нейрохирургии работы [34,35]. Изучение областей повреждения мозга, при различных динамических воздействиях на него, путем математического моделирования в рамках МСС, может дать важную информацию о последствиях травм. В газете “Известия” от 13.08.2003 г. опубликована статья “Бомба для мозгов”, в которой утверждается, что в настоящее время составлена трехмерная компьютерная карта мозга человека, над которой в течение нескольких лет в шести странах мира работали 20 тыс. нейрохирургов и невропатологов. Данные для этой карты получены в результате сканирования мозга 7 тыс. человек. Идея такой деятельности крайне важна, поскольку теперь с высокой степенью точности известно, какие области мозга управляют теми или иными его функциями. Понятно, что определив области повреждения мозга при различных динамических нагрузках на него (или при черепно-мозговых травмах), мы сможем определить те области функции организма, которые могут оказаться поврежденными или измененными. В нейрохирургической практике хорошо известен следующий факт: области поражения мозга при черепно-мозговой травмах не всегда совпадают с областями, прилежащими к месту удара. Примером этому является феномен “противоудара”: при ударе затылком область повреждения мозга локализуется в лобной части головного мозга. Объяснение этому явлению можно дать только путем проведения численного исследования сложнейших волновых процессов, образующихся в неоднородной механической конструкции, которую представляет собой систему череп-мозг (разумеется, речь идет только о механической системе, о принципах функционирования мозга, как об управляющей системе, речь не идет). Кроме того, механико-математические модели реакции рассматриваемой механической системы на динамические нагрузки можно использовать как инструмент для изучения различных параметров травмы (место приложения удара, сила удара, геометрия головы, атрофия мозга и др.) на вероятную степень риска. Эти задачи ставились в работах [36,37]. В [38] численно изучали процесс распространения ударного импульса через 5-слойную биоконструкцию черепа человека, которая, как показали расчеты, оказалась идеальной противоударной защитой для мозгового вещества при естественных динамических воздействиях (разумеется, при действии боевых поражающих средств необходимо использовать дополнительные средства защиты головного мозга).


В работе [39] приводятся примеры расчетов многослойных защитных конструкций, которые могут использоваться для защиты как техники, в которой находится человек, так и тела (бронежилеты, каски, др.). Разработанные авторами вычислительные методы и алгоритмы позволяют изучать также динамические процессы, происходящие в теле человека, при воздействии динамических нагрузок. Такая задача оказывается весьма актуальной, поскольку появляется возможность априорного предсказания областей и степеней поражения тела при каких-либо опасных нагрузках. Области приложения подобных задач очевидны (кроме медицинской): спортивные травмы, травмы при производстве опасных видов работ в экстремальных условиях (например, горные или взрывные работы), боевых и полицейских операциях, автомобильных, авиационных, иных техногенных катастрофах.

К проблеме математического моделирования травматологических процессов относится и задача о залечивании ран. Численному изучению этого процесса посвящены работы [40,41], в которых получено количественное описание динамики залечивания резаной раны кожного покрова человека.

Важнейшей областью в травматологии является проблема математического моделирования движения ног человека при ходьбе, с целью построения ортопедических протезов, имитирующих их движение. Авторы работы [42] не только строят такие модели, но и реализуют их. Для этих целей они использовали аппарат ОДУ. Однако для моделирования распределения динамических нагрузок и деформаций при перемещении полной стопы, уже требуется использование аппарата дифференциальных уравнений в частных производных, в частности системы уравнений МДТТ. Создание подобных моделей для нужд травматологии и ортопедии представляется новой и актуальной задачей вычислительной медицины.

Перспективным направлением вычислительной медицины представляется компьютерная реализация виртуальных хирургических операций и предсказания их последствий. Разумеется, это очень сложное направление, которое только начинает появляться, а постановка некоторых задач до конца не ясна (например, моделирование работы хирурга с помощью скальпеля). Однако реализация некоторых виртуальных операций представляется реальной задачей. Так, в работе [43] представлено численное моделирование операций литотрипсии (дробление почечных камней акустическим волнами, инициируемыми искровым разрядом или лазерным импульсом). Цель таких исследований – найти режимы работы литотриптора (длительность и интенсивность импульса, количество импульсов), при которых фрагменты разрушенного камня были бы достаточно малыми для их выведения из организма естественным путем. Для этого численно исследовалась картина распространения акустического импульса в теле и в камне, а также решалась задача его разрушения.


Другой пример – моделирование офтальмологической операции экстракции (удаления) катаракты. Суть операции в том, чтобы с помощью лазера или ультразвукового факоэмульсификатора разрушить помутневший хрусталик (точнее его плотное ядро) так, чтобы не повредить сетчатку и роговицу глаза и вывести мутные хрусталиковые массы. Эта задача условно подразделяется на три части: первая – расчет импульсного воздействия на хрусталик, вторая – распространение акустического импульса в стекловидном теле до сетчатки и расчет динамического воздействия на нее (поскольку в результате взаимодействия импульса с сетчаткой последняя может расслаиваться), третья – вымывание мутных хрусталиковых масс из передней камеры глаза. Численное решение последней задачи, которое сводится к решению уравнения Пуассона, в двухмерном случае рассматривалось в работе [44], причем с использованием как прямоугольных, так и треугольных расчетных сеток. Расчеты сводились с целью определения застойных зон в передней камере глаза при вымывании мутных масс и оптимизации рабочих режимов хирургических инструментов. Проблема распространения импульса через хрусталик и стекловидное тело к сетчатке представляет сложную динамическую задачу, поскольку глаз представляет собой неоднородную механическую систему с рядом поверхностей раздела сред. Численное ее решение представлено в работе [45]. Моделирование этого процесса сводится к решению динамической системы уравнений гиперболического типа (уравнений МДТТ). Эти работы выполнялись в тесном контакте с сотрудниками Центра микрохирургии глаза им. академика С.Н. Федорова, который является соавтором этих работ.

Важным приложением вычислительной медицины представляются проблемы предсказания динамики развития онкологических заболеваний, т.е. развития опухолей, в том числе с учетом кровообращения. Для их численного решения используются уравнения гидродинамики, жесткие системы ОДУ, уравнения типа реакция-диффузия. Численное моделирование этих процессов проводится в работах [46,71].


Отметим, что с помощью нелинейных уравнений параболического типа (реакция-диффузия) проводится также и численное моделирование процессов структурообразования в активных биосредах, колониях бактерий, микроорганизмов (например, Esherichia coli, Distyostelium discoicleum). Этим задачам посвящены работы [47-50], в которых численно решаются двух и трехмерные динамические задачи об образовании таких структур. Задачи структурообразования при свертывании крови и тромбообразовании рассматривались в работах [50,51].

Численно-математический аппарат для решения рассматриваемого класса биомедицинских задач является аппаратом решения уравнений в частных производных эллиптического, параболического и гиперболического типов, в основном, нелинейных. Кроме того, использовались системы уравнений МСС, в частности, гидродинамики несжимаемой идеальной или вязкой жидкости и МДТТ; чаще всего это линейно и нелинейно-упругие, а также вязкоупругие среды (модель Максвелла).

Для решения задач гидродинамики используются уравнения Эйлера или Навье – Стокса, для задач МДТТ – система уравнений Ламэ (в напряжениях или деформациях).

Эти системы можно найти, например, в монографиях [52-54].

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

Численные методы решения эллиптических уравнений (уравнения Лапласа или Пуассона) подробно рассматриваются в известных монографиях [56-59]. Среди методов, представленных в литературе для решения этого класса задач, отметим наиболее широко используемые: в вычислительной практике – верхней релаксации, трехслойный метод Чебышева, метод Дугласа–Ренфорда, попеременно-треугольный метод, многосеточный метод Р.П. Федоренко. В работе [60] предлагаются монотонные разностные схемы для численного решения этого типа уравнений на нерегулярных сетках. С использованием этого подхода рассматривалась задач о гидродинамическом течении в передней стенке глаза при экстракции катаракты.


Широкий класс процессов в биологии и медицине моделируется при помощи нелинейных уравнений параболического типа (реакция-диффузия). Исследованию свойств решений уравнений этого типа посвящены работы [4,61,62 и др.]. Для численного решения таких уравнений обычно используются двух и трехслойные разностные схемы типа Кранка – Никольсона, Нумерова, Толстых с последующим проведением итераций по нелинейности. Эти вычислительные методы описали в [12,58,59].

По-видимому, наибольшие трудности при численном решении рассматриваемых задач вызывают системы уравнений в частных производных гиперболического типа (МСС). Разработке вычислительных методов решения систем этого типа посвящены работы [63-67, 69, 70]. Наиболее эффективными для численного решения этого класса задач являются сеточно-характеристические методы [65-67,69,70]. Эти методы оказываются не только наиболее точными, но и позволяют наиболее корректно строить вычислительные алгоритмы на границах области интегрирования и границах раздела сред. Одним из наиболее важных показателей качества численного решения является его близость к точному вблизи областей с большими градиентами решений. Это вызвано немонотонным, в одномерном случае, поведением численных решений, получаемых с помощью т.н. немонотонных схем (в линейном случае это, в соответствии с теоремой С.К.Годунова, схемы порядка точности выше первого) и “размыванием” разрывов схемами первого порядка точности. По этой причине многие работы, посвященные этой теме, направлены на регуляризацию решений уравнений гиперболического типа, (см., например, [59,60,66,70]). В работе [70] приводятся определения и критерии монотонности разностных схем для этого типа уравнений в частных производны. В приложении к задачам медицины эти методы использовались при изучении процессов литотрипсии, офтальмологических операциях, черепно-мозговых травм, деформации сердечной мышцы [22,37,43,44,45].

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


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

На рис. 1 представлена схема черепно-мозговой системы, которая выбиралась для моделирования последствий черепно-мозговых травм. Расчетные сетки (в горизонтальной проекции), как четырехугольные, так и треугольные показаны на рис. 2 (соответственно, a,b,c).

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

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

Сравнение расчетной и полученной при томографических исследованиях областей поражения головного мозга видно на рис. 4 a,b (томография проведена в Главном военном клиническом госпитале имени Бурденко).

На рис. 5 показана схема глаза человека, которая использовалась при моделировании динамических процессов, происходящих при проведении офтальмологических операций по экстракции катаракты при помощи лазера или факоэмульсификатора (ультразвуковой инструмент).

Расчетные сетки, поля скоростей и области возможных поражений глаза при лазерном разрушении хрусталика приведены на рис. 6 (расчеты проводились сеточно-характеристическим методом). Расчет давления в хрусталике глаза при лазерном воздействии и расчетная сетка при использовании метода конечных элементов приведен на рис. 7

Результаты операции по разрушению почечного камня (процесс литотрипсии) показаны на рис. 8a,b. Характер разрушения камня показан на рис. 8b (вертикальный срез дискообразного камня в почке и окружающей среде).


Схематичное изображение сердца человека (в разрезе) представлено на рис. 9.

Расчетная сетка, используемая при моделировании процесса начала сжатия сердца (систолы) и начального выталкивания крови, а также поля скоростей в сердце, дана на рис. 10.

Столкновение двух импульсов в волокнах Пуркинье в автоколебательном режиме в четыре момента времени показаны на рис. 11 a,b,c,d (в расчетах были получены также режимы аннигиляции и солитоноподобные режимы).

Динамика залечивания кожной раны (пореза) представлена на рис. 12 a,b. Здесь показаны трехмерные картины распределения плотности коллагена, являющиеся “подложкой”, на которой растут клетки кожи, в два момента времени: t=0 (1 в области поражения и 1 на неповрежденной коже; -плотность коллагена); t=40 дней (плотность коллагена заметно повышается внутри раны вследствие процесса хемотоксиса). Величина плотности коллагена представлена светлой частью поверхностей, изображенных на рисунке.

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

Спиралевидное двух- и трехмерное распространение волн возбуждений в нелинейной активной среде, в частности, в колонии одноклеточных организмов Dictyostelium discoideum, представлено на рис. 14 и 15. Особенность этих микроорганизмов состоит в том, что благодаря своей подвижности и хемотоксису, в определенных внешних условиях (при недостатке субстрата) амебы объединяются в псевдоорганизм, который благодаря согласованному движению составляющих его клеток ведет себя, как одно целое. В активной среде, сформированной этими клетками, образуется спиральная волна ЦАМФ (циклический аденозинмонофосфат) с центром в точке агрегации.


Для численного решения также использовалась система уравнений параболического типа как в двух, так и в трехмерном случаях. Начальные распределения концентрации ЦАМФ и ее дальнейшее “закручивание” в спираль, в двухмерном случае показано на рис. 14. Формирование трехмерное волны в активной среде с распределенным вдоль оси коэффициентом возбудимости, показано на рис. 15. Понятно, что вихреобразная волна в передней части закручивается быстрей, чем в задней, из-за изменения этого коэффициента по координате, что приводит к ее неравномерной закрутке.

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

Автор выражает благодарность за внимание к циклу работ, о котором шла речь в статье действительному члену РАН О.М. Белоцерковскому, члену-корреспонденту РАН А.С, Холодову, профессору Г.Г. Малинецкому, профессору А.А. Полежаеву, профессору А.И. Лобанову, профессору Г.Ю. Резниченко.

СПИСОК ЛИТЕРАТУРЫ


  1. Г.И. Марчук. Математические модели в иммунологии. — М.: наука, 1985 г., С.240.

  2. А.Ю. Лоскутов, А.С. Михайлов. Введение в синергетику. — М.: Наука. 1990, С. 270.

  3. Г.Ю. Резниченко. Лекции по математическим моделям в биологии. — Ч. 1, Москва-Ижевск. Научно-издательский центр “Регулярная и хаотическая динамика”, 2002, С. 231.

  4. Т.С. Ахромеева, С.П. Курдюмов, Г.Г. Малинецкий, А.А. Самарский. — Нестационарные структуры и диффузионный хаос. М.: Наука. 1992, С. 541.

  5. Г.Г. Малинецкий, С.П. Курдюмов (отв. ред.). Новое в синергетике. Взгляд в третье тысячелетие. — М.: Наука, 2002, С. 482.

  6. Ю.В. Ракитский, С.М. Устинов, И.Г. Черноруцкий. Численные методы решения жестких систем. — М.: Наука. 1979, С. 208.
  7. Р.П. Федоренко. Жесткие системы обыкновенных дифференциальных уравнений. В кн. Вычислительные процессы и системы. Вып. 8, под ред. Г.И. Марчука. — М.: Наука, 1991, С. 381.


  8. Э. Хайер, Г. Винер. Решение обыкновенных дифференциальных уравнений. Жесткие и дифференциально-гиперболические задачи. — М.: Мир, 1999, С. 685.

  9. И.М. Макаров (отв. ред.) Информатика и медицина. — М.: Наука, 1997, С. 208.

  10. О.М. Белоцерковский, А.С. Холодов (отв. ред.). — Компьютерные модели и прогресс медицины – М.: Наука. 2001, С. 300.

  11. О.М. Белоцерковский (отв. ред.). Компьютер и мозг. Новые технологии — М.: Наука, 2005, С. 322.

  12. А.А. Самарский. Теория разностных схем. – М.: Наука. 1977, С. 654.

  13. Б.Л. Рождественский, Н.Н. Яненко. Системы квазилинейных уравнений. — М.: Наука, 1978, С. 687.

  14. С.К. Годунов, В.С. Рябенький. Разностные схемы. – М.: Наука, 1973 г., С. 400.

  15. О.М. Белоцерковский. Численное моделирование в механике сплошных сред. — М.: Физматлит, 1994, С. 442.

  16. К.М. Магомедов, А.С. Холодов. Сеточно-характеристические методы. — М.: Наука, 1988, С. 288.

  17. С.А. Регирер. Лекции по биологической механике. – М.: МГУ, 1980.

  18. В.И. Кондауров, А.В.Никитин. Конечные деформации вязкоупругих мышечных тканей.Прикладная математика и механика. 1987, т. 51, вып. 3, С. 443-452.

  19. В.И. Кондауров, Л.В. Никитин. Модель биологически активного вязкоупругого тела. Тематический сборник “Методы расчета изделий из высокоэластичных материалов”. — Рига. 1986, С. 107-108.

  20. K.D. Gosfa, P.J. Hunter, J.M. Pogers, G.M. Gussione, L.K. and A.D. Waldmen. A three-dimensional limite elements method for large elastic deformations of ventricular myocardium. Part I. ASME J.Biomech, Eng. 1996.118 (4), p.p. 452-463.
  21. S.C. Panda and R. Natarajon. Finite-element method of stress analysis in the human left ventricular layered wull structure. Med.Biol.Eng.Comp. 1977, 15 p.p. 67-71.


  22. И.Б. Петров. О численном моделировании биомеханических процессов в медицинской практике. // Информационные технологии и вычислительные системы. 2003, 1-2, С.102-111.

  23. О.В. Асланиди, О.А. Морнев. Эхо в возбудимых волокнах сердца. Математическое моделирование. 1999, т. 11, №9, С. 3-22.

  24. Р.А. Пашко, И.Б. Петров. Моделирование распространения импульсов в волокнах Пуркинье. В пер. сб. Обработка информации и моделирование. — М.: МФТИ, 2002, С. 171-181.

  25. О.В. Асланиди, О.А. Морнев. Могут ли нервные импульсы отражаться? Письма в ЖЭТФ, 1997 , т.65, С. 553-558.

  26. О.В. Асланиди, О.А. Морнев. Новые в пространственно-неоднородных возбудимых средах с рефракторностью: отражение сталкивающихся импульсов возбуждения. // Биологические мембраны. 1997 , т.14, С. 621-625.

  27. А.Н. Ремизов. Медицинская и биологическая физика. — М.: Высшая школа, 1987 г., С. 638.

  28. П.И. Бегун, П.Н. Афонин. Моделирование в биомеханике. — М.: Высшая школа. 2004, С. 389

  29. М.В. Абакумов, И.В. Ашметов, Н.Б. Ешкова, В.Б. Кошелев, С.И. Мухин, Н.В. Соснин, В.Ф. Тишкин, А.П. Фаворский, А.Б. Хруменко. Методики математического моделирования сердечно-сосудистой системы. — Математическое моделирование. 2000 т. 12, №2, С. 106-117.

  30. И.В. Ашметов, А.Я. Буничева, С.И. Мухин, Т.В. Соколова, Н.В. Соснин, А.П. Фаворский. Математическое моделирование гемодинамики в мозге и в большом круге кровообращения. В кн. Компьютер и мозг. Новые технологии. — М.: Наука, 2005, С. 321.

  31. А.С. Холодов. Некоторые динамические модели внешнего дыхания и кровообращения с учетом их связности и переноса вещества. // В кн. Компьютерные модели и прогресс медицины. — М.: Наука, 2001, С. 127-165.
  32. А.В. Евдокимов, А.С. Холодов. Квазистационарная пространственно-распределенная модель замкнутого кровообращения организма человека. В кн. Компьютерные модели и прогресс медицины. — М.: Наука, 2001, С. 164-193.


  33. T.J. Peadley, R.C. Schroter, M.F. Sudllow. Energy loses and pressure drop in models of human airways. Respir Physiol, 1970, №9, pp. 371-386.

  34. А.П. Лебедев, В.В. Крылов. Замечания к патогенезу ушибов мозга, возникающих по противоудраному механизму, в остром периоде их развития. Нейрохирургия, 1198, С. 22-25.

  35. J.H. Adams, D.I. Graham and T.A. Genmarelli. Head injury in man and experimental animals: neuropathalogy. Acta Neuro Chir, 32, p.p. 15-30.

  36. П.И. Агапов, О.М. Белоцерковский, И.Б. Петров. Численное моделирование последствий механического воздействия на мозг человека при черепно-мозговой травме. // Журнал вычислительной математики и математической физики, т. 49, №9, С. 1711-1720.

  37. П.И. Агапов, И.Б. Петров. Расчет повреждений мозга при черепно-мозговой травме. В кн. Компьютер и мозг, Новые технологии (отв. ред. О.М. Белоцерковский). — М.: Наука, 2006, С. 28-38.

  38. П.И. Агапов, А.В. Васюков, И.Б. Петров. Компьютерное моделирование волновых процессов в покровах мозга при черепно-мозговой травме. // В пер. сб. Процессы и методы обработки информации. — М.: МФТИ, 2006, С. 154-163.

  39. И.Б. Петров, Ф.И. Челноков. Численное исследование волновых процессов и процессов разрушения в многослойных преградах. Журнал вычислительной математики и математической физики, 2003 г., т. 43, №10, С. 1562-1579.

  40. L. Olsen, J.A. Sherratt, J.A. and P.K. Maini. A Mechanical Model for Adult Dermal Wound contraction and Permanence of the Contracted Tissue Displacement Profile. J. Theor. Biol. 1995, 177, p.p. 113-128.

  41. Р.А. Пашков. Численное моделирование контракции кожной раны. В сб. Процессы и методы обработки информации. – М.: МФТИ, 2005 г., С. 194-200.
  42. Б.С. Фарбер. А.С. Витензон, И.Ш. Морейнис. Теоретические основы построения протезов нижних конечностей и коррекция движения. — М.: ЦНИИПП, 1995, ч. 2, С. 574.


  43. Д.С. Жуков, И.Б. Петров, А.Г. Тормасов. Численное и экспериментальное изучение разрушения твердых тел в жидкости. — Известия АН СССР, сер. Механика твердого тела, 1991, С. 183-190.

  44. С.Н. Федоров, Э.В. Егорова, А.С. Холодов. А.В. Бубнов. О численном моделировании процессов ирригации и аспирации при экстракапсулярной экстракции катаракты. В кн. Вопросы кибернетики (под ред. О.М. Белоцерковского). — М.: ВИНИТИ, 1982, С. 99-114.

  45. Н.Н. Балановский, А.В. Бубнов, А.С. Обухов, И.Б. Петров. Расчет динамических процессов в глазу при лазерной экстракции катаракты. Математическое моделирование, 2003, т. 15, №11, С. 37-44

  46. A.V. Kolobov, A.A. Polezhaev, G.I. Solyanyk. Stability of shape in pre-angiogenic stage of growth depends on the migration capacity of cancer ulls. In: Mathematical Modeling Computing in Biology and Medicine (Ed.: V.Capasso). Bologna: Progetto Leonerdo, 2003, p.p. 603-609.

  47. E.O. Budriene, A.A. Polezhaev and M.O Ptitsyn. Mathematical modeling of intercullar regulation causing the bacterial colonies. J. Theor Biol. 1998, 135, p. 323-341.

  48. И.Б. Петров. А.А. Полежаев, А.С. Шестаков. Численное моделирование волновых процессов в нелинейных активных средах. // Математическое моделирование, 2000 г., т. 12, №1, С. 38-44.

  49. A.A. Polezhaev, R.A. Paskov, A.I. Lobanov and I.B. Petrov. Spatial patterns formed by chemotactic bacteria Escherichia coli: Internetional Journal of Development Biology, 2006, 50, p.p. 309-314.
  50. А.И. Лобанов, Т.К. Старожилова, Г.Т. Гурия. Численное исследование структурообразования при свертывании крови. Математическое моделирование, 1997, т. 9, №8, С. 83-95.


  51. Г.Т. Гурия, А.И. Лобанов, Т.К. Старожилова. Моделирование роста оторвавшегося тромба в пристеночном потоке. В кн. Компьютерные модели и прогресс медицины. — М.: Наука, 2001 г., С. 250-263.

  52. Л.И. Седов. Механика сплошной среды. – М.: Наука. 1976, т. 1-2, 3-е издание

  53. Ю.Н. Работнов. Механика твердого деформируемого тела. М.: 19179, 632 с.

  54. С.К. Годунов. Е.И. Роменский. Элементы механики сплошных сред и законы сохранения. Новосибирск: Научная книга, 1998, 267 с.

  55. А.А. Самарский, Е.С. Николаев. Методы сеточных уравнений. М.: Наука, 1978 г., С. 590.

  56. Г.И. Марчук. Методы вычислительной математики. М.: Наука, 1989, С. 608.

  57. А.А. Самарский, П.Н. Вабищевич, П.П. Матус. Разностные методы с операторными множителями. — Минск: ЦОТЖ, С. 441.

  58. И.Б. Петров, А.И. Лобанов. Лекции по вычислительной математике. — М.: Интернет-Университет, 2006, С. 523.

  59. Р.П. Федоренко. Введение в вычислительную физику. М.: МФТИ. С. 526.

  60. А.С. Холодов. Монотонные разностные схемы на нерегулярных сетках для эллиптических уравнений в области со многими несвязанными границами. // Математическое моделирование, 1991, т. 3, №9. С. 104-113.

  61. А.А. Самарский (ред.). Компьютеры, модели, вычислительный эксперимент. — М.: Наука, 1988, С. 170.

  62. А.А. Самарский, Н.В. Змитриенко, С.П. Курдюмов, А.П. Михайлов. Эффект метастабильностной локализации тепла в среде с нелинейной теплопроводностью. — ДАН СССР, 1975 г., т. 233, №6, С. 1344-1347.

  63. С.К. Годунов, А.В. Забородин, М.Я. Иванов, А.Н. Крайко, Г.П. Прокопов. Численное решение многомерных задач газовой динамики. — М.: Наука, 1976, С. 400.
  64. А.Г. Куликовский, Н.В. Погорелов, А.Ю. Семенов. Математические вопросы численного решения гиперболических систем уравнений. — М.: Физматлит, 2001, С. 607.


  65. И.Б. Петров, А.С. Холодов. Численное моделирование некоторых задач механики деформируемого твердого тела сеточно-характеристическим методом. Журнал вычислительной математики и математической физики, 1984, т. 24, №5, С. 722-739.

  66. И.Б. Петров, Холодов А.С. О регуляризации разрывных численных решений уравнений гиперболического типа. // Журнал вычислительной математики и мат. физики; т. 24, N 8, 1984, С. 1172-1188.

  67. А.И. Толстых. Компактные разностные схемы и их применение в задачах аэродинамики. — М.: Наука, 1990, С. 230.

  68. В.В Смолянов. Математические модели биологических тканей. М.: Наука, 1980, С. 368.

  69. В.М. Головизин, С.А. Карабасов. Метод прыжкового переноса для численного решения гиперболических уравнений. Точный алгоритм для моделирования конвекции на эйлеровых сетках. Препринт ИБРАЭ РАН №IBRAE-2000-04. — М.: РАН, 2000, С. 40.

  70. А.С. Холодов, А.Я. Холодов. О критериях монотонности разностных схем для уравнений гиперболического типа. // Журнал вычислительной математики и математической физики, 2006, т. 46, №9, С. 1638-1667.

  71. С.А. Астанин, А.В. Колобов, А.И. Лобанов. Влияние пространственной гетерогенной среды на рост и инвазию опухоли. Анализ методами математического моделирования. В кн. Медицина в зеркале информатики. — М.: 2006, С. 163-194.



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






а

b

c

Рис.2. Расчетные сетки в горизонтальной проекции (четырехугольные и треугольные), использовавшиеся при моделировании последствий черепно-мозговых травм.



а b

Рис.3. Области максимальных сжимающих и растягивающих напряжений.



а b

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




Рис.5. Схема глаза человека, которая использовалась при моделировании динамических процессов, происходящих в глазу при проведении офтальмологических операций.




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



Рис.7. Расчет давления в глазу человека при проведении лазерной операции и расчетная сетка.





Рис.8. Результаты операции по разрушению почечного камня (процесс литотрипсии). Представлен характер разрушения камня и изоповерхности давления в теле человека после прохождения через почечный камень ударного импульса.





Рис.9.



Рис.10. Расчетная сетка, используемая при моделировании процесса начала сжатия сердца (систолы) и начального выталкивания крови, а также поля скоростей в сердце.



a b



c d

Рис.11. Столкновение двух импульсов в волокнах Пуркинье в автоколебательном режиме в четыре момента времени.

a b

Рис.12. Динамика залечивания кожной раны (трехмерные картины распределения плотности коллагена в два момента времени).



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







Рис.14. Распределения концентрации ЦАМФ и ее дальнейшее “закручивание” в спираль, в двухмерном случае.





Рис.15. Формирование трехмерное волны в активной среде с распределенным вдоль оси коэффициентом возбудимости.