Разработка научных методов расчета нестационарного взаимодействия тонкостенных элементов с жесткими односторонними связями и математических моделей волновых передач тема диссертации и автореферата по ВАК РФ 01.02.06, доктор технических наук Люминарский, Игорь Евгеньевич

• Люминарский, Игорь Евгеньевич
• доктор технических наук
• 2009

Оглавление диссертации доктор технических наук Люминарский, Игорь Евгеньевич

Введение

Глава 1. Методы расчета упругих систем с односторонними связями

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

1.1.1. Шаговый метод последовательных нагружений

1.1.2. Метод попыток

1.1.3. Метод продолжения решения по параметру

1.1.4. Мытоды, основанные на минимизации функционала энергии

1.1.4.1. Метод приведенного градиента

1.1.4.2. Метод проекционного градиента

1.1.4.3. Метод штрафа

1.1.4.4. Классические Лагранжевы методы

1.1.5. Статический расчет силового взаимодействия звеньев волновой зубчатой передачи

1.2. О положительной определенности матриц упругих систем с односторонними связями

1.3. Метод введения восстанавливающих сил

1.4. Поперечный удар шара по балке

1.4.1. Сравнение методов численного решения интегрального уравнения, описывающего поперечный удар шара по балке

1.4.2. Динамический расчет взаимодействия шара и балки при поперечном ударе

1.5. Динамический расчет линейных упругих систем с упругими односторонними связями

1.5.1. Расчет систем с неинерционными односторонними связями

1.5.2. Расчет систем с инерционными односторонними связями

1.6. Выводы

Глава 2. Изгибные колебания балок с односторонними связями

2.1. Колебания тонких балок

2.1.1. Математическая модель колебаний

2.1.2. Определение собственных частот и форм изгибных колебаний балки

2.1.3. Оценка достоверности результатов расчета

2.1.4. Численные и экспериментальные исследования

2.1.4.1. Свободные колебания

2.1.4.2. Вынужденные колебания

2.2. Влияние деформации поперечного сдвига и инерции вращения на колебания балок с односторонними связями большой жесткости

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

2.4. Выводы

Глава 3. Изгибные колебания прямоугольных пластин с односторонними связями 3.1. Колебания тонких пластин

3.1.1. Математическая модель колебаний

3.1.2. Определение собственных частот и форм колебаний тонкой прямоугольной пластины

3.1.2.1. Метод последовательных нагружений

3.1.2.2. Метод Ритца

3.2.2.3. Метод конечных элементов

3.2.2.4. Метод конечных разностей

3.1.3. Сравнение методов определения собственных форм и частот прямоугольной пластины

3.1.4. Численные исследования сходимости получаемого решения

3.1.5. Расчет контактной системы

3.1.6. Вынужденные колебания пластин с односторонней связью

3.1.7. Экспериментальные исследования колебаний пластины с односторонней связью 180 3.2. Влияние деформацию сдвига и инерцию вращения на колебания пластины с односторонними связями большой жесткости

3.2.1. Математическая модель

3.2.2. Определение собственных частот и форм колебаний

3.2.2.1. Метод последовательных нагружений

3.2.2.2. Метод Ритца

3.2.3. Численные исследования колебаний пластины типа

С.П. Тимошенко с односторонней связью

Введение диссертации (часть автореферата) на тему «Разработка научных методов расчета нестационарного взаимодействия тонкостенных элементов с жесткими односторонними связями и математических моделей волновых передач»

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

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

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

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

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

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

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

Одним из первых расчеты конструкций с односторонними связями в общей постановке выполнил М. Геллер [184]. Он исследовал статически неопределимые системы, у которых некоторые лишние связи являлись односторонними. Геллер вполне правильно видел основную трудность расчета систем с односторонними связями в том, чтобы отыскать среди возможных систем действительную, отвечающую заданной нагрузке. Критерии правильности нахождения действительной системы, которые ввел Геллер, следующие: 1) усилия во всех односторонних связях, найденные из уравнений упругости, должны быть не отрицательными; 2) зазоры по направлению односторонних связей должны быть положительными.

Однако Геллер не решил основную проблему систем с односторонними связями — нахождение действительной системы. Указания Геллера сводятся к тому, что нужно путем подбора найти систему, которая удовлетворяет указанным выше двум условиям. Эта система и будет действительной.

Развитию методов статического расчета упругих систем с односторонними связями посвящены работы И.М. Рабиновича [142-145]. Рабинович указывает на два метода поиска действительных рабочих систем. Согласно первому методу поиск начинается с приближенного задания рабочей системы. Затем производится проверка усилий и зазоров в односторонних связях. Согласно ее результатам рабочая система корректируется, т.е. включаются или выключаются некоторые односторонние связи. Описанный процесс повторяется до тех пор, пока очередная проверка покажет, что зазоры во всех односторонних неработающих связях и усилия во всех работающих связях будут неотрицательными.

Описанный метод можно назвать методом последовательных попыток. Поскольку в указанном методе поиск действительной рабочей системы не носит целенаправленного характера, то рост числа односторонних связей в системе приводит к резкому возрастанию количества вариантов «ложных» рабочих систем необходимых для пересчета, что значительно увеличивает время расчета. В [142] показано, что число таких вариантов для системы, в состав которой входит L односторонних связей, будет равно 2Ь. Следует отметить также, что описанный метод может не сходится. При численной реализации этого метода есть возможность зацикливания расчета.

Идея второго метода предложенного И. М. Рабиновичем, состоит в том, чтобы рабочую систему, найденную для какой-либо нагрузки О, использовать для определения рабочей системы, отвечающей нагрузке Р. При этой методике разыскиваются последовательно сменяющие одна другую действительные рабочие системы для суммарной нагрузки F = (l - v)0 + vP в процессе постепенного возрастания множителя v от нуля до единицы. Увеличение коэффициента v производится по шагам. Каждый шаг заканчивается сменой рабочей системы, т.е. включением в работу или выключением из работы одной или нескольких односторонних связей. Поэтому этот метод можно назвать шаговым методом последовательных нагружений. И.М. Рабинович рассматривает этот метод на примере шарнирно опертой балки с тремя односторонними связями. К недостатку метода относится тот факт, что приходится перебирать все промежуточные рабочие системы, не пропуская ни одной, которые образуются в процессе увеличения коэффициента v от нуля до единицы. Это обстоятельство может приводить к значительному увеличению машинного времени, в случае большого числа односторонних связей.

Новый этап в развитии теории расчета статически неопределимых систем с односторонними связями начался со времени использования методов математического программирования.

В работах [45,133] А.В. Перельмутер и В.Н. Гордеев рассматривают задачу расчета упругой системы с односторонними связями как задачу квадратичного программирования. Вводится функционал потенциальной энергии вида

J(R) - ^RTCR + (р-v;T + д0т )R , (1) np где R — вектор-столбец реакций в односторонних связях (R — векторстрока реакций); С — матрица податливости упругой системы, освобожден* ной от односторонних связей; w — вектор перемещений свободной упругой системы по направлениям односторонних связей от воздействия внешней нагрузки с единичной интенсивностью; р* — интенсивность внешней нагрузки; А0— вектор зазоров между односторонними связями и упругой недеформированной системой.

Действительная система удовлетворяет минимуму функционала (1) при ограничениях вида

- 10

Rj> 0 , j = \,.,L , (2) где Z - число односторонних связей.

Таким образом, задача расчета упругой системы сводится к задаче квадратичного программирования вида min{j(R), R > 0}. (3)

Применив условия Куна—Таккера [73] для задачи (3), Перельмутер приходит к следующей системе равенств и неравенств:

A = CR + Jpw;+A0;

А,. >0, i = l,.,L;

Д,. >0, i = 1,.,Z; (4)

ДД=0, i = 1,.,Z.

Перельмутер в [133] предлагает решать систему (4) методом последовательных нагружений. Суть этого метода заключается в следующем. Пусть состояние системы полностью определено при р = 0, тогда состояние системы при р = р* предлагается получать путем последовательного увеличения р от нуля до р*. Таким образом, система (4) заменяется задачей вида min{F = -р,р* - р > о} A = CR + ^w;+A0;

А, > 0, i = l,.,L; (5)

Rt> 0, i = l,.,L; АД=0, i = \,.,L.

Для решения (5) в [133] предлагается использовать алгоритм Вульфа [115], решающий задачу квадратичного программирования.

Метод, предложенный Перельмутером в [133], имеет тот же недостаток, что и метод Рабиновича. Для того чтобы получить рабочую систему при р = р\ необходимо перебрать все рабочие системы в процессе увеличения р от нуля до р*.

Алфутовым Н.А. и Клениковым С.С. [7] разработан шаговый метод расчета упругих систем с односторонними связями, где вся «история нагружения» системы представляется последовательной совокупностью отдельных этапов ее нагружения. При этом в качестве основных неизвестных задачи на каждом этапе нагружения системы принимаются не сами силы и зазоры между линейно-упругой системы и односторонними связями, а их приращения. Для определения границ смены этапов нагружения авторами разработана процедура линейного нормирования этапных решений.

Разработке алгоритмов статического расчета упругих систем с односторонними связями при действии возрастающей нагрузки посвящены работы [9,140,170,169]. В [9] А.Х. Астрахан для нахождения действительной рабочей системы на каждом шаге нагружения решает задачу квадратичного программирования методом сопряженных градиентов [24]. В работе [140] Л.П. Портаев рассматривает поставленную задачу как задачу линейного программирования.

