CC BY
30
УДК 531/534: [57+61]
МОДЕЛЬ ТРАНСВЕРСАЛЬНО-ИЗОТРОПНОГО СФЕРИЧЕСКОГО СЛОЯ ДЛЯ РАСЧЕТА ИЗМЕНЕНИЯ ВНУТРИГЛАЗНОГО ДАВЛЕНИЯ ПРИ ИНТРАСКЛЕРАЛЬНЫХ ИНЪЕКЦИЯХ
С.М. Бауэр*, Л.А. Замураев*, К.Е. Котляр**
* Кафедра теоретической и прикладной механики, Санкт-Петербургский государственный университет, Россия, 198504, Санкт-Петербург, Университетский проспект, 28, e-mail: [email protected] ** Department of Ophthalmology, Munich University of Technology, 22 Ismaninger Strasse, 81675, Munich, Germany, e-mail: [email protected]
Кафедра офтальмологии, Мюнхенский технический университет, Мюнхен, Германия
Аннотация. Рассмотрена задача трехмерной теории упругости о деформации трансверсально-изотропного сферического слоя, находящегося под действием внутреннего давления. Получена зависимость, характеризующая изменение внутреннего давления в оболочке глаза в зависимости от введенного объема несжимаемой жидкости, а также получено изменение толщины слоя склеры под действием давления.
Ключевые слова: биомеханика глаза, внутриглазное давление, инъекции, модели трансверсально-изотропного сферического слоя.
Введение
Одним из современных способов лечения некоторых глазных заболеваний является интрасклеральная инъекция небольшой (до 0,2 мл) дозы лечебного препарата (например, триамцинолона, авастина и т.д.) [1, 2]. За счет кратковременного увеличения внутреннего объема глазного яблока при введении таких инъекций в первый момент происходит резкое увеличение внутриглазного давления. Даже кратковременное увеличение внутриглазного давления выше определенного индивидуального уровня может привести к нарушению кровообращения на сетчатке и в диске зрительного нерва [3]. Важно в каждом конкретном случае оценить возможный уровень изменения внутриглазного давления в результате инъекции и риск для отдельного пациента, а также, возможно, на основе этих данных рекомендовать уменьшенную дозу препарата для отдельных больных (в настоящее время инъекции имеют стандартный объем: 0,05; 0,1; 0,2 мл, определяемый фирмой-производителем).
С точки зрения механики - это задача определения изменения внутреннего давления в оболочке, заполненной несжимаемой жидкостью, при введении дополнительного объема несжимаемой жидкости.
© Бауэр С.М., Замураев Л.А., Котляр К.Е., 2006
Постановка задачи. Математическая модель
Внешняя оболочка глаза - корнеосклеральная оболочка - состоит из роговицы и склеры, причем склера занимает 93% этой оболочки, поэтому в задачах, связанных с изменением объема глазного яблока под действием внутреннего давления биомеханические свойства склеры играют решающую роль. Роговица в таких задачах обычно не включается в модель, а вся внешняя оболочка глаза рассматривается как состоящая целиком из склеры [4].
При постановке данной задачи и проведении предварительных расчетов авторы пришли к выводу, что расчет изменения внутриглазного давления с использованием различных теорий оболочек, полагающих, что толщина склеральной оболочки остается неизменной, приводит к результатам, которые не очень хорошо согласуются с экспериментальными данными [5]. Вероятно, это связано с тем, что склера является анизотропной (в первом приближении - трансверсально-изотропной [6, 7]) оболочкой, причем модуль упругости склеры в направлении ее толщины на два порядка меньше модуля упругости в тангенциальном направлении [6, 7]. В связи с этим рассмотрим задачу трехмерной теории упругости о деформации сферического слоя с внутренним радиусом Яг и внешним радиусом Я2 (Яг < р < Щ), находящегося под действием
внутреннего давления Р. Для изотропной сферы эта задача описана, например, в [8, 9].
Как и в [8, 9], из соображений симметрии положим, что в сферической системе координат (р, ф, 9 ) справедливы равенства:
аее (р) = СФФ(рХ °рр=срр (р)
РР
РР'
и
°р9 °рф °9ф 0
В этом случае уравнения равновесия приводятся к одному уравнению:
РР + _ (2°рр - °фф - °99 ) = 0 .
Ф Р
(1)
Закон упругости для трансверсально-изотропного тела в сферических координатах имеет вид [10]:
1
єфф е ^°фф уаее) е ' °РР ’ Єфр '
1
а
е.
= - — (а
РР е ' 0
+ а
) 2.
фф ) + Е' аРР ’
е=
афЄ
~в
Число упругих коэффициентов в этом случае равно пяти [10]: Е и Е' -соответственно, модули Юнга при растяжении - сжатии в поверхности изотропии и в направлении, перпендикулярном к ней; V и V' коэффициенты Пуассона; О' - модуль сдвига для любой плоскости, перпендикулярной к поверхности изотропии; модуль сдвига для поверхности изотропии О определяется соотношением
е
о = . Е
2(1 + V)
Упругие коэффициенты удовлетворяют системе неравенств [10]:
Ч Ї)"-' - ¥ ■
Предполагается также, что компоненты вектора перемещения удовлетворяют соотношениям
иФ = ие = 0, иР = иР (Р).
В силу этого компоненты тензора деформаций в сферической системе координат примут следующий вид:
йир иР п
е =—-, е = еаа = —, еа = е = е = 0.
рр йр фф ее р фЄ єр фр
Выразим компоненты тензора напряжений через компоненты тензора деформаций и компоненты вектора перемещений:
ир йир
°ФФ=°»=АеФФ+Верр=VВ-йр •
ир йир
п = Се + Бе = С-^ + Б- р
рр фф рр р йр
где
Е Е1
А =---------—--------—, В = IV', С = 2IV', Б = А(1 - V) Е'/ Е.
Е '(1 -V) - 2(v ')2 Е
При подстановке этих компонент в уравнение равновесия (1) получается дифференциальное уравнение:
й2ир 2 йир 2Е(V'-1)
и„
2 +----р +—^------’--Рт = 0 (2)
йр р йр Е '(1 -V) р
с граничными условиями:
арр (Яг) = -Р, арр (Л2) = 0. (3)
Решение уравнения (2) имеет вид:
. = С№+с
где
ир (Р, Р) = -^ + С2(Р)Р”
1 ф + 4р 2 Е(1 -V')
да = — + ----------, В = — ---------------.
2 2 Е '(1 -V)
Отметим, что в изотропном случае т = 1, так как Е = Е', v = v', Р = 2 Из граничных условий (3) можно получить
_ р [(1 -у)Е2у'2 Е] Я^2Ят+1
1 _ Е [(т +1)(1 -у)Е'- 2у ' Е] (Я2^1 - Я12т+1) ’
_ р [(1 -у)Е'- 2у '2 Е] Я1т+2
2 _ Е' [т(1 -у)Е'+ 2у'Е] (Ят+1 -Я12т+1) '
Изменение объема исходной сферы, находящейся под действием давления Р, можно определить по следующему соотношению:
4 , 4
АУ _ 3<Я + пр(Я1,Р))3 -3<Я1)3,
а отсюда можно вывести и обратную зависимость Р _ /(АУ), характеризующую изменение внутреннего давления от дополнительно введенного объема жидкости.
Результаты расчета и выводы
Для сферического слоя с параметрами Я1 _ 11,75 мм, Я2 _ 12,25 мм, Е _ 14,3 МПа и у _ 0,4 - реальными параметрами склеральной оболочки эмметропического глаза (глаза с нормальным зрением), полученными в экспериментальных исследованиях [4, 6], и при начальном значении внутриглазного давления равного 15 мм рт. ст., при различных соотношениях Е/ Е' для инъекции 0,1 мл получены значения коэффициента Пуассона у' и изменения внутриглазного давления, приведенные в таблице.
Таблица
Изменение внутриглазного давления при введении инъекции 0,1 мл
Соотношение Е / Е1 Коэффициент Пуассона у' Р (мм рт. ст.)
1 0,40 70,2
2 0,30 68,8
10 0,15 60,6
100 0,01 56,3
200 0,0075 47,4
Данные таблицы показывают, какое существенное значение имеет учет соотношения трансверсального и нормального модулей упругости склеры. Отметим, что результаты, представленные в предпоследней строке таблицы, наиболее соответствуют экспериментальным данным, полученным при интраоперационном измерении внутриглазного давления тонометром Шиотца сразу после инъекции 0,1 мл препарата тримцинолон [5].
Паликарис (Pa.llika.ri^) с соавторами [11] на основе прижизненных измерений изменения внутриглазного давления от увеличения объема глазного яблока в группе здоровых лиц приводит экспериментальную линейную зависимость роста внутриглазного давления от объема добавляемой в глаз малыми порциями жидкости. Авторы приводят данные для дополнительных объемов до 0,08 мл. Экстраполяция этих
А
Г
♦
0 0,02 0,04 0,06 0,08 0,1
Изменение объема полости сферы V, мл
♦ Е/£=1 □ Е1Е=2 А £/£'=10 • £/£=100 □ £/£=200
Рис. 1. Зависимость изменения толщины склеральной оболочки от объема введенной
несжимаемой жидкости
данных для дополнительного объема 0,1 мл при начальном внутриглазном давлении, равном 15 мм рт. ст., дает значение соответствующего изменения внутриглазного давления в 55 мм рт. ст., что хорошо согласуется с данными расчетов, проведенных в статье, при Е / Е' = 100.
Данные о коэффициенте ригидности корнеосклеральной оболочки глаза, приводимые авторами других прямых манометрических измерений изменения внутриглазного давления от увеличения объема глазного яблока [12-14], хорошо соотносятся с данными работы [11]. Хотя понятие ригидности тканей в офтальмологии и в биомеханике различаются, математическая аппроксимация изменения внутриглазного давления при изменении внутриглазного объема на 0,1 мл по этим экспериментальным данным, видимо, приведет к похожим результатам.
С возрастом ригидность склеральной оболочки возрастает (увеличиваются модули упругости склеры, причем модуль упругости склеры в направлении ее толщины увеличивается быстрее) [6]. Это приводит, видимо, к уменьшению отношения Е / Е', и тот же самый объем инъекции в более пожилом возрасте может вызывать более высокий рост внутриглазного давления, чем в молодом возрасте. Кроме того, поскольку упомянутые клинические эксперименты [5, 11] проводились в основном на пожилых пациентах, можно предположить, что отношение Е / Е' в глазах молодых людей еще больше, чем 100.
Отметим также, что по соотношению для нормального перемещения можно получить изменение толщины слоя склеры под действием давления Р:
На рис. 1 показана рассчитанная по представленной математической модели зависимость относительного изменения толщины склеральной оболочки от объема введенной несжимаемой жидкости.
Изменение толщины склеральной оболочки при увеличении внутриглазного давления как один из результатов представленного моделирования, согласуется с имеющимися данными о том, что при открытоугольной глаукоме, глазном заболевании, сопровождающемся, как правило, повышенным уровнем внутриглазного давления, также наблюдается уменьшение толщины склеры [15].
М = ир (^ р) - ир р).
При этом при глаукоме также увеличивается ригидность корнеосклеральной оболочки [16, 17]. Возможно, наблюдавшееся в эксперименте увеличение ригидности склеры при глаукоме [16] также сопровождается уменьшением отношения Е / Е', отличающегося в динамике от физиологического.
Заключение
Необходимость учета трансверсальной анизотропии корнеосклеральной оболочки глаза для конкретных задач моделирования физиологических и патологических процессов в глазу, а также возрастные зависимости значения Е / Е' в норме и при патологиях глаза предстоит еще изучить. Первые результаты моделирования изменения внутриглазного давления при интрасклеральных инъекциях показали, что представленный математический аппарат может быть использован для адекватной оценки роста внутриглазного давления в результате инъекции для глаз, имеющих форму, близкую к сферической форме, по известным биометрическим данным.
Представленные результаты для сферического слоя являются основой для решения более общей задачи оценки изменения внутреннего давления при изменении внутреннего объема оболочки эллипсоидальной формы. Такая задача позволит оценивать внутриглазное давление после инъекций в гиперметропических (дальнозорких) и миопических (близоруких) глазах.
Благодарности
Работа выполнена при поддержке РФФИ, грант № 04-01-00258.
Список литературы
1. Jonas, J.B. Intravitreal use of triamcinolone acetonide: an experimental clinical procedure / J.B. Jonas // Acta Ophthalmol. Scand. - 2005. - Vol. 83, No. 5. - P. 526-527.
2. Avery, R.L Intravitreal bevacizumab (Avastin) for neovascular age-related macular degeneration / R.L. Avery, D.J. Pieramici, M.D. Rabena, A.A. Castellarin, M.A. Nasir, M.J. Giust // Ophthalmology. -2006. - Vol. 113, No. 3. - P. 363-372.
3. Nagel, E. Retinal vessel reaction to short-term IOP elevation in ocular hypertensive and glaucoma patients / E. Nagel, W. Vilser, I.M. Lanzl // Eur. J. Ophthalmol. - 2001. - Vol. 11, No. 4. - P. 338-344.
4. Бауэр, С.М. Простейшие модели теории оболочек и пластин в офтальмологии / С.М. Бауэр, Б.А. Зимин, П.Е. Товстик. - СПб.: Изд-во СПбГУ, 2000. - 92 с.
5. Maier, M. Effect of Intravitreal Injections and Volume Change on Intraocular Pressure - a Biomechanical model / M. Maier, S.M. Bauer, I.M. Lanzl, N. Feucht, E. Fisher, C. Winkler von Mohrenfels, C. Lohmann, K.E. Kotliar // Ophthalmic Research. - 2005. - Vol. 37, Suppl. 1. - P. 197.
6. Иомдина, Е.Н. Биомеханические свойства склеры и возможности ее укрепления при миопии / Е.Н. Иомдина: автореф. дис. ... канд. биол. наук. - М., 1984. - 24 c.
7. Battaglioli, J.L. Measurements of the compressive properties of scleral tissue / J.L. Battaglioli, R.D. Kamm // Invest. Ophthal. Vis. Sci. - 1984. -Vol. 25. - P. 59-65.
8. Лурье, А.И. Пространственные задачи теории упругости / А.И. Лурье. - М.: Гостехиздат, 1955. - 492 c.
9. Атанацкович, Т. Лекции по теории упругости / Т. Атанацкович, А. Гуран. - СПб., 2003. - 400 c.
10. Родионова, В.А., Пpикладная теоpия анизотpопных пластин и оболочек / В.А. Родионова, Б.Ф. Титаев, К.Ф. Чеpных // СПб.: Изд-во С.-Петерб. ун-та, 1996. - 278 c.
11. Pallikaris, I.G. Ocular rigidity in living human eyes / I.G. Pallikaris, G.D. Kymionis, H.S. Ginis, G.A. Kounis, M.K. Tsilimbaris // Invest. Ophthalmol. Vis. Sci. - 2005. - Vol. 46, No. 2. - P. 409-414.
12. Ytteborg, J. The effect of intraocular pressure on rigidity coefficient in the human eye / J. Yttenborg // Acta Ophthalmol. - 1960. - Vol. 8. - P. 548-561.
13. Eisenlohr, J.E. Manometric studies of the pressure-volume relationship in living and enucleated eyes of individual human subjects / J.E. Eisenlohr, M.E. Langham, A.E. Maumenee // Br. J. Ophthalmol. - 1962. -Vol. 46. - P. 536-548.
14. Silver, D.M. Pressure-volume relation for the living human eye / D.M. Silver, O. Geyer // Curr. Eye Res. -2000. - Vol. 20. - P. 115-120.
15. Нестеров, А.П. Глаукома / А.П. Нестеров. - М.: Медицина, 1995. - 256 с.
16. Козлов, В.И. Новый метод изучения растяжимости и эластичности оболочек глаза при иэменении офтальмотонуса / В.И. Козлов // Вестник офтальмологии. - 1967. - № 2. - С. 5-9.
17. Светлова, О.В. Развитие теории оттока водянистой влаги и перспективные гипотензивные воздействия / О.В. Светлова, М.В. Засеева, А.В. Суржиков, И.Н. Кошиц // Глаукома. - 2003. - № 1. -С. 51-59.
MODEL OF THE TRANSVERSELY ISOTROPIC SPHERICAL LAYER FOR ESTIMATION OF INTRAOCULAR PRESSURE CHANGES AFTER INTRAVITREAL INJECTIONS
S.M. Bauer, L.A. Zamuraev (St.-Petersburg, Russia),
K.E. Kotliar (Munich, Germany)
A problem of 3D-theory of elasticity on the deformation of transversely isotropic spherical layer under the internal pressure is considered. This problem is developed for the prediction of the intraocular pressure increase after intravitreal injections in the eye. Dependence of internal pressure change in the closed shell cavity on the introduced additional volume of incompressible liquid is obtained. Change of the layer thickness is calculated as a function of the internal pressure.
Key words: eye biomechanics, intraocular pressure, intravitreal injection, transversal-isotropic spherical layer.
Получено 4 апреля 2006
