Mathematical model of the movement of a tracked train for off-road container transportation

Cover Page

Cite item

Abstract

To determine the indicators of the mobility properties of unmanned tracked trains and the formation of traffic control laws, it is advisable to use mathematical simulation. On the example of a tracked semitrailer train, a mathematical model is being developed. It consists of a large number of rigid bodies connected by kinematic and power connections; therefore, it is proposed to use software systems that make it possible to synthesize a system of equations for the dynamics of a tracked train in an automated mode. The article discusses general approaches to the formation of such systems of equations in modern software packages for automated analysis of the dynamics of body systems. There are given a structural diagram of the developed mathematical model of a tracked train and a description of its main structural elements (blocks), which are responsible for the operation of various systems: a power plant, transmission and running system.

The features of the interaction of active sections of a caterpillar propeller with a low-deformable support base of the “dense soil” type are considered. On the example of the maneuver “turn with a given radius”, an imitation mathematical modeling of the movement of a caterpillar train was carried out, and on the basis of an analysis of the results obtained, the operability of the developed mathematical model was confirmed. The developed mathematical model of the movement of a tracked semitrailer train makes it possible to study its movement in space, taking into account the mutual influence of the tractor and the semitrailer, taking into account the interaction of the active sections of the tracked propeller with the sup-port base and taking into account the modeling of the working processes of its systems.

Based on the study, it can be concluded that the use of the method of computer modeling of the dynamics of body systems is an effective tool for predicting the performance indicators of tracked trains, determining the laws of traffic control, and also for assessing energy efficiency when performing virtual maneuvers. The developed mathematical model of the dynamics of a tracked train allows a wide range of studies to assess the mobility of tracked trains of various design and layout.

Full Text

Введение

В последние десятилетия наблюдается существенный рост объема контейнерных перевозок по сравнению с другими видами грузов [1]. Для наземной транспортировки контейнеров вне дорог могут быть использованы автопоезда с колесным и гусеничным движителем. Однако использование многоосных колесных транспортных средств и автопоездов все же требует подготовленного опорного основания [2].

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

 

Рис. 1. Конструктивно-компоновочная схема беспилотного гусеничного поезда

 

Цель исследований

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

Материалы и методы

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

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

Примером таких программных комплексов автоматизированного анализа динамики систем тел являются: ADAMS [3], Euler [4], «Универсальный механизм» [5], ФРУНД [6], MATLAB Simscape/Multibody [7]. В различных программных пакетах подход к автоматизированному составлению уравнений динамики системы может отличаться: например, синтез может осуществляться не в абсолютных координатах, а в относительных. Кроме того, уравнения динамики могут составляться в форме второго закона Ньютона, а не уравнений Лагранжа первого рода. Однако с позиции пользователя программных продуктов существенной разницы нет.

Уравнение пространственного движения твердого тела

Рассмотрим уравнения пространственного движения твердого тела, не обремененного связями (рис. 2); в общем случае их можно представить в следующей матричной форме [8, 9]:

ms¨=i=1kfFi;Jω˙=i=1kfRTFiri+j=1kMMj+Kω, (1)

где m – диагональная матрица масс тела; s – вектор координат цента масс тела в неподвижной (глобальной) системе координат ; Fi – вектор i-й внешней силы, приложенный к некоторой точке тела, заданный в глобальной системе координат ; kf – количество внешних сил, приложенных к телу; J – тензор инерции тела; ω – вектор угловых скоростей тела в подвижной системе координат xyz, связанной с телом в центре масс; R – ортогональная матрица поворота, переводящая вектор, заданный в подвижной системе координат xyz, связанной с телом в центре масс, в неподвижную XYZ; ri – радиус-вектор точки приложения i-й силы в подвижной системе координат xyz, связанной с телом в центре масс; Mj – вектор j-го внешнего момента, действующего на тело, заданный в подвижной системе координат xyz; kM – количество внешних моментов, действующих на тело; Kω – вектор моментов гироскопических сил (K – вектор кинетического момента тела).

 

Рис. 2. Расчетная схема движения твердого тела в пространстве

 

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

Уравнения кинематических связей

