DOI: 10.15593/RZhBiomeh/2017.4.08 УДК 531/534:[57+61]
МОДЕЛИРОВАНИЕ КОНТАКТНОГО ВЗАИМОДЕЙСТВИЯ ЗАХВАТЫВАЮЩЕГО ИНСТРУМЕНТА С БИОЛОГИЧЕСКОЙ ТКАНЬЮ
А.А. Яковенко1' 2
1 Лаборатория трибологии, Институт проблем механики им. А.Ю. Ишлинского Российской академии наук, Россия, 119526, Москва, проспект Вернадского, 101, корп. 1, e-mail: anastasiya.yakovenko@phystech.edu
2 Московский физико-технический институт, 141701, Московская область, Долгопрудный, Институтский переулок, 9
Аннотация. Захватывающие инструменты используют для временного сдавливания биологических тканей. Чтобы оценить значения контактных напряжений, возникающих в биологической ткани, и оптимизировать прилагаемые нагрузки, в данной работе была решена контактная задача о взаимодействии зажимного устройства с биологической тканью. Рассмотрен дискретный контакт с биологической тканью. Изучено контактное взаимодействие зажимного инструмента, имеющего на контактирующей поверхности несколько выступов вытянутой формы, с биологической тканью, для описания механических свойств которой используется модель линейно упругого материала. Разработана модель для исследования распределения нагрузок между выступами медицинского инструмента, используемого для захвата мягких биологических тканей. Модель также позволяет произвести анализ зависимости силы, приложенной к зажиму, от его внедрения при разных расположениях выступов. На основании проведенного анализа установлено, что нагрузка на выступы распределена неравномерно, степень неравномерности зависит от расстояния между выступами. Рассмотрены случаи расположения выступов, при которых нагрузка на некоторые выступы стремится к нулевым значениям. Установлено, что разница между нагрузками на выступы возрастает с уменьшением отношения длины выступов к их ширине. Предложен алгоритм расчета высотного распределения выступов, которое обеспечивает их равномерную нагруженность. Используемая в работе форма выступа в виде узкого прямоугольного в плане штампа хоть и является упрощенной моделью, но позволяет на основании полученных аналитических зависимостей дать рекомендации относительно расположения выступов и их геометрических размеров с целью обеспечения их равной нагруженности, а также уменьшения общей прикладываемой силы для достижения нужной величины внедрения (захвата ткани).
Ключевые слова: захватывающий инструмент, биологическая ткань, система выступов, упругий материал, распределение нагрузок, равнонагруженность, высотное распределение.
Введение
Захватывающие инструменты, или по-другому зажимы, - это медицинские инструменты, которые используют для временного сдавливания биологических тканей с целью остановки кровотечения, перекрытия просвета полых органов или для фиксации тканей и органов. Основное требование для этих инструментов - прочное удерживание тканей при лапароскопических операциях. Для этого они имеют рабочие кремальеру и губки, у которых, в зависимости от области применения, различаются
© Яковенко А.А., 2017
Яковенко Анастасия Александровна, инженер, Москва
профили и размеры. По силе действия на ткани хирургические зажимы бывают эластичные, не вызывающие разрушения органов и тканей, жесткие, вызывающие обратимые или необратимые изменения ткани, и раздавливающие, вызывающие разрушение органов и тканей. Воздействие инструмента на биологическую ткань будет зависеть от свойств ткани, приложенных нагрузок, а также механических характеристик и формы рельефа поверхности инструмента. Для того чтобы оценить значения контактных напряжений, возникающих в биологической ткани, и оптимизировать прилагаемые нагрузки, необходимо решить контактную задачу о взаимодействии зажимного устройства с биологической тканью. Подробный класс задач контактного взаимодействия с биологической тканью рассмотрен в работе [7], где представлены методы решения задач Герца о сжатии упругих тел, внедрении штампа в многослойные основания, а также контактные задачи для тел с включениями. Также изложен подход, применяемый при решении контактных задач для линейных вязкоупругих тел. Данные методы позволяют изучать изменение напряженного состояния биологической ткани при взаимодействии с различными инструментами. Кроме того, в данной книге изложены подходы к решению динамических задач, например, задача о движении жесткого тела по границе упругого полупространства с включением. В работе [1] разработана математическая модель контактного взаимодействия головки сенсора с мягкой биологической тканью, которая позволяет исследовать напряженно-деформированное состояние, в том числе и контактные характеристики. В работе [10] решена задача о внедрении с постоянной скоростью индентора в двухслойное основание с учетом релаксационных характеристик слоев. Также в данной работе получены аналитические зависимости нагрузки и радиуса контакта от времени для различных значений скорости вдавливания индентора и выражение для распределения давлений под индентором в разные моменты времени. Дискретный характер взаимодействия тел рассматривается в работах [2, 4-6, 9]. В них изучается контактное взаимодействие системы N связанных между собой сферических штампов с упругим полупространством. На основании решения контактной задачи изучено распределение усилий между штампами при заданном внедрении системы штампов в упругое тело. В данной работе объединены две эти задачи и рассмотрен дискретный контакт с биологической тканью. Изучено контактное взаимодействие зажимного инструмента, имеющего на контактирующей поверхности несколько выступов вытянутой формы, с биологической тканью, для описания механических свойств которой используется модель линейно-упругого материала [8].
Постановка задачи
Рассмотрим контактную задачу о внедрении жесткого штампа, имеющего на поверхности систему N выступов, в упругое полупространство с заданными характеристиками (Е - модуль Юнга, V - коэффициент Пуассона). Выступы моделируются вытянутыми вдоль оси Ох штампами, где 2аг - длина г -го выступа, 25г. - его ширина (рис. 1). Все выступы зажимного инструмента имеют одинаковую
высоту, плоское основание и прямоугольную в плане форму. Расстояние между выступами и размеры площадок контакта будем считать заданными (1^ - расстояние
между г -м и ] -м выступами). Инструмент вдавливается без перекоса в упругое основание заданной силой Р.
Введем систему координат Охух. Ось Ох совпадает с осью, направленной вертикально вниз через центр фиксированного г -го штампа, а плоскость Охуг с
недеформированной поверхностью полупространства. Оси Ох и Оу направлены вдоль длинной и короткой граней г -го штампа соответственно.
Запишем условие контакта для г -го выступа. На основании принципа суперпозиции смещение границы упругого полупространства под данным выступом можно представить как сумму вертикальных смещений и^(1)(х, у) и и^(2)(х, у), где
<(2)( х, у) - смещение, обусловленное распределенным по рассматриваемому пятну
контакта давлением рг(х,у), а и'г(1)(х,у) - смещение, обусловленное влиянием
остальных выступов. Таким образом, условие контакта для произвольного г -го выступа будет иметь вид
б = и;(1)(х, у) + и;(2)( х, у), (1)
где Б - внедрение инструмента в основание, которое не зависит от координат вследствие одинаковой высоты выступов. Смещение границы выступа, вызванное давлением р.(х,у), распределенным по другим выступам ( j = 1,2,...,N,. Ф г),
вычисляется по формуле
<(1)( x, у) =
1 -v
2 N
ХЯ-
Pj (x', у ')dx' dy у
(2)
и
ПЕ у!(Ху + х'- х)2 + ( + у'- у)2'
где (х., У.) - координаты центральной точки области контакта .-го выступа, интеграл берется по области контакта .-го выступа. В данном случае X. = 0, У. = I. .
Чтобы определить зависимость упругих смещений от приложенных давлений для единичного выступа, воспользуемся решением Галина [3] для узкого прямоугольного в плане штампа длиной а и шириной 8, согласно которому выражение для контактных давлений можно представить в виде
Р( х)
P( X У) =
2 2 2 У
где р(х) - распределение вдоль оси Ох силы, действующей по нормали к границе полупространства у = 0.
Рис. 1. Схема внедрения инструмента
При этом упругие смещения определяются соотношением
и (х, 0) = ЕЮ, к
где коэффициент постели вычисляется по формуле
пЕ 1
к =
2 (1 -v2 )lg (а/ 5)
Отсюда смещение в среднем сечении /-го выступа, вызванное приложенным к нему давлением, вычисляется по формуле
2(1 -V2) (а/О .) О, 2(1 -V2) (а/О .)
К (х) = Л-'Е I Р'(х,У)Ф = --'Е Р(х), (3)
пЕ • пЕ
( . 5,. . -5 >
/г 4.7 Рг (х) / 4у pi (х) . о,. . -О,.
где I р1 (х, у)ау = —— I , = = —- агс81п—1 - агс81п—- = рг (х).
-о, п -О, VО,.2 - у2 п I О/ О, )
Условие равновесия в данном случае будет иметь вид
N
Р = 1Р, (4)
,=1
где р =|| р, (х, у)ахау - нагрузка, действующая на i-й выступ (интеграл берется по области контакта ,'-го выступа).
Метод определения распределения нагрузок между выступами
Определим дополнительные смещения под произвольным ,-м выступом от действия давлений р. (х, у), (, Ф.), приложенных к другим выступам. В силу
предположения об узости выступа (О. << а.) распределение давления р.(х, у), (, Ф .)
можно заменить распределением погонной силы р(х) на отрезке (-а., а.).
Рассмотрим задачу, в которой действие соседних выступов заменяется силами, распределенными по линии, проходящей через центр выступа, т.е. погонными силами р(х). Тогда из соотношений (2) и (3) получаем
и'(2) =2(1 >'8^О)р.(х), и«" = (5)
г ПЕ гл ь г _Е ^ J г— "2 — V >
ПЕ ПЕ 1 =1,}*, -а^(х '- х) + I.
Подставив (5) в условие контакта (1), будем иметь
2(1 -V2)'в(а7О) 1 -V2 " } р}(х*)ах•
° =-ПЕ-р,(х)I I
nE nE j=1.j-, -L (x'- x)2 +1 2
\l(xx )2
Отсюда получаем
D = 2 (1 -v2) lg (aj 5,) p + Х f f PJ (x') dx' dx
' ^j=1, j* -f, ^(x x)2 + j
a. %E
(6)
Данная система вместе с условием равновесия (4) позволяет найти внедрение и распределение нагрузок, если известны взаимное расположение выступов, вид функции распределения давления под каждым штампом и величина силы, приложенной к инструменту.
Будем считать, что все выступы имеют одинаковую форму и длину (2а). Тогда формула (6) приобретает вид
D = 2(1 -v2)lg(a/5)p + * yrJjixWdxx_
' anE j=1, j-a f^(x '- x)2 + l y2
a%E
(7)
P
Рассмотрим равномерное распределение давления, т.е. pi (х) = —- . Тогда
2а
} } Pj (x ')dx' dx Pj a a J J 2a ^
-a -a W (x - x) + l j -a -a
dx' dx
V(x'- x )2 +1 j2
=P
ln
2 + 4a2 + 2a
+ 4a2 - 2a a
aa -1
И окончательно формула (7) приобретает вид
D =
2(1 -v2 )lg (a/5)p 1
-v
2 N
a%E
Х Pj
anE , j
ln
- 1 -11
Jl„2 + 4a2 -2a aVVJ J>
Y j
(8)
Поскольку при выводе (8) использован ряд упрощающих предположений (в частности, предположение, что длина 2а сохраняется для каждого выступа с учетом их взаимного влияния), то при исследовании системы уравнений (8) необходимо принимать во внимание ограничение на минимальное расстояние между выступами .
Для оценки минимального расстояния /тт достаточно рассмотреть влияние
ближайшего штампа к фиксированному и оценить, как погрешность замены зависит от расстояния между двумя рассматриваемыми штампами. Тогда запишем условия контакта в первом и во втором случаях:
d = 2 (1 -v2) lg (a/ 5) р + ff
a%E
p,.-1( x ')dx' dx anE -a J^(x '- x)2 + l2
D=2(1-У2)lg(a/52p+хХ ((-1+д)
a%E
a%E
j=1, j * i
2 ln#+g+2a -1 -1)
Vl2 + 4a2 -2a a\ >
Различаются только вторые слагаемые, поэтому погрешность запишется как
Д =
2lg(a/ 5,)
P
\ Ф2 + 4a2 + 2a 1
ln
Vl2 + 4a2 - 2a a
. P-
(Vl2 +4a2 -l) -f p,-1(x')ln^
) -a v(
2lg(a/ 5)
^ln^l2 + 4a2 + 2a 1/ ПГПТл Л
-1 (л/l2 + 4a2 -1 a
a - x ')2 +12 + a - x' a + x ')2 +12 - a - x'
)
dx'
<
Здесь для двойного интеграла использовалась теорема о среднем, а также то, что
^(а - х ')2 +/2 + а - х'
минимальное значение функция 1п
yj(a + х ')2 +12 - a - x'
достигает на границе области
интегрирования.
Перейдем к безразмерной длине / = I/2а . Тогда
1
P-1 21g(a/5,) А
1n-
1 +1
l
(Л
- 2 NI2 +1 -1
Найдем /тп, оценив отношение как
P
< 0,1. Например, для ^ (а/5 ,) = 1, /тп « 0,5.
Заметим, что минимальное расстояние зависит от количества выступов N. Погрешности в определении силы на произвольном фиксированном выступе от замены давлений на других выступах суммируются и убывают с увеличением расстояния.
// ~ Б « Р (1 -V2)
Введем для удобства безразмерные величины: /. = —, Б = —, Р =
P =
P (1 -V2)
2a
2a
4Ea2
и
. На расстояние между штампами накладывается условие / . > 8г + 8..
4Еа2
Тогда система уравнений (8) приобретает вид
' N ^ Ц
D = 1 п
_ 1 N _
1g(a/S)P + - I P
2 j=1, j&
1n
+1 +1
VF
- 2 (Л
+1 -1
+1 -1
(9)
J J
или в матричном виде АР = Б, где А - симметрическая матрица размера N х N,
' ■ Л
. (а/ 8,) _ + 2 элементы которой а// =-8,,. +—
п
п
1n
¿1+1+1 - 2 (#77 -1) /2 +1 -1 ^ '
#
(S.. - символ
Кронекера), а Р, Б - векторы размерности N. Данная система позволяет исследовать внедрение и распределение нагрузок при варьировании ширины выступов и их взаимного расположения.
Анализ распределения давления между выступами
при их равномерном расположении
Пусть все выступы располагаются на равном расстоянии / друг от друга (I/ = /|/ — / ). Тогда условие контакта (9) примет вид
D = 4
п
1 N
1g(a/S)р +1 I P
2 А
j=1, j *i
1n
ß2 (i - j )2 +1 +1 V'2 (i - j )2+1 -1
- 2
(2 (i - j)2 +1 -1 (i - j
(10)
VJ
Проанализируем зависимость нагрузки от расстояния при стремлении последнего к бесконечности / ^да . Воспользуемся следующей оценкой, справедливой при больших / :
ln
+1 -2()2 +,-/.)2)|п[Ы-2(-.-цг-л) = 0.
{12(-]) +1 -1 ^ ^ ^ '' '|г-. и 1 1 "
Тогда из (10) следует, что величина внедрения стремится к некоторому
41§( а 8г)
значению, не зависящему от расстояния, т.е. Б ^
. Таким образом, при
больших расстояниях между выступами распределение нагрузок становится равномерным.
Просуммировав формулу (9), получаем
( ( РТ,-^— ^
4 , . , _ 1 N А ~ .// (- 7) +1 +1
DD = ■
nN
1 N N
lg(«A2P+2х Х P
2 ,=1 j=1, j*,
ln
>/'2 (- j 22
--2
+1 -1
(2 (- j )2+1 -1 (- j)
которая отображает зависимость внедрения от прилагаемой к инструменту нагрузки.
Для анализа распределения нагрузок на выступы возьмем определенное количество выступов, а также будем считать, что выступы имеют одинаковую ширину (28). Ниже
приведены графики распределения нагрузок между выступами (рис. 2, а, б) и зависимости отношения нагрузок (рис. 3, а, б), действующих на средний и крайний выступы от расстояния / между выступами для количества выступов, равного трем (рис. 2, а) и пяти (рис. 2, б).
0.4 0:3 ОД ОД о
I 2а = 1
/ 2а= 1,75
0.7 0 6 0,5 0,4 0J
од од о
im2=I.7d
Г2а = 2
1:2а--
б
¡:2а = 2.5
Рис. 2. Диаграммы распределения нагрузок между выступами: а - N = 3 ; б - N = 5 для разных расстояний между ними (а/8 = 10)
а
Результаты расчетов показывают, что при малых расстояниях между выступами нагрузки на крайние выступы больше, чем на средний, а при больших расстояниях, как и ожидалось, эти величины стремятся к одинаковому значению. Также из полученных зависимостей следует, что при уменьшении соотношения длины и ширины выступов увеличивается разница между нагрузками на средний и крайние выступы.
Ръ!Р\
1 2 3 4 > 7 3 9 ю Ц
б
Рис. 3. График зависимости N = 3 отношения нагрузки на средний выступ к нагрузке на крайний (Р2/Р1) (а), (Р3/Р1) (б) и от расстояния между выступами (/ ) для четырех значений отношения длины выступа к его ширине (а/8). Желтая линия - а/ 8 = 100, красная линия - а/ 8 = 30, синяя линия - а/ 8 = 10, зеленая
линия - а 8 = 7
Определение высотного распределения равнонагруженной системы штампов
Рассмотрим теперь контактную задачу о внедрении жесткого штампа с системой N выступов на поверхности, имеющих разную высоту, но одинаковую форму. Тогда условие контакта (9) будет иметь вид
' 1 N ' + , _ ( ^ ^
, = ±
П
lg (/s)p +1 Z P
2 ^ 1
21=1,1 *i
ln
#
^ - 2 U +1 -1 v
l2 +1 -1
(11)
JJ
где 3i = к - 30, к = к/2а - высотное распределение выступов, 30 = 30/2а -сближение тел под действием нагрузки Р.
Найдем такое высотное распределение, чтобы нагрузки на все выступы были одинаковы, т.е. р = Р/N. Также пусть все выступы располагаются на равном расстоянии I друг от друга (!р = ¡^ - ). Из формулы (11) высота каждого выступа определяется как
h =
4 P
lg (/«) +1 Z
' 1=11 *i
Jl2 ((- j )2+1+1
ln\ V ;--2
Ф2 (i -1 )2 +1 -1
(2 (i -1 )2 +1 - l(i -1))
+ D0. (12)
J J
Из (12) следует, что высотное распределение равнонагруженных выступов существенно зависит от расстояния между ними; с увеличением расстояния второе слагаемое в скобках выражения (12) стремится к нулю, и высоты всех выступов выравниваются.
Из (12) также следует, что внедрение каждого выступа Di линейно зависит от
нагрузки, прилагаемой к штампу. На рис. 4 представлен график зависимости соотношения высот от расстояния между штампами для случая пяти выступов.
Рис. 4. Графики зависимости (N = 5), высотного распределения ((к - 30Р)
от расстояния между выступами (I), а!5 = 10: желтая линия - первый и пятый, красная линия - второй и четвертый, синяя линия - третий)
На малых расстояниях при / ^ 0 логарифмические члены стремятся к члену
же расстоя
4Р (а/ 8)
1п С~, где С - некая константа. На больших же расстояниях подтверждается /2
равномерное распределение штампов по высоте, Di ^
nN
Выводы
Разработана модель для исследования распределения усилий между выступами зажима, используемого для захвата мягких биологических тканей, а также для анализа зависимости приложенной к зажиму силы от его внедрения при разных расположениях выступов. На основании проведенного анализа установлено, что нагрузка на выступы распределена неравномерно, и степень неравномерности зависит от расстояния между выступами. Также могут возникать ситуации, когда нагрузка на некоторые выступы будет стремиться к нулевым значениям. Разница между нагрузками на выступы возрастает с уменьшением отношения длины выступа к ширине. Показано, что при выбранной форме прямоугольных в плане выступов внедрение возрастает при уменьшении расстояния между выступами. Предложен алгоритм расчета распределения выступов по высоте, обеспечивающего их равную нагруженность.
Используемая в работе форма выступа в виде узкого прямоугольного в плане штампа является упрощенной моделью выступов зажима. Однако она позволяет на основании полученных аналитических зависимостей дать рекомендации относительно расположения выступов и их геометрических размеров с целью обеспечения их равной нагруженности, а также уменьшения общей прикладываемой силы для создания нужной величины внедрения (захвата ткани).
Благодарности
Автор выражает благодарность своему научному руководителю, доктору физико-математических наук, профессору Московского государственного университета имени М.В. Ломоносова и Московского физико-технического института (государственного университета), академику Российской академии наук, Горячевой Ирине Георгиевне, за ценные советы и замечания при постановке и решении задачи, за необходимые исправления и внесенные правки, а также за рекомендации по оформлению статьи.
Работа выполнена при финансовой поддержке РФФИ (грант 16-58-52033).
Список литературы
1. Антонов Ф.К., Горячева И.Г., Мартыненко Ю.Г., Рязанцева М.Ю., Су Фончин. Анализ контактного взаимодействия искусственного видеотактильного механорецептора с мягкими биологическими тканями // Российский журнал биомеханики. - 2013. - Т. 17, № 1. - С. 8-21.
2. Аргатов И.И., Дмитриев Н.Н. Основы теории упругого дискретного контакта. - СПб.: Политехника, 2003. - 234 с.
3. Галин Л.А. Контактные задачи теории упругости и вязкоупругости. - М.: Наука, 1980. - 304 с.
4. Горячева И.Г. Механика фрикционного взаимодействия. - М.: Наука, 2001. - 478 с.
5. Горячева И.Г. Периодическая контактная задача для упругого полупространства // Прикладная математика и механика. - 1998. - Т. 62, № 6. - С. 136-144.
6. Горячева И.Г. Расчет контактных характеристик с учетом макро- и микрогеометрии взаимодействующих тел // Трение и износ. - 1999. - Т. 20, № 3. - С. 239-248.
7. Горячева И.Г., Садовничий В.А. Применение методов механики контактных взаимодействий при диагностике патологических состояний мягких биологических тканей - М.: Изд-во Моск. ун-та, 2009. - 306 с.
8. Джонсон К. Механика контактного взаимодействия - М.: Мир, 1989. - 510 с.
9. Любичева А.Н. Контактное взаимодействие и изнашивание неоднородных тел: дис. ... канд. физ.-мат. наук. - М., 2008. - 88 с.
10. Yakovenko A., Goryacheva I., Dosaev M. Estimating characteristics of a contact between sensing element of medical robot and soft tissue. - Springer, 2017. - P. 561-569.
SIMULATION OF CONTACT INTERACTION OF A GRIPPING TOOL
WITH A BIOLOGICAL TISSUE
A.A. Yakovenko (Moscow, Russia)
Gripping tools used for the temporary compression of biological tissues. To estimate the values of the contact stresses in biological tissue, and to optimize the applied force, in this paper, we solved the contact problem of interaction of a gripping device with biological tissue. We considered a discrete contact with biological tissue. We studied the contact interaction of the gripping tool with the contact surfaces of several projections of elongated shape, with biological tissue, for describing the mechanical properties using a model of linearly elastic material. We developed a model for the study of the distribution of forces between the projections of the medical tool used to grip soft biological tissues. The model also allows analyzing the dependence of the force applied to the tool from its introduction at different locations of the projections. Based on the conducted analysis we found that the force on the projections is irregularly distributed, and the extent of the variation depends on the distance between the projections. In addition, there may be situations when the force on some of the projections will tend to zero values. We established that the difference between the forces on the projections increases with a decrease in the ratio of the length of the projections to their width. The proposed algorithm calculates the height distribution of the projections, which ensures uniform loading. The used form of a narrow rectangular stamp of a projection is a simplified model but it allows based on obtained analytical dependencies make recommendations relative to the location of the projections and their geometric dimensions to ensure equal loading and reduce the total applied force to provide the desired value of the introduction (capture the tissue).
Key words: gripping instrument, biological tissue, system of the projections, elastic material, distribution of forces, equal load, height distribution.
Получено 05 декабря 2017