В работах [112,113] Э.Н. Кузнецов положил в основу исследования свойств систем, содержащих односторонние связи, систему неравенств, которая выражает условия двухсторонних связей и ограничивающие условия, налагаемые на односторонние связи. Совокупность этих равенств и неравенств дает полную информацию о статико-кинематических свойствах системы.

Практическое применение и развитие методов расчета упругих систем с односторонними связями получило в работах [77,83-96,1,117, 118,151,165,168]. В этих работах рассматриваются задачи статического расчета силового взаимодействия элементов волновых передач. Для определения контактных усилий между элементами передачи использовались итерационный [77,86,87,165,168] и шаговый методы [1,7,83,84,8891,94-96,117,118,152].

Отличительной особенностью задач, рассматриваемых в работах [1,90,91], является большое количество односторонних связей. Так, при расчете сдвоенной волновой передачи [1,90] число односторонних связей достигает 500. В связи с этим разработанные алгоритмы предусматривали возможность автоматизации уточнения действительной рабочей системы.

На рис. 1 представлен алгоритм, используемый в работах [1,90], который определяет действительную рабочую систему, вектора сил R и зазоров А между упругой системой нагруженной внешними силами и односторонними связями.

Главным недостатком шаговых методов последовательных нагружений является тот факт, что решение задачи начинается с известного начального положения (положения, при котором известны силы и зазоры между упругой системой и односторонними связями). Таким состоянием, как правило, является положение, в котором реакции в односторонних связях равны нулю. Поэтому в случае расчета подобных систем при фиксированном параметре нагрузки р необходимо построить весь спектр решений от р=0 до р~р*. В задачах с большим числом односторонних связей это приводит к значительному увеличению числа этапов решения, и следовательно, к большим затратам машинного времени. В значительной степени этот недостаток проявляется в том случае, когда необходимо производить множество расчетов, в которых зоны контакта между упругой системой и односторонними связями меняются незначительно. Такие расчеты производятся, например, при оптимизации конструкции. В этом случае выгодно применять итерационные методы поиска рабочей системы, т.к. они начинаются с приближенного задания предварительной зоны контакта, которая потом уточняется по определенному алгоритму. Однако недостатком итерационных методов является возможность зацикливания при расчете. к-к+l нет

Определение реакций и зазоров между упругой системой и односторонними связями к концу решения.

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

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

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

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

Широкое распространение при решении задач теории упругости с односторонними краевыми условиями получили вариационные методы [191,185,41]. Эти методы основаны на сведении исходной задачи к эллиптическим вариационным неравенствам или к задаче минимизации некоторого функционала, отвечающим некоторым выпуклым односторонним ограничениям. Этими методами решаются многие задачи пластичности [188,141].

К. Байокки и А. Капело в [20] рассматривают следующие три задачи.

Задача 1

Найти и0 > 0, такое, что

А(и0)-р>0 и (А(и0)-р,и0) = 0.

Задача 2 Найти u0<e К, такое, что F{u0)<F{v), для Vv е К, где F(v) = ^(A(u),v)-(p,v)

Задача 3 Найти и0е К, такое, что {a(u0),u0-v)<(p,uo-v), для VveK.

При формулировке задач 1-3 введены следующие обозначения: А : Н —> Н' линейный непрерывный оператор, где Н' - пространство, сопряженное к гильбертову пространству Н \ К а II - непустое замкнутое выпуклое множество; (a(u),v) = Jw(x)- v{x)dx - скалярное произведение двух D функций, где D — область в n-мерном евклидовом пространстве Rn.

Сведение задач с односторонними ограничениями (задача 1) к вариационным неравенствам основано на следующем утверждении, доказанном в [20].

Утверждение

Пусть оператор^ положительный и симметричный,- множество неотрицательных элементов, тогда задачи 1-3 эквивалентны и имеют единственное решение.

Сущность вариационных методов рассмотрим на примере изгиба балки под действием распределенной нагрузки р{х). На некотором расстоянии а от балки находится жесткое препятствие. При достаточно большой нагрузке р{х) будет существовать область контакта балки с препятствием. Для случая жестко заделанной по краям балки такая задача записывается в следующем виде: 0 <*</; ах и(х)= 0,^ = 0, х = 0,/; (7) ах dx = 0, где и — и — a. и(х) >а, 0 < х < /, где D — изгибная жесткость балки; 1 — длина балки. dAu u)=D——- симметричный и положительный, ах то задача (7) эквивалентна задачами 5-6.

Задача 5 Найти и > 0, такое, что

