УДК 531/534: [57+61]
РОССИЙСКИЙ
Журнал Биомеханики
www.biomech.ru
ОЦЕНКА КОРРЕКТНОСТИ СТАНДАРТНОЙ МЕТОДИКИ РАСЧЕТА ВНУТРИГЛАЗНОГО ДАВЛЕНИЯ ПРИ ТОНОМЕТРИИ ПО ШИОТЦУ
И.Н. Моисеева, А.А. Штейн
Институт механики Московского государственного университета им. М.В. Ломоносова, Россия, 119192, Москва, Мичуринский проспект, 1, e-mail: [email protected], [email protected]
Аннотация. Проанализированы зависимости, применяемые при стандартной обработке данных процедуры измерения внутриглазного давления по Шиотцу. В тонометре Шиотца осуществляется совместное нагружение роговицы глаза вогнутой тяжелой стопой и тонким, более легким, вогнутым стержнем (плунжером). Нагружающее устройство имеет с роговицей две зоны контакта - всю вогнутую поверхность стержня и область стопы, определяемую внутренним давлением. При математическом моделировании использована предложенная ранее авторами модель, в которой роговица представлена безмоментной (мягкой) упругой поверхностью, а склеральная область заменяется упругим элементом, откликающимся изменением объема на изменение давления. Стандартная обработка данных тонометрии по Шиотцу базируется на двух эмпирически обоснованных алгоритмах: сначала по измеренному заглублению плунжера под стопу и известному весу плунжера определяется тонометрическое давление, а затем по последнему находится внутриглазное давление в ненагруженном глазу («истинное»). Показано, что использование стандартной методики вполне правомерно для тонометра с минимальным весом стержня 5,5 г, если упругие характеристики глазного яблока близки к средним. При большем весе стержня или значениях упругих констант, сильно отличающихся от нормы, необходимо учитывать поправки. Влияние геометрических параметров роговицы не слишком велико, однако при сильно увеличенном радиусе ее кривизны (например, в результате операции) поправки к результатам стандартного расчета следует учитывать.
Ключевые слова: глаз, математические модели, оболочка глаза, внутриглазное давление, тонометрия, тонометр Шиотца.
Введение
Прямое измерение внутриглазного давления в условиях стандартного офтальмологического исследования невозможно. Используется неинвазивная измерительная процедура - тонометрия, при выполнении которой давление обычно оценивается по деформациям глазного яблока под нагрузкой, в частности при статической тонометрии - в результате приложения груза к роговице [5]. Основная цель тонометрии - определение истинного внутриглазного давления, т.е. такого давления, которое имело место до измерительного воздействия. Попутно определяется также тонометрическое давление (давление в нагруженном глазу), которое может иметь и самостоятельную диагностическую ценность. При стандартной клинической
© Моисеева И.Н., Штейн А.А., 2015
Моисеева Ирина Никитична, к.ф.-м.н., старший научный сотрудник лаборатории общей гидромеханики, Москва
Штейн Александр Александрович, к.ф.-м.н., ведущий научный сотрудник лаборатории общей гидромеханики, Москва
тонометрии плоским штампом по Маклакову, например, несмотря на существование таблиц пересчета от тонометрического давления к истинному [5], в карту больного вносится именно тонометрическое давление. При тонометрии по Шиотцу, когда роговица испытывает вдавливающее (импрессионное) воздействие, наложенное на напряженное состояние, вызванное воздействием вогнутой тяжелой стопы, разница между тонометрическим и истинным давлением велика и расчет истинного давления (калибровка тонометра) - обязательная составляющая обследования.
Методика расчета истинного давления по показаниям тонометра Шиотца разрабатывалась на основании эмпирических и полуэмпирических исследований в 30-70 гг. прошлого века [6, 8, 9], и для нее принят международный стандарт. Цель настоящей работы - оценить уровень возможных ошибок при определении истинного и тонометрического внутриглазного давления на основе стандартной методики в зависимости от индивидуальных упругих и геометрических характеристик глазного яблока. Такая оценка необходима в свете нескольких обстоятельств. Известно, что упругие характеристики глаза могут у разных людей различаться на порядок и более [10]. Кроме того, упругие свойства могут по-разному распределяться между роговичной и склеральной областями и изменяться вследствие заболеваний, например глаукомы. А в результате хирургических операций значительно изменяются и геометрические параметры глаза, в первую очередь кривизна роговицы.
В качестве аппарата для исследования будет использована предложенная в работе [3] модель глазного яблока, в которой упругие свойства задаются тремя упругими константами, из которых две оказываются существенными, характеризуя отдельно упругие свойства роговичной и склеральной областей. Модель [3] достаточно проста и не требует знания множества параметров, определение которых для индивидуального обследуемого глаза в клинике невозможно. Вместе с тем она позволяет учесть неоднородность глазного яблока как упругой системы. Роговица представляется мягкой (безмоментной) двумерной однородной изотропной линейно упругой поверхностью, а склера и прилегающие к ней ткани - упругим элементом, откликающимся изменением объема на изменение давления. Пренебрежение сопротивлением роговицы изгибу обосновано эмпирически и впервые было реализовано в механически корректной постановке в работе [1].
Разработанный в статье [3] метод был эффективно применен к анализу механического смысла данных, получаемых с помощью тонометров, в которых к роговице приложен единственный груз. Форма груза может быть разной - это плоский штамп (тонометр Маклакова) [3], выпуклый штамп (конвекс-тонометр Краснова) [2], вогнутый штамп [8], тонкий стержень [4]. Нагружение глазного яблока в тонометре Шиотца отличается от всех этих ситуаций тем, что в этом тонометре осуществляется комбинированное нагружение вогнутым штампом и стержнем.
Тонометр Шиотца: устройство и стандартная методика интерпретации данных
Схема нагружения в импрессионном тонометре Шиотца представлена на рис. 1. Роговица лежащего пациента первоначально нагружается так называемой стопой -вогнутым грузом с характеристиками, одинаковыми или близкими для всех тонометров рассматриваемого типа: вес стопы равен 11,5 г, а радиус кривизны - 15 мм, т.е. приблизительно вдвое больше, чем у нормальной роговицы человека. В центральной части стопы имеется отверстие, в которое свободно опускается грузик в виде тонкого стержня (плунжера). Диаметр стержня 3 мм, его вдавливающаяся в роговицу поверхность так же, как и стопа, вогнута и имеет тот же радиус кривизны. Вес стержня варьируется. Стандартный вес 5,5 (наиболее употребителен); 7,5; 10; 15 г.
О внутриглазном давлении судят по глубине погружения стержня под контактную поверхность стопы.
В отличие от тонометра Маклакова, калибровка тонометра Шиотца, т.е. разработка алгоритма пересчета измеряемой геометрической характеристики в давление, невозможна по элементарным формулам даже приблизительно и потребовала разработки специальной теории.
Стандартная обработка данных тонометрии по Шиотцу базируется на двух алгоритмах: сначала по измеренной глубине погружения И и известному весу плунжера О определяется тонометрическое давление р, а затем по р находится истинное внутриглазное давление р0.
Зависимость р (И) принимается в виде
Рг = 0 . , (1)
С: + С 2 И
где С и Сг - раз и навсегда определенные эмпирические константы, С = 0,107, С2 = 0,0138. Формула (1) предполагает фиксированные внесистемные единицы измерения: О измеряется в граммах силы; р - в миллиметрах ртутного столба; И -
в единицах Шиотца (1 ед. Шиотца = 50 мкм). Формула чисто эмпирическая, она основывается на экспериментальных исследованиях, предшествовавших тарировке тонометра Шиотца и показавших, в частности, слабую зависимость функции р (О, И) от упругих свойств оболочки глаза и линейных размеров, характеризующих глазное яблоко, причем не только для человека, но и для нескольких изученных видов животных [7]. Формулу (1) легко интерпретировать на языке теории размерности: из предположения, что существует связь вида р = р( (О, И) (с дополнительным присутствием в числе аргументов геометрических параметров - размерности длины и безразмерных), общий вид такой связи И = / (О / р), а конкретная форма (1) означает
эмпирически установленную линейную зависимость И от величины обратной р . Функция р( (И) в этом случае представляется гиперболой.
Рис. 1. К постановке задачи: а - глазное яблоко до нагружения; б - роговица, нагруженная тонометром Шиотца: 1 - плунжер; 2 - стопа; 3 - роговица. Остальные
обозначения в тексте
Расчет р по р, вычисленному по И в соответствии с формулой (1), базируется на полуэмпирических представлениях, разработанных Фриденвальдом. Для этого вводится так называемый смещенный объем V - геометрическая характеристика, однозначно (по специальному правилу) определяемая по искривленной поверхности роговицы после наложения груза [6]. Для данного обследуемого глаза (с некоторыми фиксированными геометрическими и механическими свойствами, включая давление до нагружения р) и фиксированного нагружающего устройства связь между И и Vd
однозначна; поэтому должна существовать связь р / р =ф (V). Эта связь предполагается универсальной (т.е. не зависящей от индивидуальных характеристик глаза) и экспоненциальной по V :
р / р = 10а^ , (2)
-3
где а = 0,0215 мм - константа Фриденвальда, найденная им как средняя по множеству обследованных глаз.
Неявно считается, что при одинаковых глубине погружения плунжера И и весе плунжера О (а следовательно, с учетом принятой зависимости (1), при одинаковом тонометрическом давлении р ) геометрия искривленной роговицы одинакова. Отсюда
следует однозначная связь V (Рг, О), которая находится эмпирически и в логарифмических координатах аппроксимируется линейными функциями = \gVrf (1§ р) с разными коэффициентами для разных весов плунжера. В частности,
при О = 5,5 и 7,5 г результат такой аппроксимации имеет вид (V, мм3; р, мм рт. ст.)
О = 5,5 г: 1вV, = 2,016(2,029- 1вр), О = 7,5 г : = 2,174(2,092-^р).
Константы, присутствующие в формуле (3), как и константа Фриденвальда в зависимости (2), найдены как средние для некоторой обширной популяции испытуемых.
Зависимости (1)-(3) непосредственно зашиты в программное обеспечение для используемого в отечественной медицинской практике тонографа G/aw7esí-60.
Постановка задачи и метод решения
Общая постановка задачи подробно описана и обоснована в работе [3]. Она включает в себя систему обыкновенных дифференциальных уравнений для описания статических деформаций нагруженной извне оболочки (роговицы), которая дополняется конечным соотношением для склеральной части внутриглазного объема, находящейся вне роговицы (т.е. ниже плоскости AB на рис. 1).
Рассматриваем осесимметричную деформацию поверхности (роговицы), распертой внутренним давлением p, причем предполагается, что при некотором
давлении p она имела форму сферического сегмента радиуса R (рис. 1, а) с радиусом основания R и углом полураствора фс (sin фс = R / R). В области свободной роговицы (вне контакта с приложенной нагрузкой) уравнения равновесия и состояния для роговицы как мягкой двумерной упругой поверхности в предположении малых деформаций с учетом конечности смещений имеют вид
г
(r0T1) " T2C0S Ф= 0 ,
—T sin ф — r0ф'^ + Xr0p = 0,
r' = X cos ф, (4)
r0 = Re sin Фo, фо = 4" , X = 1 +S1, S2 = — — 1 ,
Re Г0
S1 = Tí [T1 — VT — T0(1 — V )] , S2 = "ir [T2 — VeT1 — T0(1 — V )] .
ee
Здесь независимым параметром считается длина s дуги образующей в ненагруженном состоянии, отсчитываемая от ее пересечения с осью симметрии ( —фсRc < s <фсRc ); T и T - касательные усилия в меридиональном и перпендикулярном к нему направлениях; r - расстояние точки поверхности, представляющей роговицу, от оси; ф - угол между нормалью к этой поверхности и
осью симметрии (рис. 1, б); S и s2 - деформации в направлении главных осей; E* и V - модуль растяжения и коэффициент Пуассона роговицы как двумерной упругой поверхности [3]. Угол ф считается положительным, если нормаль направлена во внешнюю относительно оси сторону. Функции r0 С s) и ф0 ( s) - значения функций r С s ) и ф^) при начальной сферической конфигурации (см. рис. 1, а). В этой конфигурации усилия T и Г2 равны между собой и T = T = T = PoR / 2.
В случае тонометра Шиотца постановка задачи должна учитывать две области, в которых роговица свободна и ее поведение управляется уравнениями (4), а также две области, в которых она прилегает к жестким поверхностям. Эти четыре области показаны на рис. 1 , б. В центральной области I (в окрестности точки s = о ) роговица контактирует с вогнутой поверхностью плунжера, в области III - с вогнутой поверхностью стопы, в областях II и IV роговица свободна. В предположении свободного проскальзывания роговицы по поверхностям контакта система (4) остается справедливой в областях I и III за исключением проекции уравнения равновесия на нормальное к роговице направление (второе уравнение (4)): это соотношение теперь должно включать силу нормальной реакции, действующую со стороны груза. Если не интересоваться распределением силы реакции, уравнение (4) можно не рассматривать и заменить его геометрическими связями: под плунжером (область I)
r
sin ф = —, (5)
гр
где r - радиус кривизны плунжера, и под стопой (область III)
r
sin ф = —, (6)
rf
где r - радиус кривизны стопы.
В силу симметрии решение задачи ищется в области 0 < s < фсR. Система уравнений - третьего порядка для свободной роговицы и второго - под плунжером и под стопой. Поскольку три значения s, соответствующие границам между областями, заранее неизвестны, требуются 13 дополнительных условий. Эти условия выставляются
как на внешних границах области решения (s = 0 и s = фсR), так и на неизвестных внутренних границах. В качестве граничных условий на концах области принимаем (см. [3]) условие r(фсR) = R, соответствующее закреплению роговицы на склере, деформации которой считаются малыми, и условие ограниченности усилия Т в точке пересечения с осью симметрии. Последнее принимает форму асимптотического соотношения [3]. Остальные условия формулируются на границах отрыва роговицы от поверхностей. С учетом неразрывности поверхности роговицы, а также условий свободного проскальзывания (сила реакции нормальна к поверхности контакта) и исчезновения силы реакции в точках отрыва получаем следующие условия:
r = r+ = r*, Т~ = Т+,
I—II: 1 1 ; (7) G = лг» • p - Т si^-ф ) • 2лг»
II—III: r = r+, ф =ф+, Т- = Т+; (8)
r = r+, ф- =ф+, Т1- = т:+,
III—IV: 1 1 (9) G + Gf = лг • p - Т sin ф • 2лг.
Здесь G и Gf — вес стержня и стопы соответственно; r — радиус плунжера, а
индексами «—» и «+» обозначены значения переменных по разные стороны линии отрыва, причем «—» соответствует области, лежащей от оси дальше. В соотношениях, содержащих вес (условиях равновесия), функции, непрерывные на рассматриваемой поверхности разрыва, индексами «—» и «+» не помечены. Каждая из систем соотношений (7)—(9) соответствует границе между указанными областями. Граничные условия (7) сформулированы ранее в работе [4] при решении задачи о чисто стержневом нагружении роговицы.
Система уравнений (4)—(5) с указанными граничными условиями дает возможность вычислять геометрические характеристики деформированной под грузом роговицы по известным p, G и Gf. В числе таких характеристик приращение AVC
подроговичного объема (часть внутриглазного объема над плоскостью AB), а также расстояния от точек деформированной роговицы до плоскости AB, необходимые для расчета измеряемого погружения плунжера. Таким образом, в силу постановки задачи геометрия деформированной роговицы определяется, помимо внутреннего давления, только геометрическими и механическими характеристиками самой роговицы. Чтобы связать AVC с приращением полного заполненного жидкостью внутриглазного объема
AV, вводится склеральный объем (под AB) с приращением AV, подчиняющимся соотношениям
AVs = K (p - Po), AV = AV +AVs, (10)
где K — интегральная константа, определяющая упругое поведение всей склеральной области и зависящая не только от механических, но и от геометрических характеристик, в частности от ее объема. При сохранении механических свойств склеры и ее окружения больший объем внутрисклерального пространства (например, при миопии) приведет к увеличению K, т.е. к уменьшению эффективной жесткости склеральной области.
В рамках модели свойства объекта (глазного яблока) определяются пятью параметрами. В их числе три, характеризующие его упругие свойства (но включающие в
себя зависимость и от геометрических характеристик): Е*, ус и К. Два других - чисто геометрические, характеризующие форму роговицы: ее радиус Лс и радиус основания роговицы Яь). Как продемонстрировано в работе [3], вместо линейной жесткости роговицы Е* удобно ввести ее поверхностную жесткость Е*^ = Е* / (1 — ус ) . В задачах, где главную роль играет изотропное растяжение роговицы (а именно таковы задачи тонометрии), влияние коэффициента Пуассона ус отдельно от Е* несущественно.
Связь этих констант с традиционно используемыми обсуждалась ранее [3]. Для удобства сравнения с данными других авторов были формально определены роговичная и склеральная жесткости Ес и Е, измеряемые в мегапаскалях и отличающиеся от Е* и К —1 фиксированными размерными множителями: Е* =8ЕС, К —1 =СЕ, где 5 = 0,5 мм - средняя в норме толщина роговицы, а £ = 1,44-105 мм-3. Эти величины имеют простой физический смысл в рамках сильно упрощенных моделей (см., например, [3]). Расчеты в основном проводились при стандартных для человеческого глаза значениях Кс = 8 мм, Яь = 5,8 мм. Изменение коэффициента
Пуассона роговицы приводит лишь к переопределению Е , по той же причине (удобство сравнения с другими данными) он выбирался равным 0,45, что соответствует в случае трехмерной изотропной оболочки почти несжимаемому материалу. Выполнялись еще расчеты при других значениях геометрических параметров: Яс=7...9 мм, Яь =5,4... 6,2 мм, а также при коэффициентах Пуассона ус =0,25... 0,65. В последнем случае сравнение проводилось с расчетами для
измененного значения линейного модуля Ес' = Ес(1 — ус)/(1 — 0,45) , т.е. при том же поверхностном упругом модуле Есж = Ес / (1 — ус ). Подтверждена слабая зависимость результатов от коэффициента Пуассона при таком выборе параметров.
В работе [3] оценены «средние» по множеству относительно здоровых (неглаукомных) глаз значения упругих констант: Е3 = 10 МПа, Ес = 0,3 МПа. Всюду
далее, когда говорится о средних (нормальных) значениях, имеются в виду именно вышеуказанные. С учетом высокого индивидуального (в том числе вследствие патологий) разброса влияние этих констант исследовалось в диапазонах Е6, =2...20 МПа, Ес =0,1...0,7 МПа.
Единственный параметр нагружающего устройства, который варьировался, - вес плунжера О. Расчеты проводились для О = 5,5 и 7,5 г. Остальные параметры тонометра, как и в медицинской практике, выбирались следующими: радиус контактной поверхности стопы ^ = 15 мм, радиус вогнутого основания плунжера
г = ^ = 15 мм, вес стопы О = 11,5 г.
В основе рассмотрения было изучение функции полного объема ДV = ДV(р,О; Ес,Е), состоявшее из трех этапов: определение А^ из решения краевой задачи для роговицы, а затем А^ и ДV по формулам (10). Формально в числе аргументов функции АV должно присутствовать давление р0, соответствующее условно выбранной начальной сферической конфигурации, однако его влияние на исследуемые зависимости мало [6]; поэтому этот аргумент явно не указан. При численном решении применялись метод стрельбы и схема Рунге-Кутты четвертого порядка.
Связь глубины погружения плунжера с тонометрическим давлением
В случае тонометра Шиотца измеряемая геометрическая характеристика -глубина погружения плунжера И, отсчитываемая от внутренней поверхности стопы. При стандартной обработке по формуле (1) считается, что на связь ее с тонометрическим давлением влияет только вес плунжера.
В рамках используемой модели тонометрическое давление есть функция не только О, И и геометрических параметров роговицы Яс и Яь, но и упругости
роговицы, а также веса стопы: р = р (О, И, Яс, ^, Ес, О/ ) . Зависимость от склеральной
жесткости отсутствует в силу постановки задачи для роговицы (в которой деформация склеры не учитывается). Поскольку вес стопы в тонометре Шиотца всегда одинаков, далее исследуется влияние на результат восстановления значения тонометрического давления по величине И только параметров £с, Яс и Яь.
На рис. 2 представлена зависимость тонометрического давления р от глубины погружения плунжера И для различных жесткостей роговицы £с и двух стандартных весов плунжера в сравнении с эмпирической формулой (1). Рассчитанные кривые подтверждают как весьма слабую зависимость кривых р (И) от £с, так и пригодность в первом приближении формулы (1). Во всем практически важном диапазоне тонометрическое давление, соответствующее данной фиксированной И, меньше, чем определяемое формулой (1) (кроме давлений ниже 20 мм рт. ст.) и растет с ростом жесткости роговицы. В области низких тонометрических давлений (ниже 25 мм рт. ст.) совпадение всех кривых (при разных £с и эмпирической) практически полное, однако
с ростом давления отклонение нарастает и для О = 5,5 г при р ~30 мм рт. ст. имеет порядок 3 мм рт. ст., а для О = 7,5 г прир ~ 35 .. .50 мм порядок 5 мм рт. ст. Для груза О = 5,5 г с ростом р выше 30 мм рт. ст. наблюдается сокращение разницы результатов, рассчитанных по модели и полученных по эмпирической формуле, и при 40-50 мм рт. ст. они вновь практически совпадают.
Влияние радиуса роговицы на зависимость тонометрического давления от заглубления представлено на рис. 3, где приведены результаты расчетов зависимостей р (И) при разных радиусах роговицы Яс для О = 5,5 г и Ес = 0,3 МПа. Эта величина может отличаться от стандартного для человека значения в результате операций на роговице. Также она иная и у детей. Возрастание радиуса кривизны роговицы приводит к уменьшению реального тонометрического давления при той же И. Таким образом, при увеличенном радиусе кривизны использование зависимостей, полученных при стандартном значении радиуса Яс, дает завышенный результат, что указывает на возможность ошибок измерения на прооперированных глазах. Это расхождение особенно заметно при достаточно малых И , когда возможна ошибка до 3 мм рт. ст. Изменение радиуса кривизны роговицы на 1 мм практически не сказывается на кривых р(И) при р ~ 35 мм рт. ст. и выше и достигает 2-3 мм рт. ст. при р ~ 20 мм рт. ст.
Вычисление зависимостей р (И) при разных радиусах основания роговицы Яь =5,4... 6,2 мм продемонстрировало практически полное совпадение соответствующих кривых. Таким образом, влияние параметра Яь на вид зависимости р (И) можно не учитывать.
Рис. 2. Зависимость тонометрического давления pt от веса плунжера и жесткости роговицы. Семейство кривых 1 соответствует весу плунжера G = 5,5 г, семейство 2 -весу G = 7,5 г. Сплошные кривые: нижняя в каждом семействе - Ec = 0,2 МПа, верхняя - E = 0,6 МПа; пунктирные - рассчитанные по стандартной формуле (1)
для соответствующего веса G
Рис. 3. Зависимость тонометрического давления pt от радиуса кривизны роговицы R при весе плунжера G = 5,5 г и жесткости роговицы Ec = 0,3 МПа: 1 -R = 7 мм; 2 - R = 8 мм; 4 - R = 9 мм; пунктирная кривая 3 построена по
стандартной формуле (1)
В целом расчеты показывают, что формула (1) вполне пригодна для оценки тонометрического давления, но в некоторых случаях все же дает заметные отклонения от реальных значений р . Эти отклонения наиболее заметны на интервале умеренно больших р (разном для разных весов плунжера) и увеличиваются при возрастании веса плунжера и радиуса кривизны роговицы.
Определение истинного давления по тонометрическому
Задача тонометрии решалась численно подбором такого значения давления р ,
при котором приложение нагрузки не меняет объем: АV = 0. В этом случае р0
соответствует внутриглазному давлению до нагружения. Целью исследования было изучить влияние индивидуальных механических и геометрических характеристик глаза на зависимость между р и р (тонометрическую зависимость) и оценить область
приложимости формул (2), (3).
На рис. 4, 5 приведены результаты расчетов тонометрических зависимостей при различных роговичной и склеральной жесткостях для О = 5,5 г. Представлены зависимость истинного давления р от тонометрического р (см. рис. 4) и тонометрической разности Ар = р — р0 от р0 (см. рис. 5). Тонометрическая разность зависит как от уровня давления, так и от обеих жесткостей. С ростом давления эта величина уменьшается (см. рис. 5).
Увеличение склеральной жесткости Е (рис. 4, а) приводит к уменьшению истинного давления р , соответствующего измеренному тонометрическому давлению р. Различие значительно: при изменении Е от 6 до 20 МПа (Ес = 0,3 МПа) оно достигает 6 мм рт. ст. Практически важно, что Е слабо сказывается на форме тонометрических зависимостей: кривые р0(р) при разных Е визуально подобны и различаются лишь сдвигом по оси ординат.
Совсем по-иному сказывается жесткость роговицы Ес (рис. 4, б). Возрастание этой жесткости приводит уже к возрастанию р , соответствующего некоторому измеренному р. Как функция от р0 тонометрическая разность Ар и в этом случае убывает, но существенно неравномерно для разных р0: при Е3 = 10 МПа изменение Ес от 0,2 до 0,7 МПа практически не влияет на Ар при р0 = 10 мм рт. ст. (различие не более 0,3 мм рт. ст.), но с возрастанием р0 это влияние становится все более заметно, и при р = 30 мм рт. ст. различие достигает 3,5 мм рт. ст. Визуально это проявляется в том, что семейство кривых р0 (р) при разных Ес расходится веером, расширяющимся в сторону возрастания давлений (см. рис. 4, б).
Результаты расчета связи р0 (р) по нашей модели при средних значениях
упругих констант, т.е. близких к Е3 = 10 МПа, Е3 = 0,3 МПа, и весе плунжера О = 5,5 г практически совпадают с результатом применения стандартных формул (2), (3) (рис. 6). Однако совпадение имеет место лишь для нормальных глаз: при повышенной жесткости роговицы стандартный расчет по тонометрическому давлению занижает, а при повышенной склеральной жесткости завышает истинное давление.
a
б
Рис. 4. Влияние упругих свойств роговицы и склерального сегмента на зависимость истинного давления p0 от тонометрического pt : a - Ec = 0,3 МПа, Es = 6 МПа
(кривая 1), E = 8 МПа (кривая 2), Es = 10 МПа (кривая 3), Es = 12 МПа
(кривая 4), E = 16 МПа (кривая 5), Es = 20 МПа (кривая 6); б -Es = 10 МПа,
E = 0,7 МПа (кривая 1), Ec = 0,6 МПа (кривая 2), Ec = 0,5 МПа (кривая 3),
E = 0,4 МПа (кривая 4), Ec = 0,3 МПа (кривая 5), Ec = 0,2 МПа (кривая 6)
Рис. 5. Влияние упругих свойств роговицы и склерального сегмента на зависимость тонометрической разности Ар от истинного давления р0: кривая 1 -
Е = 0,3 МПа, Е = 12 МПа; 2 - Е = 0,3 МПа, Е = 10 МПа; 3 - Е = 0,3 МПа,
с" " с ~ ' э ' с ' '
Е = 12 МПа; 4 - Е = 0,3 МПа, Е = 10 МПа
При весе плунжера О = 7,5 г расчетная кривая р0 (р,) дает при нормальных упругих константах значения истинного давления приблизительно на 3 мм рт. ст. выше, чем стандартная методика. При значениях упругих констант, отличных от нормальных, наблюдаются те же закономерности, что и при О = 5, 5 г. Таким образом, для О = 7, 5 г
Рис. 6. Сопоставление зависимостей истинного давления р0 от тонометрического р, рассчитанных по принятой модели, с зависимостями, получаемыми по стандартным формулам (2), (3) при разных весах плунжера. Кривые 1-4 соответствуют весу О = 5,5 г: 1 - Ес = 0,6 МПа, Е = 10 МПа; 2 - расчет
по формулам (2), (3); 3 - Ес = 0,2 МПа, Е = 10 МПа; 4 - Ес = 0,2 МПа, Е = 12 МПа; кривые 5-7 соответствуют весу О = 7,5 г: 5 - Ес = 0,2 МПа, Е = 10 МПа; 6 - Е = 0,2 МПа, Е = 12 МПа; 7 - расчет по формулам (2), (3)
повышенная жесткость роговицы и (или) пониженная склеральная жесткость увеличивают отклонение от результатов стандартного расчета, а повышенная склеральная жесткость его уменьшает.
Тонометрические кривые претерпевают некоторые изменения при изменении геометрических параметров роговицы (форма и размеры склеры учитываются условной склеральной жесткостью Е ). На рис. 7 представлена зависимость р(р) для разных радиусов кривизны роговицы Ес при средних значениях упругих параметров. Увеличение радиуса кривизны приводит к некоторому повышению р , соответствующего тому же р . При малых давлениях ( р0 ~ 10 мм рт. ст.) изменение Ес на 1 мм соответствует изменению р порядка 2 мм рт. ст. С увеличением р (или р ) влияние Ес ослабевает и при р0 > 20 мм рт. ст. практически пропадает.
Влияние радиуса основания роговицы более равномерно (рис. 8): при возрастании кривые р0 (р) остаются визуально подобны, сдвигаясь кверху вдоль оси ординат. При увеличении на 0,8 мм от 5,4 до 6,2 мм давление р0, соответствующее заданному р, возрастает на величину порядка 2 мм рт. ст. Следует, впрочем, иметь в виду, что радиус у человека очень стабилен. Различия при его
определении связаны не с индивидуальным разбросом, а с анатомической привязкой границы роговицы. В рамках рассматриваемого огрубленного моделирования выбор значения в известной мере условен и его небольшое изменение фактически сведется к некоторому переопределению упругих констант.
30
25
н
^ 20 ■ о. s s
¿15-
10
А А
// /
Ж/ У
3 21
20 24 28 32 36 40 44 pt, мм рт. ст.
Рис. 7. Влияние радиуса кривизны роговицы на зависимость истинного давления p0 от тонометрического pt при G = 5,5 г, Ec = 0,3 МПа, Es = 10 МПа. Сплошные
кривые: 1 - R = 7 ; 2 - R = 8; 3 - R = 9 мм; пунктирная кривая - расчет по
стандартным формулам (2), (3)
32 28
S324'
¡20
о
16 12
/уу* v/s г
///У 7/
ft
21
20 24 28 32 36 40 44 р{, мм рт. ст.
Рис. 8. Влияние радиуса основания роговицы на зависимость истинного давления р0 от тонометрического р при О = 5,5 г, Ес = 0,3 МПа, Ех = 10 МПа. Сплошные
кривые: 1 - Яь = 5,4; 2 -Яь = 5,8; 3 -Яь = 6,2 мм; пунктирная кривая - расчет
по стандартным формулам (2), (3)
Таким образом, формулы, используемые при стандартной обработке данных, вполне применимы, если упругие характеристики глазного яблока близки к средним. При наличии существенных отклонений этих величин от средних значений оценка по
стандартным формулам может давать существенную ошибку. Кроме того, ошибка возрастает для плунжера большего веса. Возможно, применение большего веса плунжера правомерно для глаз с большой склеральной жесткостью.
Выводы
Разработанная авторами ранее простая механическая модель глазного яблока позволила механически корректно исследовать вопрос о применимости стандартных формул, используемых при обработке данных, получаемых при тонометрии по Шиотцу. Расчеты показали, что оценка тонометрического давления для тонометров с весами плунжера 5,5 и 7,5 г по стандартной методике в целом правомерна, хотя в некоторых диапазонах давлений ошибка может и увеличиваться. Оценка истинного давления для тонометра с минимальным весом плунжера 5,5 г по стандартным формулам соответствует расчетным данным при средних значениях упругих констант, но заметно отклоняется от них, если эти константы сильно отличаются от средних значений. При большем весе плунжера (7,5 г) расчетные кривые заметнее отличаются от полученных по стандартной методике даже и для средних значений констант. Влияние геометрических параметров роговицы не слишком велико, однако при сильно увеличенном радиусе ее кривизны (например, в результате операции) поправки к результатам стандартного расчета следует учитывать.
Вопрос о диагностической ценности различных оцениваемых при тонометрии параметров (глубины погружения плунжера, истинного и тонометрического давления и т.д. ) должен служить предметом медицинских исследований.
Благодарности
Авторы благодарны Г.А. Любимову за ценные замечания. Работа выполнена при поддержке гранта РФФИ № 14_01_00475.
Список литературы
1. Бауэр С.М., Любимов Г.А., Товстик П.Е. Математическое моделирование метода Маклакова измерения внутриглазного давления // Изв. РАН. Мех. жидк. и газа. - 2005. - № 1. - С. 24-39.
2. Любимов Г.А., Моисеева И.Н., Штейн А.А. Исследование свойств двухкомпонентной механической модели глазного яблока и возможности ее использования при практической оценке механических свойств глаза человека // Изв. РАН. Мех. жидк. и газа. - 2014. - № 6. - С. 5-16.
3. Моисеева И.Н., Штейн А.А. Анализ зависимости давление - объем для глазного яблока, нагруженного плоским штампом, на основе двухсегментной упругой модели // Изв. РАН. Мех. жидк. и газа. - 2011. - № 5. - С. 3-15.
4. Моисеева И.Н., Штейн А.А. Исследование зависимости давление - объем для глазного яблока при нагружении тонким стержнем // Изв. РАН. Мех. жидк. и газа. - 2013. - № 4. - С. 3-14.
5. Нестеров А.П., Бунин А.Я., Кацнельсон Л.А. Внутриглазное давление. Физиология и патология. - М.: Наука, 1974. - 381 с.
6. Friedenwald J.S. Contribution to the theory and practice of tonometry // Am. J. Ophthalmol. - 1937. -Vol. 20, № 10. - P. 985-1024.
7. Friedenwald J.S. Some problems in the calibration of tonometers // Am. J. Ophthalmol. - 1948. - Vol. 31. - P. 935-944.
8. Friedenwald J.S. Standartization of tonometers; decennial report by the committee on standartization of tonometers. - Omaha: Am. Acad. Ophthalmol., 1954. - 177 c.
9. Moses R.A. Theory of the Schi0tz tonometer and its empirical calibration // Trans. Am. Ophthalmol. Soc. -1971. - Vol. 69. - P. 494-562.
10. Pallikaris I.G., Kymionis G.D., Ginis H.S., Kounis G.A., Tsilimbaris M.K. Ocular rigidity in living human eyes // Invest. Ophthalm. Vis. Sci. - 2005. - Vol. 46, № 2. - P. 409-414.
ESTIMATING THE CORRECTNESS OF THE STANDARD METHOD USED IN CALCULATING THE INTRAOCULAR PRESSURE FROM SCHI0TZ TONOMETRY
I.N. Moiseeva, A.A. Stein (Moscow, Russia)
The dependences used in the standard processing of the Schi0tz tonometry data for estimating the intraocular pressure are analyzed. In the Schi0tz tonometer, the eye cornea is loaded together by a heavy concave footplate and a thin lighter concave rod (plunger). The loading device contacts the cornea in two zones: the entire concave rod surface and a footplate region determined by the internal pressure. The mathematical modelling uses the model proposed by the authors earlier, in which the cornea is represented by a momentless (soft) shell and the scleral region is modeled by an elastic element that responses to changes in pressure by changes in volume. The standard processing of the Schi0tz tonometry data is based on two algorithms justified empirically: first, from the measured depth of plunger indentation under the footplate and the known plunger weight the tonometric pressure is determined and then, from the latter the intraocular pressure in the unloaded eye ("true" pressure) is calculated. It is shown that the standard method is quite acceptable for the tonometer with a minimum rod weight of 5.5 g if the elastic characteristics of the eyeball are close to average values. For a greater rod weight or elastic constants substantially different from the normal values, it is necessary to take corrections into account. The influence of the geometric parameters of the cornea is moderate but for corneas with strongly increased radius of curvature, for example, owing to surgery, corrections to the standard calculation results should be taken into account.
Key words: eye, mathematical models, eye shell, intraocular pressure, tonometry, Schi0tz tonometer.
Получено 27 июня 2015