Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ МЕХАТРОННОЙ СИСТЕМЫ ОРТОПЕДИЧЕСКОГО АППАРАТА АВТОМАТИЗИРОВАННОГО ОСТЕОСИНТЕЗА

Мамаев И.М. 1 Федотов О.В. 1 Филимонов В.Н. 1
1 ФГБОУ ВПО «Владимирский государственный университет имени Александра Григорьевича и Николая Григорьевича Столетовых»
Испытания аппарата выявили ряд конструктивных недостатков, устранение которых потребовало наряду со структурной модернизацией, направленной на улучшение его качественных показателей – массогабаритных, жесткостных и точностных, решения математической задачи оптимального проектирования. С этой целью разработана математическая модель мехатронной системы аппарата как совокупная модель его функциональных частей. Механическая система аппарата представлена в виде пятимассовой дискретной динамической модели. Зависимость момента шагового двигателя от параметров движения задана динамической характеристикой. Модель микропроцессорной системы управления построена как для дифференцирующего звена с замедлением. Динамический анализ мехатронной системы аппарата позволит обоснованно задать значения фиксированных конструктивных параметров, назначить диапазоны изменения варьируемых параметров и сформулировать задачу оптимального проектирования в окончательном виде.
ортопедический аппарат
роликовинтовая передача
динамическая модель
оптимизация
1. Автономный ортопедический аппарат автоматизированного остеосинтеза для реабилитации инвалидов / В.В. Козырев, О.И. Алборов, О.В. Федотов, С.М. Воробьев, С.Л. Тихомиров, И.М. Мамаев // Медицинская техника. – М., 2011. – № 4 – С. 30–33. – ISSN 0025-8075.
2. Мамаев И.М. Обоснование выбора планетарной роликовинтовой передачи для модернизированного привода ортопедического аппарата / И.М. Мамаев, О.В. Федотов // Информационно-измерительные и управляющие системы военной техники: материалы третьей Всероссийской научно-технической конференции, Владимир, 14–16 ноября 2012 г. – М.: Изд-во РАРАН, 2012. – С. 132–134 – ISBN 978-5-8311-0706-7.
3. Модернизация конструкции ортопедического аппарата автоматизированного остеосинтеза. Физика и радиоэлектроника в медицине и экологии / О.В. Федотов, В.Н. Филимонов, И.М. Мамаев, А.В. Крылов, А.В. Кириллов // Труды 10-й межд. науч.-техн. конф. с элем. научн. молод. школы. Кн 2. – Владимир, 2012. – С. 98–101. – ISBN 978-5-9055-2706-7.
4. Свидетельство о государственной регистрации программы для ЭВМ № 2011610101 Российская Федерация. Управление приводами ортопедического аппарата автоматизированного остеосинтеза / Мишулин Ю.Е., Федотов О.В., Филимонов В.Н., Мамаев И.М.; заявитель и патентообладатель: Общество с ограниченной ответственностью «Малое инновационное предприятие «Трикарол»; поступл. 8.10.10; зарегистр. 11.01.12, Бюл. № 1(П ч.).
5. Система измерений осевых усилий автоматизированного ортопедического аппарата. Физика и радиоэлектроника в медицине и экологии / В.В. Козырев, О.В. Федотов, А.В. Крылов, И.М. Мамаев // Труды 9-й межд. науч.-техн. конф. с элем. научн. молод. школы. – Владимир – Суздаль: Влад. гос. ун-т, 2010. – С. 274–278. – ISBN 978-5-904875-03-9.
6. Кенио Т. Шаговые двигатели и их микропроцессорные системы управления: пер. с англ. – М.: Энергоатомиздат, 1987. – 200 с.

1. Конструкция ортопедического аппарата автоматизированного остеосинтеза

Ортопедический аппарат автоматизированного остеосинтеза (далее – аппарат) с тремя независимыми прецизионными электромеханическими приводами (ЭМП) разработан в ВлГУ [1] взамен компрессионно-дистракционного аппарата с ручным приводом проф. Г.А. Илизарова.

