Содержание к диссертации
Введение
1. Современные проблемы механики контактного взаимодействия 17
1.1. Классические гипотезы, применяемые при решении контактньгх задач ддя гладких тел 17
1.2. Влияние ползучести твердых тел на их формоизменение в области контакта 18
1.3. Оценка сближения шероховатых поверхностей 20
1.4. Анализ контактного взаимодействия многослойных конструкций 27
1.5. Взаимосвязь механики и проблем трения и изнашивания 30
1.6. Особенности применения моделирования в трибологии 31
Выводы по первой главе 35
2. Контактное взаимодействие гладких цилиндрических тел 37
2.1. Решение контактной задачи для гладких изотропных диска и пластины с цилиндрической полостью 37
2.1.1. Общие формулы 38
2.1.2. Вывод краевого условия для перемещений в области контакта 39
2.1.3. Интегральное уравнение и его решение 42
2.1.3.1. Исследование полученного уравнения 4 5
2,1.3.1.1. Приведение сингулярного интегро-дифференциального уравнения к интегральному уравнению с ядром, имеющим логарифмическую особенность 46
2.1.3.1.2. Оценка нормы линейного оператора 49
2.1.3.2. Приближенное решение уравнения 51
2.2. Расчет неподвижного соединения гладких цилиндрических тел 58
2.3. Определение перемещения в подвижном соединении цилиндрических тел 59
2.3.1. Решение вспомогательной задачи для упругой плоскости 62
2.3.2. Решение вспомогательной задачи для упругого диска 63
2.3.3. Определение максимального нормального радиального перемещения 64
2.4. Сопоставление теоретических и экспериментальных данных исследования контактных напряжений при внутреннем касании цилиндров близких радиусов 68
2.5. Моделирование пространственного контактного взаимодействия системы соосных цилиндров конечных размеров 72
2.5.1. Постановка задачи 73
2.5.2. Решение вспомогательных двумерных задач 74
2.5.3. Решение исходной задачи 75
Выводы и ос новные результаты второй главы 7 8
3. Контактные задачи для шероховатых тел и их решение с помощью корректировки кривизны деформированной поверхности 80
3.1. Пространственная нелокальная теория. Геометрические предположения 83
3.2. Относительное сближение двух параллельных кругов, определяемое деформацией шероховатости 86
3.3. Метод аналитической оценки влияния деформирования шероховатости 88
3.4. Определение перемещений в области контакта 89
3.5. Определение вспомогательных коэффициентов 91
3.6. Определение размеров эллиптической области контакта 96
3.7. Уравнения для определения области контакта близкой к круговой 100
3.8. Уравнения для определения области контакта близкой к линии 102
3.9. Приближенное определение коэффициента а в случае области контакта в виде круга или полосы
3.10. Особенности усреднения давлений и деформаций при решении двумерной задачи внутреннего контакта шероховатых цилиндров близких радиусов 1и5
3.10.1. Вывод интегро-дифференциального уравнения и его решение в случае внутреннего контакта шероховатых цилиндров 10"
3.10.2. Определение вспомагательных коэффициентов
3.10.3. Напряженная посадка шероховатых цилиндров
Выводы и основные результаты третьей главы
4. Решение контактных задач вязкоупругости для гладких тел
4.1. Основные положения
4.2. Анализ принципов соответствия
4.2.1. Принцип Вольтерра
4.2.2. Постоянный коэффициент поперечного расширения при деформации ползучести 123
4.3. Приближенное решение двумерной контактной задачи линейной ползучести для гладких цилиндрических тел
4.3.1. Общий случай операторов вязкоупругости
4.3.2. Решение для монотонно возрастающей области контакта 128
4.3.3. Решение для неподвижного соединения 129
4.3.4. Моделирование контактного взаимодействия в случае
однородно стареющей изотропной пластины 130
Выводы и основные результаты четвертой главы 135
5. Ползучесть поверхности 136
5.1. Особенности контактного взаимодействия тел с низким пределом текучести 137
5.2. Построение модели деформирования поверхности с учетом ползучести в случае эллиптической области контакта 139
5.2.1. Геометрические предположения 140
5.2.2. Модель ползучести поверхности 141
5.2.3. Определение средних деформаций шероховатого слоя и средних давлений 144
5.2.4. Определение вспомогательных коэффициентов 146
5.2.5. Определение размеров эллиптической области контакта 149
5.2.6. Определение размеров круговой области контакта 152
5.2.7. Определение ширины области контакта в виде полосы 154
5.3. Решение двумерной контактной задачи для внутреннего касания
шероховатых цилиндров с учетом ползучести поверхности 154
5.3.1. Постановка задачи для цилиндрических тел. Интегро-
дифференциальное уравнение 156
5.3.2. Определение вспомагательных коэффициентов 160
Выводы и основные результаты пятой главы 167
6. Механика взаимодействия цилиндрических тел с учетом наличия покрытий 168
6.1. Вычисление эффективных модулей в теории композитов 169
6.2. Построение самосогласованного метода вычисления эффективных коэффициентов неоднородных сред с учетом разброса физико-механических свойств 173
6.3. Решение контактной задачи для диска и плоскости с упругим композиционным покрытием на контуре отверстия 178
6.3. 1 Постановка задачи и основные формулы 179
6.3.2. Вывод краевого условия для перемещений в области контакта 183
6.3.3. Интегральное уравнение и его решение 184
6.4. Решение задачи в случае ортотропного упругого покрытия с цилиндрической анизотропией 190
6.5. Определение влияния вязкоупругого стареющего покрытия на изменение параметров контакта 191
6.6. Анализ особенностей контактного взаимодействия многокомпонентного покрытия и шероховатости диска 194
6.7. Моделирование контактного взаимодействия с учетом тонких металлических покрытий 196
6.7.1. Контакт шара с пластическим покрытием и шероховатого полупространства 197
6.7.1.1. Основные гипотезы и модель взаимодействия твердых тел 197
6.7.1.2. Приближенное решение задачи 200
6.7.1.3. Определение максимального контактного сближения 204
6.7.2. Решение контактной задачи для шероховатого цилиндра и тонкого металлического покрытия на контуре отверстия 206
6.7.3. Определение контактной жесткости при внутреннем контакте цилиндров 214
Выводы и основные результаты шестой главы 217
7. Решение смешанных краевых задач с учетом износа поверхностей взаим одействующих тел 218
7.1. Особенности решения контактной задачи с учетом изнашивания поверхностей 219
7.2. Постановка и решение задачи в случае упругого деформирования шероховатости 223
7.3. Метод теоретической оценки износа с учетом ползучести поверхности 229
7.4. Метод оценки износа с учетом влияния покрытия 233
7.5. Заключительные замечания по постановке плоских задач с учетом износа 237
Выводы и основные результаты седьмой главы 241
Заключение 242
Список использованных источников
- Влияние ползучести твердых тел на их формоизменение в области контакта
- Интегральное уравнение и его решение
- Относительное сближение двух параллельных кругов, определяемое деформацией шероховатости
- Приближенное решение двумерной контактной задачи линейной ползучести для гладких цилиндрических тел
Введение к работе
Актуальность темы диссертации. В настоящее время значительные усилия инженеров в нашей стране и за рубежом направлены на поиск путей определения контактных напряжений взаимодействующих тел, так как для перехода от расчета изнашивания материалов к задачам конструкционной износостойкости решающую роль имеют контактные задачи механики деформируемого твердого тела.
Следует отметить, что наиболее широкие исследования контактного взаимодействия выполнены с помощью аналитических методов. При этом применение численных методов значительно расширяет возможности анализа напряженного состояния в области контакта с учетом свойств поверхностей шероховатых тел.
Необходимость учета структуры поверхности объясняется тем, что выступы, образующиеся при технологической обработке имеют различное распределение высот и касание микронеровностей происходит только на отдельных площадках, образующих фактическую площадь контакта. Поэтому при моделировании сближения поверхностей необходимо использовать параметры, характеризующие реальную поверхность.
Громоздкость математического аппарата, применяемого при решении контактных задач для шероховатых тел, необходимость использования мощных вычислительных средств существенно сдерживают применение имеющихся теоретических разработок при решении прикладных задач. И, несмотря на достигнутые успехи, пока трудно получить удовлетворительные результаты с учетом особенностей макро- и микрогеометрии поверхностей взаимодействующих тел, когда элемент поверхности, на котором устанавливаются характеристики шероховатости твердых тел, соизмерим с областью контакта.
Все это требует разработки единого подхода к решению контактных задач, наиболее полно учитывающего как геометрию взаимодействующих тел, микрогеометрические и реологические характеристики поверхностей, характеристики их износостойкости, так и возможность получения приближенного решения поставленной задачи с наименьшим количеством независимых параметров.
Контактные задачи для тел с круговыми границами составляют теоретическую основу расчета таких элементов машин как подшипники, шарнирные соединения, соединения с натягом. Поэтому данные задачи обычно выбираются в качестве модельных при проведении подобных исследований.
Интенсивные работы, проводившиеся в последние годы в Белорусском национальном техническом унивеуиіі сі е. Дьшк іиіікишеньї
на решение этой проблемы и составляют осної у настдзддодоод^ы.
Связь работы с круппыми научными программами, темами.
Исследования выполнены в соответствии со следующими темами: "Разработать метод расчета контактных напряжений при упругом контактном взаимодействии цилиндрических тел, не описываемом теорией Герца" (Министерство образования РБ, 1997 г., № ГР 19981103); "Влияние микронеровностей соприкасающихся поверхностей на распределение контактных напряжений при взаимодействии цилиндрических тел, имеющих близкие по величине радиусы" (Белорусский республиканский фонд фундаментальных исследований, 1996 г., № ГР 19981496); "Разработать метод прогнозирования износа опор скольжения с учетом топографических и реологических характеристик поверхностей взаимодействующих деталей, а также наличия антифрикционных покрытий" (Министерство образования РБ, 1998 г., № ГР 1999929); "Моделирование контактного взаимодействия деталей машин с учетом случайности реологических и геометрических свойств поверхностного слоя" (Министерство образования РБ, 1999 г. №ГР2000Г251)
Цель и задачи исследования. Разработка единого метода теоретического прогнозирования влияния геометрических, реологических характеристик шероховатости поверхности твердых тел и наличия покрытий на напряженное состояние в области контакта, а также установление на этой основе закономерностей изменения контактной жесткости и износостойкости сопряжений на примере взаимодействия тел с круговыми границами.
Для достижения поставленной цели требуется решить следующие проблемы:
Разработать метод приближенного решения задач теории упругости и вязкоупругости о контактном взаимодействии цилиндра и цилиндрической полости в пластине с использованием мипимального количества независимых параметров.
- Разработать нелокальную модель контактного взаимодействия тел
с учетом микрогеометрических, реологических характеристик
поверхностей, а также наличия пластических покрытий.
- Обосновать подход, позволяющий корректировать кривизну
взаимодействующих поверхностей за счет деформации шероховатости.
Разработать метод приближенного решения контактных задач для диска и изотропного, ортотропного с цилиндрической анизотропией и вязкоупругого стареющего покрытий на отверстии в пластине с учетом их поперечной деформируемости.
Построить модель и определить влияние микрогеометрических особенностей поверхности твердого тела на контактное взаимодействие с пластическим покрытием на контртеле.
- Разработать метод решения задач с учетом износа цилиндрических тел, качества их поверхностей, а также наличия антифрикционных покрытий.
Объектом и предметом исследования являются неклассические смешанные задачи теории упругости и вязкоупругости для тел с круговыми границами с учетом нелокальности топографических и реологических характеристик их поверхностей и покрытий, на примере которых в настоящей работе разработан комплексный метод анализа изменения напряженного состояния в области контакта в зависимости от показателей качества их поверхностей.
Гипотеза. При решении поставленных граничных задач с учетом качества поверхности тел используется феноменологический подход, согласно которому деформация шероховатости рассматривается как деформация промежуточного слоя.
Задачи с изменяющимися во времени краевыми условиями рассматриваются как квазистатические.
Методология и методы проведенного исследования. При проведении исследований использовались основные уравнения механики деформируемого твердого тела, трибологии, функционального анализа. Разработан и обоснован метод, позволяющий корректировать кривизну нагруженных поверхностей за счет деформаций микронеровностей, что существенно упрощает проводимые аналитические преобразования и позволяет получить аналитические зависимости для размера площади контакта и контактных напряжений с учетом указанных параметров без использования предположения о малости величины базовой длины измерения характеристик шероховатости относительно размеров области контакта.
При разработке метода теоретического прогнозирования износа поверхностей наблюдаемые макроскопические явления рассматривались как результат проявления статистически усредненных связей.
Достоверность полученных в работе результатов подтверждается сравнениями полученных теоретических решений и результатов экспериментальных исследований, а также сравнением с результатами некоторых решений, найденных другими методами.
Научная новизна и значимость полученных результатов. Впервые на примере контактного взаимодействия тел с круговыми границами проведено обобщение исследований и разработан единый метод комплексного теоретического прогнозирования влияния нелокальных геометрических, реологических характеристик шероховатых поверхностей взаимодействующих тел и наличия покрытий на напряженное состояние, контактную жесткость и износостойкость сопряжений.
Комплекс проведенных исследований позволил представить в диссертации теоретически обоснованный метод решения задач механики твердого тела, основанный на последовательном рассмотрении макроскопически наблюдаемых явлений, как результата проявления статистически усредненных по значительному участку контактной поверхности микроскопических связей.
В рамках решения поставленной проблемы:
- Предложена пространственная нелокальная модель контактного
взаимодействия твердых тел с изотропной шероховатостью поверхности.
Разработан метод определения влияния характеристик поверхности твердых тел на распределение напряжений.
Исследовано интегро-дифференциальное уравнение, получаемое в контактных задачах для цилиндрических тел, что позволило определить условия существования и единственности его решения, а также точность построенных приближений.
Практическая (экономическая, социальная) значимость полученных результатов. Результаты теоретического исследования доведены до приемлемых для практического использования методик и могут быть непосредственно применены при проведении инженерных расчетов подшипников, опор скольжения, зубчатых передач. Использование предлагаемых решений позволит сократить время создания новых машиностроительных конструкций, а также с большой точностью прогнозировать их служебные характеристики.
Некоторые результаты выполненных исследований были внедрены на Н П П «Циклопривод», НПО «Алтех».
Основные положения диссертации, выносимые на защиту:
- Приближенное решите задачи механики деформированного
твердого тела о контактном взаимодействии гладких цилиндра и
цилиндрической полости в пластине, с достаточной точностью
описывающих исследуемое явление при использовании минимального
количества независимых параметров.
Решение нелокальных краевых задач механики деформируемого твердого тела с учетом геометрических и реологических характеристик их поверхностей на основе метода, позволяющего корректировать кривизну взаимодействующих поверхностей за счет деформации шероховатости. Отсутствие предположения о малости геометрических размеров базовых длин измерения шероховатости по сравнению с размерами области контакта позволяет переходить к разработке многоуровневых моделей деформирования поверхности твердых тел.
Построение и обоснование метода расчета перемещений границы цилиндрческих тел, обусловленных деформацией поверхпостных слоев. Полученные результаты позволяют разработать теоретический подход,
определяющий контактную жесткость сопряжений с учетом совместного влияния всех особенностей состояния поверхностей реальных тел.
- Моделирование вязкоупругого взаимодействия диска и полости в
пластине из стареющего материала, простота реализации результатов
которого позволяет использовать их для широкого круга прикладных
задач.
Приближенное решение контактных задач для диска и изотропного, ортотропного с цилиндрической анизотропией, а также вязкоупругого стареющего покрытий на отверстии в пластине с учетом их поперечной деформируемости. Это дает возможность оценить влияние композиционных покрытий с низким модулем упругости на нагруженность сопряжений.
Построение нелокальной модели и определение влияния характеристик шероховатости поверхности твердого тела на контактное взаимодействие с пластическим покрытием на контртеле.
Разработка метода решения краевых задач с учетом износа цилиндрических тел, качества их поверхностей, а также наличия антифрикционных покрытий. На этой основе предложена методология, сосредотачивающая математические и физические методы при исследовании износостойкости, что дает возможность вместо исследований реальных узлов трения делать основной упор на исследовании явлений, происходящих в области контакта.
Личный вклад соискателя. Все результаты, выносимые на защиту, получены автором лично.
Апробация результатов диссертации. Результаты исследований, приведенных в диссертации были представлены на 22 международных конференциях и конгрессах, а также конференциях стран СНГ и республиканских, среди них: "Понтрягинские чтения - 5" (Воронеж, 1994, Россия), "Математические модели физических процессов и их свойства" (Таганрог, 1997, Россия), Nordtrib'98 (Ebeltoft, 1998, Дания), Numerical mathematics and computational mechanics - "NMCM'98" (Miskolc, 1998, Венгрия), "Modelling'98" (Praha, 1998, Чехия), 6th International Symposium on Creep and Coupled Processes (Bialowieza, 1998, Польша), "Вычислительные методы и производство: реальность, проблемы, перспективы" (Гомель, 1998, Беларусь), "Полимерные композиты 98" (Гомель, 1998, Беларусь), "Mechanika'99" (Kaunas, 1999, Литва), П Белорусский конгресс по теоретической и прикладной механике (Минск, 1999, Беларусь), Internat. Conf. On Engineering Rheology, ICER'99 (Zielona Gora, 1999, Польша), "Проблемы прочности материалов и сооружений на транспорте" (Санкт-Петербург, 1999, Россия), International Conference on Multifield Problems (Stuttgart, 1999, Германия).
Опубликованность результатов. По материалам диссертации
опубликована 40 печатных работ, среди них: 1 монография, 19 статей в журналах и сборниках, в том числе 15 статей под личным авторством. Общее количество страниц опубликованных материалов составляет 370.
Структура и объем диссертации. Диссертация состоит из введения, семи глав, заключения, списка использованных источников и приложения. Полный объем диссертации составляет 2-М" страниц, в том числе объем, занимаемый иллюстрациями - 14 страниц, таблицами - 1 страницу. Количество использованных источников включает 310 наименований.
Влияние ползучести твердых тел на их формоизменение в области контакта
Практическое получение аналитических зависимостей для напряжений и перемещений в замкнутой форме для реальных объектов даже в простейших случаях сопряжено с существенными трудностями. Вследствие этого при рассмотрении контактных задач принято прибегать к идеализации. Так, считается, что если размеры самих тел достаточно велики по сравнению с размерами области контакта, то напряжения в этой зоне слабо зависят от конфигурации тел вдали от области контакта, а также способа их закрепления. При этом напряжения с достаточно хорошей степенью достоверности можно вычислить, рассматривая каждое тело как бесконечную упругую среду, ограниченную плоской поверхностью, т.е. как упругое полупространство [41, 65,142,146,248].
Поверхность каждого из тел предполагается топографически гладкой на микро- и макроуровне. На микроуровне это означает отсутствие или неучет микронеровностей контактирующих поверхностей, которые обусловили бы неполное прилегание поверхностей контакта. Поэтому реальная область контакта, которая образуется на вершинах выступов, значительно меньше теоретической. На макроуровне профили поверхностей считаются непрерывными в зоне контакта вместе со вторыми производными [65, 297].
Указанные предположения впервые были использованы Герцем при решении контактной задачи. Получаемые на основе его теории результаты удовлетворительно описывают деформированное состояние идеально упругих тел в отсутствии трения по поверхности контакта, однако неприменимы, в частности, к низкомодульным материалам. Кроме того, условия, в которых используется теория Герца, нарушаются при рассмотрении контакта согласованных поверхностей. Это объясняется тем, что вследствие приложения нагрузки размеры области контакта быстро растут и могут достигать величин, сравнимых с характерными размерами контактирующих тел, так что тела не могут рассматриваться как упругие полупространства [65, 210].
Особый интерес при решении контактных задач вызывает учет сил трения. Вместе с тем последнее на поверхности раздела двух тел согласованной формы, находящихся в условиях нормального контакта;, играет роль только при относительно высоких значениях коэффициента трения [69, 210].
Развитие теории контактного взаимодействия твердых тел связано с отказом от перечисленных выше гипотез. Оно осуществлялось по следующим основным направлениям: усложнение физической модели деформирования твердых тел и (или) отказом от гипотез гладкости и однородности их поверхностей.
Интерес к ползучести резко возрос в связи с развитием техники. В числе первых исследователей, обнаруживших явление деформирования материалов во времени при постоянной нагрузке, были Вика, Вебер, Кольрауш [181]. Максвелл впервые представил закон деформирования во времени в виде дифференциального уравнения. Несколько позднее Болыгман [257, 258] создал общий аппарат для описания явлений линейной ползучести. Этот аппарат, значительно развитый впоследствии Вольтерра [305, 306] является в настоящее время классическим разделом теории интегральных уравнений.
До середины прошлого столетия элементы теории деформирования материалов во времени находили малое применение в практике расчетов инженерных конструкций [174 - 178]. Однако с развитием энергетических установок, химико-технологических аппаратов, работающих при более высоких температурах и давлениях, стал необходим учет явления ползучести. Запросы машиностроения привели к огромному размаху экспериментальных и теоретических исследований в области ползучести [95, 182, 220]. Вследствие появившейся необходимости в точных расчетах явление ползучести стали учитывать даже в таких материалах, как древесина и грунты [12, 96],
Изучение ползучести при контактном взамодействии твердых тел важно по ряду причин прикладного и принципиального характера. Так, даже при постоянных нагрузках форма взаимодействующих тел и их напряженное состояние, как правило, изменяется [65, 243, 244], что необходимо учитывать при проектировании машин.
Качественное объяснение происходящих при ползучести процессов можно дать, опираясь на основные представления теории дислокаций [95, 197, 295, 296]. Так, в строении кристаллической решетки могут встречаться различные местные дефекты. Эти дефекты называются дислокациями. Они перемещаются, взаимодействуют друг с другом и вызывают различного типа скольжения в металле. Результатом движения дислокации является сдвиг на одно межатомное расстояние [95, 197, 299]. Напряженное состояние тела облегчает движение дислокаций, снижая потенциальные барьеры [95].
Временные законы ползучести зависят от структуры материала, которая меняется с течением ползучести. Экспериментально получена экспоненциальная зависимость скоростей установившейся ползучести от напряжений при относительно высоких напряжениях (-10" и более от модуля упругости). В значительном интервале напряжений экспериментальные точки на логарифмической сетке обычно группируются около некоторой прямой линии. Это означает, что в рассматриваемом интервале напряжений (-10" -10" от модуля упругости) имеется степенная зависимость скоростей деформаций от напряжения [95]. Следует отметить, что при низких напряжениях (10" и менее от модуля упругости) эта зависимость линейная. В ряде работ приведены различные экспериментальные данные по механическим свойствам различных материалов в широком интервале температур и скоростей деформирования [25, 27, 99, 127, 170, 201, 218, 253, 255, 277, 299].
Интегральное уравнение и его решение
Отметим, что если упругие постоянные диска и пластины равны, то ух= О и данное уравнение становится интегральным уравнением первого рода. Особенности теории аналитических функций позволяют в этом случае, используя дополнительные условия, получить единственное решение [46, 145, 146]. Это так называемые формулы обращения сингулярных интегральных уравнений, позволяющие получить решение поставленной задачи в явном виде. Особенность состоит в том, что в теории краевых задач обычно рассматриваются три случая (когда V составляет часть границы тел): решение имеет особенность на обоих концах области интегрирования; решение имеет особенность на одном из концов области интегрирования, а на втором обращается в ноль; решение обращается в ноль на обоих концах. В зависимости от выбора того или иного варианта строится общий вид решения, в состав которого в первом случае входит общее решение однородного уравнения. Задаваясь поведением решения на бесконечности и угловых точках области контакта, исходя из физически обоснованных предположений, строится единственное решение, удовлетворяющее указанным ограничениям [46,145].
Таким образом, единственность решения указанной задачи понимается в смысле принятых ограничений. Следует отметить, что при решении контактных задач теории упругости наиболее распространенными ограничениями являются требования обращения в ноль решения на концах области контакта и предположение об исчезновении напряжений и вращений на бесконечности [145, 147, 245]. В случае, когда область интегрирования составляет всю границу области (тела), то единственность решения гарантируется формулами Коши [145]. При этом наиболее простым и распространенным методом решения прикладных задач в этом случае является представление интеграла Коши в виде ряда [145].
Следует отметить, что в приведенных выше общих сведениях из теории сингулярных интегральных уравнений никак не оговариваются свойства контуров исследуемых областей, т.к. в данном случае известно, что дуга окружности (кривая, вдоль которой выполняется интегрирование) удовлетворяет условию Ляпунова [146]. Обобщение теории двумерных краевых задач в случае более общих предположений на гладкость границы областей можно найти в монографии ИИ. Данилюка [59].
Наибольший интерес представляет общий случай уравнения, когда 7i 0. Отсутствие методов построения точного решения в этом случае приводит к необходимости применения методов численного анализа и теории приближений. Фактически, как это уже отмечалось, численные методы решения интегральных уравнений обычно основаны на аппроксимации решения уравнения функционалом определенного вида. Объем накопленных результатов в этой области позволяет выделить основные критерии, по которым эти методы обычно сравниваются при их использовании в прикладных задачах [158]. Прежде всего простота физической аналогии предлагаемого подхода (обычно это в том или ином виде метод суперпозиции системы определенных решений); объем необходимых подготовительных аналитических вычислений, используемых для получения соответствующей системы линейных уравнений; необходимый размер системы линейных уравнений для достижения требуемой точности решения; использование численного метода решения системы линейных уравнений, максимально учитывающего особенности ее структуры и, соответственно, позволяющего с наибольшей скоростью получить численный результат [158]. Следует отметить, что последний критерий играет существенную роль лишь в случае систем линейных уравнений большого порядка. Все это определяет эффективность используемого подхода. Вместе с тем следует констатировать, что к настоящему времени существуют лишь отдельные исследования, посвященные сравнительному анализу и возможным упрощениям при решении практических задач с помощью различных аппроксимаций [210].
Отметим, что интегро-дифферешщальное уравнение может быть приведено к виду: V дуга окружности единичного радиуса, заключенная между двумя точками с угловыми координатами -сс0 и а0, а0 є(0,л/2); у1 - вещественный коэффициент, определяемый упругими характеристиками взаимодействующих тел (2.6); f(t)- известная функция, определяемая приложенными нагрузками (2.6). Кроме того, напомним, что стг(т) обращается в нуль на концах отрезка интегрирования.
Относительное сближение двух параллельных кругов, определяемое деформацией шероховатости
Задача о внутреннем сжатии круговых цилиндров близких радиусов впервые была рассмотрена И.Я. Штаерманом. При решении поставленной им задачи [248] принято, что внешняя нагрузка, действующая на внутренний и внешний цилиндры по их поверхностям, осуществляется в виде нормального давления, диаметрально противоположного давлению контакта. При выводе уравнения задачи использовано решение о сжатии цилиндра двумя противоположными силами и решение аналогичной задачи для внешности кругового отверстия в упругой среде [210, 248]. Им было получено явное выражение для перемещений точек контура цилиндра и отверстия через интегральный оператор от функции напряжений. Это выражение использовалось рядом авторов для оценки контактной жесткости [210, 138].
Используя эвристическую аппроксимацию для распределения контактных напряжений для схемы И.Я. Штаермана, А.Б. Милов получил упрощенную зависимость для максимальных контактных перемещений. Однако им было установлено, что полученная теоретическая оценка существенно отличается от экспериментальных данных. Так, перемещение, определенное из эксперимента, оказалось меньше теоретического в 3 раза. Этот факт объясняется автором существенным влиянием особенностей пространственной схемы нагружения и предлагается коэффициент перехода от трехмерной задачи к плоской [138].
Аналогичный подход использовал М.И. Теплый, задавшись приближенным решением несколько иного вида [210]. Следует отметить, что в этой работе, кроме того, получено линейное дифференциальное уравнение второго порядка для определения контактных перемещений в случае схемы, приведенной на Рисунке 2.1. Указанное уравнение следует непосредственно из способа получения интегро-дифференциального уравнения для определения нормальных радиальных напряжений. При этом сложность правой части определяет громоздкость результирующего выражения для перемещений. Кроме того, в этом случае остаются неизвестными величины коэффициентов в решении соответствующего однородного уравнения. Вместе с тем отмечается, что, не устанавливая значений постоянных, можно определить сумму радиальных перемещений диаметрально противоположных точек контуров отверстия и вала [210].
Таким образом, несмотря на актуальность задачи определения контактной жесткости анализ литературных источников не позволил выявить метода ее решения, позволяющего обоснованно установить величины наибольших нормальных контактных перемещений, обусловленных деформацией поверхностных слоев без учета деформаций взаимодействующих тел в целом, что объясняется отсутствием формализованного определения понятия "контактная жесткость".
При решении поставленной задачи будем исходить из следующих определений: перемещения под действием главного вектора сил (без учета особенностей контактного взаимодействия) будем называть сближение (удаление) центра диска (отверстия) и его поверхности, не приводящее к изменению формы его границы. Т.е. это жесткость тела в целом. Тогда контактная жесткость это максимальные перемещения центра диска (отверстия) без учета перемещения упругого тела под действием главного вектора сил. Данная система понятий позволяет разделить перемещения;, полученные из решения задачи теории упругости, и показывает, что оценка контактной жесткости цилиндрических тел, полученная А.Б. Миловьш из решения ИЛ. Штаермана [248], верна только для данной схемы нагружения.
Рассмотрим задачу, поставленную в п. 2.1. (Рисунок 2.1) с краевым условием (2.3). Учитывая свойства аналитических функций, из (2.2) имеем, что [145]:
Важно подчеркнуть, что первые слагаемые (2.30) и (2.32) определяются решением задачи о сосредоточенной силе в бесконечной области. Это объясняет наличие логарифмической особенности. Вторые слагаемые (2.30), (2,32) определяются отсутствием касательных напряжений на контуре диска и отверстия;, а также условием аналитического поведения соответствующих слагаемых комплексного потенциала в нуле и на бесконечности. С другой стороны суперпозиция (2.26) и (2.29) ((2.27) и (2.31)) дает нулевой главный вектор сил, действующих на контур отверстия (или диска). Все это позволяет выразить через третье слагаемое величину радиальных перемещений в произвольном фиксированном направлении С, в пластине и в диске. Для этого найдем разность Фпд(г), (z) и Фп 2(2), 4V2(z):
Приближенное решение двумерной контактной задачи линейной ползучести для гладких цилиндрических тел
Идея о необходимости учета микроструктуры поверхности сжимаемых тел принадлежит И.Я. Штаерману [248]. Им введена модель комбинированного основания, согласно которой в упругом теле, кроме перемещений, вызванных действием нормального давления и определяемых решением соответствующих задач теории упругости, возникают дополнительные нормальные перемещения, обусловленные чисто местными деформациями, зависящими от микроструктуры контактирующих поверхностей. И.Я.Штаерман предположил, что дополнительное перемещение пропорционально нормальному давлению, причем коэффициент пропорциональности является для данного материала величиной постоянной. В рамках этого подхода им впервые было получено уравнение плоской контактной задачи для упругого шероховатого тела, т.е. тела, имеющего слой повышенной податливости.
В ряде работ предполагается, что дополнительные нормальные перемещения за счет деформации микровыступов контактирующих тел пропорциональны макронапряжению в некоторой степени [130, 168, 173, 210]. Это основано на приравнивании усредненных значений перемещений и напряжений в пределах базовой длины измерения шероховатости поверхности. Однако, несмотря на достаточно хорошо разработанный аппарат решения задач подобного класса, ряд трудностей методического характера не преодолен. Так, используемая гипотеза о степенной связи напряжений и перемещений поверхностного слоя с учетом реальных характеристик микрогеометрии верна при малых базовых длинах, т.е. высокой чистоте поверхности, а, следовательно, при справедливости гипотезы о топографической гладкости на микро и макроуровне [60]. Следует также отметить существенное усложнение уравнения при использовании подобного подхода и невозможность описания с его помощью влияния волнистости.
Несмотря на достаточно хорошо разработанный аппарат решения контактных задач с учетом слоя повышенной податливости [3, 50, 130, 173, 205, 210, 272, 301], остался ряд вопросов методического характера, затрудняющих его применение в инженерной практике расчетов. Как уже отмечалось, шероховатость поверхности имеет вероятностное распределение высот. Соизмеримость размеров элемента поверхности, на котором определяются характеристики шероховатости, с размерами области контакта является главной трудностью при решении поставленной задачи и определяет некорректность применения некоторыми авторами непосредственной связи между макродавлениями и деформациями шероховатости в виде: где s - точка поверхности.
Следует отметить также решение поставленной задачи с использованием предположения о трансформации вида распределения давления в параболический, если деформациями упругого полупространства в сравнении с деформациями шероховатого слоя можно пренебречь. Этот подход приводит к существенному усложнению интегрального уравнения и позволяет получать только численные результаты. Кроме того, авторами использовалась уже упомянутая гипотеза (3.1).
Необходимо упомянуть, попытку разработки инженерного метода учета влияния шероховатости при внутреннем касании цилиндрических тел [235], основанного на предположении о том, что упругие радиальные перемещения в области контакта, обусловленные деформацией микро-неровности, постоянны и пропорциональны среднему контактному напряжению т в некоторой степени к. Однако, несмотря на свою очевидную простоту, недостатком этого подхода является то, что при таком способе учета шероховатости ее влияние постепенно возрастает с возрастанием нагрузки, что не наблюдается на практике (Рисунок 3 Л ,[65, ПО]).
