Полный текст

Введение Роговица глаза обладает весьма сложными, недостаточно изученными и сильно индивидуальными механическими свойствами [3, 9, 13]. Если даже ограничиться только упругим поведением, роговица анизотропна и пространственно неоднородна, а ее упругие свойства нелинейны [10, 11, 14, 15]. Кроме того, роговица неоднородна геометрически - ее толщина и кривизна существенно меняются от центральной точки (апекса) к периферии. Тем не менее при моделировании многих механических процессов, в том числе статической тонометрии, эффективными оказываются относительно простые модели, рассматривающие роговицу как сферическую, изотропную и в тангенциальном направлении однородную линейно упругую оболочку [2, 4, 6]. Учет влияния пространственной неоднородности становится существенным, когда рассматривается роговица, подвергшаяся значительным механическим изменениям в результате корректирующих зрение механических операций. Такого рода операции в последние десятилетия стали массовыми. В частности, широко применяются операции, при которых удаляется материал роговицы из некоторой области, прилегающей к ее внешней поверхности [1]. Локализация этой области может быть различной - в центре или на периферии - в зависимости от исправляемого дефекта зрения. Поскольку толщина удаляемого слоя составляет заметную долю всей толщины роговицы, а именно ближние к поверхности слои в наибольшей степени отвечают за ее прочностные свойства [10], происходит значительное механическое ослабление роговицы в оперированной зоне - ее эффективная жесткость на растяжение сильно падает. Такое ослабление неизбежно влияет на связь характеристик, определяемых при непрямых офтальмологических измерениях, с физически существенными параметрами механического состояния глаза: методы оценки, разработанные для интактного глаза, могут оказаться недостоверными. К числу диагностических процедур, базирующихся на оценке деформации роговицы под механическим воздействием, относятся собственно тонометрия, направленная на оценку внутриглазного давления, и дифференциальная тонометрия, дающая возможность комплексной оценки механических свойств и механического состояния глаза. Если для нагружения используется тонометр Маклакова (широкий плоский штамп), эту процедуру принято называть эластотонометрией. Диагностические возможности дифференциальной тонометрии (и в частности эластотонометрии) велики; поэтому возможность ее использования для оперированного глаза имеет принципиальное значение. Разработанная авторами максимально простая механическая модель глазного яблока, рассматривающая роговицу как двумерную упругую поверхность [6], а склеру как нульмерный упругий элемент, была эффективно применена к анализу влияния пространственной неоднородности роговицы на механические свойства глаза в работе [7], где, в частности, рассмотрена тонометрия по Маклакову для случая неоднородной роговицы, в том числе при ослабленной апикальной области. Оказалось, что учет неоднородности заметно влияет на существенные для диагностики результаты расчетов лишь при очень высоком ее уровне, не достигаемом при естественной неоднородности. Однако неоднородность, возникающая в результате операций, может оказывать на них значительное влияние. В настоящей работе анализируется влияние присутствия искусственно ослабленных областей на механическое поведение глаза при нагружении широким плоским штампом (тонометром Маклакова). В центре внимания процедура эластотонометрии и влияние неоднородностей такого рода на определяемую при этой процедуре характеристику - коэффициент эластоподъема. Сравниваются два качественно различных случая локализации ослабленной области - апикальной и периферической - соответствующих операциям, направленным на коррекцию миопии и гиперметропии соответственно. Применяемый метод позволяет наглядно оценивать следствия различных гипотез, что отличает его от громоздких вычислительных методик, для которых такая наглядность отсутствует. Постановка задачи В соответствии с развиваемым авторами подходом будем рассматривать глазное яблоко как заполненный несжимаемой жидкостью с давлением объем, состоящий из двух частей - роговицы (на рис. 1, а выше плоскости AB) и склеральной области, включающей склеру и прилегающие к ней ткани (ниже плоскости AB). Роговица рассматривается как материальная поверхность, т.е. слой нулевой толщины, а склеральная область - как упругий элемент, деформация которого характеризуется только объемом находящейся в нем жидкости. Поверхность предполагается безмоментной, т.е. не сопротивляющейся изгибу, изотропной и линейно упругой, однако неоднородной пространственно, т.е. ее упругие константы зависят от положения на поверхности. Распределение значений этих констант, тем не менее, симметрично относительно оси, что позволяет рассматривать задачу как осесимметричную. В [7] показано, что в предположении малых деформаций (но не обязательно малых перемещений) для неоднородной в тангенциальном направлении поверхности сохраняется форма уравнений, совпадающая с уравнениями для пространственно однородной роговицы. Имеется единственное отличие: упругие модули являются теперь функциями расстояния соответствующей точки от апекса по большой дуге. Система обыкновенных дифференциальных уравнений имеет, таким образом, вне контакта со штампом (тонометром) вид (1) где - расстояние точки поверхности, представляющей роговицу, от оси; - угол между нормалью к этой поверхности и осью симметрии; - длина дуги образующей роговицы, отмеряемая от точки пересечения роговицы с осью симметрии (апекса). p G A B а б Рис. 1. К постановке задачи: а - глазное яблоко до нагружения; б - роговица, нагруженная тонометром Маклакова. Обозначения в тексте Величина (длина дуги в начальном состоянии) рассматривается как сопутствующая координата и играет роль независимой переменной. Другие обозначения: и - касательные усилия в меридиональном и перпендикулярном к нему (окружном) направлениях; и - деформации в направлениях главных осей, отсчитываемые от начального состояния ненагруженной извне роговицы. Последнее условно считается соответствующем некоторому сферическому пространственно однородному сегменту радиуса (радиус роговицы) с внутренним давлением В силу малости деформаций это, вообще говоря, не реализуемое состояние мало отличается от любого реализуемого (и, вообще говоря, неоднородного) состояния ненагруженной извне роговицы при том же давлении. В этом состоянии При симметричном относительно оси наложении широкого штампа (тонометра) в предположении свободного проскальзывания система уравнений для роговицы имеет тот же вид (1) с исключением второго уравнения этой системы, заменяемого геометрической связью Приращение «подроговичного» объема (над плоскостью AB) связано с приращением всего заполненного жидкостью внутриглазного объема соотношениями (2) где - приращение «склерального» объема (под AB). Переобозначим упругие константы, присутствующие в (1), по формулам где мм - средняя в норме толщина роговицы, а мм-3. Физический смысл такого переобозначения подробно обсуждается в [6]. Будем называть и роговичной и склеральной жесткостями соответственно, а - коэффициентом Пуассона представляющей роговицу поверхности. В силу симметрии решение задачи ищется в области где - значение координаты на границе роговицы и склеры (в лимбе). Граничные условия подробно рассмотрены в [6] и включают условие (жесткое закрепление роговицы на склере), условие ограниченности усилия в апексе, а также выполнение соотношения (3) при неизвестном заранее значении , соответствующем отрыву роговицы от штампа, где - вес груза, и непрерывность функций и Оценка пределов приложимости развиваемого подхода к задаче о плоском штампе для безмоментной оболочки с заметно меняющимся в сравнении с ее толщиной профилем внешней поверхности (а значит, и с сильно меняющейся толщиной) нуждается в специальном исследовании. Необходимость такого исследования, которое в настоящее время авторами проводится, возникает в связи с тем, что на практике значение в условии (3) определяется не для оболочки в среднем, а именно для внешней поверхности. Уже полученные результаты показывают, что в практически важных случаях использование рассмотренной выше модели допустимо. В настоящей работе решается задача о тонометрии по Маклакову, т.е. при симметричном наложении широкого плоского штампа на роговицу и постоянном суммарном объеме внутриглазной жидкости При расчетах будем полагать мм, мм, что соответствует стандартным для человеческого глаза значениям. Далее независимую переменную (длина дуги образующей) будем обезразмеривать по радиусу роговицы сохраняя для этой величины прежнее обозначение. При принятых значениях геометрических параметров роговицы Для коэффициента Пуассона принималось значение Склеральная жесткость во всех расчетах считалась равной ее среднему для нормальных глаз человека значению МПа [6]. Жесткость роговицы варьировались в диапазоне 0,1-0,8 МПа, и рассматривались два типа ее изменения по дуге большого круга. Первый вариант моделирует резкое ослабление прочности роговицы в области, прилегающей к апексу (рис. 2, а). Это распределение соответствует удалению в некоторой центральной зоне части материала роговицы, прилегающей к ее внешней поверхности. Такое удаление имеет место при операциях, корректирующих миопию [1]. Второй вариант (см. рис. 2, б) моделирует резкое ослабление роговицы в ограниченной области, лежащей между апексом и лимбом. Этот вариант соответствует удалению прилегающего к внешней поверхности материала из такой области, что выполняется при коррекции гиперметропии. Предполагается, что до операции упругие свойства роговицы были пространственно однородны: на рис. 2 жесткость интактной роговицы соответствует участкам, примыкающим к правому концу (т.е. контакту со склерой), и равна Вес штампа (тонометра) варьировался от 5 до 15 г. Границы этого интервала соответствуют минимальному и максимальному весу в стандартном наборе грузов, применяемых при тонометрии по Маклакову. В числе физически существенных характеристик распределения жесткость интактной роговицы и среднее значение жесткости. Последнее можно вводить различно. Если не оговорено, под средним понимается среднее интегральное значение по длине дуги Именно это значение использовалось в работе [7]. а б Рис. 2. Общий вид зависимостей эффективной жесткости роговицы от характеризующего положение на роговице безразмерного параметра при радикальном ослаблении апикальной зоны (а) и кольцевой области на периферии (б) в первоначально однородной роговице Когда далее говорится о среднем значении жесткости без уточнений, имеется в виду Возможно рассматривать и, например, среднее интегральное по площади поверхности роговицы Как будет видно из дальнейшего, выбор определения среднего достаточно заметно влияет на некоторые результаты. Математически упомянутые средние определены формулами (4) где - площадь поверхности сферического сегмента, представляющего роговицу в начальном состоянии, а координата определяется как функция для этого состояния в соответствии с (1). В качестве основной характеристики, отражающей степень неоднородности роговицы, выберем амплитуду жесткости т.е. разницу между ее максимальным при (на границе со склерой) и минимальным значениями. В случае ослабленной центральной зоны (см. рис. 2, а) минимальное значение жесткости достигается при (в апексе) и При периферическом ослаблении (см. рис. 2, б) минимальное значение жесткости имеет место в некоторой внутренней точке отрезка Важна еще локализация перехода между характерными жесткостями, в первую очередь между уменьшенной и интактной. Для апикального ослабления роговицы (см. рис. 2, а) эта локализация характеризуется одним пограничным значением ( - измеренный по дуге диаметр ослабленной зоны). В случае ослабления нецентральной зоны (рис. 2, б) локализация ослабления определяется двумя параметрами: и где и - внутренний и внешний диаметры ослабленной зоны соответственно. Типичные радиусы зоны ослабления жесткости равны при выполнении операций по коррекции миопии (4,3 мм), при коррекции гиперметропии (3,2 мм) и (4,6 мм). Расчеты показывают, что наклон переходных между зонами кривых несуществен. При моделировании ослабления в центральной зоне любые три из названных четырех характеристик ( ) независимы. В случае ослабления кольцевой периферической зоны их четыре из набора Методику численного решения см. в [6]. Сравнительное моделирование эластотонометрии по Маклакову для роговицы, ослабленной в центральной и периферической зонах Процедура эластотонометрии состоит в последовательном нагружении роговицы глаза широкими плоскими штампами (тонометрами Маклакова) разного веса и измерении формирующихся при этом давлений (называемых тонометрическими, в отличие от давления до нагружения - «истинного»). Для роговицы как безмоментной поверхности тонометрическое давление определяется в соответствии с соотношением (3) (формула Имберта-Фика). Рассчитанные кривые давление - вес штампа оказались, как и в случае однородной роговицы, близки к прямым, по крайней мере в практически существенной области значений. Заметное отклонение от прямолинейности, обнаруженное при расчетах на трехмерной конечно-элементной модели в [8], связано с достаточно произвольным выбором этой модели, что делает сопоставление затруднительным. Практическая прямолинейность зависимости позволяет ввести [6] в качестве характеристики конкретного глаза коэффициент наклона этой прямой названный коэффициентом эластоподъема. Коэффициент эластоподъема зависит от упругих характеристик роговицы и склеры, а также от истинного давления На рис. 3 представлены зависимости коэффициента эластоподъема в случае ослабленной апикальной зоны от средней жесткости роговицы (рис. 3, а) и от радиуса этой зоны (рис. 3, б) при двух различных амплитудах неоднородности. Кривые приближены ортогональными полиномами, чтобы исключить мелкомасштабные по аргументу вариации функций. Разные диапазоны изменения аргумента на рис. 3, а связаны с чисто арифметическими ограничениями на параметры распределения жесткости. Видно, что зависимость от средней жесткости ведет себя иначе, чем для однородной роговицы с той же жесткостью. В физиологически существенной области значений эта зависимость оказывается существенно немонотонной, возрастая при малых значениях и убывая при достаточно больших, тогда как при однородной жесткости расчет дает при малых почти не различимое возрастание коэффициента эластоподъема, за которым следует убывание, а при достаточно больших истинных давлениях возрастающий участок вообще исчезает [6]. Аналогичное немонотонное поведение имеет место и для зависимости Максимум кривых тем более выражен, чем больше амплитуда жесткости а при фиксированной становится заметней с убыванием интактной жесткости Максимум кривой там, где он выражен, достигается при близких значениях радиуса ослабленной зоны D / 2 ~ 0,45-0,5 (для разных значений параметров). Для сравнения, минимальный а б Рис. 3. Зависимость коэффициента эластоподъема g для неоднородной, апикально ослабленной роговицы от ее средней жесткости (а) и радиуса зоны ослабления (б) при (пунктирные кривые, интактная жесткость ) и (кривые 1 - 2 - 3 - ). Зависимость g от жесткости для однородной роговицы показана штрих-пунктиром (а). Все жесткости даны в МПа. Истинное давление мм рт. ст. безразмерный радиус области контакта груза с роговицей, определяемый по формуле (3), в рассматриваемом диапазоне варьирования параметров равен 0,28 (при весе тонометра г), а максимальный - 0,43 ( г). Таким образом, убывание кривых соответствует ослабленному участку, перекрывающему область контакта роговицы со штампом. Отклонение от значений γ, соответствующих однородной жесткости, возрастает при увеличении амплитуды, а при одной и той же амплитуде растет с ростом достигая максимума при некотором значении близком к тому, в котором достигается максимум функции а затем начинает падать. На рис. 4 кривые сопоставлены для различных значений истинного давления и двух значений склеральной жесткости . Коэффициент эластоподъема растет с уменьшением при приблизительном сохранении формы кривых и . Он также растет с увеличением , причем при возрастании склеральной жесткости точка максимума функции сдвигается вправо (в сторону бо́льших значений ), так что при МПа и мм рт. ст. функция становится в исследуемой области значений практически монотонной (теперь возрастающей). Монотонной (убывающей) при этих значениях становится и функция . Разница между значениями g для неоднородной при апикальном ослаблении и однородной роговиц не слишком велика: она не превышает 0,11 (11%) абсолютной величины g при истинном давлении 12 мм рт. ст. и 0,15 (13%) при истинном давлении а б Рис. 4. Зависимость коэффициента эластоподъема g от средней жесткости апикально ослабленной роговицы (а) и радиуса зоны ослабления (б) при МПа, склеральной жесткости МПа и истинных давлениях (кривая 1), 16 (2), 20 (3) мм рт. ст., а также при склеральной жесткости МПа и истинном давлении мм рт. ст. (кривая 4). Зависимости для однородной роговицы от ее жесткости при тех же значениях и обозначены индексом а 20 мм рт. ст. в случае ΔEc = 0,78, а в случае ΔEc = 0,34 при истинном давлении 12 мм рт.ст. - 0,06 (8%). Эта разница сопоставима с вариациями коэффициента эластоподъема вследствие изменений склеральной жесткости и истинного давления, а порой и значительно меньше их. Таким образом, качественно отличающаяся от случая однородной роговицы немонотонная зависимость в этом случае не слишком существенна. При умеренной амплитуде вполне можно использовать зависимость полученную для однородной роговицы, если понимать под среднее значение жесткости . Более того, при относительно малых значениях истинного давления (~12 мм рт. ст.) и близких к среднему значениях склеральной жесткости (~10 МПа) сама зависимость γ от , вообще, несущественна и коэффициент эластоподъема определяется практически только и . На рис. 5 представлены зависимости коэффициента эластоподъема от и полуширины оставшегося интактным апикального участка при различных других параметрах распределения жесткости роговицы для случая периферического выреза. При не слишком узкой зоне ослабления (что при равенстве остальных параметров соответствует умеренным значениям ) функция ведет себя в известном смысле противоположно случаю апикального ослабления: она теперь выпукла не кверху, а книзу: сначала падает с ростом , а затем начинает расти. Точка минимума этой функции сдвигается в сторону бо́льших при увеличении а б Рис. 5. Зависимость коэффициента эластоподъема роговицы, ослабленной в зоне близкой к склере, от ее средней жесткости (а) и внутренней (ближней к апексу) границы зоны ослабления (б) при различных значениях интактной жесткости Ec0 = 0,6 (кривые 1) и 0,7 (кривые 2 и 3) МПа и разной локализации внешней (удаленной от апекса) границы ослабленной зоны D2/2 = 0,6 (кривые 1 и 2) и 0,7 (кривые 3); пунктирная кривая 4 - зависимость для однородной роговицы (в этом случае ). Амплитуда неоднородности, истинное давление и склеральная жесткость во всех случаях фиксированы и равны 0,46 МПа, мм рт. ст. и МПа интактной жесткости и отдалении внешней границы зоны ослабления от границы со склерой (лимба). Однако при дальнейшем росте достигается максимум и функция возобновляет падение. Функция выглядит аналогично. Сначала падает с ростом (т.е. с отдалением ослабленной зоны от апекса и уменьшением ее ширины), а затем начинает расти. При этом отклонение коэффициента эластоподъема от его соответствующего значения для однородной роговицы с тем же заметно меньше, чем для апикального выреза, и не превышает 2%. Из рис. 6 видно, что эти отклонения значительно меньше вариаций функции из-за вполне реалистичных вариаций склеральной жесткости и истинного давления . Таким образом, неоднородностью жесткости роговицы при анализе зависимости коэффициента от этой величины можно свободно пренебрегать даже при больших ее амплитудах и пользоваться расчетной кривой для однородной роговицы с жесткостью, равной среднему значению ее распределения. Рис. 7 демонстрирует практически полное совпадение эластотонометрических кривых (зависимостей тонометрического давления от веса тонометра) в случае однородной роговицы и ослабленной в кольцевой зоне роговицы с тем же средним значением жесткости. Отличие от аналогичной зависимости для интактной роговицы Рис. 6. Зависимость коэффициента эластоподъема роговицы, ослабленной в зоне близкой к склере, от ее средней жесткости при различных значения параметров (сплошные кривые). Интактная жесткость Ec0 = 0,7 (кривые 1, 3) и 0,6 МПа (2); внешний радиус ослабленной зоны D2/2 = 0,6 (1, 2) и 0,7 (3). Кривые группируются в три семейства - а, б и в, различающиеся значениями склеральной жесткости и истинного давления а - МПа, мм рт. ст.; б - МПа, мм рт. ст.; в - МПа, мм рт. ст. Пунктирные линии - зависимость для однородной роговицы () при значениях и соответствующего семейства. Амплитуда во всех случаях равна 0,46 МПа Рис. 7. Зависимость тонометрического давления от веса груза: 1 - роговица с апикально ослабленной зоной; 2 - роговица с кольцевой зоной ослабления; 3 - однородная роговица с жесткостью, равной средней; 4 - однородная роговица с жесткостью, равной интактной жесткости (до операции). Во всех случаях МПа, мм рт. ст. Для кривых 1 и 2 средняя жесткость одинакова и равна 0,48 МПа также незначительно. Таким образом, при расчете истинного давления использование соотношений, полученных для однородной роговицы с жесткостью равной интактной, дает лишь незначительную ошибку, а с жесткостью равной средней - ничтожную. В то же время эластотонометрическая кривая для апикально ослабленной роговицы располагается выше, хотя и не слишком существенно: ее подъем над линией, соответствующей среднему значению порядка 1 мм рт. ст., а над линией, соответствующей интактному - порядка 1,5 мм рт. ст. Это указывает на то, что в случае апикального ослабления при оценке истинного давления поправку, связанную с неоднородностью, все же следует учитывать. Что касается наклона (равного коэффициенту ) линий, соответствующих неоднородной роговице, то он равен для апикального ослабления (линия 1) 0,92 мм рт. ст./г, а для периферийного (2) - 0,86 мм рт. ст./г. Для однородной роговицы с жесткостью, равной средней жесткости обоих неоднородных распределений = 0,84, а для интактной - 0,8 мм рт. ст./г. Таким образом, при апикальной неоднородности коэффициент эластоподъема отклоняется от случая «средней» роговицы незначительно, а при периферийной этим отклонением, безусловно, можно пренебречь. Во всех расчетах, рассмотренных выше, среднее значение жесткости роговицы получено как среднее по длине дуги образующей (первая формула (4)). Использование осреднения по площади роговичного сегмента (вторая формула (4)) несколько изменяет результаты, хотя и не слишком существенно. На рис. 8 показано влияние способа осреднения. Видно, что в случае апикального ослабления роговицы кривые, полученные при осреднении по поверхности, располагаются правее кривых с осреднением по дуге, что приводит к некоторому изменению разности между коэффициентом эластоподъема а б Рис. 8. Зависимость коэффициента эластоподъема от средней жесткости при разных способах осреднения: а - апикально ослабленная роговица; б - роговица, ослабленная в кольцевой области вблизи склеры. Кривые 1 - осреднение по длине дуги, 2 - по площади поверхности роговицы. Штрих-пунктир соответствует однородной роговице с той же средней жесткостью для неоднородной и однородной роговицы с тем же средним значением. Имеется пороговое значение жесткости роговицы, ниже которого эта разность меньше, чем при осреднении по дуге, а выше оказывается больше. В случае ослабления в кольцевой зоне взаимное расположение кривых противоположно, но в зоне вогнутости кривой (именно она показана на рис. 8, б) поведение упомянутой разности качественно такое же. Способ осреднения не влияет ни на форму кривых ни на их взаимное расположение при тех же значениях параметров, например, на рис. 6. Заключение Численно изучено влияние искусственно вызванного ослабления жесткости роговицы глаза на физический смысл результатов эластотонометрии по Маклакову. Проведенное исследование свидетельствует о незначительном отклонении значений коэффициента эластоподъема для неоднородной роговицы от значения этого коэффициента для однородной роговицы при ее умеренной жесткости. Это отклонение невелико при сравнении с измерением на интактной роговице и еще заметно меньше, если для сравнения выбирается однородная роговица с той же средней жесткостью. При этом способ осреднения (по длине дуги образующей или по площади поверхности) влияет на результат незначительно. При ослаблении роговицы в апикальной зоне такое отличие не превышало 13%, а при ослаблении в кольцевой зоне, близкой к склеральному кольцу, - 2%. Эта разница сопоставима с вариациями коэффициента эластоподъема вследствие изменений склеральной жесткости и истинного давления, а порой и значительно меньше их. Таким образом, к прооперированным глазам с исходно умеренно жесткой роговицей можно применять как известные ранее, так и предложенные нами [5] методики клинической оценки результатов эластотонометрии. Наши предварительные расчеты, однако, показывают, что использование этих методик при обследовании глаз с очень жесткой роговицей может привести к значительным ошибкам. Между тем, по современным данным, жесткость корнеосклеральной оболочки, включая роговицу, при глаукоме может существенно возрастать [5, 12]. Изучение особенностей тонометрии и эластотонометрии для таких глаз представляет собой отдельную задачу, которая требует подробного исследования.

Об авторах

И. Н Моисеева

Московский государственный университет им. М.В. Ломоносова

А. А Штейн

Московский государственный университет им. М.В. Ломоносова

Список литературы

Статистика

Просмотры

Аннотация - 76

PDF (Russian) - 37

Ссылки

• Ссылки не определены.