D^M.p{x)> о и \[*Щй dxA V J {{ dx

Задача 6 Найти и0 е К, такое, что F(u0)<F(v), для Vv е К, 1 ' dAv где F(v) = -D^—vdx - Jp(x)v(x)dx; А о

К- множество достаточно гладких функций таких, что К = jv/v(x)>a,0<x< /,v(x) = ^ = 0,х = 0,/| .

С помощью интегрирования по частям функционал F(v), используемый в задачи 6, преобразуется к виду

F(v) — Z)J ^ dx - J p{x)v(x)dx. 2 о \dx )

Существует несколько подходов к решению задачи 6. Первый подход основан на непосредственном применении методов бесконечномерной оптимизации [124]. Для примера опишем решение задачи 6 методом штрафа [116,149], приведенным в работах [32,33]. Метод основан на замене исходного функционала F(v) функционалом вида

F£(v) = F(v) + ijpE(vXv-«)2^, о где коэффициент штрафа взят в виде

- 17-[в ,v<a.

Решение поставленной задачи, которое обозначается иЕ(х), дает приближенное решение исходной задачи. Уравнение Эйлера, соответствующее задаче минимизации функционала Fe(v), имеет вид да dx D

Уравнение (8) описывает изгиб балки в условиях, когда на расстоянии а от нее находится упругая среда с коэффициентом упру2 гости, пропорциональным S . Случай е —> О соответствует жесткому препятствию.

Для решения нелинейной краевой задачи (8) использовался итерационный процесс вида + 0<х<1, (9) dx D dx

Аналогичный подход используется в работах [97,98,189] для решения задачи о контакте упругой пластины с препятствием. В работе [97] доказывается сходимость предложенных итерационных методов решения вариационного неравенства с использованием оператора штрафа.

Второй подход к решению вариационных задач заключается в замене бесконечномерных функционалов конечномерными функционалами с дальнейшим применением методов конечномерной оптимизации. Так, например, задача 2 заменяется задачей 7.

Задача 7

Найти uh g Kh, такое, что Для где uh,vh — сеточные функции; Fh{uh\Kh — сеточные аналоги функционала F(u) и выпуклого множества К с параметром дискретизации (например, шагом сетки) h.

В зависимости от природы изучаемой задачи этот подход может принимать различные формы: методы конечных разностей, методы конечных элементов, методы Ритца—Галеркина.

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

Третий подход к решению задач с односторонними ограничениями заключается в представлении решения в виде интегрального уравнения. Так, например, решение задачи (7) представляется в виде и(х)~ Ji?(x)G!(jc,x]yx + M;)(x), (Ю) о где ир{х) — функция перемещений свободной балки от внешней нагрузки р(х); Z?(x) - функция реакций от одностороннего ограничения; Cj(x,x) -функция Грина.

Интегральное уравнение (10) должно удовлетворять условию и(х)>а, 0 <х<1. (11)

В некоторых простейших задачах [53] можно найти точное аналитическое решение задачи (10) — (11). В большинстве же задач решение интегрального уравнения (10) с ограничением (11) ведется численными методами. Опишем один из них.

Функция реакции R(x) представляется в виде ряда

ЯМ = 2>,5(х,ХД (12) i где ) — функция Дирака, Xj =—-—j (в точках с координатами х=0 и

7 L +1 х-а реакция основания принималась равной нулю). Записав уравнение (10) для точек с координатами , задача (10) - (11) сводится к системе неравенств вида

А = u(x)-a = Gr + up -а, А, >0,/ = 1,.Х, ri>0,i = l,.L, (13)

A,r = 0,i = 1 ,.L.

Неравенства (13) аналогичны неравенствам (5) , которые записаны для системы с односторонними связями. Таким образом, задачи деформирования упругих систем, ограниченных жесткой преградой и односторонними связями, описываются одинаковыми по структуре неравенствами.

Описанный метод был применен в работах [1,90,91] при расчетах упругого взаимодействия элементов волновой передачи. К недостаткам описанного метода можно отнести необходимость вычисления функций Грина (j(.r(,xJ для узловых точек сетки, покрывающей возможную область контакта взаимодействующих тел. Это обстоятельство может привести к значительным затратам машинного времени. Однако следует заметить, что для различных методов вычисления функций Грина (метод конечных разностей, методы Ритца-Галеркина, метод конечных элементов и т.п.) можно найти способы экономии машинного времени. Так, например, большинство методов определения функций Грина сводятся к решению систем линейных алгебраических уравнений вида

AGj = 5J , j = \,.,L, (14) где 8/ — дельта Кронеккера.

Поскольку в процессе вычислений векторов GJ матрица А не меняется, то ее приведение к треугольному виду осуществляется один раз, что значительно сокращает затраты машинного времени.

Начало динамического расчета систем с односторонними связями было положено М.В. Остроградским. Он впервые получил уравнения движения систем с односторонними ограничениями [130], которые сейчас в аналитической механике называются системами с неудерживающими связями. М.В. Остроградский указал алгоритм отбора тех связей, которые должны быть приняты во внимание в данный момент времени при интегрировании уравнений движения. Алгоритм М.В. Остроградского основан на слежении за знаками множителей Лагранжа — как только лагранжев множитель в процессе движения меняет знак, то соответствующая связь должна быть отброшена. Поскольку лагранжев множитель является реакцией связи, то условием схода со связи является смена знака реакции в односторонней связи.

Вопросам движения твердых тел с односторонними связями посвящены работы [63,158,74,57-59].

Строгое определение движения системы с односторонней связью дано Дерябиным М.В. в работе [59]. Пусть дана натуральная механическая система с конфигурированным пространством R" =х, на которую наложена голономная идеальная односторонняя связь f(x) > 0. Под движением с односторонней связью понимают отображение х: [0, т] —» R", которое удовлетворяет уравнениям dUdTл dt дху d (дТ\ дТ = F + f(x) = 0, (15) дх ох dt jr = /М>0, (16) ydxj дх где Т — кинетические энергия системы, X — неотрицательный множитель.

В работе [74] А.П. Иванов указывает типы движения систем с односторонними связями. Основные из них следующие: перелеты; "скользящий" режим; сход со связи; удар; касание ограничения.

Перелеты характерны тем, что существует зазор между системой и односторонней связью (Л>0), т.е. связь не нагружена (R=0). При перелетах наличие односторонней связи никак не сказывается на движении системы.

Скользящий режим движения соответствует нагруженному состоянию односторонней связи (R > 0; А = 0). На этом режиме в уравнения движения необходимо добавить реакцию односторонней связи.

Сход со связи характеризуется уменьшением реакции в односторонней связи до нуля (i? = 0,A = 0) с дальнейшим увеличением зазора между системой и связью. Уменьшение реакции до нуля может происходить плавно и скачком.

Классические условия схода со связи даны в работе [59]. Система с односторонней связью движется по связи, если в процессе движения выполняется условие f(x) = 0. Система сходит со связи в момент времени т, если до этого момента времени система двигалась по связи и для любого достаточно малого е>0 найдется такой момент /Ее(т,т + в], что /(х(тс))> 0. Данное определение схода со связи не представляет конструктивного условия схода со связи.

Фаза удара характеризуется переходом от фазы полета к фазе «скользящего» режима. В начальный момент удара зазор между системой и связью и реакция связи равны нулю. Условия удара в момент времени t = t* следующие:

-22В этом случае, для того чтобы определить движение системы при t>t* необходимо использовать дополнительные уравнения, описывающие природу удара.

Литература, посвященная соударению твердых тел, весьма обширная. Имеются несколько широко известных монографий непосредственно по теории соударения твердых тел [3,5,6,23,29,51,52,60,80,81,131,146].

Первым определить то, как протекает соударение шаров во времени, т.е. найти силу соударения как функцию времени, удалось Г. Герцу [25,60,187]. Его решение основано на квазистатической теории соударения упругих тел, согласно которой соударение считается продолжающимся сколь угодно долго по отношению к периоду первой основной формы колебаний каждого из соударяющихся тел.

Теория Г. Герца относится не только к соударению шаров, но и к случаю прямого центрального удара двух тел, ограниченных в окрестности точки контакта поверхностями второго порядка.

Учет и местной деформации и возникающих при ударе волн для случая продольного соударения стержней со сферическими концами впервые выполнены Дж. Сирсом [194].

Идея Сирса совместного учета местных деформаций и упругих колебаний в соударяющихся телах получила широкое распространение. Применительно к поперечному удару шара по балке она была использована С.П. Тимошенко [156]. Местные деформации С.П. Тимошенко определяет по квазистатической теории Герца, т.е. a = KhrK где (18)

R — сила соударения, Кн — некоторая постоянная, <2— местная деформация при ударе.

Обобщая теорию С.П. Тимошенко на случай одномерной или двумерной упругих систем, B.JI. Бидерман записывает разрешающие уравнение движения относительно силы соударения. Разработке численных методов решения этого уравнения посвящены работы [11,12,26,128,147, 154,159,192,193,44]. Актуальным этот вопрос остается и сегодня.

Вопросам аналитического решения задач об ударе шара по балке и пластине посвящены работы С.А. Зегжды [64-70 ].

Ударному взаимодействию объекта с одномерной средой посвящены работы [119,120]. Используя вариационный принцип Гамильтона, в [119] получены условия контакта упругой среды и сосредоточенного объекта, соотношения справедливые в момент начала и конца контакта. На основании полученных соотношений рассмотрены аналитические решения следующих задач: центральный удар материальной точки по ограниченной покоящейся струне; косой удар материальной точки по бесконечной струне на упругом основании.

Движению твердых тел по шероховатой вибрирующей поверхности посвящена работа Блехмана И.И. и Джанилидзе Г.Ю [31]. В этой работе рассматривается движение частиц по гармонически колеблющейся плоскости без учета и с учетом отрыва частицы.

Ударному взаимодействию твердого тела с ограничителем, которое сопровождается повторными затухающими соударениями посвящена работа Нагаева Р.Ф. [126].

Вынужденным колебаниям упругих тел с одной или несколькими ограничителями посвящена работа [10] Бабицкого В.И. и Крупенина B.JI. В работе рассматриваются приближенные и точные аналитические методы, основанные на методе гармонической линеаризации [34] и на припасовывании решений [131]. Трудоемкость отыскания решений точными методами ограничивает возможность их применения для анализа сложных систем большой размерности. Такие системы исследуются приближенными методами.

Точная постановка задачи об ударе упругих тел рассматривается в работах [105-110,181,182]. В этих работах данная задача рассматривается как смешанная задача динамической теории упругости, которая решается численными методами. Подходы, описанные в этих работах, учитывают переменность площадки контакта в процессе соударения, возможность пластических деформаций в зоне контакта. При переменной площадке контакта на начальном этапе удара скорость расширения области контакта выше скорости продольной волны в упругом теле. В работах А.Г. Горшкова и Д.В. Тарлаковского [48-50] эти особенности учитываются при ударе затупленным твердым телом о поверхность упругого полупространства.

Следует отметить следующее. При ударе шара по балке, пластине, оболочке время начального этапа взаимодействия (этапа на котором скорость расширения области контакта выше скорости продольной волны в упругом теле) на несколько порядков меньше, чем время удара [70]. Однако его учет значительно затрудняет расчеты.

Расчету динамических задач, относящихся к проблеме односторонних связей в деформированных механических системах, посвящены работы [13-16,55,56,103,104,171-179,72,111,122,180].

В работах [13-16,55,56] рассмотрены вопросы вынужденных периодических колебаний стержневых, пластинчатых и оболочечных элементов конструкций с односторонними связями и упругим основанием.

При расчете упругих тонкостенных элементов в этих работах используются уравнения движения, основанные на гипотезах Кирхгофа-Лява. Дифференциальное уравнение в частных производных, описывающее движение упругого тела методом Бубнова-Галеркина сводилось к системе обыкновенных дифференциальных уравнений. В работах [13-16,55,56] рассматривались колебания упругих тел с простейшими граничными условиями, соответствующими шарнирному опиранию. Такие граничные условия позволяют решение представлять в виде тригонометрического ряда v{x,y,t)= ^W^^ ' (19) m=1 n=l \ a J \ & ) где a,b- размеры в плане.

Для учета влияния односторонних связей в дифференциальное уравнение вводилось слагаемое iб(х,у,х,y)ks(1 - sign(w(x,y,t)-A0))|w(x,y,t) - A01, (20) где x,y — координаты срединной поверхности упругого тонкостенного тела; w(x,y,i) - функция прогибов; х,у— координата односторонней связи; — коэффициент жесткости односторонней связи; Д0 — зазор между связью и недеформированной пластиной; б(х,>',х,>') — функция Дирака.

Физический смысл слагаемого (20) — сосредоточенная реакция односторонней связи, которая является линейной неинерционной пружиной.

Решение системы нелинейных уравнений относительно коэффициентов Amn(t) базируется на совместном использовании алгоритма продолжения решения по параметру и метода Ньютона-Канторовича. Согласно этому подходу параметр нагружения системы р последовательно увеличивался с шагом Ар. При этом разрешающая система дифференциальных уравнений записывается относительно приращений коэффициентов разложения в тригонометрический ряд AAmn(t). Решение системы дифференциальных уравнений ведется методом конечных разностей. Значения коэффициентов разложения, соответствующих интенсивности нагружения рк+] = рк + Арк, находятся по формуле

Ak::(t) = Akmn(t)+AAlXt)- (21) /t+1 При каждом вновь наиденном значении интенсивности р анализируется функция прогибов упругой системы и устанавливаются моменты включения в работу и выключения из работы односторонних связей. За начальное состояние принимается положение упругой системы с нулевым значением амплитуды р внешнего воздействия, что соответствует

Отметим следующие недостатки метода, предложенного в [16 ]:

1) Точность полученного решения зависит от двух параметров дискретизации: шага последовательного нагружения Ар и шага по времени

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

2) На каждом шаге нагружения необходимо производить суммирования рядов (19) для анализа зазоров. Общее количество этих суммирований в процессе решения задачи определяется формулой

K = KrKuL, (22) где L — число односторонних связей; Кр— количество шагов нагружения; К&1 — количество интервалов, на которые разбивается период движения.

3) Методика описана на примере расчета пластин и оболочек с простейшими граничными условиями, позволяющими представлять решение в виде ряда (19).

Необходимость получения более точного решения приводит к уменьшению шагов Ар и At, что в свою очередь может значительно увеличить число суммирований по формуле (19). Это обстоятельство может затруднить расчет по предложенной в работе [16] методике в связи с большими затратами машинного времени.

При заданной точности решения задачи количество суммирований К в первую очередь зависит от коэффициента к = max k' Л кК; = 1,.Х, где к'ос жесткость односторонней связи с номером I; к' — жесткость упругой системы под связью с номером i. В работах Баженова В.А. и др. [16,55,56] рассматриваются системы, в которых коэффициент к меняется в пределах от 0,1 до 3. Если же жесткость односторонних связей на несколько порядков выше жесткости упругой системы (к = 102 -Ю5), то расчет по методике, предложенной Баженовым В.А., затруднен из-за больших затрат машинного времени и большой размерности разрешающей системы линейных уравнений, которая зависит от шага At.

Значительное увеличение времени расчета с ростом коэффициента к также связано с необходимостью учета большого числа членов ряда (19).

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

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

Влияние конечности скорости волновых процессов в материале на напряженно-деформированное состояние упругих элементов позволяет учесть теория С.П. Тимошенко [156]. Дифференциальные уравнения колебаний балочных систем при этом содержат члены, учитывающие инерцию вращения и искажение поперечных сечений. Подобные уравнения для пластин и оболочек можно найти в работах [54,148]. Решению уравнений колебаний балок, пластин и оболочек, учитывающих инерцию вращения сечений и искажение поперечных сечений, посвящены работы [18,19,27, 28,8,46,47] и др.

При импульсных нагрузках расчет с использованием приближенной теории Кирхгофа-Лява может давать большую погрешность. В связи с этим в работе [103] используются уточненные уравнения движения упругих элементов, которые называются уравнениями типа С.П. Тимошенко. Искомые функции разлагаются в тригонометрические ряды Фурье, которые удовлетворяют условиям шарнирного опирания. Реакция между взаимодействующими балками представляется в виде ступенчатой функции координаты jc

R(x,t) = ±Rp(tiHixP-Х)~Н(Х' (23> р=1 где хр = рАх.

Входящие в выражение (23) неизвестные функции Rp(t) представляются также в виде ступенчатых функций * v - 0 - # (с, - 0L (24)

V=1 где tv = vAt, Rpv — неизвестные постоянные величины.

Представление реакции односторонней связи в виде рядов (23) и (24) позволяет на каждом шаге по времени разрешающую систему интегральных уравнений заменить системами алгебраических линейных уравнений относительно коэффициентов R . Полученные системы отражают совпадение нормальных перемещений взаимодействующих балок в ряде точек, расположенных равномерно вдоль длины в различные моменты времени. Из этих систем последовательно, с использованием итерационного уточнения, для значений v = \,2,.,т находятся коэффициенты разложения R .

Итерационное уточнение производилось для каждого момента времени tv следующим образом. Полученные из решения системы линейных уравнений величины R являются первым приближением. На второй итерации отрицательные значения R приравниваются нулю. На участках, где на предыдущем приближении были нулевые значения R , проверялись зазоры между взаимодействующими балками. Если на некоторых участках были отрицательные зазоры, то на них вводились контактные реакции. После этого переформировывалась система алгебраических уравнений и пересчитывались коэффициенты R . Третье приближение проводилось аналогично второму и т.д. Если величины контактных давлений Rpv на всех участках, отвечающие двум последним приближениям, совпадают, итерационный процесс считается законченным. С помощью описанного итерационного процесса определяются в момент времени tv = vA t реакции и зазоры между балками.

Следует отметить, что итерационный процесс, предложенный в [103], может не сходиться. При численной реализации такого процесса есть возможность зацикливания программы.

В работах Янютина Е.Г. и Ермоленко JI.B. [175,178] рассматривается нестационарное деформирование балки и цилиндрической оболочки на неинерционном упругом одностороннем основании. Распределенная нагрузка r{x,t)со стороны основания описывается функцией где - функция, описывающая изменение жесткости основания вдоль длины балки.

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

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

В работах [175,178] рассматриваются колебания шарнирно опертой балки на инерционном основании. Расчет этой задачи основан на совместном решение дифференциальных уравнений описывающих колебания балки типа С.П. Тимошенко и уравнений, моделирующих деформирование инерционного основания [150]. Методика решения этих уравнений аналогична методу, описанному в [103].

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

В работе [111] рассматривается задача динамического поведения многослойных оболочек с учетом особенностей контакта слоев. В основу решения положен стандартный метод характеристик [111], который базируется на замене дифференциальных соотношений конечноразностными уравнениями. В уравнениях учитывается местное сближение слоев в точках их контакта, связанное с шероховатостью материала. Это сближение описывается функцией вида X = (хд, где а — коэффициент пропорциональности; 8 — сближение слоев.

В работе Ковтуненко В.А. [99] рассматривается динамическая задача внедрения жесткого штампа в упругопластическую балку, шарнирно закрепленную по краям. Для описания поведения материала балки используется модель пластического течения. Решение поставленной задачи ведется итерационными методами, основанными на методе штрафа [116,42]. Итерационные схемы используют приближенную замену производных конечными разностями. Решение разрешающей системы дифференциальных уравнений каждого временного слоя ведется методом конечных разностей.

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

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

Наиболее эффективными методами динамического расчета упругих систем с односторонними связями, как представляется автору, являются методы, изложенные в работах Кохманюка С.С., Янютина Е.Г., Баженова В.А., Гуляева В.И. [103,104,174-178]. Однако разрешающая система интегральных уравнений в этих работах получается относительно коэффициентов разложения в тригонометрический ряд функции перемещений упругой системы. Такой метод приводит к необходимости на каждом временном шаге решения выполнять суммирование гармоник решения для получения функции перемещений упругой системы. Это обстоятельство ввиду значительного увеличения времени расчета, затрудняет решение задач, в которых жесткость односторонних связей на несколько порядков выше жесткости упругой системы, т.к. в этом случае:

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

2) необходимо суммировать большое число гармоник для получения решения с достаточной точностью.

Методы, предложенные в работах [103,104,174-178], реализованы на примерах, в которых упругие элементы (балки, пластины, оболочки) шарнирно оперты. Математические модели не учитывают внешнее трение и внутреннее трение в упругих элементах. В этих работах не исследуется сходимость решения.

Упругой системой с большим числом односторонних связей является волновая зубчатая передача. Статическому расчету сил взаимодействия звеньев волновой передачи, как упругой системы с односторонними связями посвящены работы Ковалева Н.А., Кленикова С.С., Цейтлина Н.И. и др. Однако для расчета кинематической точности волновых передач методики, предложенные в этих работах, не применялись. Это связано с необходимость усложнения расчетной схемы.

Исследованию кинематической точности волновых зубчатых передач посвящены работы Е.Г. Гинзбурга, И.И. Васильевой, П.К. Попова, И.Л. Смирновой, Г.А. Тимофеева, А.Ф. Фирсанова, Ф.И. Фурсяка, Л.С, Черновой, Л.О. Штриплинга и др.

Влияние многопарности и многозонности зацепления на кинематическую погрешность рассматривалась в работах [37,61,75,121, 136,137,138,160,162,163,164]. Несмотря на противоречивость количественных оценок, во всех работах отмечается, что многопарность зацепления ведет к значительному повышению кинематической точности волновых передач. Степень точности волновой передачи на 2, 3 степени выше точности составляющих передачу колес.

Одним из важных достоинств волновых зубчатых передач, так же как и передач с зацеплением Новикова М. Л. [127], является малая чувствительность передачи к монтажным погрешностям.

Анализу влияния точности изготовления зубчатых колес и основных деталей на кинематическую погрешность волновых передач посвящены работы [40,61,62,79,93,121,136,155,160,161, 162,163,164, 166,167]. В этих работах показано, что основное влияние на кинематическую точность оказывают погрешности генератора волн и деталей, поддерживающих зубчатые колеса. Волновая передача становится чувствительной к точности изготовления зубчатых колес лишь при высокой точности других деталей передачи. Так в работе [79] отмечается, что:

1) повышение точности изготовления зубчатых колес с 9-й по 6-ю при 7-й степени точности прочих основных деталей не приводит к заметному уменьшению кинематической погрешности волновой зубчатой передачи;

2) повышение точности изготовления зубчатых колес с 9-й по 6-ю при 4-й степени точности прочих основных деталей снизило кинематическую погрешность волновой зубчатой передачи в 2 раза.

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

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

