MATHEMATICAL MODELLING OF EYE CORNEA DEFORMATION BY EXTERNAL PRESSURE
- Authors: Moiseeva I.N1, Stein A.A1
- Affiliations:
- Issue: Vol 23, No 4 (2019)
- Pages: 511-525
- Section: Articles
- URL: https://ered.pstu.ru/index.php/rjb/article/view/2334
- DOI: https://doi.org/10.15593/RZhBiomeh/2019.4.02
- Cite item
Abstract
The deformation of the eyeball under the effect of external pressure applied to the cornea in a limited area, for example, as a result of loading with an air jet, is considered. The cornea is modeled by a homogeneous, isotropic, momentless surface and, in the general case, is considered to be linearly Voigt-type viscoelastic. Its elastic properties are determined by the quadratic dependence of the strain energy on the components of the strain tensor, and viscous by the quadratic dissipative function of the strain rates. The scleral region is described by the differential relationship between the scleral volume and intraocular pressure.The process of slow deformation is considered in detail, when the derivatives of the deformation characteristics with respect to time can be neglected in the constitutive relations. It is shown that, at a certain external pressure, a section is formed in the central region of the zone of its influence, which can be considered to be flat with high accuracy. From the external pressure at which this zone occurs, the true intraocular pressure,i.e., pressure in an unloaded eye, can be determined with good accuracy. The relationship between these pressures is weakly affected by the elastic properties of both the cornea and the scleral region. In the area of external pressures corresponding to “flattening”, the rate of subsidence of the corneal apex at external pressure uniformly increasing changes. Possible generalizations of them ode land the problem formulation are discussed.
Full Text
Введение В настоящее время разработано и применяется в клинике немало устройств, основанных на воздействии на роговицу глаза струей воздуха и предназначенных как для измерения внутриглазного давления, так и для более детального исследования механических свойств глазного яблока. Непосредственно измеряются различные деформационные характеристики роговицы в зависимости от задаваемых параметров струи. В частности, может определяться смещение центральной точки роговицы (апекса) под действием приложенного давления. Во многих приборах оцениваются давление в струе или длительность ее воздействия, при которых достигается уплощение некоторой области роговицы, окружающей апекс. При более тонких исследованиях (например, в популярном в отечественных клиниках приборе ORA [11]) оценивается разность давлений, соответствующих уплощению при быстрых нагружении и разгрузке (корнеальный гистерезис). Как и при других принятых в клинике непрямых офтальмологических измерениях, возникают две проблемы. Основной для медицины является возможность использования таких измерений для диагностики патологий, прежде всего глаукомы, и объективной оценки результатов лечения. Эта проблема решается в первую очередь накоплением статистических данных. Однако большой индивидуальный разброс механических свойств глазных тканей [12] и геометрических характеристик глазного яблока делает актуальным получение надежных результатов при обследовании конкретного пациента, что невозможно без решения другой задачи - отчетливой оценки механического смысла получаемых результатов. Для решения этой задачи оказывается необходимым математическое моделирование происходящих при выполнении измерения процессов. В зависимости от характера нагружения и непосредственно измеряемых характеристик при исследовании механики струйного нагружения глаза может оказаться достаточным моделирование процесса либо в рамках чисто упругой модели, либо ее дополнение учетом вязкоупругого поведения глазных тканей. Последнее существенно, если измерения (как в ORA) предполагают получение результатов в быстрой развертке во времени. Авторами ранее был разработан [1] максимально упрощенный подход к механике глазного яблока, который характеризуется небольшим количеством подлежащих определению параметров и оказался полезен при моделировании различных типов тонометрии. В рамках этого подхода роговица представляется безмоментной (мягкой) изотропной однородной двумерной поверхностью, а склеральная область нульмерным элементом, откликающимся изменением объема на изменение внутриглазного давления. В базовом варианте зависимость деформации роговицы от напряжений и объема склеральной области от изменений давления считалась линейно-упругой. В дальнейшем были рассмотрены обобщения базовой модели на случаи пространственно неоднородной и нелинейно упругой роговицы [2, 3]. В предлагаемой работе авторы ограничиваются линейными определяющими соотношениями для роговицы и склеральной области, допуская в общем случае учет их вязкоупругого поведения фойгтовского типа. Однако многие вопросы могут исследоваться на основе чисто упругой модели. Механические, в частности упругие, свойства глазных тканей не только сильно индивидуальны - им свойственны нелинейность, анизотропия и пространственная неоднородность. Особенно обширный материал накоплен по детальному исследованию механических свойств роговицы [4, 6, 8, 14]. В литературе нагружение глазного яблока струей воздуха обычно моделируется, как и для других способов его нагружения, на основе громоздких конечно-элементных моделей, содержащих много констант [5, 10]. Эти константы определяются, как правило, в экспериментах на трупном материале с осреднением результатов по многим глазам. При моделировании быстрых нестационарных процессов, имеющих место при струйном воздействии на роговицу, зачастую используются данные, полученные в опытах с существенно бо́льшими характерными временами. Все эти обстоятельства делают затруднительным эффективное решение обратной задачи непосредственно для обследуемого глаза. В последние годы получило распространение использование на базе конечно-элементных моделей методов оптимизации с целью определения неизвестных параметров. Этот подход применялся и к струйному нагружению глаза [9]. Однако на таком пути плохо просматривается физический смысл используемых соотношений, причем оказывается, что выбор недоступных при прямых измерениях параметров может осуществляться в рамках моделей, некорректно учитывающих основные особенности механического поведения объекта. Исследование на основе элементарной модели с небольшим набором определяющих параметров позволяет понять механизмы наблюдаемых явлений, выяснить физический смысл измеряемых величин и наметить пути минимально необходимых обобщений для описания более тонких, не описываемых этой моделью, эффектов. В этом преимущество такого подхода в сравнении с расчетными методиками, требующими задания множества характеристик, не могущих быть надежно оцененными, тем более для конкретного обследуемого глаза. Постановка задачи: уравнения равновесия роговицы и геометрические соотношения На рис. 1 схематически представлены глазное яблоко, не нагруженное извне (а), а также нагруженное внешним давлением, приложенным в ограниченной центральной области роговицы (б). Поскольку в общем случае предполагается рассмотрение достаточно быстрых процессов, необходимо оценить, возможно ли в этом случае использовать статические соотношения, т.е., во-первых, пренебрегать инерционными членами в динамических уравнениях для глазных тканей и, во-вторых, не рассматривать нестационарные процессы в жидкости. Первое заведомо возможно, если выполнено условие >> а второе - если >> где и - характерное время процесса и характерный размер объекта соответственно; - наименьший из характерных упругих модулей; - характерная плотность тканей (практически равная плотности воды) и - скорость звука в воде. Принимая см, м/с, МПа [1], получаем, что нестационарность процессов в жидкости становится существенной для глазного яблока при временах порядка микросекунд, а нестационарность для тканей - при временах, меньших миллисекунды. Поскольку рассматриваемые процессы имеют характерное время не менее десятков миллисекунд, оба предположения выполняются. Допущение о стационарности процессов перераспределения жидкости заведомо выполняется с очень высокой точностью, тогда как для тканей оболочки нестационарные эффекты могли бы вносить поправку лишь при точных расчетах, каковые для практического исследования рассматриваемых процессов не имеют смысла. С другой стороны, на таких временах, как показывают эксперименты, становятся существенны эффекты немгновенности деформирования (фойгтовской вязкоупругости) [11]. Малость вызываемых внешним воздействием а б Рис. 1. К постановке задачи: а - глазное яблоко до нагружения; б - роговица нагружена внешним давлением. Плоскость AB отграничивает роговицу от склеральной области. Стрелки (б) показывают область действия внешнего давления. Остальные обозначения в тексте при тонометрии смещений жидкости позволяет также пренебрегать эффектами вязкости маловязкой жидкости (по существу, воды) в передней камере под роговицей и полагать, что давление устанавливается в ней мгновенно. Поэтому будем пользоваться для роговицы статическими уравнениями, которые с сохранением предположений, сформулированных в [1], приобретают в случае присутствия распределенного внешнего давления, например, создаваемого струей, следующий вид: (1) (2) (3) (4) Здесь - расстояние точки поверхности, представляющей роговицу, от оси; - угол между нормалью к этой поверхности и осью симметрии; - длина дуги образующей роговицы в начальном состоянии, отмеряемая от точки пересечения роговицы с осью симметрии (апекса) и играющая роль независимой переменной; - текущее внутриглазное давление; и - касательные усилия в меридиональном и перпендикулярном к нему (окружном) направлениях; и - деформации в направлениях главных осей, отсчитываемые от начального состояния ненагруженной извне роговицы. Последнее считается условно соответствующим некоторому сферическому пространственно однородному сегменту радиусом (радиус роговицы) с внутренним давлением , которое в дальнейшем будет отождествляться с внутриглазным давлением до нагружения. До начала деформирования () усилия и равны между собой и определяются формулой Функция радиальной координаты и времени определяет приложенное к роговице внешнее давление. В случае струйного воздействия эта функция отлична от нуля при в ограниченной области , где - радиус области, в которой струя воздуха действует на роговицу. Тангенциальной составляющей воздействия струи на роговицу будем пренебрегать в сравнении с нормальной силой, обеспечиваемой внешним давлением. Определяющие соотношения и краевые условия Рассмотрим теперь линейные определяющие соотношения для усилий и с учетом вязкоупругости фойгтовского типа. Остаемся в рамках базового для нашего подхода представления роговицы двумерной безмоментной, изотропной и однородной поверхностью. В рамках такого подхода поверхность не сопротивляется изгибу и существенны лишь деформации и напряжения в плоскости, касательной к этой поверхности, которые характеризуются тензорами второго ранга в двумерном пространстве. В задачах, для которых в процессе деформации сохраняется условие осевой симметрии, присутствуют лишь компоненты этих тензоров, направленные по главным осям, соответствующим в каждой точке меридиональному и перпендикулярному к нему направлениям. Введенные выше деформации и оказываются главными компонентами двумерного тензора деформаций, а усилия и - главными компонентами двумерного тензора напряжений. Введем энергию деформации и диссипативную функцию зависящую от скоростей деформации и В силу малости скоростей деформаций правомерно отождествлять их с производными по времени от и соответственно, безразлично частными или субстанциональными. Отсутствие сопротивления изгибу проявляется в том, что из числа аргументов этих функций исключены параметры, отвечающие за изменение кривизны поверхности. Тогда связь напряжений с деформациями определяется общими формулами (5) С учетом изотропии материала роговицы в тангенциальных направлениях будем задавать энергию деформации и диссипативную функцию в виде квадратичных функций, симметричных по своим аргументам: (6) где коэффициенты , , , - постоянные. Зависимость энергии деформации и диссипативной функции от главных значений соответствующих тензоров обеспечивает инвариантность этих зависимостей. Определяющие соотношения получаются в виде (7) Положительная определенность свободной энергии и диссипативной функции приводит к неравенствам (8) Соотношения (7) приводятся к виду (9) где , , , (10) Очевидно, , Если дополнительно потребовать, что (т.е. , что не следует из термодинамических неравенств), то . Знак не определен, даже если дополнительно потребовать (по аналогии с ). В силу (8) и (10) выполняются неравенства (11) Коэффициенты и имеют смысл упругих модулей (модуля растяжения и коэффициента Пуассона) роговицы как двумерной упругой поверхности при медленных процессах. Будем заменять склеральную область с окружающими ее тканями вязкоупругим элементом, поведение которого описывается зависимостью (12) где - приращение склерального внутриглазного объема (части внутреннего объема, ограниченной склерой и плоскостью, отграничивающей склеру от роговицы) вследствие изменения давления, а точкой обозначена производная по времени. Строго говоря, обобщенное соотношение (12) учитывает не только вязкоупругое поведение склеры и окружающих ее тканей, но и вязкие свойства заполняющей ее среды (стекловидного тела), а при детальном анализе давление в этой среде может оказаться не только отличным от внутриглазного давления (т.е. давления в передней камере), но и пространственно неоднородным. Полное приращение объема жидкости во внутриглазном пространстве определяется как сумма и приращения «подроговичного» объема (части объема, лежащей выше плоскости ). Приращение вычисляется исходя из функций и , рассчитанных при решении краевой задачи для роговицы. В рассматриваемой задаче полное приращение объема равно нулю вследствие возможности пренебречь притоком и оттоком внутриглазной жидкости за время тонометрии: (13) В силу симметрии решение краевой задачи для роговицы ищется только в одной полуплоскости (). Условие закрепления роговицы в плоскости AB внешней окружности, по которой она контактирует со склерой, и условия в точке ее пересечения с осью - коллинеарности нормали оси и ограниченности меридиональной деформации - совпадают с использованными в [1]: (14) Последнее из условий (14) является более общим, чем приведенное в [1] условие для напряжений, которое может быть получено из (14) в рамках использованной в [1] линейно-упругой модели. Упругие константы далее переобозначены по формулам ,, где мм - средняя в норме толщина роговицы, а мм-3 [1]. Будем называть и роговичной и склеральной жесткостями соответственно. Геометрические характеристики роговицы (см. рис. 1, а) принимаются равными мм, мм, что соответствует стандартным для человеческого глаза значениям. Моделирование медленного нагружения внешним давлением Характерное время установления деформаций равно , поскольку второе присутствующее в соотношениях (9) характерное время в силу неравенств (11). Надежные оценки для этого времени отсутствуют. Наличие эффектов вязкоупругости при струйном нагружении [11] позволяет предположить, что для соответствующих процессов оно имеет порядок десятков миллисекунд. Более точно параметрможет быть оценен при сопоставлении результатов моделирования с данными экспериментов. В настоящей работе далее будут рассматриваться относительно медленные процессы, для которых характерное время >>. В этом случае можно ограничиться упругой моделью. В соответствии с этим в соотношениях (9) будем полагать . Таким образом, изучаемая в этом разделе модель роговицы, по существу, совпадает с линейно-упругой моделью [1] с различием лишь в уравнении (2) (отличное от нуля давление ). Если считать, что за характерное время задачи успевает установиться также склеральный объем (в соотношении (12) << и можно положить ), отсутствие жесткого нагружающего устройства с заменой его отличным от нуля внешним давлением остается единственным отличием от реализовывавшейся ранее постановки. Далее рассматривается тонометрическая задача для такого нагружения. Изучается изменение геометрии роговицы. В центре внимания оценка внешнего давления, приводящего к уплощению примыкающей к апексу области роговицы (размер этой области заранее не известен и на практике, как правило, не отслеживается). Все линейные размеры далее обезразмерены по невозмущенному радиусу роговицы , в соответствии с чем независимая переменная меняется в диапазоне от 0 до Распределение давления вблизи поверхности роговицы определяется конструкцией устройства и эволюцией струи в процессе ее распространения. По экспериментальным данным для ORA, представленным в [5], это распределение имеет общий вид, показанный на рис. 2. Будем базироваться на распределении такого вида, используя для его моделирования функцию где - максимальное значение действующего на роговицу давления, которое далее для определенности будем называть давлением в струе. В дальнейшем давление будет варьироваться, а полуширина зоны влияния струи (т.е. области где положительно) полагается равной 0,385, если не оговорено противное. Эта величина соответствует размерному значению, приведенному в [5]. Малость деформаций и узость области непосредственного струйного воздействия позволяют по мере надобности предполагать совпадение зависимостей и в области действия струи, а также отождествлять расчетные функции от , характеризующие форму pj -0,4 -0,2 0 0,2 0,4 r Рис. 2. Принятая в расчетах зависимость приложенного извне давления от безразмерного расстояния до центральной оси . При давление достигает максимального значения роговицы в этой области, с соответствующими функциями от . Упругие характеристики глазного яблока будут варьироваться, но, если особо не оговорено, полагаются равными характеристикам «стандартного» глаза: МПа, МПа [1]. Коэффициент Пуассона во всех случаях считается равным На рис. 3 представлена зависимость угла , характеризующего наклон моделирующей роговицу поверхности к центральной оси для пяти различных давлений в струе при одних и тех же значениях упругих констант и истинного внутриглазного давления. Видно, что этот угол (а следовательно, и форма роговицы) при рассмотренных давлениях заметно меняется лишь в области непосредственного воздействия струи . При возрастании формируется область, в которой роговица становится вогнутой (т.е. область отрицательных ). Радиус этой области зависит от давления в струе и на этапе ее формирования, т.е. при давлениях, не слишком отличающихся от того давления, при котором появляются значения , и имеет порядок, не превышающий . На рис. 4 прослежен процесс вдавливания роговицы в центральной области в характерном для этого процесса масштабе. На рис. 4, а продемонстрировано изменение в зависимости от давления в струе, а на рис. 4, б показана эволюция собственно формы. Из рис. 4, б хорошо видно, что при некотором давлении в самом деле имеется практически плоский участок, который при расчетных значениях параметров приблизительно соответствует давлению 14 мм рт. ст. и имеет полуширину . На графике (см. рис. 4, а) видно, что углы в этой области близки к нулю (при выполняется неравенство ). Строго говоря, области, где у роговицы появляется плоский участок, не существует, но визуально он присутствует. В практически используемых устройствах область «уплощения» определяется по пику интенсивности светового сигнала, отраженного от роговицы. Ширина светового пучка определяется конструкцией 1 0,8 0,6 0,4 0,2 0 -0,2 0 0,2 0,4 0,6 0,8 1 s Рис. 3. Зависимость угла от максимального значения давления в струе для «стандартного» глаза: мм рт.ст.; МПа; МПа. Значения давления в струе мм рт.ст. (сверху вниз) 0,4 0,3 0,2 0,1 0 -0,1 0,05 0,04 0,03 0,02 0,01 0 -0,01 0 0,1 0,2 0,3 0,4 s 0 0,1 0,2 0,3 0,4 а б Рис. 4. Распределение угла отклонения нормали к роговице от оси (а) и форма роговицы в области воздействия струи (б) при давлениях мм рт. ст. (сверху вниз). На рис. б по осям расстояния указаны в безразмерных единицах прибора и обычно подбирается эмпирически. Расчеты, проведенные нами в широком диапазоне физиологически допустимых значений параметров, показали, что в рамках используемой модели целесообразно исследовать процесс уплощения в области с безразмерным радиусом Если выбрать эту область существенно уже, зона «уплощения», сформировавшись в этой области, продолжает при увеличении давления в струе расширяться до значительно большего размера.Тогда диапазон давлений, соответствующих «уплощению», окажется слишком широким. При существенно же большем радиусе выбранной области она будет захватывать участок, где роговица при любых внешних давлениях остается заметно выпуклой, что существенно исказит общую картину. В дальнейшем всегда полагаем, что О точном значении давления, при котором достигается «уплощение», говорить нельзя, а можно лишь о диапазоне давлений, в котором на участке выбранного радиуса отклонение нормали от оси достаточно мало. Это допущение соответствует неравенству (15) где выбирается достаточно малым, т.е. таким, чтобы диапазон давлений струи, в котором это неравенство выполняется, был бы достаточно узким и приемлемым в соответствии с требуемой точностью измерений. Можно также воспользоваться критерием (16) При соответствующем выборе констант и оба критерия эквивалентны. Заметим, что эти константы не могут быть выбраны сколь угодно малыми, поскольку точное уплощение не достигается. В дальнейшем будем выбирать в качестве давления уплощения такое давление , при котором достигается минимум функции . Этот подход удобен тем, что нет необходимости задавать критерий малости или, выбор которого представляет самостоятельную (требующую расчетов) задачу, поскольку такой критерий не может быть задан произвольно малым. Из рис. 5, а непосредственно видно, что при выборе значение, при котором достигается минимум функции , определяется достаточно отчетливо: Для такого значения угол на отрезке удовлетворяет неравенству (16) с . При этом давление уплощения оценивается достаточно грубо: зависимости при изменении в диапазоне шириной порядка 1 мм рт. ст. различаются незначительно (рис. 5, б). Названное обстоятельство практически важно: оно позволяет не требовать от процедуры измерений чрезмерной точности. В принципе можно не оговаривать радиус области, для которой исследуется уплощение, заранее, а искать при каждом давлении в струе область, в которой выполняется критерий (15) или (16), затем находить давление уплощения как такое давление , при котором эта область максимально широка. Однако использованный нами метод проще и ввиду продемонстрированной грубости определения не приведет к существенному изменению результата. Заметим, что, хотя ширина области уплощения и может быть определена в рамках используемой модели, в реальности (как и при тонометрии по Гольдману) эта величина сильно зависит от распределения упругих свойств по толщине и поэтому вряд ли может играть роль достаточно эффективного критерия для оценки внутриглазного давления. На рис. 6 показано при тех же значениях параметров, что и в предыдущих расчетах, изменение с ростом давления в струе расстояния по оси от центральной точки 12,8 13,2 13,6 14 14,4 14,8 pj 0 0,04 0,08 0,12 0,16 0,2 s а б Рис. 5. Расчетный график функции (а) и распределение в области при (сверху вниз) (б) 0,08 0,06 0,04 0,02 0 -0,02 0,6 0,5 0,4 0,3 0,2 0,1 2,04 2 1,96 1,92 1,88 1,84 10 12 14 16 18 pj, мм рт. ст. 1,96 1,95 1,94 1,93 1,92 1,91 12,8 13,2 13,6 14 14,4 14,8 pj, мм рт. ст. а б Рис. 6. Зависимость (в разных масштабах) высоты апекса над границей роговицы со склерой от давления в струе: а - в широком диапазоне давлений; б - в области давлений, близких к давлению уплощения (апекса) роговицы до ее границы со склерой, т.е. проседание апекса под действием струи. При равномерном увеличении действующего на роговицу давления апекс проседает практически равномерно, но в области давлений, где достигается уплощение, наблюдается сначала ускорение, а затем замедление проседания апекса. По выходе из этого интервала давлений равномерное проседание возобновляется. Этот эффект в принципе может быть использован для обнаружения давления уплощения. Далее будем исследовать связь давления уплощения с параметрами, характеризующими механические свойства и механическое состояние глазного яблока. В рамках применяемой модели таких параметров три: истинное внутриглазное давление , т.е. давление в ненагруженном глазу, а также упругие характеристики роговицы и склеральной области и Основная задача любой тонометрии - определение истинного внутриглазного давления . На рис. 7 представлена зависимость давления уплощения от . Зависимость эта близка к линейной. Давление уплощения всегда больше , и разность между этими величинами возрастает с давлением: мм рт. ст. при мм рт. ст. и 4 мм рт. ст. при мм рт. ст. Если оценивать внутриглазное давление по давлению в струе, полагая разность между ними средней между этими значениями величиной мм рт. ст., ошибка не превысит 1 мм рт. ст., что для практических целей приемлемо. Более корректно учитывать в соответствии с зависимостью, представленной на рис. 7, реальные значения , соответствующие разным . 35 30 25 20 15 10 12 16 20 24 28 32 p0, мм рт. ст. Рис. 7. Зависимость давления уплощения от истинного давления; МПа, МПа 1 2 3 4 5 12,8 13,2 13,6 14 14,4 14,8 pj, мм рт. ст. 1,2 1 0,8 0,6 0,4 0,2 0 5 4 3 2 1 0,45 0,4 0,35 0,3 0,25 0,2 0,15 12,8 13,2 13,6 14 14,4 14,8 pj, мм рт. ст. а б Рис. 8. Функция примм и различных значениях упругих констант: а - МПа, = 0,3; 0,4; 0,5; 0,6; 0,7 МПа (1-5); б - = 0,3 МПа, 10, 20, 30, 40, 60 МПа (1-5) Оценим теперь влияние на давление уплощения упругих характеристик глазного яблока. На рис. 8 представлены зависимости при значительных превышениях нормальных значений: на рис. 8, а роговичной жесткости и на рис. 8, б склеральной жесткости . Минимумы кривых соответствуют давлению уплощения. С увеличением роговичной жесткости давление уплощения, соответствующее тому же истинному давлению, уменьшается, а с возрастанием склеральной жесткости растет, т.е. ужесточения обеих структур действуют, как и при тонометрии по Маклакову, разнонаправленно. В обоих случаях различие со стандартным глазом невелико и не превышает в рассмотренном диапазоне значений величины порядка 0,5 мм рт. ст. Таким образом, при медленном нагружении внешним давлением внутриглазное давление с практической точностью может быть определено по внешнему воздействию, приводящему к определенной специфической деформации (давлению уплощения). Индивидуальные механические свойства роговицы и склеральной области, в отличие от статической тонометрии, влияют на связь между истинным давлением и давлением уплощения незначительно. Заключение В рамках разработанного авторами общего подхода, основанного на представлении роговицы двумерной безмоментной поверхностью, а склеральной области нульмерным элементом, откликающимся изменением объема на изменения внутриглазного давления, разработана модель, учитывающая вязкоупругие (фойгтовского типа) свойства глазных тканей. Поставлена задача о деформировании глазного яблока под действием приложенного к роговице извне в ограниченной области давления, в частности, вследствие воздействия струей воздуха. Подробно рассмотрены относительно грубые эффекты, проявляющиеся при медленном нагружении, когда возможно пренебречь вязкостью роговицы и склеральной области. Модель воспроизводит ряд важнейших особенностей струйного воздействия на роговицу глаза. Показано, что при определенном давлении в струе в центральной области зоны ее воздействия формируется участок, который с высокой точностью можно рассматривать как плоский. Именно формирование зоны уплощения отслеживают многие используемые в медицинской практике приборы. По давлению в струе, при котором возникает эта зона, с хорошей точностью можно определять истинное внутриглазное давление, т.е. давление в ненагруженном глазу. Практическую важность имеет слабое влияние на соотношение, связывающее эти давления, упругих свойств как роговицы, так и склеральной области. В области внешних давлений, соответствующих «уплощению», меняется скорость проседания апекса роговицы при равномерном увеличении внешнего давления: эта скорость, приблизительно постоянная вне этой области, сначала растет, затем падает, снова выходя на постоянное значение за ее пределами. Представленная выше модель, учитывающая вязкоупругость, позволяет получать более тонкие оценки и моделировать поведение глаза в высоком временном разрешении при его быстрых нагружении и разгрузке. Именно быстрые процессы, для которых эффект вязкоупругости существен, реализуются в большинстве современных приборов, осуществляющих струйное воздействие на роговицу. Легко показать, что модель в этом случае описывает различие давлений, соответствующих уплощению в этих процессах (т.н. корнеальный гистерезис). Соответствующие детальные расчеты будут представлены в последующих публикациях. Дальнейшая модификация модели позволит учесть неоднородность толщины роговицы, которая, как показывают многочисленные исследования [5, 7, 10, 13], имеет существенное значение для оценки реальных механических характеристик глаза, в первую очередь внутриглазного давления, по данным тонометрии. Первым этапом такого исследования может служить разработанный авторами ранее подход [2], при котором неоднородность толщины и механических свойств роговицы учитывается через интегральную характеристику - эффективную жесткость. Однако тот факт, что воздействие струи в основном сосредоточено в достаточно узкой области (ширина которой сравнима с толщиной роговицы), скорее всего, потребует разработки менее грубой модели, учитывающей неоднородность толщины отдельно от неоднородности жесткости.About the authors
I. N Moiseeva
A. A Stein
References
- Моисеева И.Н., Штейн А.А. Анализ зависимости давление-объем для глазного яблока, нагруженного плоским штампом, на основе двухсегментной упругой модели // Известия РАН. Механика жидкости и газа. - 2011. - № 5. - С. 3-15.
- Моисеева И.Н., Штейн А.А. Влияние пpоcтpанcтвенной неодноpодноcти pоговицы на дефоpмационные cвойcтва глазного яблока и pезультаты аппланационной тонометpии по Маклакову // Биофизика. - 2017. - Т. 62, № 6. - С. 1193-1203.
- Штейн А.А., Моисеева И.Н., Любимов Г.А. Математическая модель роговицы глаза с учетом экспоненциальной нелинейности ее упругих свойств при условии геометрической малости деформаций // Российский журнал биомеханики. - 2019. - Т. 23, № 3. - С. 375-390.
- Corneal biomechanics and refractive surgery. Ed. by F.A. Guarnieri. - Berlin: Springer, 2015. - 146 p.
- Elsheikh A., Alhasso D., Kotecha A., Garway-Heath D.F. Assessment of the ocular response analyzer as a tool for intraocular pressure measurement // J. Biomech. Eng. - 2009. - Vol. 131. - ArticleID 081010. - P. 1-9.
- Fernandez D.C., Niazy A.M., Kurtz R.M., Djotyan G.P., Juhasz T. Finite element analysis applied to cornea reshaping // J. Biomed. Opt. - 2005. - Vol. 10, № 6. - 064018. - P. 1-11.
- Goebels S.C., Seitz B., Langenbucher A. Precision of ocular response analyzer // Current Eye Res. - 2012. - Vol. 37, № 8. - P. 689-693.
- Hjortdal J.Ø. Regional elastic performance of the human cornea // J. Biomech. - 1996. - Vol. 29, № 7. - P. 931-942.
- Jannesari M., Kadkhodaei M., Mosaddegh P., Kasprzak H., Behrouz M.J. Assessment of corneal and fatty tissues biomechanical response in dynamic tonometry tests by using inverse models // Acta Bioeng. Biomech. - 2018. - Vol. 20, № 1.- P. 39-48.
- Joda A.A., Shervin M.M.S., Kook D., Elsheikh A. Development and validation of a correction equation for Corvis tonometry // Comp. Meth. Biomech. Biomed. Eng. - 2016. - Vol. 19, № 9. - P. 943-953.
- Luce D.A. Determining in vivo biomechanical properties of the cornea with an ocular response analyzer // J. Cataract Refract. Surg. - 2005. - Vol. 31, № 1. - P. 156-162.
- Pallikaris I.G., Kymionis G.D., Ginis H.S., Kountis G.A., Tsilimbaris M.K.Ocular rigidity in living human eyes // Invest. Ophthalm. Vis. Sci. - 2005. - Vol. 46, № 2. - P. 409-414.
- Ramm L., Herber R., Spoerl E., Raiskup F., Pillunat L.E., Terai N. Intraocular pressure measurement using Ocular response analyzer, dynamic contour tonometer, and scheimpflug analyzer Corvis ST // J. Ophthalm. - 2019. - Vol. 2019. - ArticleID 3879651. - 9p.
- Shin T.J., Vito R.P., Johnson L.W., McCarey B.E. The distribution of strain in the human cornea // J. Biomech. - 1997. - Vol. 30, № 5. - P. 497-503.
Statistics
Views
Abstract - 104
PDF (Russian) - 39
Refbacks
- There are currently no refbacks.