DIFFERENTIAL TONOMETRY WITH A SCHIØTZ TONOMETER: MATHEMATICAL MODELING WITH ACCOUNT FOR THE NONLINEARITY OF THE ELASTIC BEHAVIOR OF THE CORNEA AND COMPARISON WITH CLINICAL DATA
- Authors: Moiseeva I.N1, Stein A.A1
- Affiliations:
- Institute of Mechanics, Moscow State Lomonosov University
- Issue: Vol 26, No 4 (2022)
- Pages: 10-18
- Section: Articles
- URL: https://ered.pstu.ru/index.php/rjb/article/view/3594
- DOI: https://doi.org/10.15593/RZhBiomeh/2022.4.01
- Cite item
Abstract
Mathematical modeling of the diagnostic procedure of differential tonometry based on loading the eye with a Schiotz tonometer using two or more weights of the plunger, one of the two structural elements of this tonometer, is performed. The model of the eyeball developed earlier by the authors is used, based on the representation of the cornea as a momentless two-dimensional elastic surface and the scleral region as a zero-dimensional elastic element that responds with a change in pressure to a change in the volume of the fluid contained in it. Unlike previous works, in which a theoretical study of differential tonometry was carried out based on a physically linear model of the cornea, the nonlinearity of its elastic behavior is taken into account on the basis of an effective approach that made it possible to describe this nonlinearity by means of a single parameter. It is shown that the dependence of the tonometric pressure on the weight of the plunger remains almost linear, which makes it possible to introduce the pressure difference coefficient γ as the slope of the corresponding straight line. This coefficient increases with both scleral and corneal stiffness and decreases with true pressure. Accounting for nonlinearity leads to increase in the calculated values of γ at sufficiently high both true pressures and stiffnesses of the two components of the eyeball, as well as to an expansion of the pressure range for which these values are large. It is shown that two groups of eyes can be distinguished, in one of which, characterized by high values of the pressure difference coefficient, both components of the corneoscleral coat are very rigid, and in the other, where these values are significantly lower, the stiffness of the eye tissues can be considered normal. At the same time, an intermediate range of γ values remains, for which such definite conclusions cannot be drawn only from the result of differential tonometry. However, the ability to isolate a group of very stiff eyes that are particularly likely to have glaucoma makes this procedure useful in diagnosing this disease.
Full Text
Уже в первой половине прошлого века, на заре развития тонометрии, стало ясно, что показания тонометров сильно зависят от упругих свойств глазных тканей, которые значительно различаются для разных глаз. Тогда же была сформулирована задача: разработать методику оценки этих свойств на основе неинвазивных клинических измерений. Цель такой оценки с самого начала была двоякой: во-первых, учесть индивидуальную жесткость глазного яблока для более точного определения внутриглазного давления и, во-вторых, получить дополнительный диагностический показатель. Дело в том, что уже тогда для врачей стало ясно, что повышенная жесткость глаза коррелирует с большей вероятностью глаукомы и не менее значима для ее прогноза, чем внутриглазное давление. Для решения этой задачи были разработаны процедуры эластометрии (впоследствии переименованной в эластотонометрию) с использованием аппланационного тонометра Маклакова и дифференциальной тонометрии с использованием импрессионного тонометра Шиотца. Обе процедуры основываются на результатах нагружения роговицы несколько раз (два или более) одним и тем же тонометром, но грузами разных весов. Предполагалось, что разница между тонометрическими давлениями под фиксированной парой грузов тем больше, чем выше так называемая ригидность. Под последней подразумевалась по-разному и не слишком отчетливо определявшаяся величина, которая истолковывалась как некоторая интегральная характеристика жесткости глаза. В настоящее время обе процедуры используются редко, во многом потому, что стала понятна неоднозначная связь их результата с жесткостью. Подробное исследование их физического смысла стало возможным, когда мы применили к изучению эластотонометрии и дифференциальной тонометрии разработанный нами эффективный подход, базирующийся на удобной в расчетах - сильно упрощенной, но тем не менее содержательной - механической модели глазного яблока, основанной на представлении роговицы безмоментной двумерной поверхностью, а склеральной области - нульмерным элементом, отвечающим изменением давления на изменение объема содержащейся в нем жидкости [4; 9]. Небольшое число подлежащих определению параметров отличает развиваемый подход от исследований, основанных на учете множества параметров, характеризующих упругие свойства глазных тканей [1; 2; 13; 21; 32]. Пространственная неоднородность этих свойств хорошо известна и составляет предмет тщательного изучения [17; 18; 25] Однако, поскольку упругие свойства сильно различаются для разных глаз [33] и не могут быть заранее определены для конкретного исследуемого глаза, получаемые выводы в этом случае неизбежно сводятся лишь к общим рекомендациям. Первоначально мы считали оба компонента модели линейно упругими. В предположении изотропии и пространственной однородности роговицы как поверхности такая система характеризуется тремя упругими константами, из которых существенны две, характеризующие раздельно жесткость роговицы и склеральной области. Было показано, что и в случае тонометра Маклакова, и в случае тонометра Шиотца зависимость тонометрического давления от веса приложенного груза близка к линейной, что делает нагружение более чем двумя грузами информационно бессмысленным и полезным лишь для уточнения результатов. Наклон линии, представляющей эту зависимость, т.е. отношение разности давлений к разности весов (коэффициент эластоподъема) оказался существенно зависимым как от жесткостей роговицы и склеральной области, так и от давления в ненагруженном глазу (истинного давления), подлежащего определению при тонометрии. Коэффициент эластоподъема и для аппланационного, и для импрессионного тонометра всегда возрастает с жесткостью склерального объема и падает с давлением. Влияние роговичной жесткости не столь однозначно: как правило, для аппланационного тонометра ее возрастание приводит к уменьшению коэффициента эластоподъема, а для импрессионного - к его увеличению. Такой результат имеет место при соотношении между двумя жесткостями, не слишком отклоняющемся от нормы. Если же отклонение велико, зависимость коэффициента эластоподъема от роговичной жесткости может качественно меняться. Однако тот факт, что все кратко сформулированные выше результаты получены в рамках физически линейной модели, оставляет некоторое сомнение в их окончательности. Дело в том, что и отдельно роговица [12; 15, 20], и глазное яблоко в целом [16; 19; 29; 30; 34] демонстрируют нелинейно упругое поведение, и эта нелинейность была учтена (хотя и не слишком корректно) при разработке и последующих уточнениях применяющейся в настоящее время стандартной процедуры обработки данных тонометрии по Шиотцу [22-24; 31]. Анализ влияния реальной нелинейности оказался возможен после разработки нами эффективной модели нелинейно упругой роговицы [11]. Эта модель была тестирована, и область изменения ее параметров была оценена на основании сопоставления с экспериментальными данными [10; 11]. Также показано, что нелинейностью упругого поведения склеральной области при характерных для клинических измерений давлениях можно пренебречь. Поправки, вносимые учетом нелинейности упругого поведения роговицы в расчет внутриглазного давления по показаниям тонометра Маклакова, оценены в работе [6]. В целом они оказались не слишком значительны, как и вклад нелинейности в коэффициент эластоподъема, определяемый при эластотонометрии. Этого следовало ожидать, поскольку при тонометрии по Маклакову давление в процессе нагружения тонометром меняется в весьма узком диапазоне. В статьях [7; 8] мы рассмотрели влияние нелинейности упругих свойств роговицы на тонометрию по Шиотцу. Для этого тонометра ситуация иная, чем для тонометра Маклакова: тонометрическое давление может превосходить истинное в два раза и более. Тем не менее влияние параметра, характеризующего нелинейность упругого поведения роговицы, на расчет истинного давления по тонометрическому оказалось умеренным и качественно влияет на его оценку лишь при достаточно больших давлениях. Однако при заметно 1 1 C В 1 h 1 А С 3 2 а б Рис. 1. К постановке задачи: роговица до внешнего нагружения (а) и нагруженная тонометром Шиотца (б); 1 - свободная роговица, 2 - вогнутый штамп (стопа), 3 - стержень (плунжер); стрелки и буквы A и B указывают точки отрыва роговицы от стопы, а C - от плунжера. На схеме показана конфигурация изучаемого объекта справа от оси симметрии (оси глаза); Rc - радиус не нагруженной извне роговицы, Rb - радиус круга, на который она опирается; - заглубление плунжера отклоняющихся от нормы характеристиках механических свойств и состояния глазного яблока (упругие модули, истинное давление) зависимость между истинным и тонометрическим давлениями продемонстрировала некоторые аномалии в сравнении с линейным случаем. Влияние таких аномалий, в частности, может сказаться на результатах дифференциальной тонометрии. Последнее важно, в частности, потому, что в рамках линейной модели нам до сих пор не удавалось получить относительно высокие значения коэффициента эластоподъема (2 мм рт. ст./г и выше), которые были измерены в клинике у достаточно большого количества пациентов. Такое несоответствие теории и эксперимента затрудняло обоснование использования этого коэффициента в качестве диагностического параметра. В настоящей работе при теоретическом анализе дифференциальной тонометрии по Шиотцу, помимо учета нелинейной упругости роговицы, расширен в сравнении с нашими предыдущими исследованиями диапазон упругих параметров и истинных давлений. Это сделано для того, чтобы включить в рассмотрение область тех аномалий, о которых говорилось выше. Постановка задачи Схема нагружения в импрессионном тонометре Шиотца представлена на рис. 1. Роговица глаза лежащего пациента нагружается вогнутым штампом (стопой). Через ее центральное отверстие свободно опускается стержень (плунжер) с вогнутой торцевой поверхностью, вдавливающейся в роговицу. Непосредственно измеряемая величина - заглубление плунжера под поверхность стопы. Как показали оценки, сделанные нами ранее, в задачах тонометрии роговица демонстрирует физическую нелинейность упругого поведения при малых деформациях, а нелинейностью упругого отклика склеральной области на внутреннее давление можно пренебречь. В соответствии с этим используется система уравнений равновесия, описывающая роговицу как двумерную безмоментную поверхность в предположении геометрической малости деформаций при конечных смещениях, которая рассматривалась ранее (см., например, соотношения (1.1) в [4]). Для роговицы принимается эффективное нелинейное определяющее соотношение [11], включающее три константы, а для склеральной области линейная связь приращений ее объема с приращениями внутреннего давления, характеризуемая одной константой - склеральной жесткостью В качестве констант, характеризующих упругое поведение роговицы, модель рассматривает три упругих модуля: модуль растяжения (как поверхности) и коэффициент Пуассона для физически малых деформаций в окрестности некоторого однородного состояния, именуемого отсчетным, а также параметр нелинейности ответственный за степень физической нелинейности. Однородное состояние было бы возможно для сферической однородной поверхности, нагруженной только внутренним давлением. Однако можно с практической точностью принять, что не нагруженная извне роговица (рис. 1, а), закрепленная по границе, имеет при некотором давлении форму сферического сегмента радиуса , и тогда ее напряженное состояние можно рассматривать как однородное. При другом давлении ее сферичность утратится, деформация уже не будет пространственно однородной, и тензор жесткости перестанет быть однородным и изотропным. Однако при нагружении только внутренним давлением можно пренебречь неоднородностью деформаций и изменением формы, что позволяет выбирать отсчетное состояние произвольно, считая сегмент сферическим с тем же радиусом и полагая, что упругие константы для этих состояний связаны такими же соотношениями, как для однородной деформации [11], т.е. и не изменяются, а модуль растяжения преобразуется по формуле . (1) Таким образом, давление, при котором форма роговичного сегмента считается сферической, может практически выбираться произвольно с пересчетом по формуле (1). Поскольку растяжимость распертой давлением роговицы сильно зависит от этого давления, сравнивать роговицы по их упругим свойствам целесообразно при одном и том же отсчетном давлении, которое далее всегда принимается равным 12 мм рт. ст. Такое давление находится в области физиологически нормальных значений. Вместо модуля растяжения роговицы как поверхности далее рассматривается условный модуль ее материала где мм - характерная толщина роговицы. Значение модуля всегда указывается для отсчетного давления, и его аргумент опускается. Влияние жесткостей и исследуется в широком диапазоне. В соответствии с оценками, выполненными в [4], под средними (нормальными) значениями этих величин понимаются МПа и МПа. Коэффициент Пуассона принимается равным 0,45, что соответствует слабой сжимаемости роговицы как трехмерного объекта. В силу симметрии решение задачи ищется в области, лежащей по одну сторону от оси симметрии ( ). Граничные условия при нагружении тонометром Шиотца - те же, что обоснованы нами ранее [5], с учетом представленной в [6] несколько иной формы для асимптотического условия в апексе, соответствующего допущению об ограниченности усилий. Принималось также условие постоянства полного объема жидкости под оболочкой глаза за время тонометрии, т.е. сумма приращений объема под роговицей и склерального объема равна нулю. Полагаем геометрические параметры роговицы равными их нормальным для человеческого глаза значениям: мм, мм. Радиусы кривизны двух компонентов тонометра Шиотца - стопы и плунжера - принимаем равными 15 мм, а вес стопы - 11 г. Это стандартные характеристики прибора. Численное решение осуществлялось методом стрельбы с использованием схемы Рунге - Кутта четвертого порядка. Результаты численных расчетов Выполненное нами численное исследование зависимости тонометрического давления от веса груза при линейно упругой роговице показало, что эта зависимость близка к линейной. Таким образом, можно ввести параметр (коэффициент эластоподъема) по формуле (2) где и - тонометрические давления при весах груза (плунжера) и соответственно. Этот параметр естественно рассматривать как физически адекватную характеристику дифференциального нагружения роговицы, практически не зависящую от выбора значений и . Расчеты с учетом нелинейности упругого поведения роговицы показали, что и в этом случае линии остаются в практически важной области значений параметров почти прямыми. Это видно из рис. 2, на котором сопоставлены кривые для случаев практически линейной ( ) и нелинейной ( ) роговицы. Значение было ранее оценено нами как характерное из сравнения расчетов с экспериментальными данными. Оценка тонометрического давления по измеряемой геометрической величине - заглублению плунжера - представляет собой отдельную задачу, которая рассмотрена нами для нелинейно упругой роговицы в [7]. Далее при численных расчетах для определения применялась линейная регрессия кривой на участке от 4 до 8 г. Были выполнены расчеты коэффициента эластоподъема в широком диапазоне значений параметров , и , характеризующих механические свойства и механическое состояние глаза. На рис. 3 показано, каким образом нелинейность роговицы влияет на зависимость коэффициента эластоподъема от склеральной жесткости при нормальном давлении и средней жесткости роговицы. Видно, что это влияние практически отсутствует при значениях склеральной жесткости, не более чем вдвое превышающих среднее значение этой величины МПа. Независимо от степени нелинейности роговицы, с возрастанием склеральной жесткости коэффициент растет. При МПа этот коэффициент увеличивается при увеличении параметра нелинейности, и тем значительнее, чем выше В случае нелинейной роговицы при достаточно больших достигаются высокие значения которые не могут быть получены для тех же значений и ни при каких значениях для линейно упругой роговицы. Из рис. 4 видно, что, как и для случая линейно упругой роговицы [9], увеличение не только склеральной но и роговичной жесткости ведет к росту значений коэффициента эластоподъема γ. Этот коэффициент существенно зависит также от истинного давления , как правило, уменьшаясь с его ростом. Скорость этого уменьшения нарастает с ростом . Кривые , представленные на рис. 4, б, становятся все более пологими при уменьшении и при достаточно больших роговичной и склеральной жесткостях коэффициент γ в области низких и нормальных значений практически Рис. 2. Зависимость тонометрического давления от веса плунжера при (штриховые линии) и 60 (сплошные линии) для мм рт. ст., МПа и различных значениях склеральной жесткости , 20, и 60 МПа (чем выше линия, тем больше значение ) Рис. 3. Зависимость коэффициента эластоподъема от склеральной жесткости для различных значений коэффициента нелинейности , 30, 60, и 90 (кривые 1, 2, 3 и 4 соответственно); мм рт. ст.; МПа а б Рис. 4. Влияние обеих жесткостей и истинного давления на величину эластоподъема: а - три семейства кривых различаются склеральной жесткостью: снизу вверх , 20 и 40 МПа; в каждом семействе снизу вверх падает истинное давление: , 16 и 10 мм рт. ст.; ζ = 60; б - три семейства кривых различаются отсчетной роговичной жесткостью: снизу вверх (3), 0,6 (2) и 0,9 (1) МПа; в каждом семействе нижняя (сплошная) кривая соответствует ζ = 60, верхняя (штриховая) - ζ = 90; МПа перестает от зависеть. Этот эффект становится более выраженным с возрастанием параметра нелинейности. При большом и не слишком высокой (например, нормальной) жесткости роговицы возможно даже существование области истинных давлений , в которой с ростом коэффициент γ может возрастать. Рис. 5. Сравнение зависимостей коэффициента эластоподъема γ от истинного давления для очень жесткого и нормального глаза. Верхнее семейство кривых: МПа, МПа; нижнее семейство: МПа, МПа; ζ = 90 (кривые 1), 60 (2), 30 (3), 1 (4) На рис. 5 представлена зависимость коэффициента эластоподъема от истинного давления для двух глаз, сильно различающихся жесткостями обоих компонентов при разных параметрах нелинейности. Один глаз можно назвать нормальным: его роговичная жесткость МПа, а склеральная - МПа, что соответствует средним значениям этих величин. Другой - очень жесткий: обе жесткости в несколько раз превышают норму: МПа, МПа. Видно, что независимо от степени нелинейности роговицы и истинного давления коэффициент γ для жесткого глаза оказывается заметно больше, чем для нормального, достигая при малых давлениях величины 2,2 мм рт. ст./г. Для глаза с нормальной жесткостью компонентов значения коэффициента эластоподъема, хотя и существенно зависят от давления, не могут превышать 1,4 мм рт. ст./г. Если для нормального глаза γ уменьшается с ростом давления более или менее равномерно, то для жесткого при малых давлениях возникает область его медленного падения - практически плато. В случае нелинейной роговицы область плато оказывается заметно шире, чем в случае линейной, захватывая умеренные давления - до 20 мм рт. ст. и выше. Подчеркнем, что для достижения высоких значений коэффициента эластоподъема - порядка 2 мм рт. ст./г и выше - необходима высокая жесткость обоих компонентов глазного яблока. Это видно из рис. 4, б: уменьшение жесткости роговицы ведет к значительному уменьшению γ даже при очень большой склеральной жесткости. Такой эффект понятен: хотя в общей поверхности корнеосклеральной оболочки доля роговицы многократно меньше, чем склеры, но при очень жесткой склере растяжимость хоть и необширной, но существенно более мягкой роговицы становится определяющей для компенсации роста давления. Практически эта ситуация вряд ли встречается часто. Во всяком случае, с возрастом и при глаукоме возрастают жесткости как склеры, так и роговицы [3; 14; 18; 26-28]. Обсуждение Анализ, проведенный на основе эффективной модели, учитывающей нелинейно упругое поведение роговицы, уточнил и качественно подтвердил основные результаты, полученные нами ранее с использованием для роговицы физически линейной модели: измеряемый при дифференциальной тонометрии коэффициент эластоподъема зависит не только от жесткостей компонентов глаза, но (не менее значительно) и от истинного давления. Исследование в широком диапазоне практически возможных значений параметров (включая теперь и параметр нелинейности) показало, что коэффициент эластоподъема растет с увеличением как склеральной, так и роговичной жесткостей и падает с ростом истинного давления. Учет нелинейности приводит к некоторому увеличению расчетных значений этого коэффициента. Влияние параметра нелинейности сильнее сказывается как при больших жесткостях (и роговичной, и склеральной), так и при высоких давлениях. В области значений, где нелинейность роговицы существенна, увеличение этого параметра ведет к росту коэффициента эластоподъема. В расчетах, представленных в настоящей работе, удалось получить при определенных значениях параметров большие коэффициенты эластоподъема (до 2,2 мм рт. ст./г), чего ранее не удавалось достичь в рамках физически линейной модели. Такие значения достаточно часто встречаются при клинических измерениях. Для иллюстрации приведем диаграмму, в которой сведены данные обследования группы пациентов с подозрением на первичную открытоугольную глаукому. Два измерения, давшие значения коэффициента эластоподъема порядка 3 мм рт. ст./г и выше, скорее всего, были не вполне качественно выполнены, что не редкость при клиническом обследовании. Однако можно предположить, что и для этих измерений коэффициент был большим, а обнаружения этого факта достаточно для диагностики (рис. 6). Тот факт, что ранее нам не удавалось получать в расчетах высокие значения коэффициента эластоподъема, не связан только с неучетом физической нелинейности роговицы. В самом деле, рис. 5 демонстрирует достижение им весьма высокого значения мм рт. ст./г при параметре нелинейности , соответствующем практически линейному случаю. Однако такой коэффициент эластоподъема получается лишь при сочетании высоких жесткостей обоих компонентов глаза с достаточно низкими давлениями, а соответствующий набор параметров не попадал в их ранее исследованный диапазон. Нелинейность роговицы тем не менее сказывается существенно, но Рис. 6. Представленные в виде диаграммы данные дифференциальной тонометрии для группы пациентов: по оси абсцисс - порядковый номер обследуемого, по оси ординат - измеренный коэффициент эластоподъема (в мм рт. ст./г). Клиническое исследование проводилось в НИИ глазных болезней им. Гельмгольца с использованием тонографа GlauTest 60 при весах плунжера и г. Данные обрабатывались по формуле (2) несколько в ином отношении: в медленности убывания с возрастанием истинного давления. В широкой области давлений кривая даже практически образует плато. При линейной же роговице с ростом давления коэффициент падает довольно быстро. Таким образом, нелинейность способствует большей вероятности обнаружения высоких значений независимо от давления. Из проведенных расчетов видно, что высокие значения коэффициента эластоподъема наблюдаются лишь при очень высоких (в разы превышающих норму) значениях обеих жесткостей. Возможность теоретического объяснения больших значений дает дополнительный довод в пользу эффективности разрабатываемого нами подхода. Кроме того, сформулированное выше утверждение позволяет до известной степени реабилитировать процедуру дифференциальной тонометрии как метода оценки жесткости глаза. Действительно, первоначальное представление о том, что результат дифференциальной тонометрии сам по себе (без дополнительной информации) позволяет судить о жесткости (ригидности) глаза ошибочно, поскольку определяемый показатель (коэффициент эластоподъема) сильно зависит от внутриглазного давления. Однако при высоких значениях этого коэффициента можно с уверенностью сказать, что независимо от истинного давления жесткость обоих компонентов глаза ненормально высока. А именно обнаружение очень жестких глаз, в которых наиболее вероятно развитие глаукомы, и составляет основную задачу дифференциальной тонометрии. Таким образом, теоретически обоснована эффективность этой процедуры как метода диагностики глаукомы, что было подтверждено ранее клиническими исследованиями, в которых мы принимали участие. Можно утверждать, что при коэффициенте эластоподъема порядка 2 мм рт. ст./г и выше жесткость обоих компонентов корнеосклеральной оболочки - роговицы и склеры - заведомо аномально велика. С другой стороны, коэффициент эластоподъема порядка 1,5 мм рт. ст./г и ниже говорит о значениях жесткостей роговицы и склеры, близких к их средним значениям. Таким образом, по данным только дифференциальной тонометрии можно выделить две группы глаз, в одной из которых корнеосклеральная оболочка очень жесткая, а в другой - жесткость тканей глаза можно считать нормальной. В промежуточной области значений коэффициента эластоподъема с ростом возрастает вероятность (от нулевой до стопроцентной) повышенной жесткости корнеосклеральной оболочки. Реальная же жесткость может быть оценена только при наличии данных об истинном давлении. В полном объеме задача оценки параметров, характеризующих упругие свойства глаза, может быть решена лишь при комплексном определении всех трех основных механических характеристик глаза: двух жесткостей (роговичной и склеральной) и истинного внутриглазного давления. Такая оценка не может быть надежно выполнена с помощью лишь одного тонометра, хотя бы и с использованием разных грузов [9].About the authors
I. N Moiseeva
Institute of Mechanics, Moscow State Lomonosov University
A. A Stein
Institute of Mechanics, Moscow State Lomonosov University
References
Statistics
Views
Abstract - 213
PDF (Russian) - 251
Refbacks
- There are currently no refbacks.