CC BY
40
УДК 616-71
В.Н. Учаев, А.В. Кобзев, В.Д. Семенов, Д.О. Пахмурин, А.Ю. Хуторной, А.В. Литвинов
Согласование рабочего инструмента с ультразвуковой колебательной системой для введения иглы в биологическую ткань при реализации локальной гипертермии
Предложена методика расчета ультразвуковой колебательной системы (УЗКС) с рабочим инструментом (РИ) на одну длину волны ультразвуковых колебаний в материале, получена формула зависимости части длины УЗКС от длины РИ. Получена 3D модель УЗКС с учетом размеров РИ, рассчитанных по формуле, предложенной авторами. Проведен модальный анализ рассчитанной конструкции в среде ANSYS .
Ключевые слова: ультразвук, колебания, модальный анализ, пьезоэлемент.
Для реализации способа управляемой локальной гипертермии с применением игольчатых нагревателей, далее РИ [1], необходимо их введение в биологическую ткань под действием ультразвуковых колебаний (УК). Исходя из медицинских требований, амплитуда УК должна лежать в пределах от 3 до 30 мкм, частота - в пределах от 25 до 35 кГц. Применение УК для введения РИ необходимо, чтобы уменьшить усилия ввода при прохождении тканей, кровопотерю и болевые ощущения, а также для ускорения заживления раневых каналов. Эти положительные эффекты связаны с нагревом за счет ультразвука и коагулированием тканей раневого канала, кроме того при таком введении имеет место раздвигание крупных кровеносных сосудов, что, в свою очередь, снижает вероятность возникновения кровотечений [2, 3].
Для того чтобы вводимый РИ имел продольные колебания с заданными параметрами, необходимо согласовать его с ультразвуковой колебательной системой (УЗКС) с учетом того, что длина и диаметр РИ могут меняться в широких пределах. Согласование необходимо чтобы по всей длине УЗКС с РИ помещалась одна длина волны ультразвуковых колебаний либо ее половина, это достигается применением отражающей и излучающей накладок с концентратором, от них и РИ зависят размеры рис.1,
УЗКС РИ
e z d l c a b
IV 0
x —x 5
2
3
4
У4
Х/4 m
^-
Рис. 1. УЗКС с РИ и распределение амплитуды колебания: 1 - отражающая накладка; 2 - пьезоэлемен-ты; 3 - излучающая накладка; 4 - составной концентратор (экспоненциальный и линейный); 5 - РИ
В общем виде уравнение продольных колебаний может быть представлено дифференциальным уравнением в частных производных (1) [4]:
1
s =dU
dt2 dx V дх
где S - площадь поперечного сечения; U - перемещение РИ; с — скорость звука в материале. Рассмотрим общее решение этого уравнения для участка (a, b) (см. рис.1), (2), (3) [4]:
Ui(x,t) = | A■ cos—• x+B-sin — I-sin(ю• t+ф), I c1 c1 J
U2(x,t) =|c■ cos—■ x+D■ sin —j-sin(ю■ t + ф) ,
где A, B, C, D — неизвестные коэффициенты полученные в результате решения уравнения (1).
Граничные условия здесь записываются следующим образам, с учетом того, что Ui -это перемещение по длине a , U 2 - по длине b :
dU1
(1)
(2) (3)
При х = a, U1 = 0; при х = 0 U2 = U1, E2 ■ S2 dU2 = E1 ■ S1■dUi; при x
дх
dx
dx
• = 0 [4].
т т , ^ * Юг
Используя граничные условия и зная что а + Ь = — , о^ =— можно получить выражение зави-
4 с
симости длины участка УЗКС а от размера РИ Ь (4):
arctg
E1 ■ S1 ■ш*
a = -
E2 ■ S2 ■ш ■ tg(b■Ю2),
*
Ш1
(4)
где S2 - это площади соприкосновения участка УЗКС с РИ; Е\, Е2 - модули Юнга металла
УЗКС и РИ; с1, с2 - скорости звука в УЗКС и РИ.
Для длины участка с УЗКС можно воспользоваться уравнением, формула (5) [4]:
* = , (5)
где п - количество полуволн, помещающихся в концентраторе; / - заданная частота работы 30 кГц. Для длины участка I УЗКС можно воспользоваться формулами (6) и (7) [5]:
,=" ■ с
2 ■ /\
1 +
In (К)
V п ■n J
К =
Dl
D2
(6) (7)
где К — коэффициент трансформации; Е>\, Д2 - диаметры большей и меньшей части концентратора.
Используя методику расчета, приведенную для длины Ь РИ и длины а УЗКС, можно определить размеры оставшихся участков: е = 1 см, г = 1,5 см, ё = 0,8 см. Участок е - это толщина отражающей накладки; г - толщина пьезоэлементов; ё — толщина излучающей накладки [4].
Используя формулы (5)-(7), получим размеры участков: с = 4,15 см, I = 4,8 см. Если подставить длину РИ Ь = 4 см в формулу (4), длина линейного участка УЗКС а будет составлять 1 см.
С использованием полученных размеров в среде БоШШэгк^ была создана 3Б-модель УЗКС с РИ, рассчитанная на 30 кГц (рис. 2).
1
Рис. 2. Внешний вид модели УЗКС с РИ
Для определения собственных частот колебаний данной модели УЗКС был проведен модальный анализ в среде АКБУБ [6]. В данной модели материалом для излучающей и отражающей на-
кладок, экспоненциального и прямолинейного участков УЗКС был выбран титан - марки ВТ6. Материал РИ - медицинская сталь 12Х18Н10Т.
На рис. 3, 4 приведены результаты модального анализа:
Туре: Тс^а1 Ое^лггмЬоп Ргечиепсу: 29197 Нг
итЬ т
09,12.2014 10:43
96,13 Мах
85,449 74,768 64,087 53,406 42,724 32,043 21,362 10,681 ОМт
/\Г\
Рис. 3. Изгибные колебания УЗКС на частоте 29197 Гц
Туре: Тс^а| С'сЪитыЪсп Р^иепсу: 3 1775 Нг
|_1т£ т
09.12.2014 10:47
Рис. 4. Продольные колебания на частоте 31775 Гц
Проведенный модальный анализ в среде АКБУБ показывает, что в пределах заданной расчетной частоты 30 кГц УЗКС с РИ может иметь изгибные колебания на частоте 29197 Гц и продольные на частоте 31775 Гц. Необходим режим с продольными колебаниями, соответственно после конст -руирования и согласования УЗКС с генератором возбуждающих импульсов можно будет настроить ультразвуковое устройство на механический резонанс с продольными колебаниями РИ, их амплитуда будет зависеть от прикладываемого напряжения и суммарного коэффициента усиления УЗКС.
Полученная формула зависимости длины участка УЗКС а от длины РИ Ь, учитывающая площади соприкосновения и материалы, совместно с формулами для расчета участков УЗКС могут применяться для расчета различных видов УЗКС и РИ.
Литература
1. Управление электронными игольчатыми нагревателями при реализации метода локальной гипертермии и его экспериментальная проверка / А.В. Кобзев, В.Д. Семенов, Д.О. Пахмурин, А.В. Литвинов, А.Ю. Хуторной, В.Н. Учаев // Доклады ТУСУРа. - 2010. - № 2 (22), ч. 2. - С. 300302.
2. Миллер Э. Применение ультразвука в медицине / Э. Миллер, К. Хилл. - М.: Мир, 1989. -568 с.
3. Акопян В.Б. Основы взаимодействия ультразвука с биологическими объектами / В.Б. Акопян, Ю.А. Ершов. - М.: МГТУ, 2005. - 224 с.
4. Кумабэ Д. Вибрационное резание. - М.: Машиностроение, 1985. - 182 с.
5. Гершгал, Д. А. Ультразвуковая технологическая аппаратура / Д. А. Гершгал, В.М. Фридман. -М.: Энергия, 1976. - 500 с.
6. Басов К. А. АКБУБ в примерах и задачах / К. А. Басов, Д.Г. Красновский. - М.: Компьютер -Пресс, 2002. - 224 с.
Учаев Виктор Николаевич
Инженер ООО «ПромЭл», г. Томск
Тел.: 8-923-423-17-00
Эл. почта: [email protected]
Кобзев Анатолий Васильевич
Д-р техн. наук, профессор, зав. каф. промышленной электроники (ПрЭ) ТУСУРа
Тел.: (382-2) 51-05-30
Эл. почта: [email protected]
Семенов Валерий Дмитриевич
Канд. техн. наук, профессор, зам. зав. каф. ПрЭ по научной работе Тел.: (382-2) 413-918 Эл. почта: [email protected]
Пахмурин Денис Олегович
Канд. техн. наук, доцент, зав. лаб. каф. ПрЭ
Тел.: 8-903-913-46-38
Эл. почта: [email protected]
Хуторной Александр Юрьевич
Аспирант каф. ПрЭ, мл. науч. сотрудник
лаборатории интеллектуально-модуляционных энергетических систем
Тел.: 8-913-883-26-26
Эл. почта: [email protected]
Литвинов Александр Викторович
Инженер ООО «ПромЭл» Тел.: 8-923-414-13-33 Эл. почта: [email protected]
Uchaev V.N., Kobzev A.B., Semenov V.D., Pakhmurin D.O., Khutornoy A.Y. Litvinov A.V.
Harmonization of the working tool with the ultrasonic vibrating system for insertion of the needle into the
biological tissue in the implementation of local hyperthermia
We propose the design procedure of ultrasonic oscillatory system (UOS) with a working tool (WT) at one wavelength of ultrasonic vibrations in the material obtained according to the formula of the length of the UOS from WT. We received a 3D model of UOS with allowance for the WT size, calculated according to the formula proposed by the authors. The modal analysis of the calculated structure designed in the ANSYS. Keywords: ultrasound, vibrations, modal analysis, piezoelectric element.
