CC BY
34
К. Ю. Кузнецов, Ю. Т. Нагибин, В. А. Трофимов
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ НАПРЯЖЕННОГО СОСТОЯНИЯ РОГОВОЙ ОБОЛОЧКИ ГЛАЗА ЧЕЛОВЕКА
Рассматривается возможность использования информации о напряженно-деформированном состоянии глаза человека для диагностики офтальмологических заболеваний. Описывается математическая модель напряженно-деформированного состояния глаза. Приведены результаты математического моделирования картины интерференции поляризованного света.
Ключевые слова: офтальмология, диагностика, напряженно-деформированное состояние глаза.
Лечение широко распространенного заболевания глаз — глаукомы — остается актуальной проблемой и в настоящее время. Глаукома опасна тем, что проявляется на ранних стадиях заболевания лишь повышением внутриглазного давления.
Для диагностики глаукомы используются „бесконтактные" оптико-электронные тонометры [1], в основе действия которых лежит пневматическая аппланация (от англ. applanation — уплощение) роговицы глаза, что позволяет установить значение внутриглазного давления. Данный метод можно назвать бесконтактным лишь условно, так как на глаз действует пневматический импульс, а вследствие инерционности воздействия происходит нежелательная избыточная деформация роговицы глаза.
Известны также проводимые с помощью офтальмологического поляриметра [2] исследования картины интерференции поляризованного света, представляющей собой результат двойного прохождения лучей через роговицу глаза при их отражении от радужной оболочки. Наблюдаемая картина интерференции, отображающая напряженно-деформированное состояние роговицы глаза, позволяет выявить различного рода патологии и, следовательно, обеспечить возможность их диагностики [3].
Картина интерференции формируется в результате одновременного интегрального действия различных факторов, определяющих напряженное состояние тканей роговицы. Раздельное действие каждого из этих факторов невозможно определить экспериментально. Эта проблема может быть решена численными исследованиями на основе математической модели процессов взаимодействия поляризованного излучения с тканями роговицы глаза.
Математическое моделирование целесообразно осуществлять, поэтапно исследуя влияние силовых факторов на напряженно-деформированное состояние роговицы и уточняя геометрию оболочек глаза [4].
Распределения по полю роговицы нормальных усилий, вызванных действием глазодвигательных мышц и внутриглазным давлением, зависят от обеих сферических координат 9 и ф (рис. 1, здесь а и б — усилия, действующие вдоль меридианов и вдоль параллелей соответственно):
УДК 535.551
(9, ф) = Т22 +Tk +Тс = Т2 +Tk -(Sat° sin 9 + Sat1 cos 4ф-sin 9
где
16
К. Ю. Кузнецов, Ю. Т. Нагибин, В. А. Трофимов
üíq =
s sin2 е:
ati =
s sin2 е
S=-
- Eh
Я (1+у)
В приведенных формулах: Е — модуль Юнга, V — коэффициент Пуассона, Я — радиус роговицы, И — толщина роговицы; символом „*" обозначены величины, характеризующие основное напряженное состояние, „к" — напряжения при краевом эффекте, „с" — тангенциальные напряжения [5].
б)
72(9, ф), Н/м
а)
Ti(e, ф), Н/м
-13,6 -13,5 -13,4
-13,3
-13,6 -13,5 -13,4
-13,3
4
4
X, мм
Х, мм Рис. 1
При фурье-разложении функции силового воздействия глазодвигательных мышц можно допустить, что силовое воздействие определено зависимостью:
/ (ф) = (1+^4ф) г,
где г — постоянная, принимаемая в расчетах равной примерно 0,0667.
Распределение по полю роговицы уровней относительной интенсивности расчетной картины интерференции, наблюдаемой при настройке кругового полярископа в скрещенном состоянии, содержит полосы в виде концентрических окружностей (рис. 2). При возрастании внутриглазного давления у лимба роговицы возникают новые полосы, которые смещаются к центру и сгущаются по краям поля; максимум интенсивности отображается белым фоном, минимум — черным.
а)
Y, мм 4
-2
-4
б) Y, мм 4
-4
-2
-4
-2
4 X, мм
0 2 4 Х, мм
Рис. 2
При исследовании влияния глазодвигательных мышц наблюдаются различия при их симметричном и несимметричном действиях. В первом случае (рис. 2, а) интерференционная картина роговицы содержит полосы ромбовидного очертания, во втором случае (рис. 2, б) наблюдается симметрия относительно медиальной мышцы.
2
0
Результаты математического моделирования напряженного состояния роговицы могут быть использованы для определения значения внутриглазного давления или получения информации о других частных причинах патологии. Создание библиотеки расчетных картин интерференции позволит офтальмологам более точно установить диагноз заболеваний пациента.
Авторы выражают благодарность канд. техн. наук О. А. Шеломовой за помощь в проведении расчетов.
список литературы
1. Forbes M., Piko G., Grolman B. A noncontact арр^айои tonometer// Arch. Ophthalmol. 1974. Vol. 91, N 2. P. 134—140.
2. Журавлев А. И., Малышев Л. К. Офтальмологический поляриметр // Оптич. журн. 1994. № 12. C. 71—75.
3. ТамароваР. М. Оптические приборы для исследования глаза. М.: Медицина, 1982.
4. Галилеева П. С., Трофимов В. А., Шеломова О. А. Математическое моделирование наведенного двулучепреломления фиброзной оболочки глаза человека // Изв. вузов. Приборостроение. 2003. Т. 46, № 6. C. 9—14.
5. Бегун П. И., Шукейло Ю. А. Биомеханика: Учеб. пособие. СПб: СПбГЭТУ „ЛЭТИ", 1996. Константин Юрьевич Кузнецов
Юрий Тихонович Нагибин
Владимир Анатольевич Трофимов —
Рекомендована кафедрой твердотельной оптоэлектроники
Сведения об авторах студент; Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кафедра твердотельной оптоэлектроники; E-mail: ckyznetsov@gmail.com канд. техн. наук, доцент; Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кафедра твердотельной оптоэлектроники; E-mail: nagibin77@mail.ru канд. техн. наук, доцент; Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кафедра твердотельной оптоэлектроники; E-mail: troftu@mail.ru
Поступила в редакцию 07.12.09 г.