Цель работы. В связи с выше изложенным в данной работе поставлены следующие задачи:

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

2) разработать эффективный метод статического расчета упругих систем с односторонними связями;

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

4) применить разработанный метод для расчета колебаний тонкостенных упругих элементов (балок, пластин), ограниченных односторонними связями;

5) проанализировать сходимость решения;

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

7) оценить влияние инерции вращения и деформации поперечного сдвига сечений на реакции в односторонних связях и перемещения упругих элементов;

8) с целью подтверждения предложенных методов расчета провести экспериментальные исследования;

9) разработать пакет программ, реализующий предложенные в данной работе методы;

10) на основе предложенных методов разработать методику определения кинематической погрешности волновой зубчатой передачи с кулачковым и электромагнитным генераторами волн. Исследовать кинематическую погрешность этих передач.

Научная новизна состоит в следующем:

1) предложен эффективный метод статического расчета упругих систем с односторонними связями;

2) разработан метод расчета колебаний упругих систем с неинерционными и инерционными односторонними связями повышенной жесткости;

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

-354) разработаны математические модели и методики для расчета кинематической погрешности ВЗП с кулачковым и электромагнитным генераторами волн;

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

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

Практическая ценность.

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

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

В первой главе «Методы расчета упругих систем с односторонними связями» проводится сравнение быстродействия различных методов статического расчета упругих систем с нелинейными односторонними связями. Предлагается новый метод статического расчета. Рассматривается метод расчета взаимодействия шара и балки. Описываются методы расчета колебаний упругих систем с неинерционными и инерционными односторонними связями.

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

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