Для комплексного синтеза и решения уравнений динамики системы тел их необходимо представить в матричном виде. Так, система уравнений, описывающая динамику системы тел, может быть представлена в виде уравнений Лагранжа первого рода [8, 9]:

MJq¨+DTλ=fq˙,q,t;Dq¨=hq˙,q, (2)

где MJ – матрица инерции; q – вектор обобщенных координат (размерностью n, где n – количество обобщенных координат, необходимых для описания движения системы); D – матрица переменных коэффициентов уравнений кинематических связей размерностью kn, где k – число уравнений связи; λ – вектор множителей Лагранжа; fq˙,q,t – вектор внешних сил, включающий в себя силы нагрузок, силы от упругодемпфирующих элементов, а также гироскопические силы; hq˙,q – вектор правых частей уравнений связи.

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

Уравнения упругодемпфирующих связей

В качестве примера описания силовой связи между двумя телами и формирования вектора fq˙,q,t для системы (2) рассмотрим пружину и демпфер, соединяющие два твердых тела в произвольных точках A и B (рис. 3) [8]. Данный пример предполагает, что в процессе движения системы пружина и демпфер будут иметь постоянно изменяющуюся ориентацию в пространстве.

 

Рис. 3. Расчетная схема силовой связи в системе из двух твердых тел

 

Так, для получения величины и направления упругодемпфирующей силы определим положения точек A и B в глобальной (неподвижной) системе координат:

wA=s1+R1rA;

wB=s2+R2rB, (3)

где wA,wB – векторы координат точек A и B тел 1 и 2, соответственно, в неподвижной (глобальной) системе координат XYZ; rA, rB – радиус-векторы точек A и B в локальных системах координат x1y1z1 и x2y2z2 связанных с телами 1 и 2 соответственно.

Тогда величина деформации упругодемпирующего элемента определяется по формуле:

Δ=wAwBw0, (4)

где Δ – величина деформации упругогодемпфирующего элемента; w0 – вектор начального положения упругодемпфирующего элемента в пространстве.

Таким образом, величину силы в упругом элементе, как позиционной силы, зависящей от относительного смещения, можно определить:

Fp=fpΔ, (5)

где Fp – величина силы в упругом элементе; fpΔ – закон изменения упругой силы в зависимости от деформации упругого элемента.

Для определения направления силы Fp используются направляющие косинусы вектора (wAwB) в глобальной (неподвижной) системе координат XYZ:

aX=wAXwBXwAwB, aY=wAYwBYwAwB, aZ=wAZwBZwAwB, (6)

где aX,aY,aZ – направляющие косинусы вектора wAwB в глобальной системе координат XYZ.

Умножая величину Fp на вектор направляющих косинусов, получим искомый вектор силы упругого элемента в глобальной (неподвижной системе координат):

Fp=FpaXaYaZ. (7)

В случае вычисления демпфирующей составляющей силы аналогично определяются скорости точек A и B в глобальной системе координат:

Δ˙=w˙Aw˙B, (8)

w˙A=s˙1+R1ω1rA; w˙B=s˙2+R2ω2rB, (9)

Проекция вектора Δ˙ на направление (w˙Aw˙B) вычисляется при помощи направляющего косинуса:

Δ˙w=Δ˙cosαw; cosαw=Δ˙ΔΔ˙Δ, (10)

где Δ˙w – относительная скорость точек А и В в направлении вектора относительного перемещения Δ; αw – угол между векторами относительного перемещения Δ и скорости Δ˙ точек А и В.

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

Fд=fдΔ˙w, (11)

где Fд – величина силы в демпфирующем элементе; fдΔ˙w – закон изменения демпфирующей силы в зависимости от скорости деформации Δ˙w.

При этом направление вектора диссипативной силы в неподвижной системе координат:

Fд=FдaXaYaZ. (12)

Согласно третьему закону Ньютона, при формировании вектора fq˙,q,t для каждого из тел силы Fp и Fд необходимо приложить в соответствующих точках с разными знаками.

Математическая модель движения гусеничного поезда

Рассмотрим вопрос создания математической модели динамики гусеничного поезда, выполненного по полуприцепной схеме в программном комплексе Matlab, и встроенной в него среды разработки в составе Simulink и Simscape Multibody [7].

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

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