Как и в аппарате с ручным приводом в нем используются кольца со спицами, фиксируемыми на соединяемых отломках костной ткани конечности пациента (рис. 1), но раздвигаются кольца в процессе лечения управляемыми ЭМП линейного перемещения на базе шаговых микроэлектродвигателей (ШД) и планетарных роликовинтовых передач (ПРВП). Шарнирно-стержневая система аппарата (ШСС) совместно с ЭМП обеспечивает циклические продольные и угловые смещения подвижного и неподвижного колец. Корпус основного ЭМП жестко фиксируется на неподвижном кольце с возможностью регулирования его углового положения.

pic_23.tif

Рис. 1. Ортопедический аппарат автоматизированного остеосинтеза: 1 – основной привод; 2 – дополнительные приводы аппарата; 3 – фиксатор углового положения основного привода; 4 – двухстепенной шарнир; 5,6 – неподвижное и подвижное кольца аппарата

Программное управление ЭМП аппарата осуществляется микропроцессорной системой управления (МПСУ). Параметры цикла дистракции костной ткани: период цикла, величины линейных микроперемещений и скоростей движения для каждого независимого ЭМП задаются с пульта управления. За цикл дистракции время работы ШД составляет 5…7 с, а время остановки около 40 мин. Самотормозящиеся ПРВП в составе ЭМП фиксируют относительное положение колец во время остановки ШД, что существенно экономит энергию питающего аккумулятора.

2. Результаты испытаний и модернизация аппарата

Разработанный аппарат реализует широкий круг ортопедических задач. Вместе с тем недостатки, выявленные в процессе испытаний, потребовали решения ряда задач структурной модернизации сложной мехатронной системы аппарата, направленных на улучшение его главных качественных показателей – массогабаритных, жесткостных и точностных. В этой связи была проведена модернизация всех функциональных частей аппарата – ШСС, ЭМП и МПСУ.

Так, усовершенствованная система измерения и регистрации осевых усилий в составе МПСУ позволила получить точную информацию о фактически развиваемых осевых силах ЭМП [5]. Действительная осевая нагрузка 600 Н оказалась вдвое меньше проектной, что позволило отказаться от планетарных редукторов в ЭМП и заменить их меньшими по массе и габаритам сильфонными муфтами. В результате такой замены сократилась кинематическая цепь и, как следствие, повысились точность, чувствительность и надежность ЭМП.

При практической реализации полученных преимуществ возникла необходимость повышения жесткости исполнительных механизмов ЭМП и ШСС аппарата вцелом. Для этого были разработаны ПРВП с максимально возможным числом резьбовых роликов при сохранении заданного передаточного отношения и самоторможения передачи, а также устройство фиксации углового положения основного ЭМП аппарата [2, 3]. Управляющая программа как составная часть МПСУ также была модернизирована [4].

3. Постановка задачи исследования

Технический уровень аппарата во многом определяют ПРВП, выполняющие в составе трех независимых ЭМП функции исполнительных механизмов. Широкие возможности ПРВП по выбору конструктивных параметров делают актуальной задачу выбора их оптимального сочетания, которое дает наилучшее проектное решение, соответствующее экстремальным значениям доминирующих критериев качества (минимальной массе и максимальной жесткости) механической системы аппарата, при оптимальном цикле дистракции костной ткани. В связи с очевидной противоречивостью названных критериев можно назначить обобщенный относительный критерий – массу, отнесенную к жесткости, который следует минимизировать. Это достаточно простой и эффективный критерий для несущей механической системы аппарата [2].

4. Динамическая модель мехатронной системы аппарата

Для окончательной постановки и решения математической задачи оптимального проектирования необходимо построить динамическую модель сложной мехатронной системы модернизированного аппарата, включающую в себя динамические модели его функциональных частей, с целью всестороннего динамического анализа.

4.1. Математическая модель механической системы аппарата