В четвертой главе «Практическое применение разработанных методов для расчетов кинематической погрешности волновых зубчатых передач как упругих систем с односторонними связями» предлагаются методики расчета кинематической погрешности волновых передач с кулачковым и электромагнитным генераторами волн. Методики основаны на предложенных методах статического и динамического расчета упругих систем с односторонними связями. Численно исследуется кинематическая погрешность ВЗП-80. Численно и экспериментально исследуется ВЗП с электромагнитным генератором.

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

Заключение диссертации по теме «Динамика, прочность машин, приборов и аппаратуры», Люминарский, Игорь Евгеньевич

Выводы

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

2. Разработан метод динамического расчета ударного нагружения существенно нелинейных упругих систем с односторонними связями, позволяющий эффективно решать многомерные задачи, в которых жесткость односторонних связей на несколько порядков выше жесткости упругой системы. Метод основывается на теории соударения упругих систем, учитывающей местные деформации, предложенной Тимошенко С.П. и получившей развитие в работах Бидермана B.JI.

3. Уточнено тестовое решение С.П.Тимошенко об ударном взаимодействии шара с балкой. Численными расчетами показано, что в процессе удара шар с балкой взаимодействует не два, а три раза.

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

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

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

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

8. Установлено, что: а) максимальная собственная кинематическая погрешность волновых зубчатых передач с кулачковым генератором волн имеет небольшие значения. Например, для ВЗП-80 эти значения не превышают 1,8 мкм; б) максимальная собственная кинематическая погрешность у волновых зубчатых передач с нечетным числом тел качения примерно в четыре раза меньше, чем у передач с четным числом тел качения; в) зависимость максимальной кинематической погрешности F'or от смещения кулачка ен существенно нелинейная. Максимальная кинематическая погрешность F'or до некоторого значения смещения кулачка ен увеличивается незначительно. Затем F'or резко возрастает. При больших значениях смещения кулачка ен в процессе работы волновой передачи в некоторых положениях все зубья одной из двух полуволн выходят из контакта. Переход к одноволновому зацеплению сопровождается резким увеличением максимальной кинематической погрешности; г) повышение технологической точности деталей ВЗП-80 с 9-го до 7-го квалитета позволяет уменьшить максимальную кинематическую погрешность примерно в 1,7 раза. В тоже время повышение точности деталей с 7-го до 5-го квалитета уменьшает максимальную кинематическую погрешность примерно в 1,1 раза; е) применение кулачка, профиль которого определяется по формуле w{(p) = w0 cos(2(р) -Кsin2{2ф), где w0 = 0,48; К = 0,04, снижает наибольшую кинематическую погрешность передачи ВЗП-80 примерно в два раза, если основные детали передачи изготовлены по 7-8 квалитету точности, а передача нагружена моментом до 10 Н - м.

9. Показано, что собственная кинематическая погрешность торцевой волновой зубчатой передачи с электромагнитным генератором волн имеет большие значения. Так для торцевой передачи с 8-ю толкателями и диаметром гибкого колеса 90 мм она составила 314 мкм, что существенно выше собственной кинематической погрешности радиальных волновых зубчатых передач.

3.3. Заключение

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

Местные сближения Л(Я) упругой системы и односторонних связей в монографии определяется по теории Герца. Этой теорией можно пользоваться в том случае, когда в месте контакта нет пластических деформаций.

Если в зоне контакта упругой системы и односторонней связи наибольшее контактное давление сгтах > [сг] , где [сг] — давление при котором наступают пластические деформации, то для определения местного сближения л(л) необходимо использовать зависимости отличные от формулы Герца.

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

При колебаниях тонкостенных конструкций (балок, пластин, оболочек и т.д.), как правило, рассматривают 3 класса волн [50]: длинные, средние и короткие. Волны считаются длинными, если можно пренебречь инерцией вращения и сдвигом. В случае средних волн эти факторы подлежат учету, но их влияние значительно меньше влияния поперечной инерции. Короткие волны характеризуются тем, что влияние инерции вращения и сдвига имеет порядок, одинаковый с влиянием поперечной инерцией. Из приближенных расчетов прямоугольного стержня с высотой Я h в [50] получены следующие оценки: длинные волны-->40, средние h

X X 8 < — < 40, короткие--<8 (л - длина волны). h h

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

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

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

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

Если расчеты, проведенные по теории Бернулли - Эйлера, показывают существенное влияние средних волн на решение, то необходимо уточнить модель колебаний по теории С.П. Тимошенко.

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

Список литературы диссертационного исследования доктор технических наук Люминарский, Игорь Евгеньевич, 2009 год

1. Агамиров В.Л., Клеников С.С., Люминарский И.Е. Распределение нагрузок между элементами сдвоенной волновой передачи при ее сборке // Вестник машиностроения. - 1986. - № 10. - С. 17—20.

2. Аверьянов А.В. Волновые электромеханизмы / Под ред. B.C. Рыбакова. -М.: Информэлектро, 1970. —63 с.

3. Алабужев П.М., Стахановский Б.Н., Шпигельбурд И.Я. Введение в теорию удара. Новосибирск, 1970. - 157 с.

4. Александров А.В., Потапов В.Д. Основы теории упругости и пластичности: Учебник для строит, спец. вузов. — М.: Высшая школа, 1990.-400 с.

5. Александров Е.В., Соколинский В.Б. Прикладная теория и расчеты ударных систем. М.: Наука, 1969. - 250 с.

6. Алимов О.Д., Манжосов В.К., Еремьянц В.Э. Удар. Распределение волн деформации в упругих системах. М.: Наука, 1985. — 357 с.

7. Алфутов Н.А., Клеников С.С. Расчет сил взаимодействия упругих элементов волновых передач шаговым методом // Вестник машиностроения. 1978. - № 7. - С. 26-29.

8. Аргирос Дж., Шариф Д. Теория расчета пластин и оболочек с учетом деформаций поперечного сдвига на основе метода конечных элементов // Расчет упругих конструкций с использованием ЭВМ. — Л.: Судостроение, 1974.-С. 179-210.

9. Астрахан А.Х. Расчет конструкций с односторонними связями на переменные нагрузки // Механика стержневых систем и сплошных сред: Межвузовский тематический сборник трудов (Л). — 1980. -Вып. 13. -С. 45-51.

10. Бабицкий В.И., Крупенин В.Л. Колебания в сильно нелинейных системах. М.: Наука, 1985. - 320 с.

11. Багреев В.В. О численном решении задач удара по методу Сирса -Тимошенко // Труды Моск. ин-та инж. ж-д. транспорта (М). — 1966. — Вып. 225.-С. 305-312.

12. Багреев В.В. Об одном алгоритме расчета упругих систем на действие удара // Труды Моск. ин-та инж. ж-д. транспорта. — М., 1966.- Вып. 225. С. 234-240.

13. Баженов В.А. Вынужденные колебания пластин на одностороннем упругом основании // Проблемы прочности. 1983. — № 3.- С. 101-103.

14. Баженов В.А., Гуляев В.И., Кондаков Г.С. Нелинейные колебания пологой оболочки, односторонне взаимодействующей с упругим основанием // Тезисы докл. XIII Всесоюз. конф. по теории пластин и оболочек; В 4 ч. Таллинн, 1983. — 4.1. - С. 61-71.

15. Баженов В.А., Гоцуляк Е.А. Вынужденные колебания удлиненной цилиндрической панели на одностороннем основании // Прикладная механика. 1985. - № 8. - С. 51-56.

16. Устойчивость и колебания деформируемых систем с односторонними связями / В.А. Баженов, Е.А. Гоцуляк, Г.С. Кондаков и др. Киев: Выща школа. Головное изд-во, 1989.-399с.

17. Баженов В.Г., Чекмарев Д.Т. Об индексной коммутативности численного дифференцирования // Журнал вычислительной математики и математической физики. — 1989. — № 5. — С. 662-674.

18. Баженов В.Г., Чекмарев Д.Т. Вариационно-разностные схемы в нестационарных волновых задачах динамики пластин и оболочек. -Н. Новгород: Изд-во Нижегор. ун-та, 1992. — 159 с.

19. Баженов В.Г., Чекмарев Д.Т. Решение задач нестационарной динамики пластин и оболочек вариационно-разностным методом. — Н. Новгород: Изд-во Нижегор. ун-та, 2000. 118 с.

20. Байокки К., Капело А. Вариационные и квазивариационные неравенства: Пер. с англ./ Под ред. В.И. Агошкова. М.: Наука, Физматлит, 1988. - 448 с.

21. Бартеньев О.В. Фортран для профессионалов. Математическая библиотека IMSL. М.: Диалог-МИФИ, 2000. - Ч. 1. - 448 с.

22. Бартеньев О.В. Фортран для профессионалов. Математическая библиотека IMSL. -М.: Диалог-МИФИ, 2001. -Ч. 3. -368 с.

23. Батуев Г.С., Голубко Ю.В., А.К. Ефремов А.К.,Инженерные методы исследования ударных процессов. -М.: Машиностроение, 1977. -240 с.

24. Белов Е.Г. Алгоритм решения задачи квадратичного и линейного программирования методом сопряженных градиентов // Стандартные программы решения задач математического программирования (М). — 1970.-Вып. 18. -С. 186-192.

