УДК 531/534: [57+61]
ОСОБЕННОСТИ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ АРТЕРИИ ПРИ НАЛИЧИИ ПАТОЛОГИЧЕСКОЙ ИЗВИТОСТИ (ПЕТЛИ)
С.Г. Суханов1, В.Н. Аптуков2, Л.Ю. Осоргина2
1 Пермский институт сердца, филиал Научного центра сердечно-сосудистой хирургии им. А.Н. Бакулева
Российской академии медицинских наук, Россия, 614000, Пермь, ул. Луначарского, 95б,
e-mail: [email protected]
2 Кафедра математического анализа Пермского государственного университета, Россия, 614990, Пермь, ул. Букирева, 15, e-mail: [email protected]
Аннотация. На основе экспериментальных данных по нагружению образцов артерий внутренним давлением определены параметры нелинейного потенциала Лггиба-Воусв сосудистой стенки артерии. В рамках пространственной модели патологической извитости (петли) артерии с помощью вычислительного комплекса ЛЫБУБ 11.0 выявлены особенности напряженного состояния в зависимости от относительного радиуса петли.
Ключевые слова: артерия, патологическая извитость, нелинейный потенциал, напряженно-деформированное состояние, численное моделирование.
Введение
Нарушение кровообращения ишемического характера головного мозга чаще обусловлено атеросклерозом, аортоартериитом и патологическими извитостями артерий. Например, треть летальных случаев при инсульте сопровождается аномальным развитием сонных и позвоночных артерий [3]. Патологическая извитость может быть устранена хирургическим путем, однако остается открытым вопрос о целесообразности операции для тех или иных значений параметров геометрии. Существует три вида извитости: изгиб (С-образная и ^-образная элонгация), перегиб и петля, последние патологии считаются наиболее опасными. В работе [6] выдвинута гипотеза развития С-образной извитости, постепенно переходящей к ^-образной извитости, а далее к простому и ангулярному перегибу.
Большинство исследований касается изучения особенностей потока крови (распределение давлений и скоростей) при наличии патологической извитости. Некоторые из этих моделей также учитывают деформацию сосудистой стенки в рамках линейной теории упругости, например [5]. Повышенный интерес к гемодинамике объясняется установлением зон образования вихрей и обратных течений, а также развиваемой в последнее время теорией повреждения сосудистой стенки [4], которая говорит о повышенной вероятности развития атеросклеротического процесса на участках интимы, где наблюдаются пониженное значение касательных напряжений в потоке крови, десквамация (шелушение) эндотелия, развитие фиброза и аневризм сосудистой стенки.
© Суханов С.Г., Аптуков В.Н., Осоргина Л.Ю., 2010
Суханов Сергей Германович, д.мед.н., директор Пермского филиала Научного центра сердечнососудистой хирургии им. А.Н. Бакулева РАМН, Пермь
Аптуков Валерий Нагимович, д.т.н., профессор, заведующий кафедрой математического анализа, Пермь Осоргина Людмила Юрьевна, аспирант кафедры математического анализа, Пермь
09806267
Однако опасные участки артерии с точки зрения развития склеротического процесса могут быть связаны и с особенностями напряженно-деформированного состояния сосудистой стенки. Для детального изучения механического поведения кровеносных сосудов в местах патологической извитости необходимо применять уже более сложные нелинейные физические модели.
В данной статье на основе проведенного эксперимента по нагружению артерий внутренним давлением определены параметры нелинейного потенциала Лггиёа-Боусе материала сосудистой стенки. В вычислительной среде ЛЫ8У8 11.0 исследована пространственная модель двухслойного сосуда с патологической извитостью в виде петли. Выявлены особенности напряженно-деформированного состояния сосуда в области петли в зависимости от отношения диаметра петли к диаметру сосуда. Получены значения параметров патологической извитости (петли), при которых механическое состояние сосуда значительно отличается от состояния сосуда без патологии.
Экспериментальные данные
Замкнутые цилиндрические образцы артерий, заполненные
рентгеноконтрастным веществом, нагружались внутренним давлением 40, 80, 120 и 160 мм рт.ст., при каждом из которых производилась рентгенография сосуда в двух проекциях. На полученных рентгеновских снимках измерялись изменения длины, диаметра и толщины стенки артерии, подробно методика эксперимента изложена в статье [2].
Линейные размеры образцов, нагруженных давлением 40 мм рт.ст., были приняты за начальное (недеформированное) состояние артерии. Изменение геометрии артерии при внутреннем давлении 80, 120 и 160 мм рт.ст. оценивалось по отношению к этому начальному состоянию. Всего было исследовано девять образцов несклерозированных артерий, результаты измерений двух из девяти образцов были отброшены как вызывающие сомнения в достоверности. Данные статистической обработки по оставшимся семи образцам представлены в табл. 1, где используются следующие обозначения: 10, I - длина образца (начальная и текущая); Н0, к - толщина
стенки образца; Я0, Я - наружный радиус образца.
Для начального состояния (40 мм рт.ст.) средняя длина образца составляла 51,3 мм, среднеквадратичное отклонение 8,0 мм, коэффициент вариации 15,6%; средняя толщина стенки 1,0 мм, среднеквадратичное отклонение 0,1 мм; коэффициент вариации 10,0%; средний наружный радиус образца составлял 3,5 мм, среднеквадратичное отклонение 0,4 мм, коэффициент вариации 11,4%.
Полученные величины использовались при определении параметров нескольких нелинейных потенциалов в среде ЛЖ^ГО. Отметим, что удовлетворительное описание экспериментальных данных в рамках модели однородной артерии невозможно -деформация наружного радиуса (6-13%) значительно меньше деформации толщины стенки (15-33%). Это подтверждается фактом наличия сжимаемости стенки артерии, обсуждаемым ранее [1].
Таблица 1
Главные относительные удлинения образцов артерии___________________________
Давление, мм рт.ст. 5, = 1 1 *0 5 к к =Т~ к0 5Я = Я Я Я
80 1,17 0,85 1,06
120 1,26 0,75 1,09
160 1,34 0,67 1,13
Выбор и определение параметров нелинейного потенциала
В программном комплексе ЛЫ8У8 была построена двухслойная осесимметричная модель артерии с нелинейно-упругими свойствами. Соотношение толщины внешнего и внутреннего слоев взято на основе рентгенограмм испытанных артерий и составляет 1:10. Механические свойства внешнего слоя описывались линейно-упругими соотношениями с модулем Юнга Е = 9 МПа и коэффициентом Пуассона V = 0,25, а внутреннего слоя - различными нелинейными гиперупругими потенциалами: Лггиёа-Воусе, Blatz-Ko и Моопеу-ЯМт [7].
Для каждого потенциала численным перебором определялись его параметры, при которых достигалась минимальная, по возможности, погрешность между расчетными и экспериментальными данными (табл. 2) для трех значений внутреннего давления 80, 120 и 160 мм рт.ст.
Ниже представлены выражения нелинейных потенциалов и найденные параметры.
Потенциал Аггийа-Боуее:
ж = т
2(1 3)+ 2012
1 (Ч- -9)-
11 (А” - 27)-
10501
+ 19 „ (/14 - 81) +----519 8 (/15 - 24’)
V 1 / 8 \ 1 '
1
+ —
а
V
2
700016 ^ 1 ’ 67375018
Найденные коэффициенты: т = 0,004 МПа, 1 = 1, ё = 50 МПа1. Потенциал Б1а1г-Ко:
СI - >
1 2
ж=т 2
Т + 24Т3 - 5
V 1з
У
Найденный коэффициент: т = 0,008 МПа.
(1)
(2)
Таблица 2
Потенциал Погрешность, % Внутреннее давление, мм зт.ст.
80 120 160
Лггиёа-Воусе |8 * -з* / 5* 2,5 4,5 6,1
|8а -$ / $ 1,2 3,0 7,0
Blatz-Ko з* -з* / з* 2,8 4,7 6,4
/ $ 9,6 4,8 0,4
Моопеу-Я^Ип з * -з* / з* 2,8 6,0 14,3
-8э* / $ 8,8 15,7 15,3
Л
Потенциал Моопеу-Шл>Чп:
ж = Сш (11 - 3)+С„! (12 - 3) +1 (у -1)2
т = 2 ( С01 + С10 ) , а
1 - 2у
(3)
(С01+С10)
Найденные коэффициенты: С10 = 0,04 МПа, С01 = 0,01 МПа, У = 0,35.
В выражениях (1)-(3) используются следующие обозначения: 11 = 51 + 52 + 53 -первый инвариант тензора деформаций; 12 = 6162 + 6263 + 6361 - второй инвариант тензора деформаций; 13 = 616263 - третий инвариант тензора деформаций; ё -коэффициент сжимаемости; т - начальный модуль сдвига; У - детерминант градиента упругой деформации; 1 - коэффициент сходимости; С10, С01 - параметры потенциала.
Анализ таблицы показывает, что для моделирования механических свойств артерий человека при внутреннем давлении в диапазоне от 80 до 160 мм рт.ст. наиболее применим потенциал Лггиёа-Воусе (погрешность не более 7%). Потенциал Blatz-Ko неудовлетворительно описывает поведение материала стенки артерий при низких давлениях, тогда как потенциал Моопеу-КШт дает наиболее высокую погрешность при высоких давлениях. Найденные коэффициенты потенциала Лггиёа-Воусе позволяют проводить адекватный численный анализ механического поведения пространственных моделей артерии.
Расчет напряженно-деформированного состояния патологической извитости
(петли) артерии в среде АМ5У5
В пакете ЛЫ8У8 11.0 была построена пространственная модель петли с двухслойной стенкой: линейно-упругим внешним слоем и нелинейным внутренним слоем со свойствами, заданными потенциалом Лггиёа-Воусе (рис. 1). Геометрию модели условно можно разделить на две части: собственно петля со срединным радиусом * и расстоянием между центрами концов И7; удлиненные концы артерии,
параллельные друг другу и проходящие по касательной к центральной части петли.
Параллельные концы артерии были взяты достаточно длинными для того, чтобы граничные условия на их торцах (жесткое закрепление) не оказывали влияния на напряженно-деформированное состояние в самой петле. Нагружение модели проводилось последовательно внутренней равномерно распределенной нагрузкой, возрастающей ступенчато: 80, 120 и 160 мм рт.ст.
Расчеты показали, что под действием внутреннего давления петля деформируется как в своей плоскости (в направлении оси Х), так и перпендикулярно ей (в направлении оси 7), как показано на рис. 2. Кроме того, происходит увеличение зазора И7 между концами петли (начальное значение И7 = 1 см) и увеличение радиуса петли.
В табл. 3 приведены значения относительного изменения расстояния шага петли АИ7 при различных давлениях.
Таблица 3
Внутреннее давление, мм рт. ст. 80 120 160
АИ7 / И7, % 2,3 3,0 3,5
Рис. 2. Поле перемещений ит при внутрисосудистом давлении 80 мм рт.ст.
-0,001644
-0,001225
-0,000806
-0,000388
-0,000314
-0,000450
Рис. 3. Поле перемещений их при давлении 160 мм рт.ст.
На рис. 3 изображены перемещения вдоль оси Х, лежащей в «плоскости» петли и направленной перпендикулярно к прямолинейным цилиндрическим частям артерии, при давлении 160 мм рт.ст. Конечно-элементная сетка изображена в начальном состоянии до приложения нагрузки.
Максимальные перемещения их наблюдаются в области, наиболее удаленной от
основного ствола артерии (концов петли).
Внешний слой артерии испытывает растягивающие напряжения, в то время как внутренний слой сжимается - изменяется толщина стенки артерии. Это происходит за счет различия упругих свойств внешнего и внутреннего слоев: внутренний более мягкий слой «раздавливается» о внешний слой, обладающий большей жесткостью.
Проведение численного эксперимента в описанной выше постановке требует значительных вычислительных ресурсов, поэтому для упрощения расчетной области и уменьшения количества элементов было сделано предположение о возможности замены объемной модели петли на объемную модель тора. Очевидно, что эти модели эквивалентны за исключением области сопряжения собственно петли с прямолинейной частью артерии. Кроме того, использование модели тора позволяет проводить расчеты только для его части («восьмушки»), учитывая наличие трех плоскостей симметрии (рис. 4).
На рис. 4 отмечены контрольные точки, в которых анализировались
перемещения и напряжения при численном эксперименте.
Было построено шесть моделей тора с различными срединными радиусами Я(,
кратными внешнему радиусу артерии Яа, произведен анализ влияния относительного радиуса тора Я( = Я( / Яа на особенности напряженно-деформированного состояния.
Для тора с относительным радиусом Я( = 3 проведены тестовые расчеты на пяти различных сетках, содержащих от 31 583 до 148 760 элементов. Анализ результатов показал, что получаемые во всех пяти вариантах расчета напряженно-деформированные состояния идентичны.
4
Рис. 4. Модель «восьмушки» тороидальной части петли артерии с контрольными
точками
да /
Рис. 5. Зависимость изменения относительной площади Д5 / \ внутреннего сечения артерии от относительного радиуса петли Я для давления 160 мм рт.ст.
Результаты численного эксперимента показывают, что с уменьшением радиуса петли артерии увеличивается ее жесткость относительно воздействия внутреннего давления - это обусловлено чисто геометрическими причинами. В петле площадь внутреннего сечения сосуда под действием кровяного давления изменяется: она на 25-30% меньше, чем у прямолинейной части сосуда. Этот эффект аналогичен росту жесткости сосудистой стенки для прямолинейных участков артерий, пораженных атеросклерозом. Уменьшение внутреннего сечения отрицательно сказывается на объемном кровотоке, увеличивает внутреннее давление.
На рис. 5 показана зависимость изменения относительной площади внутреннего сечения артерии от относительного радиуса петли, где очевидно, что наиболее сильное отличие от прямолинейного сосуда проявляется при значениях Я{ < 6 .
/ о
О0 / Ое
К
Рис. 6. Зависимость изменения относительного окружного напряжения О0 / интимы в точках 1 и 2 артерии от относительного радиуса петли Я{ для давления 160 мм рт.ст.
Увеличение жесткости сосудистой стенки в петле (уменьшение деформационных свойств, что характерно для длительно существующей аномалии) ведет к росту окружных напряжений на участках интимы, расположенных ближе к центру петли. Если внешняя сторона интимы разгружена, то напряжения на внутренней
стороне для относительного радиуса
петли К < 6 в 1,5-2
раза превышают
соответствующие напряжения в прямолинейной части артерии (рис. 6).
Для внешней стенки (адвентиции) превышение окружных напряжений на внутренней части тора по сравнению с прямолинейной артерией незначительно, оно равно 1,16-1,24.
Заключение
В рамках пространственной модели нелинейно-упругого двухслойного артериального сосуда в вычислительной среде ЛЫ8У8 11.0 проведено численное моделирование процесса деформирования патологической извитости (петли) артерии при действии внутреннего давления.
Показано, что при нагружении петли происходит ее «развертывание»: увеличение радиуса петли одновременно с увеличением расстояния между концами петли. Наибольшие деформации испытывает сторона стенки артерии, располагающаяся ближе к центру петли. При увеличении радиуса петли снижается влияние патологической извитости на напряженно-деформированное состояние. Существует некоторое предельное значение относительного радиуса петли К = 5^6, ниже которого отрицательные геометрические эффекты проявляются достаточно сильно. При этих значениях К площадь внутреннего сечения сосуда под действием кровяного давления изменяется, она на 25-30% меньше, чем у прямолинейной части сосуда. Кроме того, напряжения на внутренней стороне интимы в 1,5-2 раза превышают соответствующие напряжения в прямолинейной части артерии. Найденное предельное значение К является одним из критериев для показаний к хирургической коррекции данной патологии.
Список литературы
1. Аптуков В.Н., Осоргина Л.Ю. Определение параметров потенциала нелинейного сжимаемого материала сонной артерии человека при различных стадиях атеросклероза // Российский журнал биомеханики. - 2008. - Т. 12, № 3(41). - С. 23-31.
2. Вагнер Е.А., Суханов С.Г., Аптуков В.Н. Механическое поведение сосудистого анастомоза на
склерозированных артериях и его моделирование // Механика композитных материалов. - 1982. -№ 2. - С. 336-342.
3. Казанчян П.О., Попов В. А., Гапонова Е.Н., Рудакова Т.В. Диагностика и лечение патологической извитости сонных артерий // Ангиология и сосудистая хирургия. - 2001. - № 2. - С. 93-103.
4. Красавин В.А., Дружинин С.О., Суслов И.И., Морозов К.М., Кузьмин В.А. Концепция системных нарушений иммунитета в генезе аномалий ветвей дуги аорты // Региональное кровообращение и микроциркуляция. - 2007. - № 1(21). - С. 81-83.
5. Павлова О.Е., Иванов Д.В., Грамакова А.А., Морозов К.М., Суслов И.И. Гемодинамика и механическое поведение бифуркации сонной артерии с патологической извитостью // Известия Саратовского университета. Новая серия. - 2010. - Т. 10. - Серия: Математика. Механика. Информатика. - Вып. 2. - С. 66-73.
6. Родин Ю.В. Исследование потоков крови при патологической S'-образной извитости сонных артерий // Международный неврологический журнал. - 2006. - № 4(8). - С. 57-61.
7. ANSYS Basic Analysis Procedures Guide. ANSYS Release 5.6. ANSYS Inc. - 1988.
CHARACTERISTICS OF STRESS-STRAIN STATE OF ARTERY WITH PATHOLOGICAL TORTUOSITY (KINK)
S.G. Sukhanov, V.N. Aptukov, L.Yu. Osorgina (Perm, Russia)
On the basis of experimental data on the loading of arteries samples by internal pressure, parameters of the Arruda-Boyce nonlinear potential of the vascular wall were determined. Within the 3D model of pathologic tortuosity (kink) of arteries with using computing environment ANSYS 11.0, the features of stress state depending on the relative radius of the kink were identified.
Key words: artery, pathological tortuosity, nonlinear potential, stress-strain state, numerical simulation.
Получено 20 октября 2010