CC BY
30
УДК 531/534: [57+61]
Российский Журнал
www.biomech.ru
ОПРЕДЕЛЕНИЕ ПАРАМЕТРОВ ПОТЕНЦИАЛА НЕЛИНЕЙНО СЖИМАЕМОГО МАТЕРИАЛА СОННОЙ АРТЕРИИ ЧЕЛОВЕКА ПРИ РАЗЛИЧНЫХ СТАДИЯХ АТЕРОСКЛЕРОЗА
В.Н. Аптуков, Л.Ю. Осоргина
Кафедра математического анализа Пермского государственного университета, Россия, 614990, Пермь, ул. Букирева, 15, e-mail: aptukov@psu.ru
Аннотация. Ранее авторами [1] в рамках линейной модели были даны оценки модулей упругости при различных стадиях атеросклероза, проанализированы особенности напряженно-деформированного состояния стенки сосуда
в окрестности атероматозной бляшки. Был сделан вывод о необходимости изучения данного вопроса в нелинейной постановке с учетом конечных деформаций. В данной статье на основе аналитического решения, предложенного А.И. Лурье, для цилиндра из полулинейного материала Джона представлена методика определения параметров потенциала нелинейно сжимаемого материала, получены конкретные значения параметров потенциала для сонной артерии человека при различных стадиях атеросклероза.
Ключевые слова: биомеханика, склерозированная артерия, нелинейный
потенциал, сжимаемость.
Введение
Известно, что материал стенки кровеносных сосудов человека является нелинейно упругим [4], сжимаемым [5] материалом. Представляет определенные трудности определение констант потенциала нелинейного сжимаемого материала. Обычно это требует проведения специальных экспериментов и решения серии прямых задач деформирования образца. Доступными для авторов экспериментами являлись испытания отрезков артерий в форме замкнутого цилиндра под действием внутреннего давления [2]. Анализ литературы показал, что для цилиндра из однородного нелинейно упругого сжимаемого материала (полулинейный материал Джона) известно решение в конечном виде [3].
Данная статья посвящена разработке методики определения параметров потенциала нелинейно сжимаемого материала на основе решения [3]. Получены конкретные значения параметров потенциала для сонной артерии человека при различных стадиях атеросклероза, представлены особенности распределения напряжений в стенках артерии.
Общие уравнения
Удельная потенциальная энергия для полулинейного материала Джона [3] имеет
вид
э = 2 +^2 = 2 ^/12(и - Е) + ^(и - Е)2, (1)
© Аптуков В.Н., Осоргина Л.Ю., 2008
09806267
где 51 = 51 + 52 + 53, э2 = 5^ 52 + 53, 5; - главные относительные удлинения; X, ц -
параметры Ламе; и - левая мера искажения; Е - единичный тензор; 11 - первый инвариант.
Пусть цилиндрическая поверхность ограничена радиусами г = г0 и г = г1 с заданными давлениями р0, р1. Материальными координатами в отсчетной конфигурации являются цилиндрические координаты г, ф, г с координатными ортами
ег,еф,к .
Предполагая, что радиальное перемещение зависит только от координаты г, а осевое - линейная функция г, получим выражения для радиуса-вектора точки и векторного базиса в актуальной конфигурации
Я = ег/(г) + ^ ^ = ег/'(гX Кз = еФ/(гX Яз = ка (3)
Мера деформации Коши-Грина 7
^0 0Т {2( г)
с = У Я -V Я = ег ег [ / '(г )]2 + ефеф + кка2, (3)
г
а главные относительные удлинения
5, = /\г), 52 = -f (г), 5з = а. (4)
г
Уравнения равновесия в объеме и на поверхности для полулинейного материала приводятся к виду
V u = 0, Xer V• u + 2p,er -V u =
-Poer (^ Г = ^ do
( dO )
-Pier r = rl,
do
0
где u - вектор перемещения; V - набла-оператор Гамильтона в отсчетной
ф dO f(r) й
конфигурации; -----= а-------отношение ориентированных площадок в актуальной и
do r
отсчетной конфигурациях.
А.И. Лурье [3] дает следующее окончательное решение поставленной задачи:
1—[X(2Ci + а) + 2KQ - C
'2\„ r
Or = ------C- [X(2Ci + а) + 2^(Ci - ^ - (3^ + 2^)] ;
(Cl + C2)a
(C - ^2)a r
C-[^(2C! + a) + 2KC1 + ~f) - (3^ + 2ц)] , (6)
C2 4 r
o z = f [^(2Q + а) + 2ца - (3X + 2ц)]
0ф _
где оЯ, оф, ог - физические компоненты тензора Коши; / = С1г +—- - функция,
г
определяющая вид радиального перемещения (2).
Из последнего выражения (6) находим продольную силу Q, которую нужно приложить к торцам цилиндра в актуальной конфигурации, чтобы осуществлялось рассматриваемое напряженное состояние
Я1 г1
Q = 2п | о 2ЯёЯ = 2п | о2/(г) / '(г )ёг =
Я0 г0
= п(г12 - г02 )[Х(2С1 + а) + 2ца - (3Х + 2ц)]. (7)
Таким образом, для трех неизвестных констант решения С1, С2, а существует
уравнение (7) и два уравнения (8), формулируемые на основе силовых граничных
условий на цилиндрических поверхностях
о я (го) = - Po, о я (г1) = - Р1. (8)
Методика определения параметров потенциала
Перейдем теперь непосредственно к рассматриваемой задаче равновесия отрезка артерии в виде замкнутого цилиндра конечной длины под действием внутреннего давления р, тогда Q = рлЯ0. Это требование означает, что продольная сила на торцах
цилиндра обусловлена его замкнутостью и равна просто суммарному давлению на внутреннее поперечное сечение. Краевое условие на торце удовлетворяется в интегральном смысле, в соответствии с принципом Сен-Венана.
Соотношение (7) теперь преобразуется к виду
рЯ0 = (г12 - г02 )[Х(2С1 + а) + 2ца - (3Х + 2ц)]. (9)
Соотношения (8) в развернутой форме
Х(2С1 + а) + 2ц(С1 - % = - р(С1 + %а + 3Х + 2ц, г0 г0
Х(2С1 + а) + 2ц(С1 - % = 3Х + 2ц. (10)
г1
В уравнение (9) входит внутренний радиус цилиндра в актуальной конфигурации Я00, который выражается через неизвестные константы решения
Я = С1г0 + С2 ^ г0 .
Итак, хотя и получено решение в явном виде, но система уравнений для определения констант решения (9)—(10) является нелинейной.
Нашей целью является определение физических констант потенциала Джона X , ц из экспериментальных данных по одновременному измерению давления рк и соответствующих размеров цилиндрического образца артерии в актуальной конфигурации Иэк, Я1, Гк. Получается переопределенная система уравнений (три уравнения на каждом уровне нагрузки для определения двух констант). Отметим, что
не во всяком эксперименте возможно измерение всех трех величин Иэк, Щ, Гк, в данном
случае может использоваться только часть из них.
Очевидно, что задача сводится к минимизации невязки ¥ по двум параметрам
Ламе
N
¥ = £ {[{({) - Щ ({)]2 + [А (Рк) - К (А )]2 + [>„ (Р,) - П (Р, )]2}. (11)
к=1
Причем
ЩХРк ) = С1(Рк )-0 + ) , 1к (Рк ) = а/0.
-0
Ак (Рк ) = Щ1( Рк ) - Щ0( Рк ) = С1 (Рк )(Г1 - го) + (|2( Р1) . (12)
(—)
Г1 Го
Итак, при каждом заданном текущем наборе двух параметров Ламе необходимо
для каждого уровня давления Рк, к = 1, N решить нелинейную систему уравнений (9)—
(10) для того только, чтобы определить три константы решения и вычислить с помощью (12) значение невязки (11).
Чтобы минимизировать невязку по двум параметрам Ламе и определить материальные константы, необходимо провести большой объем вычислений, например сканировать допустимую область параметров Ламе с небольшим шагом. При этом нужно каждый раз получать достаточно хорошее решение нелинейной системы уравнений (9)—(10). Таким образом, даже при наличии решения задачи о больших деформациях цилиндра в конечном виде задача определения материальных констант остается нетривиальной.
Полученные выражения можно несколько упростить, если предположить, что цилиндр бесконечно длинный, тогда а = 1. В этом случае нам не нужно уравнение (7) (или, что то же — (9)). Система (10) упрощается, она становится линейной относительно неизвестных констант решения С1, С2:
Х(2С1 +1) + 2ц(С - % = -р(С1 + % + 3Х + 2ц,
Го Го
К2С1 +1) + 2ц(С - % = 3Х + 2ц. (13)
Г12
Отсюда можно получить явные выражения для С1, С
2
2
2ц(1 --^)- Р
С =-------------------------------
М 2 2 ’
2ц(1 - ■%) - Р[1 + (1 - 2V) %
- -
С2 =---------------------------- , (14)
2ц(1 - %) - р[1 + (1 - 2 V) %
- -
где V - коэффициент Пуассона.
Здесь в качестве искомых независимых материальных констант выступают коэффициент Пуассона и модуль сдвига. Теперь можно в явном виде выразить также текущий внутренний радиус и толщину стенки в актуальной конфигурации
г2
2цт0(1 -________________Г___________
2ц(1 - 4) - Рк [1 + (1 - 2 V) 4]
Г1 Г1
(г - Го)(1 + %2ц(1 - % - Рк ]
К (Рк ) =-------г------------------ • (15)
2ц(1 - % - Рк [1 + (1 - 2
Г Г
Итак, для бесконечного цилиндра получены конечные расчетные выражения (15) толщины и радиуса цилиндра в актуальной конфигурации, входящие в невязку
N
¥ = £{ [ {(рк) - {(рк )]2 + К (Р„) - К (Рк )]2}. (16)
к=1
Ситуация несколько упрощается, поскольку для минимизации невязки не нужно при каждом вычислении текущей геометрии цилиндра решать нелинейную систему трех уравнений. Данный вариант бесконечного цилиндра может использоваться при оценке материальных констант в качестве первого приближения и для предварительной оценки параметров потенциала склерозированных артериальных сосудов.
Анализ результатов
Применим предложенную методику определения параметров потенциала Джона для склерозированных сонных артерий с помощью данных эксперимента [2].
Для нахождения параметров Ламе был использован метод сканирования величин X и ц с достаточно мелким шагом. Для каждой пары параметров были решены уравнения (9)—(10) с граничными условиями (8) для нахождения неизвестных констант С1, С2, а. На основе констант вычислялись радиальные перемещения Яс, которые в
дальнейшем сравнивались с соответствующими экспериментальными данными Яе [2]. Найденные значения параметров Ламе (МПа), соответствующие им радиусы (мм) и погрешности 8 = (|Яс -Яе| / Яе) -100% для различных давлений (в мм рт. ст. и МПа) приведены в таблице.
Степень атеросклероза ц X 80 мм рт. ст. 120 мм рт. ст. 160 мм рт. ст.
0,018072 0,016308 0,021744
Яе Яс 5 Яе Яс 5 Яе Яс 5
1 0,050 0,900 5,560 4,618 16,94 5,800 5,115 11,81 6,000 6,796 13,27
2 0,054 0,965 5,240 4,560 12,98 5,440 4,982 8,42 5,600 5,779 3,20
3 0,056 1,005 4,920 4,535 7,83 5,040 4,927 2,24 5,200 5,602 7,73
4 0,057 1,034 4,720 4,554 3,52 4,920 4,902 0,37 5,000 5,530 10,60
За нулевое значение радиуса принимался радиус артерии при давлении 40 мм рт. ст. Как следует из таблицы, величина погрешности падает с ростом степени атеросклероза, то есть полулинейный материал Джона лучше описывает артерии, находящиеся на последних стадиях атеросклероза. Это связано с тем, что с ростом патологии материал стенок сосудов становится более жестким и меньше проявляет нелинейные свойства.
Установлено, что параметры X и ц изменяются в зависимости от степени патологии по закону
ц = 0,05 • £01, (17)
X = 0,9 • £0,1, (18)
где £ = 1, 4 - стадия атеросклероза.
Ниже для различных стадий атеросклероза на основе соотношений (6) и полученных значений параметров Ламе приведены распределения тангенциальных напряжений вдоль стенок сосудов при давлении 120 мм рт. ст. (рис. 1).
Построенный нелинейный потенциал позволяет оценить уровень сжимаемости стенок артерий. Относительный объем выражается как
V / V = 8,8283 = а 1 / (г), Г '(г) = а(С,2 - С; / г4) . (19)
На рис. 2 показаны кривые распределения относительного объема по толщине сосудистой стенки при различных давлениях для 4-й стадии атеросклероза. Налицо
наличие сильной сжимаемости стенки артерии, возрастающей с ростом
внутрисосудистого давления. Характерным поведением является уменьшение, а не увеличение объема материала стенки.
На рис. 3 показаны аналогичные кривые для различных стадий атеросклероза при давлении 120 мм рт. ст. Заметно уменьшение объемной сжимаемости с ростом патологии.
Факт уменьшения объема материала стенки, отмеченный ранее авторами [2],
а*, МПа
0,20 0,18 0,16 0,14 0,12 0,10 0,08
N N ' > * Ч
* *
4 '
- -
3,0
3,2
3,4
3,6
3,8
4,0
Л, мм
Рис. 1. Распределение тангенциальных напряжений (МПа) по толщине сосудистой стенки (мм) при давлении 120 мм рт. ст. для артерий, находящихся на различных стадиях атеросклероза (1-4 соответствуют 1-4-й стадиям)
свидетельствует о том, что артерия ведет себя подобно толстостенному цилиндрическому сосуду (а не подобно тонкостенной оболочке), когда эластичный внутренний и срединный слои раздавливаются под действием внутреннего давления о
более жесткую внешнюю адвентицию. Полученные результаты говорят о необходимости использования уравнений сжимаемых материалов при моделировании механического поведения кровеносных сосудов человека.
Заключение
Таким образом, на основе аналитического решения А.И. Лурье получены конечные зависимости для вычисления параметров потенциала Джона при деформировании цилиндрического сегмента артерии под действием внутреннего давления. С помощью экспериментальных данных определены параметры потенциала для сонной артерии человека при различных стадиях атеросклероза. Показано, что полулинейный материал Джона лучше описывает кривые деформирования на последних стадиях атеросклероза. Дан анализ величины и характера сжимаемости стенки сонной артерии при различной патологии. Сделан вывод о необходимости учета эффекта сжимаемости при решении задач биомеханики кровеносных сосудов человека.
Список литературы
1. Аптуков, В.Н. Оценка влияния стадий атеросклероза на напряженно-деформированное состояние артериальных сосудов / В.Н. Аптуков, Л. А. Осоргина // Российский журнал биомеханики. - 2005. -Т. 9, № 4. - С. 45-51.
2. Вагнер, ЕА. Механическое поведение сосудистого анастомоза на склерозированных артериях и его моделирование / Е.А. Вагнер, С.Г. Суханов, В.Н. Аптуков // Механика композитных материалов. -1982. - № 2. - С. 336-342.
3. Лурье, А.И. Нелинейная теория упругости / А.И. Лурье. - М.: Наука, 1980.
4. Пуриня, БА. Биомеханика крупных кровеносных сосудов человека / Б.А. Пуриня, В.А. Касьянов. -Рига: Зинатне, 1980.
5. Tickner, E.G. A theory for the static elastic behavior of blood vessels / E.G. Tickner, A.H. Sacks // Biorheology. - 1967. - Vol. 4, No. 4. - P. 151-168.
DETERMINATION OF POTENTIAL PARAMETERS FOR NONLINEAR COMPRESSED MATERIAL OF THE HUMAN CAROTID AT VARIOUS STAGES OF THE ATHEROSCLEROSIS
V.N. Aptukov, L.Yu. Osorgina (Perm, Russia)
In the previous paper [1], the authors estimated the Young’s modules of sclerotic arteries at different stages of the atherosclerosis in the framework of linear model. The stress-strain state of a vessel’s wall around the atherosclerosis plaque was analyzed. As a result, there was made a conclusion on the necessity to solve the given problem in nonlinear statement taking into account finite deformations. On the basis of analytical solution for the cylinder made from the semi-linear John’s material proposed by A.I. Lurie, the method to determine the parameters of nonlinear compressible material potential is presented. The exact values of the parameters for the human carotid at different atherosclerosis stages are obtained.
Key words: biomechanics, sclerotic artery, nonlinear potential, compressibility.
Получено 15 мая 2008