25. Беляев Н.М. Труды по теории упругости и пластичности. — М.: Машгиз, 1957. -209 с.

26. Блох М.В. О поперечном ударе по балке конечной длины

27. Строительная механика и расчет сооружений. 1968. — № 6. — С. 16-19.

28. Болдечев В.П. О решении статических и динамических задач изгиба пластин средней толщины // Известия ВНИИГ им. Б.Е. Веденеева. -1978.-Т. 120.-С. 86-93.

29. Болдечев В.П. Повышение эффективности метода конечных элементов при решении вырождающихся задач // Вопросы динамики и прочности (Рига). 1983. - Вып. 42. - С. 38-48.

30. Бидерман B.JI. Теория удара. — М.: Наука, 1952. 76 с.

31. Бидерман B.JI. Теория механических колебаний. — М.: Высшая школа, 1980.-408 с.

32. Блехман И.И., Джанелидзе Г.Ю. Вибрационное перемещение. М.: Наука, 1964.-412 с.

33. Вабищевич П.Н. Численные методы решения задач со свободной границей. -М.: Изд-во МГУ, 1987. 163 с.

34. Вабищевич П.Н. Численное решение краевых задач для эллиптических уравнений четвертого порядка // Журнал вычислительной математики и математической физики. — 1984. — Т. 24. — С. 41-45.

35. Вавилов А.А. и др. Метод гармонической линеаризации в проектировании нелинейных систем автоматического управления. — М.: Машиностроение, 1970. 567с.

36. Вибрации в технике: Справочник; В 6т. / Под ред. К.В. Фролова. — М.: Машиностроение, 1999.-Т. 1. 504 с.

37. Вибрации в технике: Справочник; В 6т. / Под ред. К.В. Фролова. М.: Машиностроение, 1995.-Т.6. - 456 с.

38. Волков Д.П., Крайнев А.Ф. Волновые зубчатые передачи. Киев: Техника, 1976. - 222 с.

39. Воробьев А.В. Респонсивный привод. — М.: Машиностроение, 1978. -160 с.

40. Гасс С.И. Линейное программирование (Методы и приложения). М.: Физматгиз, 1961.-303 с.

41. Главачек И. Решение вариационных неравенств в механике. — М.: Мир, 1986.-270 с.

42. Голоскоков Е.Г., Филиппов А.П. Нестационарные колебания деформируемых систем. — Киев: Наукова думка, 1977. 339 с.

43. Гордеев В.Н., Перельмутер А.В. Расчет упругих систем с односторонними связями как задача квадратичного программирования // Исследование по теории сооружений (М). — 1967. — Вып. XV. —1. С. 23-29.

44. Гордон J1.K. Расчет и экспериментальные исследования пластин средней толщины / JI.K. Гордон, Я.Г. Скоморовский, Е.Ш. Фридман, Б.А. Шойхет // Известия ВНИИГ им. Б.Е. Веденеева. 1971. - Т. 95. -С. 142-166.

45. Гордон JI.K. К расчету пластин и оболочек методом конечных элементов // Известия ВНИИГ им. Б.Е. Веденеева. — 1972. Т. 99. -С. 168-188.

46. Горшков А.Г., Тарлаковский Д.В. Динамические контактные задачи для абсолютно жестких тел упругого полупространства: Препринт. -М.: Изд-во Моск. авиац. ин-та, 1989. 48 с.

47. Горшков А.Г., Тарлаковский Д.В. Результирующие реакции в пространственной задаче об ударе телом по упругому полупространству

48. Изв. АН СССР. Механика твердого тела. 1987. - № 5. - С. 21-29.

49. Горшков А.Г., Тарлаковский Д.В. Динамические контактные задачи с подвижными границами. — М.: Наука. Физматлит, 1995. 352 с.

50. Гольдсмит В. Удар. Теория и физические свойства соударяемых тел. — М.: Стройиздат, 1965. 448 с.

51. Гольдсмит В. Удар и контактные явления при средних скоростях // Физика быстропротекающих процессов; В 2 т. М., 1971. - Т. 2. -С. 21-29.

52. Григолюк Э.И., Толмачев В.М. Контактные задачи теории пластин и оболочек. — М.: Машиностроение, 1980. — 411 с.

53. Григолюк Э.И., Селезов И.Т. Неклассические теории колебаний стержней, пластин и оболочек. — М.: ВИНИТИ, 1973. — 272 с.

54. Гуляев В.И., Баженов В.А., Гоцуляк Е.А. Устойчивость нелинейных механических систем. — Львов: Вища школа, 1982 — 255 с.

55. Гуляев В.И., Баженов В.А., Попов С.Л. Прикладные задачи теории нелинейных колебаний механических систем: Учеб. пособие для втузов. -М.: Высшая школа, 1989. 383 с.

56. Дерябин М.В., Козлов В.В. К теории систем с односторонними связями // Прикладная математика и механика. — 1995. — Т. 59. — Вып. 4. -С. 531-539.

57. Дерябин М.В. О реализации неудержевающих связей // Прикладная математика и механика. 1994. - Т. 58. - Вып. 6. - С. 136-140.

58. Дерябин М.В. Общие принципы динамики и теория односторонних связей // Вестник Моск. ун-та. Сер.1. Математика, механика. 1998. -№ 1. — С. 53-59.

59. Динник А.Н. Удар и сжатие упругих тел // Избранные труды; В 3 т. Киев, 1952. - Т. 1. - С. 420-426.

60. Емельянов А.Ф., Попов П.К., Фирсаев А.Ф. Расчет кинематической погрешности волновой зубчатой передачи с учетом податливости звеньев // Вестник машиностроения. — 1983. № 7. — С. 9-12.

61. Емельянов В.И. Исследование кинематической точности волновых редукторов: Дисс. . канд. техн. наук. -М ., 1972. 155 с.

62. Журавлев В.Ф., Вуфаев Н.А. Механика систем с неудерживающими связями. М.: Наука, 1993. — 239 с.

63. Зегжда С.А. Соударение колец // Вестник Ленингр. ун-та. Сер. 1. -1986.- Вып. 1.-С. 120-127.

64. Зегжда С.А., Филиппов Н.Г. О соударении цилиндров вдоль их образующих // Вестник Ленингр. ун-та. Сер.1. - 1986. - Вып. 3. -С. 71-77.

65. Зегжда С.А., Вернигор В.Н. О колебаниях балки бесконечной длины под действием сосредоточенной силы // Вестник Ленингр. ун-та. — Сер.1.- 1986.-Вып. 4. -С. 350-357.

66. Зегжда С.А., Вернигор В.Н. Исследование влияния инерции поворота и деформации сдвига при поперечном ударе // Вестник Ленингр. унта. Сер.1. - 1987. - Вып. 2(8). - С. 54-57.

67. Зегжда С.А. Соударение кольца и балки // Колебание и устойчивость механических систем. Л., 1988. - 180 с.

68. Зегжда С.А. Уточненное решение задачи о соударении шаров

69. Прикладные задачи динамики и устойчивости механических систем. — Л., 1990.-272 с.

70. Зегжда С.А. Соударение упругих тел. — СПб.: Издательство С.-Петербургского ун-та, 1997. 316 с.

71. Зенкевич О. Метод конечных элементов в технике. М.: Мир, 1975. -541 с.

72. Златин А.Н., Пупырев В.А. Динамическая реакция круговой пластины на опорах одностороннего действия // Труды Ленингр. политех, ин-та. — 1982.-Вып. 386.-С. 50-55.

73. Зуховский С.И., Авдеева Л.И. Линейное и выпуклое программирование. М.: Наука, 1967. - 410 с.

74. Иванов А.П. Динамика систем с механическими соударениями. М.: Международная программа образования, 1997. - 336 с.

75. Иванов М.Н., ШуваловС.А., Финогенов В.А. Экспериментальное определение количества одновременно зацепляющихся зубьев и величин их деформации в волновой передаче // Известия вузов. Машиностроение. 1968. - № 9. - С. 37-40.

76. Иванов М.Н. Волновые зубчатые передачи. — М.: Высшая школа, 1981.- 184 с.

77. Иванов Ю.С. Разработка методики расчета нагруженности и жесткостных параметров упругих элементов шаговых гидромоторов с волновым зубчатым зацеплением: Дис. . канд. техн. наук. — М., 1991. — 270 с.

78. Иванова Е.А. Асимптотический и численный анализ высокочастотных свободных колебаний прямоугольных пластин // Механика твердого тела. 1998. - № 2. - С. 163-174.

79. Истомин С.Н. Кинематическая точность приборных волновых передач

80. С.Н. Истомин, С.А. Шувалов, П.К. Попов и др. — М.: Машиностроение, 1987.-160 с.

81. Кильчевский НА. Теория соударения твердых тел. — Киев: Наукова думка, 1969.-310 с.

82. Кильчевский Н.А. Динамическое контактное сжатие твердых тел.- Киев: Наукова думка, 1976. 280 с.

83. Киселев В.А. Строительная механика: Спец. курс. Динамика и устойчивость сооружений: Учебник для вузов. 3-е изд., испр. и доп. -М.: Стройиздат, 1980. - 616 с.

84. Клеников С.С. Волновая передача как упругая система с односторонними связями // Известия вузов. Машиностроение. — 1978.- № 10.-С. 51-55.

85. Клеников С.С., Сергеев B.C. Взаимодействие упругих элементов нагруженной волновой передачи с кулачковым генератором

86. Известия вузов. Машиностроение. 1978. - № 11. - С. 40—44.