На рис. 4 представлена структурная схема разработанной математической модели гусеничного поезда, выполненного по полуприцепной схеме. Модель выполнена на базе программного комплекса MATLAB с использованием стандартных блоков библиотек Simulink и Simscape/MultiBody.

 

Рис. 4. Структурная схема математической модели гусеничного поезда, выполненного по полуприцепной схеме

 

Разработанная модель состоит из функциональных подсистем, которые соответствуют различным узлам гусеничного поезда: корпус тягача, корпус прицепа, четыре гусеничных модуля, силовая установка и система управления движением и контроля параметров. Корпус тягача и корпус полуприцепа представляют собой твердые тела с соответствующими массово-инерционными параметрами, которые соединены между собой сферическим шарниром с тремя вращательными степенями свободы. Для этого используется инструмент «Spherical Joint» библиотеки Simscape/MultiBody. Сферическим шарниром в рамках разработанной модели упрощенно учитывается работа седельно-сцепного устройства. Каждый из четырех гусеничных модулей представляет собой систему твердых тел, состоящую из опорных катков, шарнирно соединенных с соответствующими осями катков с помощью шарниров «Revolute Joint», которые имеют одну вращательную степень свободы и обеспечивают возможность вращения катков относительно их осей. Для катов гусеничных модулей и их осей заданы массово-инерционные параметры. Оси катков, в свою очередь, соединяются с корпусом тягача или прицепа при помощи поступательного шарнира «Prismatic Joint», который обеспечивает возможность перемещения катка в вертикальном направлении и моделирует работу системы подрессоривания при заданных характеристиках упругого и демпфирующего элементов.

Далее при моделировании гусеничного поезда будем рассматривать его движение по горизонтальному ровному основанию с «плотным» малодеформируемым грунтом. В таком случае распределение эпюры нормальных реакций по длине активной части гусеницы имеет разрывный характер [10]. Взаимодействие гусеничного движителя с опорным основанием осуществляется за счет активных участков, расположенных под катками; остальные участки гусеницы находятся практически в ненагруженном состоянии. Кроме этого, результаты, полученные в работе [11], позволяют перейти от рассмотрения взаимодействия гусеницы с опорным основанием к рассмотрению взаимодействия отдельных активных участков гусеничного обвода, расположенных под опорными катками.

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

Для определения направления нормалей к опорной поверхности рассматриваемый участок местности разбивается на Nxop точек вдоль продольной и Nyop вдоль поперечной осей (рис. 5). Тогда нормаль в каждой точке опорной поверхности может быть найдена как:

nop i,j=ri,j+1ri,jri+1,jri,jri,j+1ri,jri+1,jri,j, (13)

где nop i,j – нормаль к опорной поверхности в точке с координатами xоp i, yоp j; ri,j=xop iyop jzopxop i,yop jT – радиус-вектор рассматриваемой точки опорной поверхности с координатами xop i, yоp j в неподвижной системе координат XYZ.

 

Рис. 5. Расчетная схема определения нормали к опорной поверхности

 

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

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

Учет взаимодействия движителя с опорным основанием осуществляется с использованием подхода, основанного на представлении об «эллипсе трения», согласно которому сила взаимодействия с опорной поверхностью направлена противоположно скорости скольжения в плоскости опорного основания в точках контакта [12–14] (рис. 6, б). Подобная модель взаимодействия движителя с опорным основанием широко используется и при моделировании движения гусеничных машин и гусеничных поездов, например, в работах [15, 16].

 

Рис. 6. Расчетная схема взаимодействия активного участка гусеницы с опорным основанием

 

Расчетная схема взаимодействия активного участка гусеницы с опорным основанием представлена на рис. 6, а. Активный участок гусеницы, как отмечалось выше, разбивается на элементарные площадки, и для каждой такой площадки определяются ее ориентация, направление вектора скорости скольжения υ скi и касательная составляющая силы взаимодействия элементарной площадки активного участка гусеницы с опорным основанием Rэxy, которая определяется как произведение коэффициента взаимодействия колеса с опорной поверхностью μi на нормальную реакцию Rэz. Далее вычисляется вектор суммарной силы взаимодействия активного участка гусеницы с опорным основанием RΣ. Подробное описание реализации описанного подхода приведено в работе [17]. Подобная модель взаимодействия активных участков гусеницы с опорным основанием используется в разработанной модели динамики движения гусеничного поезда, выполненного по полуприцепной схеме. Относительная скорость активного участка гусеницы  в данном случае соответствует скорости перемотки гусеницы υотн, которая определяется следующим выражением:

