Транспортное и энергетическое машиностроение
УДК 629.3.01:539.42 DOI 10.18698/0536-1044-2018-1-37-44
Анализ трещиностойкости торсионных валов гусеничной машины при эксплуатационных нагрузках
А.М. Покровский, Д.А. Дубин
МГТУ им. Н.Э. Баумана, 105005, Москва, Российская Федерация, 2-я Бауманская ул., д. 5, стр. 1
The Estimation of Fracture Toughness of Torsion Shafts of a Tracked Vehicle under Operational Loads
A.M. Pokrovskiy, D.A. Dubin
BMSTU, 105005, Moscow, Russian Federation, 2nd Baumanskaya St., Bldg. 5, Block 1 e-mail: [email protected], [email protected]
В настоящее время системы подрессоривания корпуса, использующие в качестве упругих элементов подвесок торсионные валы, устанавливают на большинство современных гусеничных машин. Основным видом отказов таких систем является разрушение торсионов вследствие образования усталостных трещин, что приводит к выходу из строя узлов подвесок. Проведен анализ трещиностойкости торсионных валов серийного гусеничного шасси, на базе которого разработано семейство многоцелевых транспортных машин ГТ-ТМ, ГТ-ТМС и др. Исследование выполнено при эксплуатационном уровне нагружения для трещины, расположенной на цилиндрической части торсионного вала, плоскость которой ориентирована под углом к оси торсиона и совпадает с положением главных площадок напряженного состояния. В основу расчета трещиностойкости положен силовой критерий разрушения Ирвина. Вычисления максимального коэффициента интенсивности напряжений по фронту трещины проведены с применением метода конечных элементов в программном комплексе ANSYS. Представлены результаты анализа трещиностойкости в виде зависимостей критической глубины трещины от соотношения ее полудлины к глубине. Полученные данные могут быть использованы при определении остаточного ресурса торсионных валов гусеничных машин на базе рассматриваемого шасси.
Ключевые слова: система подрессоривания, трещиностойкость торсионных валов, краевая трещина, коэффициент интенсивности напряжений, метод конечных элементов, критерий разрушения Ирвина
Nowadays, the hull suspension systems that use torsion shafts as elastic suspension elements are fitted to the majority of modern tracked vehicles. The main type of failure in such systems is the fracture of the torsion shafts due to the formation of fatigue cracks, which leads to failure of the suspension assemblies. This work presents an analysis of fracture toughness of torsion shafts of a standard tracked chassis used to develop a family of multi-
purpose transport vehicles GT-TM, GT-TMS, etc. The analysis is carried out under an operating load level for a crack located on the cylindrical part of the torsion shaft, the plane of which is at an angle to the torsion shaft axis and coincides with the position of the main areas of the stress state. The calculation of fracture toughness is based on Irwin fracture criterion. The calculations of the maximum stress intensity factor along the crack front are performed using the finite element method in the ANSYS software package. The results of the analysis of fracture toughness are presented in the form of dependences of the critical depth of the crack on the ratio of the fracture half-length to its depth. The data obtained can be used to determine the residual life of torsion shafts of the tracked vehicles based on the chassis under consideration.
Keywords: suspension system, fracture toughness of torsion shafts, edge crack, stress intensity factor, finite element method, Irwin fracture criterion
Элементы ходовой части быстроходных гусеничных машин эксплуатируются в условиях циклического и ударного нагружения. При этом усталостные разрушения торсионных валов (ТВ) являются наиболее типичной причиной отказов систем подрессоривания таких машин, где ТВ используют в качестве упругих элементов подвесок. В свою очередь, высокая нагруженность ТВ делает их очень чувствительными к различного рода поверхностным микродефектам, которые приводят к ускоренной инициации усталостных явлений и, как следствие, к образованию трещин. Отказ подвесок вследствие разрушения упругого элемента напрямую влияет на подвижность машины и, следовательно, на эффективность выполнения поставленных задач, что особенно критично для гусеничных машин специального назначения.
В настоящее время на стадии проектирования изделий в машиностроительной отрасли широко применяют методы линейной механики разрушения (ЛМР), позволяющие оценить работоспособность и период эксплуатации наиболее ответственных деталей при наличии в них трещиноподобных дефектов. В рамках ЛМР напряжения в окрестности фронта трещины описываются в полярной системе координат с помощью асимптотических формул [1]
-
f (6),
где Кг — коэффициент интенсивности напряжений (КИН), равный КI, Кц и Кш для трещин первого, второго и третьего типов соответственно; г — радиальная координата; £ (9) — тригонометрические функции угловой координаты 9 для трещины г-го типа.
Согласно силовому критерию разрушения Ирвина [1], неустойчивый рост трещины,
предшествующий полному долому, начинается, когда максимальный КИН по фронту трещины достигает критического значения (вязкости разрушения). Применительно к трещине первого типа этот критерий имеет вид
К1 < К1с, (1)
где К1с — вязкость разрушения для трещины типа I.
Таким образом, первоочередной задачей при оценке трещиностойкости ТВ является определение значения КИН, который зависит от условий нагружения, размеров вала и трещины, а также от ориентации дефекта. В соответствии с критерием Иоффе [1] траектория магистральной трещины развивается в направлении плоскости действия максимальных растягивающих главных напряжений, т. е. трещина в процессе роста стремится ориентироваться в поле напряжений так, чтобы отсутствовал сдвиг образующих ее берегов [1, 2]. Практика эксплуатации быстроходных гусеничных машин показывает, что разрушение ТВ происходит вследствие распространения трещины в плоскости, расположенной под углом к его оси (рис. 1).
Данный факт может свидетельствовать о том, что трещины растут в плоскостях, совпа-
Рис. 1. Разрушение ТВ по цилиндрической части
дающих с главными площадками напряженного состояния, в которых отсутствуют касательные напряжения и которые наклонены относительно оси нагруженного изгибающим и крутящим моментами вала.
Анализ работ и исследований в этой области показал, что задача в основном решалась для валов и стержней с трещиной, расположенной в поперечном сечении [3-6], что справедливо для валов, испытывающих чисто изгибное нагружение, и в целом не характерно для ТВ, где преобладающим видом нагружения является кручение, а напряженным состоянием — чистый сдвиг, для которого главные площадки ориентированы под углом 45° к оси вала.
Цель работы — определение критических размеров трещины, находящейся на цилиндрической части ТВ, плоскость которой расположена под углом к оси торсиона и совпадает с положением главных площадок.
Практическая ценность исследования заключается в оценке трещиностойкости ТВ системы подрессоривания серийного гусеничного шасси при реальных нагрузках, а именно в построении графиков критической глубины трещины в зависимости от ее длины. Полученные данные могут быть использованы для оценки трещиностойкости ТВ на этапе проектирования, а также при определении остаточного ресурса машин, находящихся в эксплуатации.
Объект исследования. Торсионные валы гусеничных машин изготавливают из легированной стали 45ХН2МФА, обладающей высокими прочностными свойствами. Типовая конструкция и схема нагружения ТВ приведены на рис. 2. Длина ТВ Ьт в современных системах подрессо-ривания сопоставима с шириной корпуса гусеничной машины. На концах вала имеются две шлицевые головки разной длины и диаметра со шлицами треугольного профиля. От шлицевой головки к цилиндрической поверхности выполнен плавный переход с радиусом галтели
г = (1...2)&, где й — диаметр ТВ.
Рис. 2. Типовая конструкция ТВ
Как правило, ТВ рассчитывают только на кручение, однако в действительности они испытывают упрощенное плоское напряженное состояние от действия крутящего момента Мкр (как преобладающего силового фактора со стороны рычага подвески) и изгибающего момента Мизг, возникающего вследствие несоосности опор и деформации корпуса.
Опыт эксплуатации показывает, что усталостные трещины образуются в области наибольшей концентрации напряжений — галтелей и шлицевых канавок. Однако существует высокая вероятность появления механических повреждений наружных слоев рабочей части вала, которые могут возникать при технологической обработке, транспортировке и неправильном хранении, что впоследствии приведет к его разрушению по цилиндрической части (см. рис. 1). Анализ характера разрушения ТВ по цилиндрической части указывает на распространение трещин в плоскости, расположенной под углом к оси вала в диапазоне 45...55°.
Объектом исследования являются усталостные трещины, образующиеся в ТВ диаметром dx = 37 мм. Эти валы входят в состав индивидуальной системы подрессоривания серийной гусеничной машины массой 14 т с передним расположением моторно-трансмиссионного отделения и двенадцатью независимыми подвесками. Анализ характера усталостных изломов показывает, что берега трещины близки к плоской поверхности. Максимальная глубина усталостной трещины в момент долома не превышает 5,6 мм.
Методика проведения исследования. Анализ КИН проводили с использованием метода конечных элементов в программном комплексе ANSYS [7]. Расчетная схема соответствовала консольно-закрепленному цилиндрическому стержню диаметром 37 мм и длиной 200 мм, нагруженному крутящим и изгибающим моментами (рис. 3). Для создания конечно-элементной модели ТВ с трещиной использовали 10-узловые тетраэдрические конечные элементы SOLID 187, а для описания особенностей вдоль фронта трещины — специальные сингулярные 15-узловые призматические конечные элементы SOLID 186.
Модель дефекта в зависимости от соотношения lb/la = 1,0.5,0 соответствовала краевой полукруговой или полуэллиптической трещине, расположенной в плоскости, нормаль
ОМф АЛ
m 1 V- - i.
Q 1 '
200
А-А
za, L 1ь
h
a _ iг?7/ —»1
Поверхностные слои рабочей части ТВ испытывают плоское напряженное состояние. Значение и направление первого главного напряжения можно определить по следующим формулам [8]:
h =1 ( + V0f+4rL );
(2)
0
Рис. 3. Схема нагружения сплошного цилиндра с расположенной под углом к оси трещиной
к которой образует угол а с осью цилиндра. Моделирование выполняли для двух угловых положений при а = 35° (у = 55°) и а = 45° (у = 45°) путем последовательного перестроения сетки конечных элементов, образующих берега трещины.
Количество элементов для каждой модели варьировалось в зависимости от отношения lb/la. Например, для трещины с отношением lb/la = 1,0 количество узлов и элементов равнялось соответственно 132 771 и 88 003, а для соотношения lb/la = 5,0 — 107 801 и 71 998. На рис. 4 представлена конечно-элементная сетка в области трещины.
На рис. 5 приведено распределение КИН по фронту полукруговой трещины при а = 35° и lb = 4,9 мм. Полученные в ANSYS значения КИН для самой заглубленной точки составили 1 939 МПа • мм1/2, а в точке, выходящей на поверхность, — 2 587 МПа • мм1/2, что при их переводе в общепринятые единицы соответствует 61,3 и 81,8 МПа • м1/2.
Определение критической глубины трещины la для каждого соотношения lb/la осуществлялось итерационным способом, согласно которому глубина становилась критической, когда максимальный КИН по фронту трещины на основании критерия Ирвина достигал вязкости разрушения Kic с погрешностью ±0,5 %. Условия нагружения задавались так, чтобы угол нормали к плоскости трещины а (см. рис. 3) совпадал с углом наклона первого главного напряжения р.
а, =-( + + '
о 1 2тх
Р = - arctg— 2 а*
с2 — нормальное напряжение; тХ2 — касательное напряжение.
Для трещины с углом а расчет проводился только от действия крутящего момента Мкр = = 6,1 кН-м, что соответствует эксплуатационному уровню нагружения ТВ для исследуемой системы подрессоривания [9]. При этом максимальное значение касательного напряжения тХ2 в поперечном сечении составило 613 МПа. На основании выражения (2) главное растягивающее напряжение с1 равно касательному тХ2, а направление их действия ориентировано под углом Р = 45° к оси цилиндра, что соответствует углу наклона нормали к плоско-
Рис. 4. Конечно-элементная сетка в окрестности трещины (деформированное состояние)
Рис. 5. Распределение КИН, МПа• мм1/2, по фронту полукруговой трещины при lb = 4,9 мм
сти трещины. При а = 35° расчет выполнялся при совместном действии крутящего Мкр = = 6,1 кН-м и изгибающего моментов Мизг = = 2,2 кН-м. В этом случае в поперечном сечении вала дополнительно возникают нормальные
напряжения от изгиба с 2, максимальное значение которых равно 442 МПа. Значение и угол наклона главных растягивающих напряжений на поверхности вала составили соответственно с1 = 873 МПа и Р = 35°. Следует отметить, что
K^Kn,Km, МПа-м
1/2
80 60 40 20 0
-20
—
1,0
0,5
0 а
0,5 у/с
К^К^Кт, МПа-м
1/2
80 60 40 20 0
-20
- 2 ---- 3 ■ —ъ
1 1
1,0
0,5
0,5
Ус
Ki, АГТП, МПа -м
1/2
80 60 40 20 0
-20
г 2 ' Л 3
| |
1,0
0,5
0 д
0,5
Ус
Ki, МПа -м
1/2
80 60 40 20 0
-20
1,0
0,5
0
0,5
Ус
Кт; Кп; Кт, МПа-м
1/2
80 60 40 20 0
-20
- 3 2
1 |
1,0
0,5
0
б
0,5
У/с
Ki, Кц; Кт, МПа • м
1/2
80 60 40 20 0
-20
2 ____ ——^^^^ 3
I 1
1,0 0,5
0,5 1ХПС
Ki, Ktf, Кщ, МПа • м
1/2
80 60 40 20 0
-20
2 1 Л 3 , Л.
| |
1,0
0,5
0,5
у/с
Ку, Кт, МПа • м
1/2
80 60 40 20 0
-20
2 \___ " Л 3 [. „
1 |
1,0
0,5
0,5
у/с
Рис. 6. Распределение КИН K (1), K (2) и Km (3) по фронту трещины при угле а = 45 (слева) и 35° (справа) и различных значениях отношения ее полудлины к глубине: а — lb/la = 1,0; la = 8,6 мм; б — lb/la = 1,0; la = 4,9 мм; в — lb/la = 2,0; la = 8,2 мм; г — lb/la = 2,0; la = 4,3 мм; д — lb/la = 3,0; la = 6,8 мм; е — lb/la = 3,0; la = 3,5 мм; ж — lb/la = 4,0; la = 6,2 мм; з — lb/la = 4,0; la = 3,2 мм
¡а кр> мм
Рис. 7. Зависимость критической глубины трещины ¡а кр в ТВ от отношения ¡ьИа при угле = 35 (1) и 45° (2)
для любой произвольной точки, лежащей на поверхности трещины, отношение тх2 /с 2 является постоянным, поэтому угол Р будет равен 35°.
В целях наглядности представления результатов расчета КИН для трещин разной длины использовали нормированные координаты точек фронта, которые вычисляли по зависимости
— ¡х /1с,
где ¡х — текущая координата фронта трещины; ¡с — координата точки «с» фронта трещины, выходящей на поверхность (см. рис. 3).
Результаты численного исследования. На
рис. 6 приведены результаты расчета КИН для плоских краевых полукруговых и полуэллиптических трещин, расположенных под углом а = = 45 и 35°.
Анализ представленных результатов расчета КИН по фронту трещины в условиях чистого сдвига (а = 45°) и упрощенного плоского
напряженного состояния (а = 35°) показал, что при совпадении плоскости трещины с плоскостью действия первого главного напряжения с большой степенью точности происходит растрескивание только по первому типу, о чем свидетельствуют несущественные по сравнению с Ki значения КИН второго Kii и третьего Кiii типов.
Кроме того, распределение КИН по фронту трещин для углов а = 45° и а = 35° имеет схожий характер. Для точки фронта трещины «с», выходящей на поверхность вала, наибольшее значение Ki достигается при полукруговой форме трещины (ljla = 1,0). Увеличение длины трещины, начиная с отношения lb/la > 2,0, приводит к смещению максимального КИН в наиболее заглубленную точку фронта «a» с уменьшением значения КИН по краям трещины.
Идентификация характера растрескивания по фронту трещины позволила провести анализ трещиностойкости с использованием силового критерия Ирвина (1) для трещины нормального отрыва при вязкости разрушения материала 45ХН2МФА Kic = 81,7 МПа • м1/2 [10]. На рис. 7 представлены результаты оценки трещино-стойкости ТВ гусеничной машины в виде зависимостей критической глубины трещины от отношения ее полудлины к глубине.
Как видно из рис. 7, в условиях упрощенного плоского напряженного состояния при угле а = = 35° с увеличением длины трещины критическая глубина стремится к la = 3,1 мм. В условиях чистого сдвига при угле а = 45° и отношении lb / la > 5,0 критическая глубина трещины соответствует 6 мм, что близко к таковой (5,6 мм), полученной на исследуемых валах (см. рис. 1).
Критические размеры трещин в ТВ
lb/la а = 35° а = 45°
la lb lc ls la lb lc ls
мм
1,0 4,90 4,90 4,88 4,91 8,60 8,60 8,54 8,57
2,0 4,30 8,60 8,22 8,35 8,20 16,40 14,98 15,48
2,5 3,80 9,50 8,78 8,94 7,30 18,25 15,72 16,31
3,0 3,50 10,50 9,28 9,48 6,80 20,4 16,36 17,04
3,5 3,30 11,55 9,70 9,92 6,50 22,75 16,87 17,63
4,0 3,20 12,80 10,12 10,37 6,20 24,80 17,09 17,89
4,5 3,15 14,18 10,49 10,78 6,10 27,45 17,44 18,30
5,0 3,10 15,50 10,75 11,06 6,00 30,00 17,66 18,56
Кроме графического представления, геометрические параметры критических трещин, полученные для различных значений отношения полудлины к глубине, приведены в таблице.
Выводы
1. Усталостные трещины в цилиндрической части ТВ можно считать трещинами нормального отрыва. При совпадении плоскости краевой полукруговой или полуэллиптической трещины с плоскостью действия первого главного напряжения происходит растрескивание по типу I (трещина нормального отрыва). КИН второго ^ и третьего ^^ типов по фронту трещины незначительны по сравнению с KI.
Литература
2. При эксплуатационном уровне нагруже-ния для трещин с отношением полудлины к глубине Ыи = 1,0.. .5,0 их критическая глубина в условиях чистого сдвига в ТВ варьируется от 8,6 до 6,0 мм. Для упрощенного плоского напряженного состояния критическая глубина полукруговой трещины составляет 4,9 мм, а при полуэллиптической форме с указанным отношением Ы^ разрушение произойдет при достижении глубины трещины ^ = 3,1 мм.
3. Критическая глубина плоской краевой полуэллиптической трещины, наклоненной относительно оси ТВ, уменьшается по мере увеличения ее длины. Причем начиная с десятикратного превышения длины над глубиной это влияние прекращается. Наименее опасными являются полукруговые трещиноподобные дефекты.
[1] Партон В.З., Морозов Е.М. Механика упругопластического разрушения. Основы меха-
ники разрушения. Москва, Изд-во ЛКИ, 2008. 352 с.
[2] Тихомиров В.М., Суровин П.Г. Развитие усталостных трещин смешанного типа в об-
разцах из стали. Прикладная механика и техническая физика, 2004, т. 45, № 1, с. 135-142.
[3] Fonte M., Gomes E., Freitas M. Stress intensity Factors for Femi-Elliptical Surface Cracks in
Round Bars Subjected to Mode i (Bending) and Mode ill (Torsion) Loading. European Structural Integrity Society: Multiaxial Fatigue and Fracture, 1999, vol. 25, pp. 249-260.
[4] Lissenden C.J., Tissot S.P., Trethewey M.W., Maynard M.P. Torsion response of a cracked
stainless steel shaft. Fatigue and Fracture of Engineering Materials and Structures, 2007, vol. 30, no. 8, pp. 734-747.
[5] Ismail A.E., Ariffin A.K., Abdullah S., Ghazali M.J. Stress intensity factors under combined
tension and torsion loadings. Indian Journal of Engineering and Materials Sciences, 2012, vol. 19, no. 1, pp. 5-16.
[6] Ismail A.E., Ariffin A.K., Abdullah S., Ghazali M.J. Stress intensity factors for surface cracks
in round bar under single and combined loadings. Mechanics, 2012, vol. 47, no. 6, pp. 1141-1156.
[7] Морозов Е.М., Муйземнек А.Ю., Шадский А.С. ANSYS в руках инженера. Механика
разрушения. Москва, ЛЕНАНД, 2010. 456 с.
[8] Феодосьев В.И. Сопротивление материалов. Москва, Изд-во МГТУ им. Н.Э. Баумана,
2016. 543 с.
[9] Дубин Д.А., Наказной О.А., Смирнов И.А., Шлеев А.Н. Экспериментальное определе-
ние кинематических и силовых параметров нагружения элементов системы подрес-соривания быстроходной гусеничной машины. Труды НАМИ, 2016, № 3(266), с. 45-53.
[10] Ковчик С.Е., Морозов Е.М., Панасюк В.В. Механика разрушения и прочность материалов. Т. 3. Характеристики кратковременной трещиностойкости материалов и методы их определения. Киев, Наукова думка, 1988. 436 с.
References
[1] Parton V.Z., Morozov E.M. Mekhanika uprugoplasticheskogo razrusheniia. Osnovy mekhaniki razrusheniia [Mechanics of elastoplastic fracture. Fundamentals of fracture mechanics]. Moscow, LKl publ., 2008. 352 p.
[2] Tikhomirov V.M., Surovin P.G. Razvitie ustalostnykh treshchin smeshannogo tipa v obraz-
tsakh iz stali [The development of fatigue cracks of a mixed type in steel samples]. Priklad-naia mekhanika i tekhnicheskaia fizika [Journal of Applied Mechanics and Technical Physics]. 2004, vol. 45, no. 1, pp. 135-142.
[3] Fonte M., Gomes E., Freitas M. Stress Intensity Factors for Femi-Elliptical Surface Cracks in
Round Bars Subjected to Mode I (Bending) and Mode III (Torsion) Loading. European Structural Integrity Society: Multiaxial Fatigue and Fracture, 1999, vol. 25, pp. 249-260.
[4] Lissenden C.J., Tissot S.P., Trethewey M.W., Maynard M.P. Torsion response of a cracked
stainless steel shaft. Fatigue and Fracture of Engineering Materials and Structures, 2007, vol. 30, no. 8, pp. 734-747.
[5] Ismail A.E., Ariffin A.K., Abdullah S., Ghazali M.J. Stress intensity factors under combined
tension and torsion loadings. Indian Journal of Engineering and Materials Sciences, 2012, vol. 19, no.1, pp. 5-16.
[6] Ismail A.E., Ariffin A.K., Abdullah S., Ghazali M.J. Stress intensity factors for surface cracks
in round bar under single and combined loadings. Mechanics, 2012, vol. 47, no. 5, pp. 1141-1156.
[7] Morozov E.M., Muizemnek A.Iu., Shadskii A.S. ANSYS v rukakh inzhenera. Mekhanika
razrusheniia [ANSYS is in the hands of an engineer. Mechanics of destruction]. Moscow, LENAND publ., 2010. 456 p.
[8] Feodos'ev V.I. Soprotivlenie materialov [Strength of materials]. Moscow, Bauman Press, 2016.
543 p.
[9] Dubin D.A., Nakaznoi O.A., Smirnov I.A., Shleev A.N. Eksperimental'noe opredelenie kine-
maticheskikh i silovykh parametrov nagruzheniia elementov sistemy podressorivaniia by-strokhodnoi gusenichnoi mashiny [Experimental determination of kinematic and power parameters of loading the suspension system elements of a high-mobility tracked machine]. Trudy NAMI [Proceedings of NAMI]. 2016, no. 3(266), pp. 45-53.
[10] Kovchik S.E., Morozov E.M., Panasiuk V.V. Mekhanika razrusheniia i prochnost' materialov: spravochnoe posobie. T. 3. Kharakteristiki kratkovremennoi treshchinostoikosti materialov i metody ikh opredeleniia [Mechanics of fracture and strength of materials. Vol. 3. Characteristics of short-term crack resistance of materials and methods for their determination]. Kiev, Naukova dumka publ., 1988. 436 p.
Информация об авторах
ПОКРОВСКИЙ Алексей Михайлович (Москва) — доктор технических наук, профессор кафедры «Прикладная механика». МГТУ им. Н.Э. Баумана (105005, Москва, Российская Федерация, 2-я Бауманская ул., д. 5, стр. 1, e-mail: [email protected]).
ДУБИН Дмитрий Андреевич (Москва) — ассистент кафедры «Многоцелевые гусеничные машины и мобильные роботы». МГТУ им. Н.Э. Баумана (105005, Москва, Российская Федерация, 2-я Бауманская ул., д. 5, стр. 1, e-mail: [email protected]).
Статья поступила в редакцию 12.09.2017 Information about the authors
POKROVSKIY Aleksey Mikhailovich (Moscow) — Doctor of Science (Eng.), Professor, Department of Applied Mechanics. Bauman Moscow State Technical University (105005, Moscow, Russian Federation, 2nd Baumanskaya St., Bldg. 5, Block 1, e-mail: [email protected]).
DUBIN Dmitriy Andreevich (Moscow) — Assistant Lecturer, Multipurpose Tracked Vehicles and Mobile Robots Department. Bauman Moscow State Technical University (105005, Moscow, Russian Federation, 2nd Baumanskaya St., Bldg. 5, Block 1, e-mail: [email protected]).