87. Клеников С.С., Сергеев B.C. О линейном нормировании этапного решения задачи определения нагрузок на зубья и тела качения волновой передачи с кулачковым генератором // Известия вузов. Машиностроение. 1980. - № 4. - С. 45^19.

88. Клеников С.С. Разработка методов и рациональных алгоритмов решения задач статики и динамики силового взаимодействия упругих элементов волновых передач: Дис. .докт. техн. наук. — М., 1989. -387 с.

89. Клеников С.С., Иванов Ю.С. Построение рациональных алгоритмов поиска опорных систем нагруженных односторонних связей (Сообщение 2) // Известия вузов. Машиностроение. — 1986. — № 10. -С. 14-18.

90. Клеников С.С., Люминарский И.Е., Люминарский С.Е. Шаговый поиск опорных систем нагруженных односторонних связей методом введения восстанавливающих сил // Известия вузов. Машиностроение. 1987. -№7.-С. 34-40.

91. Клеников С.С., Люминарский И.Е. Силовое взаимодействие упругих элементов нагруженных сдвоенных волновых зубчатых передач

92. Вестник машиностроения. — 1988. — № 1. С. 17-20.

93. Клеников С.С., Люминарский И.Е., Семин И.И. Расчетная модель волновых передач с учетом несимметрии нагружения элементов по волнам зацепления // Вестник машиностроения. — 1993. — № 1. — С. 17-19.

94. Клеников С.С. Исследование силового взаимодействия упругих элементов волновой передачи: Дис. .канд. техн. наук. — М.: МВТУ, 1974.-158 с.

95. Клыпин А.В. и др. Расчет кинематической погрешности волновой передачи на ЭВМ /А.В. Клыков, П.К. Попов, А.Ф. Фирсаев, А.Ф. Емельянов // Вестник машиностроения. — 1985. № 11. - С. 9-12.

96. Ковалев Н.А. Некоторые вопросы теории волновых зубчатых передач // Машиностроение. 1973. - № 2. - С. 48-53.

97. Ковалев Н.А. О распределении нагрузки по зубьям в волновой передаче // Машиностроение. 1974. - № 5. - С. 44-49.

98. Ковалев Н.А. Передачи с гибкими колесами. — М.: Машиностроение, 1979.-200 с.

99. Ковтуненко В.А. Метод численного решения задачи о контакте упругой пластины с препятствием // ПМТФ. 1994. - № 5. - С. 142-146.

100. Ковтуненко В.А. Итерационный метод решения вариационных неравенств в контактной упругопластической задаче с использованием метода штрафа // Журнал вычислительной математики и математической физики. 1993. - Т. 33. - № 9. - С. 1409-1415.

101. Ковтуненко В.А. Численное решение задачи о контакте упругопластической балки для модели Тимошенко // Изв. АН СССР. Механика твердого тела. 1996. - № 5. - С. 79-84.

102. КоллатцЛ. Задачи на собственные значения с техническими приложениями: Пер. с нем. М.: Наука, 1968. - 504 с.

103. Кончковский 3. Плиты. Статические расчеты / Пер. с пол.

104. М.П. Предтеченского; Под ред. А.И. Цейтлина. — М.: Стройиздат, 1984. -480 с.

105. Корн Г., Корн Т. Справочник по математике для научных работников и инженеров. М.: Наука, 1968. - 720 с.

106. Кохманюк С.С., Янютин Е.Г., Романенко Л.Г. Колебания деформируемых систем при импульсных и подвижных нагрузках. -Киев: Наукова думка, 1980. 232 с.

107. Кохманюк С.С. Динамика конструкций при воздействии кратковременных нагрузок / С.С. Кохманюк, А.С. Дмитриев,

108. Г.А. Шелудько, А.Н. Шупиков и др. — Киев: Наукова думка, 1989. — 304 с.-320105. Кравчук А.С. К задаче Герца для линейно- и нелинейно-упругих телконечных размеров // Прикладная математика и механика. 1977.- Т. 41, Вып. 2. С. 1020-1029.

109. Кравчук А.С. Постановка задачи о контакте нескольких деформируемых тел как задачи нелинейного программирования // Прикладная математика и механика. 1978. — Т. 42, Вып. 3. -С. 620-631.

110. Кравчук А.С. К теории контактных задач с учетом трения на поверхности соприкосновения // Прикладная математика и механика. — 1980. Т. 44, Вып. 1. - С. 732-740.

111. Кравчук А.С. Решение некоторых пространственных контактных задач с учетом трения на поверхности соприкосновения // Трение и износ.-1981.-№ 4.-С. 15-19.

112. Кравчук А.С. Вариационный метод исследования контактного взаимодействия и его реализация на ЭВМ // Расчеты на прочность: Сб. статей / Под ред. Н.Д. Тарабасова (М). 1984. - Вып. 25.-С. 220-241.

113. Кравчук А.С. Вариационный метод в динамических контактных задачах // Механика деформируемых тел и конструкций: Сб. статей. -Ереван, 1985. С. 140-147.

114. Кузнецов А.Д. Динамическое поведение многослойных оболочек с учетом особенностей контакта слоев // Расчеты на прочность и жесткость. 1982. - № 4. - С. 73-84.

115. Кузнецов Э.Н. Введение в теорию вантовых систем. — М.: Стройиздат, 1969. -340 с.

116. Кузнецов Э.Н. Статико-кинематический анализ систем, содержащих односторонние связи // Исследования по теории сооружений. — 1972. -Вып. 19.-С. 320-331.

117. Курант Р., Гильберт Д. Методы математической физики: Пер. с англ.- М. Л.: ГИТГА, 1951. - 476 с.

118. Кюнци Г.П., Крелле В. Нелинейное программирование. — М.: Советское радио, 1965. 300 с.

119. Лиос Ж.Л. Некоторые методы решения нелинейных краевых задач.- М.: Мир, 1972.-587 с.

120. Люминарский И.Е. Расчет силового взаимодействия элементов сдвоенной волновой зубчатой передачи с целью определения ее рациональных параметров: Дис. .канд. техн. наук. М., 1987. - 167 с.

121. Люминарский С.Е., Люминарский И.Е. Собственная кинематическая погрешность волновых зубчатых передач с кулачковым генератором волн // Сб. науч. трудов МАСИ. М.: МАСИ, 1995. С. 47-53.

122. Маланов С.Б., Уткин Г.А. Ударное взаимодействие сосредоточенного объекта с одномерной упругой системой // Прикладная математика и механика. 1988. — Т. 52, Вып. 1. — С. 42-46.

123. Маланов С.Б., Уткин Г А. Косой удар материальной точки по бесконечной струне на упругом основании // Прикладная математика и механика. 1988. - Т. 52, Вып. 5. - С. 881-883.

124. Малышев А.П. Переходные процессы в оболочке с расслоениями

125. Изв. АН СССР. Механика твердого тела. 1978. -№ 6. - С. 101-105.

126. МарчукГ.И. Методы вычислительной математики: Учеб. пособие.- 3-е изд., перераб. и доп. М.: Наука. Физматлит, 1989. — 608 с.

127. Мину М. Математическое программирование. Теория и алгоритмы / Пер. с фр. и предисловие А.И. Штерна. М.: Наука. Гл. ред. физ.-мат. лит., 1990. - 488 с.

128. Михлин С.Г. Вариационные методы в математической физике.- М.: Наука, 1970. 476 с.-322126. Нагаев Р.Ф. Механические процессы с повторными затухающими соударениями. — М.: Наука. Физматлит, 1985. — 200 с.

129. Новиков M.JI. Зубчатые передачи с новым зацеплением. -М.: ВВИА им. Н.Е.ЖуковскогоД958. — 180 с.

130. Основы современных методов расчета на прочность в машиностроении / Под ред. С.Д. Пономарева. М., 1950. - 703 с.

131. Осташявичюс В.В., Рудгальвис Б.В., Рагульскене B.J1. Контактные системы. — Л.: Машиностроение. Ленингр. отд-ние, 1978. -279 с.

132. Остроградский М.В. Полное собрание трудов. М.: Наука, 1966. - 327 с.

133. Пановко Я.Г. Введение в теорию механического удара. — М.: Наука, 1977.-223 с.

134. Пановко Я.Г. Внутреннее трение при колебаниях упругих систем. — М.: Физматлит, 1960. 193 с.

135. Перельмутер А.В. Использование методов квадратичного программирования для расчета систем с односторонними связями // Исследование по теории сооружений (М). -1972. Вып. XIX. -С. 250-257.

136. Писаренко Г.С., Матвеев В.В., Яковлев А.П. Методы определения характеристик демпфирования колебаний упругих систем. Киев: Наукова думка, 1976. — 86 с.

137. Полянин А.Д. Справочник по линейным уравнениям математической физики. -М.: Физматлит, 2001. 576 с.

138. Попов П.К. Исследование положения выходного вала волновой зубчатой передачи: Дисс. . канд. техн. наук (05.02.02). М., 1972. -190с.

139. Попов П.К., Фирсаев А.Ф.,Шувалов С.А. Математическая модель процесса возникновения кинематической ошибки в волновой зубчатой передаче // Труды МВТУ. 1978. - № 287 - Прочность и надежность деталей машин. — С. 50-55.

140. Попов П.К., Штриплинг JI.O. Динамическая модель возникновения кинематической погрешности волновой зубчатой передачи // Известия вузов. Машиностроение. — 1986. №1. - С. 46-50.

141. Попов П.К. Расчетно-экспериментальное обеспечение точности зубчатых передач: Дис. .докт. техн. наук. — М.: МВТУ, 1996. 241 с.