υотн=υгус=ωвкrвк

где ωвк – угловая скорость ведущего колеса; rвк – радиус ведущего колеса.

Кроме этого, в модели используется функциональная подсистема, которая отвечает за моделирование работы силовой установки с учетом выбранной схемы трансмиссии. Математическая модель позволяет реализовать различные схемы трансмиссии гусеничных поездов путем изменения компонентов внутри функциональной подсистемы; например, трансмиссия с индивидуальным приводом ведущих колес, с дифференциальной связью, с использованием механизмов передачи-поворота и т.д.

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

Mдωд,h=Mдmaxh,  если ωд<NдmaxMдmax  и  ωд0Nдmaxhωд,  если ωдNдmaxMдmax  и  ωд<ωдmax0,  если  ωдωдmax, (14)

где Mдmax – максимальный тяговый или тормозной момент, создаваемый ТЭД; Nдmax – максимальная тяговая или тормозная мощность ТЭД; ωдmax – максимальная частота вращения ротора ТЭД; ωд – текущее значение частоты вращения ротора ТЭД; h – параметр управления, который одновременно пропорционально изменяет максимальный момент и максимальную мощность ТЭД.

 

Рис. 7. Механические характеристики ТЭД при различном параметре управления

 

Гусеничный поезд в зависимости от принятых конструктивно-компоновочных решений и выбранного способа поворота может иметь различные схемы трансмиссий. Так, для гусеничного поезда, оснащенного активным полуприцепом, с использованием бортового способа поворота целесообразно применять схемы трансмиссий, позволяющих перераспределять крутящие моменты между ведущими колесами. Простейшей схемой трансмиссии в таком случае является схема с индивидуальным приводом ведущих колес, изображенная на рис. 8, которая была реализована в разработанной математической модели. Такая трансмиссия не лишена недостатков, связанных с особенностью поворота гусеничной машины: при повороте вокруг отстающего борта с малым радиусом тяга на забегающем борту может оказаться недостаточной, что потребует установки переразмеренного ТЭД. Для решения этой проблемы могут использоваться другие схемы трансмиссий с использованием механизмов передачи-поворота с возможностью перераспределения потоков мощности между ведущими колесами. Разработанная математическая модель позволяет учитывать работу различных вариантов исполнений трансмиссий путем замены функционального блока математической модели на соответствующий выбранной схеме трансмиссии. Однако в рамках статьи будем рассматривать простейший вариант трансмиссии с индивидуальным приводом ведущих колес, который представляет собой два ТЭД с тормозами, соединенных через бортовые редукторы с ведущими колесами, движение которых описывается следующей системой уравнений:

Jпр1ω˙вк1=Mд1iбр1ηтрMс1Mт1iбр1ηтрJпр2ω˙вк2=Mд2iбр2ηтрMс2Mт2iбр2ηтр, (15)

где Jпр1, Jпр2 – приведенные моменты инерции вращающихся частей трансмиссии к ведущим колесам; Mд1, Mд2 – вращающие моменты, создаваемые ТЭД; Mс1, Mс2 – моменты сопротивлений на ведущих колесах; Mт1, Mт2 – моменты механических тормозов; iбр1, iбр2 – передаточные отношения бортовых редукторов; ηтр – КПД участка трансмиссии от ТЭД до ведущего колеса.

 

Рис. 8. Схема трансмиссии с индивидуальным приводом ведущих колес

 

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

Jпр1ω˙вк1=0,5MдiгпηтрMс1Mт1Jпр2ω˙вк2=0,5MдiгпηтрMс2Mт2, (16)

