CC BY
36
НАУЧНОЕ ИЗДАНИЕ МГТУ ИМ. Н. Э. БАУМАНА
НАУКА и ОБРАЗОВАНИЕ
Эл № ФС77 • 48211. Государственная регистрация №0421200025. ISSN 1994-0408
электронный научно-технический журнал
О повышении точности определения комплексных
погрешностей сателлитных узлов планетарных передач
# 10, октябрь 2012
Б01: 10.7463/1012.0479419
Насонов Д. А., Леонтьев М. Ю., Шубин А. А.
УДК 621.83.05
Россия, КФ МГТУ им. Н.Э. Баумана nasonovda@yandex.ru newell-kaluga@mail.ru shubin55@mail.ru
Работоспособность, прочность и виброактивность зубчатых редукторов, выполненных по планетарной схеме, в значительной степени зависит от равномерности распределения нагрузки по параллельным потокам передаваемой мощности (между сателлитами) [1, 2]. В свою очередь, соблюдение названного условия во многом определяется погрешностями изготовления и сборки сателлитных узлов [3].
В работе [1] предложено в процессе сборки редуктора минимизировать указанные погрешности посредством регулировки межцентровых расстояний сателлитов с помощью эксцентриковых втулок, устанавливаемых в расточки водила под опорные шейки их осей. Однако для реализации данного предложения необходимо максимально точное определение фактических значений погрешностей сателлитных узлов, что, как показывает опыт одного из ведущих отечественных производителей планетарных редукторов большой мощности [3], является весьма нетривиальной задачей.
Рассматриваемые погрешности рассчитываются по данным технологических карт контроля точности, заполняемых специалистами ОТК завода-изготовителя на основе результатов прямых и косвенных измерений первичных ошибок изготовления элементов планетарной ступени в рамках «комплексного точностного параметра» (далее «КТП») планетарного механизма.
Изначально «КТП» был разработан для планетарных редукторов с прямозубыми цилиндрическими зубчатыми колесами [1]. Впоследствии он был распространен на шевронное зацепление и неоднократно уточнялся [4] Последняя уточненная версия «КТП» была приведена в работе [5], где на основе экспериментальных данных была выявлена корреляция указанного параметра с зубцовой вибрацией планетарного редуктора.
Согласно [5], комплексные погрешности изготовления и сборки сателлитных узлов шевронного планетарного редуктора могут быть определены, как
Л, = СОБат
( \ 2 X"1 Л V , Л „о-Ъ х л „о-а , о л „о(а-Ъ) \ и СОБ Рй
,-^ЛХкМ1 +Лв°-Ь +ЛеГ + 2Ле°(а-Ъ) -
СОБ У к=1 СОБ
V Ла У
(1)
Здесь Лг - погрешности суммарных зазоров в зацеплениях полушевронов сателлитов с
центральными колесами для случая, когда нумерация сателлитов от г до 2а совпадает с направлением окружных сил на их осях (с направлением крутящего момента на солнечной шестерне), причем погрешности вычисляются отдельно для каждого полушеврона; ат - угол зацепления в торцовом сечении (делительный угол профиля в торцовом сечении), определяемый из соотношения
ат = arаtg-— ,
СОБ рй
где ап - угол профиля исходного контура рейки в нормальном сечении, Рй - угол наклона зуба на делительном цилиндре;
г-1 П
^ ЛХк к+1 - накопленные ошибки хордальных расстояний между первой и г-й осями
к=1
расточек в щеках водила;
Ле°~а, Ле°-Ъ - погрешности диаметральных зазоров в соединениях осей с сателлитами и водилом, отнесенные к полушевронам сателлитов:
Лео-Ъ = ЛёЬ - Лё0-Ъ; Лео-а = Лёа - Лёо-а,
где ЛёЪ - погрешности диаметров расточек под оси в щеках водила,
ЛёС - погрешности средних по длине полушевронов диаметров расточек в сателлитах под оси,
Лё°-Ъ - погрешности диаметров посадочных шеек осей под установку в щеки водила, Лё°~а - погрешности средних по длине полушевронов диаметров посадочных поверхностей осей под сателлиты;
Л3°(а-Ъ) - погрешности отнесенных к полушевронам тангенциальных составляющих
эксцентриситетов посадочной поверхности г-й оси под сателлит относительно ее посадочных поверхностей под установку в водило;
ЛЪг - постоянные составляющие погрешности толщины зубьев г-го сателлита по полушевронам.
При использовании выражения (1) для расчета «КТП» наибольшие трудности возникают с определением хордальных расстояний между расточками под сателлиты в водиле ( ЛХк к+1), что связано с отсутствием у большинства отечественных производителей
редукторов большой мощности современного высокоточного оборудования, позволяющего получать координаты расточек столь крупногабаритных и имеющих сложную геометрию узлов, как водило (рис. 1), путем прямых измерений.
На практике ошибки взаимного расположения расточек в водиле определяются по результатам измерения межцентровых расстояний и скрещивания, специальных контрольных валов, устанавливаемых с заданной точностью в указанные расточки. В общем случае измерениям на каждой из щек водила подлежат расстояния между осями расточек
под сателлиты и центральной осью водила (далее - радиальные межцентровые расстояния) и между осями соседних расточек под сателлиты (далее - хордальные межцентровые расстояния). Скрещивания осей определяются для расточек водила под сателлиты относительно центральной оси и между осями соседних расточек под сателлиты.
В работе [1] из всех названных параметров для определения искомых погрешностей изготовления водила рекомендуется использовать только хордальные межцентровые расстояния, пренебрегая всеми остальными измерениями. Представляется, что такое искусственное сужение используемых в расчетах исходных данных едва ли оправдано, поскольку может давать значительное отклонение получаемых результатов от фактических значений искомых ошибок.
Для повышения точности их определения в работе [6] были предложены математическая модель и алгоритм численного решения, позволившие получать искомые координаты расточек под сателлиты в водиле с учетом всех контролируемых параметров. При этом были преодолены вычислительные трудности, связанные с необходимостью решения переопределенной (из-за значительно возросшего числа исходных данных) и несовместной (из-за неизбежных погрешностей измерения) системы уравнений.
Рассмотрим указанную модель и соответствующий ей алгоритм численного решения на примере определения координат расточек под оси для водила планетарной ступени с пятью сателлитами (рис. 1).
Для составления модели введем декартову систему координат с центром в плоскости носовой щеки водила таким образом, чтобы центральные расточки водила определяли координатную ось Х, а центр расточки под ось первого сателлита определял направление координатной оси Y (рис. 1). Обозначения остальных параметров иллюстрируют рис. 2-4.
На рис. 2 точками P0 и P'0 обозначены центры левой и правой расточек водила под опорные шейки, точками Р1 и Р'1 - центры расточек под ось первого сателлита. Из-за
у
Рис. 1 Модель водила планетарного редуктора с пятью сателлитами
погрешностей изготовления точка Р'1 может отклоняться от плоскости, определяемой точками Р1; Ро и Р'о, а величина этого отклонения характеризует скрещивание оси первого сателлита относительно центральной оси водила, принимаемой за базу (ось Х, рис. 1).
В результате измерений радиальных межцентровых расстояний водила определяются величины 81^85 (рис 3). За положительное значение Б; на
рис.3 принято смещение точки Р'; относительно Р; вокруг оси X (ось Р0-Р'0) против часовой стрелки, если смотреть с конца оси X, штрихами обозначены центры расточек на правой щеке, без штрихов - на левой.
У
Р1
Ро
1 1 1 \ \ Б X 1 р° *
р!
1
/
Рис. 2. Схема измерения взаимного расположения расточек под ось первого сателлита
относительно центральной оси водила
У
Я' '
1
-Я.
Р
г
Я
Ях
Р'
Я1
'я5
*ЛРо, Р0
Р5
// / /
Я
\ \ \ \ \ % \ \ \ N
ЧЛЯ4
Я'
Р' С±
ж
Р4 ^4
Рис 3. Вид на центры расточек (Р^ Р'1 ) водила с конца оси Х
Измерением хордальных межцентровых расстояний водила и скрещивания осей соседних расточек под сателлиты определяются величины ^=|Р; Pi.il (рис. 4), Ь'{=|Р'{ Р'м| и 5; (геометрическая интерпретация параметра 5i (рис. 5) рассмотрена ниже), где ¡=1^5. При этом, в силу цикличности индексации, имеем Ь1=|Р1 Р5|, Ь'1=|Р'1 Р'5|. Нумерация центров расточек Р; осуществляется против часовой стрелки, если смотреть с конца оси Х.
Если в качестве допущения считать, что обе щеки водила параллельны, т.е. центры всех расточек лежат либо в плоскости У2 (левая щека), либо в плоскости, параллельной У2 (правая щека), то можно не сводя задачу к двухмерной, исключить из расчетов координату X.
Учитывая, что точки Р0, Р'о приняты за базовые, а направление оси У определяется центром первой расточки (координата ъ точки Р1 равна нулю), в качестве неизвестных имеем у-ъ координаты центров расточек под оси сателлитов на левой и правой щеках, т.е. 19 величин.
Рис. 4. Схема расположения центров расточек в плоскости щеки водила У2 (Р) и угловых величин для тригонометрических соотношений (а)
Обозначив угловые величины в соответствии с рис. 4, можно записать группу из 19 тригонометрических соотношений связывающих результаты измерений и координаты центров расточек на левой щеке водила (9 неизвестных):
' * = К1, г 2 =
1 (2) ¿2 = И2 соБ(а1),
- + = ^2 ),
Аналогичные 20 уравнений, которые можно записать относительно у-ъ координат центров расточек на правой щеке (10 неизвестных), дополняют систему (2). Угловые величины в системе (2) взяты в предположении, что точки Р1 - Р5 (рис. 4) и Р'1 - Р'5 образуют правильные пятиугольники. На самом же деле, при наличии погрешностей изготовления водила, правильность этих пятиугольников нарушается. Следовательно,
зависимость расчетных координат от измеряемых параметров носит нелинейный характер, а система (2) представляет собой линеаризованную математическую модель, позволяющую, как показано ниже, получать приемлемые для практического использования приближенные решения.
Для учета всех имеющихся результатов измерений, запишем еще две группы уравнений. Первую группу получим из рассмотрения скрещивания осей расточек под сателлиты относительно центральной оси водила. Так, например, для расточки под ось сателлита №2 (1=2), параметр э2 (рис. 3) будет связывать координаты только центров расточек Р2 и Р'2, поскольку центральная ось является базой:
- ^2 + ^ = ^2 эт^),
У2 + У' ^^СОв^^ (3)
С учетом параметров скрещивания - система (3) содержит 9 уравнений относительно все тех же 19 неизвестных, названных выше.
Вторая группа уравнений получается при анализе скрещивания 5 осей расточек под соседние сателлиты. В этом случае, в отличие от предыдущего, ни одна из осей не является базовой, и параметр 5 связывает не две, а четыре точки. На рис. 5, поясняющем геометрический смысл параметра 51, показан вид на 1-ю и (1-1)-ю расточки со стороны конца оси Х.
Рис 5. Вид на центры расточек под оси соседних сателлитов со стороны конца оси Х
(рис. соответствует г = 2)
Величина 51 характеризует скрещивание 1-й и (1-1)-й осей расточек и состоит из двух компонент: 51=5^-5^, где 5ц и 512 - расстояния от точек Р'1-1 и Р'1 до линии Рм-Рь в проекции на плоскость У2. Причем за положительное значение принимается величина 5, если Р'1 (Р'1-1) и Р'о (центр водила) располагаются по разные стороны от линии Р1-1-Р1. Ту же величину имеем при измерении расстояния от точки Р"1 до линии Р1-1-Р1 (точка Р"1 получается путем параллельного переноса Р'1 на вектор е, совмещающий точки Р'1-1 и Р1-1). Именно эта величина и подлежит непосредственному контролю в процессе измерений, когда, согласно принятой технологии измерений, за базу принимается (1-1)-я ось.
Параметр 52, например, позволяет записать следующие, входящие в систему четыре, уравнения:
- + ¿2 + - = ) , У + У + У\ - У1 = $2 С08(а! ) , (4)
Система (4), учитывающая параметры 6;^55, содержит десять уравнений.
После объединения систем (2), (3) и (4), получается система из пятидесяти восьми уравнений относительно девятнадцати неизвестных, представляющая собой математическую модель, пригодную для вычисления погрешностей изготовления водила по результатам их косвенных измерений, и учитывающую, в соответствии с поставленной задачей, все имеющиеся результаты.
Как уже отмечалось, данная система является не только переопределенной и несовместной (из-за погрешностей измерений), но и преднамеренно линеаризованной с целью максимального упрощения. Для ее решения был разработан итерационный алгоритм [6], позволяющий преодолеть последствия использованных упрощений.
Согласно данному алгоритму первое приближенное решение полученной системы уравнений ищется с использованием геометрических параметров а; (рис. 4), определяемых в предположении, что искомые координаты расточек водила под сателлиты образуют правильные пятиугольники. Далее угловые величины а; уточняются с использованием результатов первой итерации, и расчет повторяется. Последующие приближения рассчитываются аналогичным образом до достижения требуемой точности.
Критерием точности получаемых результатов на каждой итерации служит сумма квадратов невязок (разностей между левыми и правыми частями уравнений уточненной обобщенной системы, при подстановке в неё вычисленных на последней итерации неизвестных), свидетельствующая о степени отклонения расчетных величин от измеряемых. Стабильное уменьшение этой величины, отмеченное при решении тестовых примеров, свидетельствует о хорошей сходимости итерационного процесса [6].
В связи с повышенными требованиями к точности расчетов описанный алгоритм реализован в пакете символьных вычислений Стивена Вольфрама «МаШешайса», позволяющем производить расчеты с любой наперед заданной точностью, независимо от разрядности используемой машины [7].
Анализ результатов выполненных тестовых расчетов показал, что одной из особенностей описанной математической модели является то, что при равноценности всех используемых в качестве исходных данных результатов измерений, их влияние на искомые координаты расточек под сателлиты в водиле существенно различается. Это происходит из-за того, что тригонометрические функции, находящиеся в правой части уравнений (2)-(4), играют роль весовых коэффициентов, и получаемое решение неизбежно «смещается» в сторону тех параметров, при которых стоят большие весовые коэффициенты.
При попытке перенести эти коэффициенты в левую часть, обусловленность системы, к решению которой приводит метод наименьших квадратов, резко ухудшается и сходимость решения недопустимо падает. Дело в том, что синусы некоторых углов в окрестностях решения столь малы, что перенос их в левую часть подобен делению на ноль. С другой стороны, вклад исходных данных (результатов измерений), при которых стоят близкие к нулю коэффициенты, в конечный результат вычислений столь ничтожен, что эти уравнения можно исключить из системы практически без потери точности.
На основании изложенного разработанная ранее модель была упрощена путем исключения из нее четырех уравнений, содержащих множитель э1п(а6), и переноса тригонометрических функций в левую часть, что исправляет допущенную ранее некорректность. Далее по тексту в ссылках на используемую математическую модель будет иметься в виду исправленная система уравнений.
Дальнейший анализ модели выявил еще одно важное свойство, благодаря которому можно существенно повысить достоверность расчетов. Оказалось, что в силу переопределенности системы и характера взаимосвязи результатов косвенных измерений, из нее можно исключить все уравнения, связанные с любым (но одним) из результатов измерений. При этом ранг матрицы коэффициентов системы не становится меньше числа неизвестных, что дает возможность ее приближенного решения.
Подставляя полученное решение «усеченной» таким образом системы уравнений в исключенные из нее уравнения, и решая их относительно исключенного параметра, можно получить его ожидаемое значение. Степень отклонения ожидаемого значения от измеренного, вполне может играть роль критерия достоверности данного параметра.
Если принять допущение, что в исходных данных (а это 30 результатов измерений) присутствует не более одной грубой ошибки, то процесс диагностирования и исключения такой ошибки легко автоматизируется и решается без человеческого участия.
Задача сводится к поочередному исключению из исходных данных каждого из результатов измерений и вычисления его ожидаемого значения. Если величина отклонения, какого-либо из результатов измерений от его ожидаемого значения больше некоторой заданной величины, то такой замер считается ошибочным, а в качестве наилучшего приближенного решения выбирается решение, полученное без учета этого результата измерения. Если величина максимального отклонения по всем измерениям не превышает заданной величины, то наилучшим считается решение, полученное с помощью полной системы (без исключения уравнений).
На рис.6 показаны среднеквадратичные отклонения исключаемых результатов измерений от их ожидаемых значений, вычисляемые по формуле
где к - номер исключаемого измерения;
а, Ь, х - коэффициенты системы [А]{Х}={В}, составленной из уравнений, исключенных из модели;
т, п - размерность матрицы [А]. Правые части этих уравнений содержат значение исключенного измерения, а вектор {X} - решение «усеченной» системы уравнений, т.е. вычисленные координаты центров расточек.
Из рис. 6 видно, что 51 - отклонение первого параметра от ожидаемого значения, вычисленного по приведенной выше методике, максимально, что вполне соответствует тестовому набору исходных данных, где все замеры, кроме первого, соответствуют идеально изготовленной конструкции, а параметр № 1 (это результат измерения Я1 согласно рис. 3) отклоняется на 0,1 мм от номинального размера (600 мм). Наличие незначительных отклонений остальных параметров объясняются влиянием первого параметра, присутствующего в уравнениях при определении этих отклонений.
Рис. 6. Среднеквадратичные отклонения результатов измерений от их ожидаемых значений: ось абсцисс - порядковый номер исключаемого измерения, ось ординат -
величина отклонения (мм)
При наличии более одной грубой ошибки в исходных данных, их диагностирование не так очевидно, особенно если величина ошибок одного порядка. Тем не менее, это возможно, но не в автоматизированном, а в ручном режиме путем анализа визуализированных результатов расчетов.
На рис. 7 приведен пример визуализации результатов тестовых расчетов с ошибочным результатом измерения Ях, причем на рис.7а отображены результаты решения при исключении из модели измерения Ях, а на рис.7б - без исключения Ях. Доверительные интервалы и отклонения центров расточек отображаются с соответствующим масштабным коэффициентом. Область 1 - доверительные интервалы положения центров расточек осей сателлитов (определяются расчетной точностью измерений). Центры областей 2 и 3 указывают положения центров расточек под оси сателлитов в правой и левой щеках водила.
Рис. 7. Визуализация результатов тестовых расчетов с ошибочным измерением R1.
Использование найденных описанным выше способом (с учетом всех контролируемых заводом-изготовителем точностных характеристик водила) координат расточек под оси сателлитов в водиле для расчета параметров АХк к+1 позволяет повысить
достоверность определения искомых комплексных погрешностей изготовления и сборки сателлитных узлов в соответствии с (1), что имеет важное практическое значение.
Вместе с тем, имеющиеся экспериментальные и расчетные данные [5, 8-10] показывают, что для дальнейшего повышения корректности определения результативной точности зубчатых зацеплений и ее влияния на нагруженность и виброактивность планетарного редуктора необходимо учитывать не только входящие в (1) первичные погрешности изготовления его деталей и узлов, но и нагрузочно-скоростные характеристики конкретных режимов эксплуатации.
Последние влияют на фактические погрешности зубчатых зацеплений в процессе работы редуктора вследствие упругих деформаций его элементов под нагрузкой, а также смещений зубчатых колес в пределах зазоров и на жесткостях опор под действием усилий в зацеплениях, весовых нагрузок и центробежных сил.
Работы в указанном направлении ведутся.
Список литературы
1. Айрапетов Э.Л., Генкин М.Д. Статика планетарных механизмов. М.: Наука, 1976. 264 с.
2. Айрапетов Э.Л., Генкин М.Д. Динамика планетарных механизмов. М.: Наука, 1980. 256 с.
3. Леонтьев Ю.А., Ямпольский И.Д., Хомяков В.П., Леонтьев М.Ю. Опыт создания на ОАО «КТЗ» судовых редукторов большой мощности // Юбилейный сборник трудов Научно-исследовательского Центра ОАО «Калужский Турбинный Завод». Калуга: Манускрипт, 2002. С. 134-144.
4. Леонтьев М.Ю. Исследование статической нагруженности мощных судовых планетарных редукторов: дисс. ...канд. техн. наук. М.: ИМАШ РАН, 2001. 229 с.
5. Леонтьев М.Ю., Полунин И.В., Лысенков В.С. Анализ влияния точности изготовления редукторов ГТЗА на виброактивность // Материалы V научно-технической конференции «Взгляд в будущее». СПб.: ФГУП «ЦКБ МТ «Рубин», 2007. С. 319-328.
6. Леонтьев М.Ю., Насонов Д.А. Об алгоритме численного решения одной из некорректно заданных задач, представляющей практический интерес // Материалы V Общероссийской научно-практической конференции с международным участием «Актуальные вопросы современной науки и образования». Вып. 2. Красноярск: Научно-инновационный центр, 2010. С. 280-287.
7. Дьяконов В.П. Mathematica 4 с пакетами расширений. М.: «Нолидж», 2000. 608 с.
8. Айрапетов Э.Л., Апархов В.И., Бедный И.А., Леонтьев М.Ю. Расчетные исследования малонагруженных планетарных передач // Теория и практика зубчатых передач: труды международной конференции. Ижевск: Изд-во ИжГТУ, 1996. С. 141-146.
9. Айрапетов Э.Л., Апархов В.И., Бедный И.А., Леонтьев М.Ю. Влияние весовой нагрузки зацеплений на динамику крупногабаритных планетарных передач // Ргосееёт§8 of the 5-th International Conference «Dynamics of machine aggregates». Gabcikovo, Slovak Republic, 2000.
10. Леонтьев М.Ю., Насонов Д.А. К вопросу о влиянии погрешностей изготовления и деформаций сателлитного узла на динамику планетарного редуктора // Материалы всероссийской научно-технической конференции «Наукоемкие технологии в приборо- и машиностроении и развитие инновационной деятельности в ВУЗе». Т. 1. М.: Изд. МГТУ им. Н.Э. Баумана, 2008. С. 76-81.
SCIENTIFIC PERIODICAL OF THE RAIJMAN MS TU
SCIENCE and EDUCATION
EL № FS77 - 48211. №0421200025. ISSN 1994-040S
electronic scientific and technical journal
On increasing the accuracy of combined errors estimation
in satellite planetary gear units
# 10, October 2012
DOI: 10.7463/1012.0479419
Nasonov D.A., Leont'ev M.Yu., Shubin A.A.
Russia, Bauman Moscow State Technical University, Kaluga Branch
nasonovda@yandex.ru newell-kaluga@mail.ru shubin55@mail.ru
The authors consider practical ways to control combined errors in satellite planetary gear units used by manufacturing plants. The authors adjusted the mathematical model for determination of coordinates of bores for satellites in the planetary gear unit according to the results of indirect measurements, developed and implemented a software algorithm to identify and reject invalid input data. Directions for further research were mapped out.
Publications with keywords:mathematical modeling, planetary gear, ill-condition tasks Publications with words:mathematical modeling, planetary gear, ill-condition tasks
References
1. Airapetov E.L., Genkin M.D. Statika planetarnykh mekhanizmov [Statics of planetary mechanisms]. Moscow, Nauka, 1976. 264 p.
2. Airapetov E.L., Genkin M.D. Dinamika planetarnykh mekhanizmov [Dynamics of of planetary mechanisms]. Moscow, Nauka, 1980. 256 p.
3. Leont'ev Iu.A., Iampol'skii I.D., Khomiakov V.P., Leont'ev M.Iu. Opyt sozdaniia na OAO «KTZ» sudovykh reduktorov bol'shoi moshchnosti [Experience of creating at JSC "KTZ" marine reduction gears of high power]. Iubileinyi sbornik trudov Nauchno-issledovatel'skogo Tsentra OAO «Kaluzhskii Turbinnyi Zavod» [Jubilee collection of the works of the Scientific Research Center of JSC "Kaluga Turbine Plant"]. Kaluga, Manuskript, 2002, pp. 134-144.
4. Leont'ev M.Iu. Issledovanie staticheskoi nagruzhennosti moshchnykh sudovykh planetarnykh reduktorov. Kand. diss. [The study of the static loading of powerful marine planetary reduction gears. Cand. diss.]. Moscow, Institute of engineering science of RAN, 2001. 229 p.
5. Leont'ev M.Iu., Polunin I.V., Lysenkov V.S. Analiz vliianiia tochnosti izgotovleniia reduktorov GTZA na vibroaktivnost' [Analysis of influence of precision of manufacturing
reducers of the turbine-geared propulsion unit at vibration activity]. Materialy 5 nauchno-tekhnicheskoi konferentsii «Vzgliad v budushchee» [Proc. of the 5th scientific-technical conference «Looking into the Future»].St. Petersburg, TsKB MT «Rubin», 2007, pp. 319-328.
6. Leont'ev M.Iu., Nasonov D.A. Ob algoritme chislennogo resheniia odnoi iz nekorrektno zadannykh zadach, predstavliaiushchei prakticheskii interes [About the algorithm of numerical solution of one of the incorrectly specified tasks, representing practical interest]. Materialy 5 Obshcherossiiskoi nauchno-prakticheskoi konferentsii s mezhdunarodnym uchastiem «Aktual'nye voprosy sovremennoi nauki i obrazovaniia» [Proc. of the 5th all-Russian scientific-practical conference with international participation «Current issues of modern science and education».]. Iss. 2. Krasnoiarsk, Nauchno-innovatsionnyi tsentr, 2010, pp. 280-287.
7. D'iakonov V.P. Mathematica 4 spaketami rasshirenii [Mathematica 4 with bump packs]. Moscow, Nolidzh, 2000. 608 p.
8. Airapetov E.L., Aparkhov V.I., Bednyi I.A., Leont'ev M.Iu. Raschetnye issledovaniia malonagruzhennykh planetarnykh peredach [Computational studies of lightly loaded planetary gears]. Teoriia ipraktika zubchatykhperedach: trudy mezhdunarodnoi konferentsii [Theory and practice of gears: Proc. of the International Conference]. Izhevsk, IzhGTU Publ., 1996, pp. 141146.
9. Airapetov E.L., Aparkhov V.I., Bednyi I.A., Leont'ev M.Iu. Vliianie vesovoi nagruzki zatseplenii na dinamiku krupnogabaritnykh planetarnykh peredach [Effect of load bearing of catching on the dynamics of large planetary gears]. Proc. of the 5-th International Conference «Dynamics of machine aggregates». Gabcikovo, Slovak Republic, 2000.
10. Leont'ev M.Iu., Nasonov D.A. K voprosu o vliianii pogreshnostei izgotovleniia i deformatsii satellitnogo uzla na dinamiku planetarnogo reduktora [To the question about the influence of the errors of manufacturing and deformations of satellite unit at the dynamics of planetary reduction gear]. Materialy vserossiiskoi nauchno-tekhnicheskoi konferentsii «Naukoemkie tekhnologii v priboro- i mashinostroenii i razvitie innovatsionnoi deiatel'nosti v VUZe» [Proc. of all-Russian scientific-technical conference «Science-Intensive technologies in instrument - and mechanical engineering and the development of the innovative activities in higher educational institutions»]. Vol. 1. Moscow, Bauman MSTU Publ., 2008. pp. 76-81.