142. Портаев Л.П. Расчет систем с односторонними связями на возрастающую нагрузку // Изв. АН СССР. Механика твердого тела. — 1978.-№ 1.-С. 21-28.

143. Прагер В. Введение в механику сплошных сред. — М.: Изд-во иностр. лит., 1963. 340 с.

144. Рабинович И.М. К теории вантовых ферм// Техника и экономика путей сообщения. 1924. — № 1. — С. 12-15.

145. Рабинович И.М. Некоторые вопросы теории сооружений, содержащих односторонние связи // Инженерный сборник. АН СССР. 1950.-Т. VI.-С. 43-49.

146. Рабинович И.М. Курс строительной механики стержневых систем. Ч. II. М.: Госстройиздат, 1954. - 390 с.

147. Рабинович И.М. Вопросы теории статического расчета сооружений с односторонними связями. — М.: Стройиздат, 1975. 220 с.

148. Размадзе Г.Н. Инженерные вопросы теории удара. Расчет элементов конструкций на удар в пределах упругих деформаций. М.: Наука, 1959.-150 с.

149. Расчеты на прочность в машиностроении; В 3 т. / Под ред. С.Д. Пономарева. М., 1959. - Т. 3. - 380 с.

150. Сагомонян А.Я. Волны напряжения в сплошных средах. — М.: Изд-во МГУ, 1985. 416 с.

151. Семин И.И., Люминарский И.Е., Люминарский С.Е. О положительной определенности матриц упругих систем с односторонними связями // Сб. науч. трудов МАСИ. -М.: МАСИ, 1995. С. 41-48.

152. Сергеев В.В. Силовое взаимодействие элементов волновых передач с кулачковым и роликовым генераторами: Дис. . канд. техн. наук. — М., 1975.- 137 с.

153. Сергеев B.C. Силовое взаимодействие элементов стандартных волновых зубчатых передач и напряженно-деформированное состояние гибкого колеса с учетом податливости звеньев и начальных зазоров: Дис. .канд. техн. наук. — М., 1985. 191 с.

154. Скляр В.А., Филиппов А.П. Колебания стержней и плит при ударе // Динамика и прочность машин: Респ. межвед. научно-техн. сб. (Киев). 1967. - Вып. 6. - С. 156-168.

155. Тимофеев Г.А. Разработка методов расчета и проектирования волновых зубчатых передач для приводов следящих систем: Дис. . .докт. техн. наук. М.: МВТУ, 1997. - 351 с.

156. Тимошенко С.П. Колебания в инженерном деле. М. - Л.: Машиностроение, 1967. — 444 с.

157. Уилкинсон Дж. Алгебраическая проблема собственных значений. -М.: Наука, 1970.-401 с.

158. Фейгин М.И. Вынужденные колебания систем с разрывными нелинейностями. — М.: Наука, 1994. — 285 с.

159. Филиппов А.П. Колебания деформируемых систем. — М.: Машиностроение, 1970. 734 с.

160. Фирсаев А.Ф. Разработка и исследование приводов подач высокоточных станков с ЧПУ на базе волновых зубчатых передач: Дисс. . канд. техн. наук. -М., 1983. 168 с.

161. Фурсяк Ф.И., Скворцова Н.А., Тимофеев Г.А. Степень влияния ошибок изготовления волновой зубчатой передачи на ее кинематическую точность // Труды МВТУ. 1976. - № 227, вып.7 -Теория механизмов. — С. 8—13.

162. Фурсяк Ф.И., Скворцова Н.А., Тимофеев Г.А. Степень влияния ошибок изготовления волновой зубчатой передачи на ее кинематическую точность // Теория механизмов. — 1976. Вып. 3. -С. 14-17.

163. Фурсяк Ф.И. Анализ действующих погрешностей волновой зубчатой передачи // Труды МВТУ. — 1981. № 352, вып.9 - Теория механизмов. - С. 86-94.

164. Чернова JI.C., Гинзбург Е.Г. О влиянии накопленных погрешностей окружного шага зубчатых колес на кинематическую точность двухволновой передачи при наличии многопарного зацепления

165. Волновые зубчатые передачи: Тез. докл. научн.-техн. конф. ЛВИКА им. А.Ф. Можайского, Л., 1969. - С. 45-49.

166. Цейтлин Н.И., Рафолович Л.Б. Алгоритмы упруго-силового взаимодействия деталей волновой зубчатой передачи // Расчеты на прочность (М). -1980. Вып. 21. - С. 243-254.

167. Шамирян-Пахлеванян Р.И. Кинематическая точность малонагруженных волновых зубчатых передач: Дисс. . канд. техн. наук. М., 1971.-174 с.

168. Шамирян-Пахлеванян Р.И. Влияние погрешностей изготовления основных деталей и нагрузки на кинематическую точность волновых зубчатых передач // Волновые зубчатые передачи: Тез. докл. Всес. симпозиума.-М., 1973.-С. 105-108.

169. Шульгин Ю.Б. Геометрическая интерпретация расчета конструкции содержащей односторонние связи // Механика стержневых систем и сплошных сред: Межвуз. темат. сб. тр. — Л.: ЛИСИ, 1979. — С. 23—28.

170. Янютин Е.Г. Нестационарное деформирование упругого пространства с расширяющейся сферической полостью // Изв. АН СССР. Механика твердого тела. 1983. - № 6. - С. 86-89.

171. Янютин Е.Г. Нестационарное деформирование балки на упругом основании переменной жесткости // Проблемы машиностроения. — 1984.-Вып. 22.-С. 11-17.

172. Янютин Е.Г. Нестационарное деформирование упругого пространства с полостью, подкрепленного замкнутой сферической оболочкой // Прикладная механика. 1984. -№ 4. - С. 23-27.

173. Янютин Е.Г. Балка на упругом одностороннем однослойном инерционном основании // Динамика и прочность машин. — 1985. Вып. 42.-С. 37—43.

174. Янютин Е.Г. Нестационарное деформирование балки конечной длины на инерционном одностороннем основании / АН УССР. Ин-т проблем машиностроения. — Харьков, 1985. 13 с. (Деп. в ВИНИТИ 31.10.85, №7603).

175. Янютин Е.Г. Нестационарное деформирование цилиндрической оболочки, односторонне контактирующей со средой // Проблемы машиностроения. 1985.-Вып. 23. —С. 6-11.

176. Янютин Е.Г. Нестационарное деформирование замкнутой сферической оболочки в среде, сопротивляющейся только вдавливанию // Проблемы машиностроения. — 1986. — Вып. 25. — С. 31—37.

177. Янютин Е.Г., Ермоленко Л.В. К исследованию нестационарных продольных колебаний низа бурильной колонны / АН УССР. Ин-тпроблем машиностроения. — Харьков, 1989. 11 с. (Деп. в ВИНИТИ 25.04.89, № 2686).

178. Янютин Е.Г., Светличная С.Д. Нестационарное упругое плоское неосесимметричное деформирование многослойного цилиндра / АН УССР. Ин-т проблем машиностроения. Харьков, 1990. - 16 с. — Деп. в ВИНИТИ 06.02.90, № 658.

179. Adams G.G. Moving loads on elastic strips with one-sided constraints // J. Eng. Sci. 1970. — V. 14. -№ 12.-P. 1071-1083.

180. Asano Naoki. A hybrid type of virtual work principle for impact contact problems of two bodies // Bull. JSME. 1986. - Vol. 29, № 252.-P. 45-51.

181. Asano Naoki. A penalty function type of virtual work principle for impact contact problems of two bodies // Bull. JSME. 1986. - Vol. 29. -№ 257. - P. 15 -24.

182. Ascione Luigi, Bruno Domenico, Olivito Renato S. On dynamical behavior of plates in unilateral contact with an elastic foundation: a finite element approach // Atti Accad. naz. Lincei. Rend.Cl. sci. fis., mat. e. natur. -1984. V. 76, № 2. - P. 93-106.

183. Geller J.M. Beitrag zur theorie veranderlich gegliederten und gestutiten // Sisteme. Der. Eisenbau. 1922. - № 8, 9. - S. 37 - 46.

184. Fichera G. Problemi elastostatici con vincoli unilaterali: in problema di Signorini con ambigue condizioni al contorno // Mem. Accad. Naz. Lei. Lincei. 1964. -№ 8. - P. 91-140.

185. Goldfard D., Idnani A. A numerically stable dual method for solving strictly convex quadratic programs // Mathematical Programming. 1983. — №27.-P. 1-33.

186. Hetz H. Uber die Beriihrung fester elastischer Korper // Zs. f. Math. (Grelle). — 1881. Bd 92.-S. 31-39.

187. Powell M.J. D. TOLMIN: A fortran package for linearly constraihed optimization calculations: DAMTP Report NA2, Cambridge, University of Cambridge, (England), 1989. P. 65-79.

188. Signorini A. Questioni di elasticity non linearizzata о semilinearizzata //Rend. Mat. 1959.-V. 18, №4.-P. 95-131.

189. Schwieger H. A simple calculation of the transverse impact on beams and its experimental verification // Exptl. Mech. 1965. — Vol. 5, №11.- P. 73-82.

190. Schwieger H. Die maximale mechanische Beanspruchung beim Zentralen. Balkenquerstob // Z. Angew. Math, und Mech. 1965. - Bd 45, № 7, 8. -S. 45-58.

191. Sears J. E. On the longitudinal impact of metal rods with rounded ends // Proc. Cambrindge Phil. Soc. 1908. - Vol. 14. -P. 11-19.-329