Введение. В современных отраслях производства все чаще применяются механизмы, использующие параллельную кинематику. Особенностью этих механизмов является подвижное выходное звено, соединенное с неподвижным основанием несколькими кинематическими цепями, параллельно передающими движение. Одним из этапов развития механизмов параллельной структуры является создание платформы Гью-Стюарта [1], позволяющей осуществлять перемещение платформы относительно основания по шести независимым координатам с помощью приводов линейных перемещений, что послужило предпосылкой для появления нового термина “Гексапод”. Гексаподы получили в последнее время достаточно широкое распространение, в первую очередь это роботизированные подвижные платформы (РПП) с кинематическими цепями по параллельной схеме, используемые в системах позиционирования инструмента в машиностроении и медицине, в тренажерах летательных аппаратов и боевой техники, для ориентации телескопов и антенн. Перспективность применения таких механизмов определяется широким спектром возможностей для решения технических задач. РПП, построенные на основе механизмов с параллельной кинематикой, могут обеспечивать следующие типы эксплуатационных задач: испытания высокотехнологичного оборудования и отдельных его компонентов; использование в качестве подвижных оснований в симуляторах для управления техникой (автомобили, самолеты, специализированная техника); создание специализированных устройств для удержания в заданном положении при движении различных приборов и исполнительных устройств. Применение механизмов с параллельной структурой достаточно эффективно реализовано в некоторых тренажерах, имитирующих процессы управления мобильными техническими средствами, где в качестве главных приводных устройств реализованы 6-ти степенные системы позиционирования (гексапод, платформа Стюарта, Hexapod, D. Stewart Platform).Широкое применение РПП параллельной структуры обусловлено конструктивными особенностями, обеспечивающими ряд эксплуатационных преимуществ: простота элементов механизма и их сборки, высокая жесткость конструкции, малый вес и низкая металлоёмкость при сравнимой жесткости, способность манипулировать большими нагрузками, повышенная точность перемещения и позиционирования, хорошие динамические характеристики, большие скорости и ускорения, возможность простого масштабирования схемы устройства. Параллельные механизмы имеют также и недостатки: сравнительно небольшая рабочая область, наличие особых положений, в которых происходит потеря управляемости, нелинейная зависимость кинематики и динамики от точки рабочей области, сложность системы управления, возможность интерференции опор.Эффективная реализация технических задач при использовании РПП параллельной структуры в различных областях послужила толчком для множества исследований, направленных на сохранение преимуществ при снижении или нивелировании недостатков этих механизмов [2–18]. Можно выделить следующие основные задачи, решаемые исследователями в этой области: - проработка и оптимизация структурной компоновки, определение возможностей и условий применения определенного типоразмера для различных эксплуатационных задач при проектировании новых конструктивных исполнений [2–5]; - определение возможностей применения определенного типоразмера платформы в соответствии с граничными значениями параметров её приводных звеньев для отработки заданных геометрических параметров траектории, например, предельно допустимых перемещений исполнительного органа в заданных направлениях [6, 7];- разработка и изучение законов управления перемещениями в приводных звеньях (электроцилиндрах) в соответствии с заданной траекторией перемещения конечного звена (обратная задача кинематики) [8–12];- тестирование системы управления с целью проверки требуемых траекторий перемещения выходного звена (исполнительного устройства или перемещаемого объекта), определение точности позиционирования [13];- кинематический, динамический и силовой анализ механизма, разработка алгоритмов определения множества достижимых положений платформы при формировании ее рабочей зоны [14–18].Для решения указанных научных задач исследователи применяют математические модели и различные аналитические методы [2, 6–9, 14, 17, 18], физические модели и лабораторные стенды [12, 14], а также прикладные программные средства [3, 5, 10–12, 16, 19, 20], которые открывают широкие возможности исследования и моделирования кинематических и динамических параметров механизмов. Данная работа является продолжением этого направления, а в качестве прикладных программных средств предлагается использовать Python и Adams.Основная часть. Объектом исследования является платформа Гью-Стюрта, представляющая собой шестиосевую конструкцию, имеющую шесть степеней свободы. Движение платформы реализуется с помощью шести телескопически связанных пар звеньев, соединенных посредством сферических шарниров с подвижной платформой и неподвижным основанием. Типовое строение платформы приведено на рис. 1, конструкция включает следующие основные компоненты: основание 1, подвижная платформа 2, электроцилиндры (ЭЦ) 3, шарниры 4.Рис. 1. Конструкция платформы Гью-Стюарта, основные механические компоненты:1 – основание; 2 – подвижная платформа; 3 – электроцилиндры, 4 – шарнирыИсследуемая роботизированная платформа подвижности предназначена для тренажеров и симуляторов в авиационной и ракетно-космической отрасли. Технические характеристики и кинематические параметры, соответствующие функциональному назначению, представлены в таблицах 1 и 2. Требуемые эксплуатационные параметры положены в основу при создании твердотельной модели: габаритные размеры, размеры базовых элементов, величина хода ЭЦ, расположение опор (шарниров), начальное угловое положение оси ЭЦ, пространственное положение центра масс. Таблица 1Технические характеристики РПП Наименование параметраЕд. изм.ЗначениеШиринамм1855Длинамм1654Высотамм1385Вескг550Ход штока ЭЦмм200Полезная нагрузкакг987Таблица 2Параметры кинематики РПП Наименование параметраПеремещениеСкоростьУскорениеПродольное перемещение-110 .. +120 мм0,66 м/с6,5 м/с2Боковое перемещение-110 .. +110 мм0,66 м/с6,5 м/с2Вертикальное перемещение-110 .. +110 мм0,66 м/с6,5 м/с2Крен-10°..+10°38°/с300°/с2Тангаж-10°..+10°38°/с300°/с2Рыскание-10°..+10°38°/с300°/с2 Для имитационного моделирования кинематических и динамических параметров, характеризующих эксплуатационные условия под действием рабочих нагрузок разработан виртуальный прототип роботизированной платформы, представленный на рис. 2. Цифровой макет, имеющий свойства параметризованной имитационной модели создан в программном комплексе MSC Adams и состоит из 14 основных идеализированных элементов (деталей): основание, подвижная платформа, шесть штоков и гильз. Данные элементы соответствуют приведенным в таблицах 1 и 2 геометрическим и функциональным параметрам рассматриваемой конструкции РПП.  Рис. 2. Цифровая имитационная модель РПП, выполненная в MSC Adams Для управления параметризованной моделью с возможностью ее автоматизированного перестроения в MSC Adams использовался интерфейс высокоуровневого объектно-ориентированного языка программирования общего назначения Python. Язык программирования Python применен в качестве альтернативы внутреннему командному языку Adams View для создания и итерационного изменения объектов моделирования в Adams. Интерфейс Adams Python представляет собой интерфейс прикладного программирования, который позволяет взаимодействовать с Adams, используя объектно-ориентированную структуру, в которой каждая сущность в Adams сопоставляется с классом в Python, имеющим свойства и методы.Выбор Python обусловлен возможностью решать поставленную задачу оптимизации конструктивных элементов и нахождения оптимального геометрического исполнения гексапода в соответствии с выбранными критериями оптимальности посредством проведения серий экспериментов, заключающихся в последовательном многократном изменении геометрических параметров, последующей симуляции модели и анализа результатов с целью нахождения вариантов конструкции, удовлетворяющих заданным критериям. Использование Python позволяет автоматизировать каждый из описанных шагов вычислительного эксперимента с помощью использования специальных процедур и функций, при этом программная реализация исключает необходимость трудоемкого выполнения перечисленных манипуляций вручную. Таким образом, применение Python позволяет существенно повысить результативность проведения вычислительного эксперимента за счет замены ручных манипуляций создания и перестроения твердотельных элементов в Adams автоматизированными программными процедурами, сокращая трудоемкость, повышая произвольность процесса и количество возможных вычислительных симуляций.Для реализации вышеописанной задачи были разработаны следующие специальные процедуры и функции:– процедура FindPoint для определения координат точек сопряжения электроцилиндров с основанием и столом, которая задает координаты соответствующих точек в полярной системе координат;– процедура CreateModel для создания модели гексапода в AdamsView, в которой с помощью интерфейса Adams Python вызываются необходимые процедуры для создания компонентов модели;– процедура ChangeModel для изменения геометрических параметров модели, в которой выполняется пересчет геометрических параметров, таких как расположение (Location) и ориентация (Orientation) для всех компонентов модели, имеющих данные свойства;– функция Simulations (qt, number_of_steps, end_time) для выполнения вычислительного эксперимента, связанного с множественными изменениями конструкции и последующими симуляциями, которая принимает на вход необходимые параметры:qt ­­– требуемое количество симуляций, number_of_steps – количество шагов в симуляции, влияющее на точность и быстродействие процесса симуляции, end_time – конечное время симуляции. В данной функции в цикле размером qt выполняются необходимые для эксперимента изменения геометрических параметров, перестроение модели (ChangeModel) и выполнение непосредственно симуляции с помощью соответствующего вызова команды интерфейса Adams Python. Функция возвращает кортеж, состоящий из трех массивов результатов, включающий амплитудные и средние значения параметра, относительно которого выполняется оптимизация, а также массив геометрических значений;– процедура AnalisesOfRes использует выходные параметры функции Simulations для анализа результатов вычислительного эксперимента и выполнения поиска симуляций, соответствующих исследуемым геометрическим параметрам модели по заданным критериям оптимальности. Результаты анализа выводятся в командную строку Adams.Для имитации реальных конструктивных сопряжений РПП использованы следующие специальные программные операторы приложения Adams: основание (Base) зафиксировано на Ground при помощи FixedJoint, крепление гильзы к основанию выполнено посредством SphericalJoint, аналогично сопряжены штоки ЭЦ с рабочей поверхностью платформы. Для имитации движения электроцилиндра гильза и шток связаны между собой в точках контакта сопряжением TranslationalJoint. В качестве внешних сил заданы сила тяжести, направленная вертикально вниз (по направлению оси – ОY) и полезная нагрузка, приложенная в точке центра масс, который задан в соответствии с требованиями проектируемого промышленного образца РПП.В качестве имитации движения РПП в шести цилиндрических шарнирах созданы поступательные движения при помощи программного оператора TranslationJointMotion, реализующего в приложении MSC Adams установленные законы линейного перемещения. Зависимость перемещения от времени может быть задана при помощи математической функции или специальной функции STEP. Функция STEP аппроксимирует идеальную математическую кусочно-заданную функцию, но без разрывов. Синтаксис функции STEP: STEP (q, q1, f1, q2, f2), где: q – независимая переменная q1 – начальное значение для q; f1 – начальное значение для f; q2 – конечное значение для q; f2 – конечное значение для f. steptime, t0, l0, tn,lk∙l+steptime, tn, lk, tm,lp∙l+…  ,                         (1) где time – независимая переменная (время); t0 – начальное значение времени для 1 участка; tn  – конечное значение времени для 1-го участка и начальное значение времени для 2-го участка; tm  – конечное значение времени для 2-го участка;   l0,lk, lp - начальные и конечные значения перемещения    штока электроцилиндра,  -1≥l0,lk, lp≤1 ; l – максимальное значение перемещения штока электроцилиндра. Результаты. Для проверки работоспособности методики моделирования и правильности взаимодействия программных модулей первоначально был проведен тестовый вычислительный эксперимент. В качестве тестируемого рабочего движения принято вертикальное движение подвижной платформы вдоль оси OY из исходного положения. Для осуществления движения по данной траектории для каждого цилиндрического шарнира созданы поступательные движения с функцией Motion = 0,2∙Sin(time) , где коэффициент равный 0,2 м определяет ход штока. На рисунке 3 представлен график изменения величины выдвижения одного из штоков во времени. В качестве оптимизируемого параметра принят диаметр стола с учетом конструктивного ограничения максимального габаритного размера при постоянном диаметре основания. В качестве критерия оптимальности выбрана минимизация затрачиваемых усилий на осуществление движения. Для выбора вариантов, удовлетворяющих заданному критерию оптимальности, исследовались амплитудные и средние значения силы, возникающие в шарнирах при осуществлении движения. Для объективности анализа влияния соотношения диаметра стола к диаметру основания на развиваемые усилия в шарнирах сопряжения электроцилиндров, приняты попарно одинаковые углы пространственного расположения шарниров основания и стола, заданные в полярных координатах.  Рис. 3. График изменения величины выдвижения штока от времени, м Таким образом, исследовалось влияние соотношения диаметра стола к диаметру основания на развиваемые усилия. Вычислительный эксперимент заключается в инициализации заданного количества симуляций, отличающихся геометрическими параметрами модели, и анализе влияния изменения этих параметров на заданный критерий оптимальности. Вычислительный эксперимент проводился при следующих условиях: диапазон изменения радиуса стола 0,25∙827–1,75∙827 мм, количество симуляций qt100 , шаг изменения радиуса стола 8,27 мм. На рис. 4 изображены твердотельные модели первой и последней симуляций. Время проведения вычислительного эксперимента составило 15 секунд, что является достаточно хорошим результатом, учитывая количество моделируемых случаев, следовательно, компоненты программных модулей взаимодействуют эффективно. На рис.5 представлен график изменения во времени возникающей при движении штока реакции силы в шарнире гидроцилиндра, соответствующий величине выдвижения штока рис. 3. Рис. 4. Твердотельные модели вычислительного эксперимента:a – первая симуляция, б – последняя симуляция Рис. 5. График изменения возникающей реакции силы в шарнире гидроцилиндра, Н В таблице 3 представлены результаты анализа параметров симуляций по критерию оптимальности: шесть симуляций с наименьшими значениями амплитудных сил, шесть симуляций с наименьшими значениями средних сил и для сравнения одна симуляция с наихудшими значениями амплитудных и средних сил. Как видно, максимальные силовые параметры возникают при минимальном размере стола (диаметр 219, симуляция № 1) и достигают минимальных значений по мере приближения к размеру стола, равному основанию (диаметр 827, симуляции №50).  Таблица 3Результаты анализа параметров симуляций по критерию оптимальности № симуляцииРадиус стола(подвижной платформы), ммАмплитудное значение силы, НСреднее значение силы, Н 12195347177449815–1613,30508273226,231613,20518393226,471613,24528523226,651613,30538643227,181613,50548773228,071613,90558893229,31– Полученные результаты соответствуют известным физическим законам и традиционным представлениям в области теоретической механики, кинематики и динамики механизмов. Следовательно, предложенная методика может использоваться для более сложных задач, например, исследования динамических параметров (сил, моментов сил) в элементах конструкции при движении платформы по сложной траектории с приложенной полезной нагрузкой, определения нежелательных положений, в которых возникают потери устойчивости или управляемости, пиковые нагрузки, значения которых превышают максимальные по конструктивным ограничениям. На рис. 6–8 представлены результаты определения одного из таких нежелательных положений: рисунок 6 демонстрирует состояние платформы в одном из крайних конструктивных положений, что подтверждается возникновением пиковых нагрузок в конструктивных элементах на 8 секунде отработки траектории (рис. 6, 7). Вычислительный эксперимент выполнялся по описанной выше методике взаимодействия Python и Adams, а для нахождения рабочей зоны использовался итерационный метод. Из исходного положения платформы симулировалась сложная траектория движения в пространстве, охватывающем все конструктивно возможные положения рабочего органа. Нежелательные положения находились программно в случае возникновения ошибки симуляции, связанной с нарушением наложенных ограничений. Рабочая зона формировалась также программно путем исключения выявленных нежелательных положений из исходного пространства. Рис. 6. Состояние платформы в момент возникновения программной ошибки (нежелательного положения)Рис. 7. Сила, возникающая в сферическом шарниреРис. 8. Силы, возникающие в приводах актуаторов Выводы. Имитационное моделирование кинематических и динамических параметров РПП посредством применения программного обеспечения MSC Adams в комбинации взаимодействия с объектно-ориентированным языком программирования общего назначения Python, имеет следующие преимущества: позволяет анализировать различные конструктивные исполнения, выдвигать и проверять гипотезы относительно геометрических параметров исследуемой модели, выявлять оптимальные варианты. Предложенная методика будет эффективна также при решении более сложных задач, например, определения рабочей зоны, выявления особых положений, сопровождающихся неустойчивой кинематикой и пиковыми нагрузками. Представленный подход применим в рамках разработки, модернизации, оптимизации конструктивных элементов РПП, а также разработки систем управления движением.