Так как работа трех независимых ЭМП осуществляется поочередно, т.е. в любой момент времени работает и создает осевое усилие только один из них, то механическую систему аппарата представим в виде трех, шарнирно соединенных между собой стержней (рис. 2). Стержень Г-образной формы схематизирует неподвижное кольцо с закрепленными на нем двумя выключенными в данный момент времени ЭМП. Прямолинейный стержень – подвижное кольцо. Функции третьего стержня (переменной длины) выполняет ПРВП, работающего в данный момент времени ЭМП. Центральный винт ПРВП через упругую муфту связан с ротором ШД, к которому приложен вращающий момент, а выходной шток шарнирно соединен с подвижным кольцом.

Все кинематические пары ШСС аппарата считаем идеальными и голономными. Наиболее податливыми звеньями этой механической системы являются упругая муфта, хрящевая ткань, расположенная между подвижным и неподвижным участками кости, а также связи между этими участками и соответствующими кольцами. Учтем также диссипативные силы внутреннего трения в материале упругих звеньев и потери на трение скольжения и качения в сопряжениях звеньев ПРВП.

В указанной постановке динамическая модель ШСС аппарата представляет собой пятимассовую систему с сосредоточенными параметрами, которая состоит из абсолютно твердых тел, соединенных между собой безынерционными упругодиссипативными элементами, подчинена стационарным, голономным связям и находится под действием потенциальных сил упругости и непотенциальных диссипативных сил.

pic_24.tif

Рис. 2. Динамическая модель механической системы аппарата

Дифференциальные уравнения движения механической системы в общем виде можно записать в форме уравнений Лагранжа в обобщенных координатах

Eqn201.wmf

где k – число степеней свободы системы; Т – кинетическая энергия системы; qj – соответствующая обобщенная координата; Eqn202.wmf и Eqn203.wmf – обобщенные силы соответственно потенциальных и диссипативных сил.

За обобщенные координаты механической системы примем: φ1 – угол поворота ротора ШД; φ2 – угол поворота блока винта с роликами; φ3 – угол поворота подвижного кольца; x4 – линейная координата подвижного участка кости относительно неподвижного; x5 – линейная координата (деформационная) неподвижного кольца относительно неподвижного участка кости; x6 – линейная координата (деформационная) блока винта с роликами относительно неподвижного участка кости.

Кинетическая энергия механической системы, найденная как сумма кинетических энергий подвижных масс, является квадратичной формой обобщенных скоростей

Eqn204.wmf

где J1, J2, J3 – моменты инерции, соответственно ротора ШД, блока винта с роликами и подвижного кольца, приведенные к оси винта; m2 – масса блока винта с роликами; m3 – масса подвижного кольца; m4 – масса подвижного участка кости, мышц и кожи; m5 – суммарная масса неподвижного кольца и связанных с ним корпусных деталей ЭМП.

Обобщенные силы для потенциальных сил упругости, соответствующих каждой из обобщенных координат, вычислим по формуле

Eqn205.wmf

запишем потенциальную энергию системы как функцию обобщенных координат

Eqn206.wmf

где см – крутильная жесткость сильфонной муфты; с11 и с12 – линейные жесткости в сопряжениях роликов с опорной и ходовой гайкой соответственно; с21 и с22 – линейные жесткости связи колец соответственно с неподвижным и подвижным участками кости; сх – линейная жесткость хрящевой ткани; as = S/2π; S – перемещение ходовой гайки за оборот винта.

Обобщенные силы диссипативных сил найдем как коэффициенты в сумме возможных работ диссипативных сил, соответствующих положительным приращениям обобщенных координат δqj

Eqn207.wmf

Тогда

Eqn208.wmf

где МД – момент ШД; μм – коэффициент демпфирования в сильфонной муфте.

Eqn209.wmf

где μ12 – коэффициент демпфирования в сопряжениях роликов с ходовой гайкой; η – КПД прямого хода ПРВП; N1 = с12(φ2as – φ3as + x4 – x5) – осевая сила в резьбовом сопряжении опорной гайки с роликами; N2 = с11(x6 – x5) – осевая сила в резьбовом сопряжении ходовой гайки с роликами; fk – приведенный коэффициент трения качения; d3 – средний диаметр резьбы гайки.

Eqn210.wmf

Eqn211.wmf

где μ11 – коэффициент демпфирования в сопряжении роликов с опорной гайкой.

Eqn212.wmf

где μ21 – коэффициент демпфирования связи кольца с неподвижным участком кости.

Eqn213.wmf

где μ22 – коэффициент демпфирования связи кольца с подвижным участком кости; μx – коэффициент демпфирования хрящевой ткани.

После соответствующих преобразований получим математическую модель механической части аппарата в виде системы линейных дифференциальных уравнений с постоянными коэффициентами

4.2. Математическая модель шагового электродвигателя

Зависимость МД от параметров движения зададим динамической характеристикой шагового электродвигателя, определяющей характер его взаимодействия с механической системой аппарата и МПСУ. С учетом дифференциального уравнения вращения ротора ШД [6] выражение для МД можно представить в виде:

Eqn214.wmf

где μД – коэффициент вязкого трения ШД; N – число витков в обмотке статора; Фм – максимальный магнитный поток; р – число пар магнитных полюсов ротора двигателя; λ – угол между обмотками ШД; I1, I2 – токи, поступающие на обмотки статора ШД с МПСУ.

Токи в обмотках статора находятся из дифференциальных уравнений:

Eqn215.wmf

где U – напряжение источника питания электродвигателя; L – собственная индуктивность каждой фазы; R – сопротивление цепи обмотки статора.

4.3. Математическая модель микропроцессорной системы управления

Управление работой ШД происходит путем чередования управляющих импульсов прямоугольной формы, формируемых МПСУ в соответствии с заданной программой. Эти импульсы преобразуются в импульсы напряжения и поступают на обмотки статора двигателя. Полагая МПСУ дифференцирующим звеном с замедлением с постоянной времени Тэл и коэффициентом усиления сигнала Ку, запишем дифференциальные уравнения изменения напряжений на обмотках статора в виде:

Eqn216.wmf

Эти уравнения представляют собой математическую модель МПСУ.

5. Направления дальнейших исследований

Таким образом, на основе построенных динамических моделей функциональных частей модернизированного аппарата сформирована математическая модель мехатронной системы аппарата в целом. Динамический анализ аппарата сводится к нахождению периодических решений системы дифференциальных уравнений математической модели при различных значениях массогабаритных и упругодиссипативных параметров для оценки их влияния на основные показатели качества аппарата. В связи с математическими трудностями, обусловленными нелинейностью и высоким порядком системы, эту задачу целесообразно решать численными методами на ЭВМ.

По результатам динамического анализа можно обоснованно задать значения фиксированных конструктивных параметров и назначить диапазоны изменения варьируемых параметров, что позволит сформировать задачу оптимального проектирования в окончательном виде.

Eqn217.wmf

Рецензенты:

Кульчицкий А.Р., д.т.н., профессор, заместитель главного конструктора по испытаниям ООО «ВМТЗ», г. Владимир;

Гоц А.Н., д.т.н., профессор кафедры тепловых двигателей и энергетических установок ФГБОУ ВПО «Владимирский государственный университет имени Александра Григорьевича и Николая Григорьевича Столетовых», г. Владимир.

Работа поступила в редакцию 04.04.2012.


Библиографическая ссылка

Мамаев И.М., Федотов О.В., Филимонов В.Н. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ МЕХАТРОННОЙ СИСТЕМЫ ОРТОПЕДИЧЕСКОГО АППАРАТА АВТОМАТИЗИРОВАННОГО ОСТЕОСИНТЕЗА // Фундаментальные исследования. – 2013. – № 6-1. – С. 46-50;
URL: https://fundamental-research.ru/ru/article/view?id=31411 (дата обращения: 29.03.2024).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1,674