где Jпр1, Jпр2 – приведенные моменты инерции вращающихся частей трансмиссии к ведущим колесам; Mд – вращающий момент, создаваемый ТЭД; Mс1, Mс2 – моменты сопротивлений на ведущих колесах; Mт1, Mт2 – моменты механических тормозов; iгп – передаточное отношение главной передачи; ηтр – КПД участка трансмиссии от ТЭД до ведущего колеса.

 

Рис. 9. Схема трансмиссии с дифференциальным приводом ведущих колес

 

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

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

Результаты и обсуждение

Для проверки работоспособности разработанной математической модели проанализируем результаты моделирования гусеничного поезда полной массой 108 тонн. В качестве дорожно-грунтовых используются характеристики опорного основания типа «плотный грунт». Примем следующие его характеристики [19]: μx max=μy max=0,8 – коэффициенты сцепления в продольном и поперечном направлениях; fгр=0,06 – коэффициент сопротивления движению.

Рассмотрим расчетный случай криволинейного движения гусеничного поезда с активным полуприцепом с использованием кинематического способа поворота. До момента времени t = 5 c гусеничный поезд движется прямолинейно, затем начинает поворачивать до заданного угла «складывания» секций, равного 27,5°, и в момент времени t = 10 c тягач движется по окружности радиуса Rп = 30 м с установившейся скоростью 8 м/с. На рис. 10 представлены траектории движении тягача и полуприцепа и изменение скорости их движения.

 

Рис. 10. Траектория движения гусеничного поезда и изменение его скорости движения

 

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

Mсrвк=Rx+fгрRzηгус, (17)

где fгр – коэффициент сопротивления качению колеса в ведомом режиме; Rx, Rz – проекции продольной и вертикальной составляющих суммарной силы взаимодействия в контакте активных участков гусеницы с опорной поверхностью; rвк – радиус ведущего колеса; ηгус – КПД гусеничного движителя.

 

Рис. 11. Изменение моментов сопротивления на ведущих колесах тягача

 

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

Заключение

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

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

×

About the authors

K. B. Yevseyev

Bauman Moscow State Technical University

Author for correspondence.
Email: kb_evseev@bmstu.ru

PhD in Engineering

Russian Federation, Moscow

References

  1. Kato H., Shibasaki R., Ducruet C. Global logistics network modelling and policy. Quantification and analysis for international freight. Amsterdam: Elsevier, 2020. 372 p.
  2. Aksenov P.V. Mnogoosn·yye avtomobili [Multi-axle vehicles]. 2 izd. Moscow: Mashinostroyeniye Publ., 1989. 280 p.
  3. ADAMS – The Multibody Dynamics Simulation Solution. URL: https://www.mscsoftware.com/product/adams (accessed: 14.07.2021).
  4. EULER – Software Complex for Automated Dynamic Analysis of Multibody Mechanical Systems. URL: http://www.euler.ru/index.php/euler (accessed: 14.07.2021).
  5. Universal’nyy mekhanizm – programmnyy kompleks dlya modelirovaniya dinamiki mekhanicheskikh system [Universal mechanism - a software package for modeling the dynamics of mechanical systems]. URL: http://www.umlab.ru/pages/index.php (accessed: 14.07.2021).
  6. FRUND – programmnaya sistema formirovaniya resheniy uravneniy nelineynoy dinamiki [FRUND - software system for generating solutions to equations of nonlinear dynamics]. URL: http://frund.vstu.ru/ (accessed: 14.07.2021).
  7. MATLAB Simscape/Multibody – Model and simulate multibody mechanical systems. URL: https://www.mathworks.com/products/simscape-multibody.html (accessed: 14.07.2021).
  8. Gorobtsov A.S., Kartsov S.K., Polyakov Yu.A. Features of the construction of spatial dynamic models of vehicles, taking into account large movements of rigid bodies. Izvestiya TulGU. Tekhnicheskiye nauki. 2013. No 6-1, pp. 102–115 (in Russ.).
  9. Gorobtsov A.S., Kartsov S.K., Pletnev A.E., Polyakov Yu.A. Komp’yutern·yye metody postroyeniya i issledovaniya matematicheskikh modeley dinamiki konstruktsiy avtomobiley. Monografiya [Computer methods for constructing and researching mathematical models of the dynamics of vehicle structures. Monograph]. Moscow: Mashinostroyeniye Publ., 2011. 463 p.
  10. Zabavnikov N.A. Osnovy teorii transportnykh gusenichnykh mashin [Fundamentals of the theory of transport tracked vehicles]. Moscow: Mashinostroyeniye Publ., 1975. 448 p.
  11. Krasnen’kov V.I., Kharitonov S.A. Dynamics of curved movement of a transport tracked vehicle. Trudy MVTU im. N.E. Baumana. 1980. No 339, pp. 3–67 (in Russ.).
  12. Rozhdestvenskiy Yu.L., Mashkov K.Yu. The formation of reactions when an elastic wheel rolls on a non-deformable base. Trudy MVTU. 1982. No 390, pp. 56–64 (in Russ.).
  13. Ellis D.R. Upravlyayemost’ avtomobilya [Vehicle handling]. Moscow: Mashinostroyeniye Publ., 1975. 216 p.
  14. Dik A.B. Raschet statsionarnykh i nestatsionarnykh kharakteristik tormozyashchego kolesa pri dvizhenii s uvodom. Dissertatsiya … kandidata tekhnicheskikh nauk [Calculation of stationary and non-stationary characteristics of the braking wheel when driving with pull. Dissertation for degree of DrSc in Engineering]: 05.05.03 / SADI. Omsk, 1988. 224 p.
  15. Gorelov V.A., Kositsyn B.B., Miroshnichenko A.V. [i dr.] Method for determining the characteristics of an individual traction electric drive of a two-link tracked vehicle at the design stage. Trudy NGTU im. R.E. Alekseyeva. 2019. No 3 (126), pp. 120–134 (in Russ.).
  16. Kotiyev G.O., Padalkin B.V., Miroshnichenko A.V. [i dr.] Theoretical studies of the mobility of high-speed tracked vehicles with electric transmissions. Materialy mezhdunarodnoy nauchno-prakticheskoy konferentsii. Pod redaktsiyey I.A. Kalyayeva, F.L. Chernous’ko, V.M. Prikhod’ko. 2018, pp. 27–36 (in Russ.).
  17. Gorelov V.A., Yevseyev K.B., Chudakov O.I. [i dr.] Evaluation of indicators of curvilinear movement of a road train using mathematical simulation. Izvestiya MGTU «MAMI». 2020. No 4 (46), pp. 2–15 (in Russ.).
  18. Kotiyev G.O., Gorelov V.A., Miroshnichenko A.V. Development of a control law for the individual drive of the propellers of a multi-axle wheeled vehicle. Izvestiya vysshikh uchebnykh zavedeniy. Mashinostroyeniye. 2012. No 1, pp. 49–59 (in Russ.).
  19. Isakov P.P. Teoriya i konstruktsiya tanka. Transmissii voennyh gusenichnyh mashin [Theory and design of the tank. Military tracked vehicle transmissions]. Vol. 5. Moscow: Mashinostroyeniye Publ., 1985. 367 p.

Supplementary files

Supplementary Files
Action
1. Fig. 1. Structural and layout diagram of an unmanned tracked train

Download (153KB)
2. Fig. 2. Design scheme for the motion of a rigid body in space

Download (53KB)
3. Fig. 3. Design diagram of the force connection in a system of two solids

Download (45KB)
4. Fig. 4. Block diagram of the mathematical model of a tracked train, made according to a semitrailer scheme

Download (153KB)
5. Fig. 5. Calculation scheme for determining the normal to the reference surface

Download (53KB)
6. Fig. 6. Design scheme of interaction of the active section of the track with the support base

Download (72KB)
7. Fig. 7. Mechanical characteristics of traction electric motor with different control parameter

Download (78KB)
8. Fig. 8. Transmission scheme with individual drive of driving wheels

Download (110KB)
9. Fig. 9. Diagram of a transmission with a differential drive of the driving wheels

Download (96KB)
10. Fig. 10. Trajectory of the tracked train and the change in its speed

Download (52KB)
11. Fig. 11. Changing the moments of resistance on the driving wheels of the tractor

Download (56KB)

Copyright (c) 2021 Yevseyev K.B.

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies