какая зависимость является исходной для решения контактной задачи
Контактные задачи железнодорожного транспорта
Содержание материала
Нормальная задача. Решение Герца.
Первое решение задачи о контакте двух твердых упругих тел при их статическом сжатии дано Герцем в 1882 г., что не было непосредственной целью его исследований [24]. Наблюдая картины интерференции, возникающие между поверхностями стеклянных линз, он пришел к выводу, что, если поверхности тел описываются уравнениями второго порядка, пятно, возникающего между ними контакта, имеет эллиптическую форму. Опираясь на знание теории электростатического потенциала, Герц воспользовался аналогией и пришел к выводу, что закон распределения контактных давлений, соответствующий такой форме пятна, должен быть эллипсоидальным.
Герц получил зависимости, позволяющие определять размеры пятен контактов и максимальные контактные давления. Они нашли исключительно широкое применение при оценке контактной прочности деталей подшипников качения, зубьев зубчатых колес, железнодорожных колес и рельсов, ходовых колес кранов, звеньев цепей и многих деталей машин. Следует отметить, что секрет такого успеха его решения не в том, что расчетные схемы этих объектов удовлетворяют допущениям, принятым при решении задачи Герцем, а в необходимости получения представления об уровне контактных напряжений с использованием простейших расчетных схем и простых инженерных методик расчета. Решения контактных задач для тел неправильных геометрических форм, клиньев, плит со скругленными ребрами и краями, с изношенными поверхностями, случаев, когда размеры контакта сопоставимы с размерами тел, когда контактирующие тела не могут быть представлены полупространствами, полученные позже, не приводят к столь простым расчетным зависимостям, в связи с чем не всегда можно воспользоваться ими.
Напряженное состояние в области контакта.
Напряженное состояние в области кругового контакта тел и случая первоначального касания цилиндров по линии исследовано в работах А. Н. Динника [25]. Позже было получено решение для эллиптического контакта в работах Н. М. Беляева [8], который сделал попытку получить решение в замкнутой форме для напряжений в любой точке области, прилегающей к эллиптическому контакту. Однако ему удалось лишь выразить напряжения через весьма сложные интегралы. Более простые зависимости, нашедшие применение в инженерных расчетах, получены для точек, лежащих на поверхности контакта и на оси z, проходящей через центр контакта внутрь тела.
Круг плоских контактных задач, поддающихся решению, значительно расширился с появлением работ Н. И. Мусхелишвили и И. Я. Штаермана[104].
Плоская тангенциальная контактная задача.
Тангенциальная задача для кругового и эллиптического контактов.
Статическая тангенциальная задача для эллиптического контакта была решена К. Коттанео [120, 121]. Пожалуй, именно в этот период обозначаются три главные проблемы, рассматриваемые при решении контактных задач: определение законов распределения нормальных и касательных сил на поверхности контакта; анализ напряженно-деформированного состояния в области, прилегающей к контакту; определение начального смещения тел при приложении касательных сил для статического контакта и крипов, спина, сил крипа для контакта качения.
Напряженное состояние в точках эллиптической площадки контакта для случая скольжения одного тела по другому в предположении закона трения Кулона, т.е. от касательных сил, распределенных по эллипсоидальному закону, впервые исследовано Б. С. Ковальским [37, 38]. Полученное решение открывало хорошие возможности для решения статических задач и задач качения с неполным скольжением. Однако оно осталось малоизвестным, и эти задачи были решены позже другими исследователями. Прежде всего нужно назвать работу Р. Д. Миндлина [163], в которой дано решение для начальных смещений тел при круговом и эллиптическом контактах при приложении к телам сдвигающих сил и момента верчения. Эта работа примечательна тем, что в ней рассмотрено распределение касательных сил в контакте для случая, когда микропроскальзывания контактирующих точек имеют место на части контакта и когда они отсутствуют. Идея рассмотреть последний случай оказалась плодотворной. На ее основе затем были построены линейные теории крипа и спина для контакта качения.
Линейная и нелинейная теории крипа.
Сначала К. Л. Джонсоном была разработана линейная теория крипа для кругового контакта качения [141], а затем линейная теория спина [142]. Позже К. Л. Джонсон с П. Д. Вермеленом получили решение для напряжений в точках эллиптического пятна контакта от касательных сил, распределенных по эллипсоидальному закону, и на основе этого решения разработали нелинейную теорию крипа [143]. Она базируется на предположении, что участок сцепления имеет форму эллипса, подобного эллипсу, ограничивающему пятно контакта, и они имеют одну точку касания на набегающем крае. Но в еще более ранней работе [141] К. Л. Джонсон обратил внимание на то, что при такой схеме деления контакта на участки сцепления и скольжения на части площадки скольжения направления проскальзываний и касательных сил совпадают, что противоречит физическому смыслу. Кроме того, не удалось создать нелинейную теорию для общего случая качения с крипами и спином.
Теория полос. Упрощенная и точная теории.
Д. Хайнес и Е. Олдертон, исследовав контакт качения методом фотоупругости, пришли к выводу, что при перекатывании в условиях, когда на поверхности контакта действуют продольные касательные силы, участок сцепления прилегает к набегающему краю, причем скольжение отсутствует не в одной точке, а на дуге набегающего края [131]. Участок сцепления имеет лимонообразную форму, сепаратриса, разделяющая участки сцепления и скольжения, симметрична относительно поперечной прямой линии дуге набегающего края. Одновременно Д. Д. Калкером была предложена теория полос. Завершение теория полос получила в работе [145]. С. В. Пуун методом фотоупругости, а К. Л. Джонсон на резиновых моделях исследовали распределения касательных напряжений в контакте, имеющем форму круга, при перекатывании со спином. С. В. Пуун пришел к выводу, что область сцепления прилегает к набегающему краю контакта и имеет форму клина [168]. Аналогичный результат получен К. Л. Джонсоном [144].
Очевидно, что с получением этих результатов для контакта качения пришло и понимание того, что эту задачу нельзя решить на основе простых аналитических подходов. Тогда Д. Д. Калкером была разработана точная теория контакта качения, разумеется, точная в рамках принятых им допущений, которая была реализована в двух программных продуктах DUVOROL и CONTACT [146]. Теоретические предпосылки точной теории построены на использовании вариационных принципов с численной реализацией на сетке элементов, покрывающей пятно контакта.
Экспериментальные исследования показали, что разработанные программы обеспечивают достаточно высокую точность решения стационарной и нестационарной задач качения. Однако они требуют больших затрат машинного времени, что не позволяет использовать их при решении задач динамики движения железнодорожного экипажа по рельсам. При интегрировании дифференциальных уравнений движения программа на каждом шаге по времени должна с использованием данных о положениях колес на рельсах и величинах крипов и спинов, по крайней мере, в доли секунды вычислить значения сил крипа.
С этой целью Д. Д. Калкером была разработана упрощенная теория, в основу которой была положена прямопропорциональная зависимость между перемещениями точек в контакте, вычисляемыми через крипы и спин, и касательными силами [147]. Наконец, на основе теории полос и упрощенной теории была предложена теория, реализованная в программе FASTSIM [148], нашедшая широкое применение в расчетах динамики движения железнодорожных экипажей.
Предложенные Д. Д. Калкером теории разрабатывались для эллиптического контакта, хотя это ограничение и не оговаривалось. При контактировании колеса и рельса с изношенными профилями пятно контакта может иметь другую форму. В этом случае оно, как правило, заменяется эквивалентным эллиптическим, как это предложено в работах [133, 154].
Быстрые алгоритмы решения нормальной и тангенциальной задач.
Разумеется, в любом случае решению тангенциальной задачи предшествует решение нормальной контактной задачи, для чего используются быстрые алгоритмы [152, 158]. Как показано в работе [152], все эти методы, в том числе и решения тангенциальной задачи с помощью программы FASTSIM, могут быть использованы и для многоточечного контакта.
Напряженное состояние от касательных нагрузок.
Напряженное состояние в точках области, прилегающей к эллиптическому контакту, вызванное действием нормальных и касательных нагрузок, исследовано М. Д. Браянтом и Л. М. Киром [112].
Законы распределения касательных сил и характер деления поверхности контакта на участки сцепления и скольжения при качении со спином исследованы С. В. Пууном и К. Л. Джонсоном экспериментальными методами [144, 168].
Д. Л. Любкиным решена задача определения закона распределения касательных напряжений в точках поверхности кругового контакта при приложении к телам момента верчения и неполном скольжении и получена зависимость угла взаимного поворота тел от момента в работе [160].
Решение контактных задач для тел с шероховатыми поверхностями.
Во всех упомянутых решениях задач и с трением, и без трения предполагалось, что поверхности контактирующих тел являются гладкими. И это допущение оказалось существенным. При более пристальном изучении поверхностей, обработанных с использованием любых технологических приемов, на них обнаруживаются неровности: волнистость, шероховатость, субшероховатость. Реальный контакт тел происходит по выступам неровностей. Фактическая площадь контакта оказывается меньше номинальной. Одной из первых работ, в которой дано решение для контакта шероховатых номинально плоских поверхностей, является работа Д. А. Гринвуда и Д. Б. П. Вильямсона [129]. Затем Д. А. Гринвуд и Д. Х. Трипп [130] рассмотрели контакт двух шероховатых шаров. Вскоре работы по исследованию контактной жесткости образовали целое научное направление. В России оно развивалось в работах И. В. Крагельского [42, 43], Н. Б. Демкина [22, 23], З. М. Левиной и Д. Н. Решетова [48], Э. В. Рыжова [76], А. Г. Суслова [91] и др. [86].
При разработке приближенных методов решения контактных задач с учетом шероховатости поверхностей предложено множество расчетных схем, моделирующих неровности: цилиндров, сфер, эллипсоидов, призм, пирамид и конусов со сферическими и усеченными вершинами, а также эмпирических зависимостей для определения жесткости стыков [48]. Для описания статистических распределений неровностей принят ряд параметров.
При использовании для решения контактных задач с учетом шероховатости конечно-элементных расчетных схем задача описания параметров неровностей упрощается. Решения контактной задачи для колеса и рельса с учетом микронеровностей их поверхностей, выполненные МКЭ, показали, что шероховатость оказывает существенное влияние на некоторые характеристики контакта. Так в работе [133] показано, что в зависимости от характеристик микронеровностей изменяется зависимость между крипами и силами крипа.
Авторы не ставят перед собой целью проследить историю развития методов решения контактных задач: упомянуты лишь немногие работы, на которые сделаны ссылки и результаты которых использованы в монографии. Более обширные сведения из этой области можно найти в литературных источниках [8, 15, 17, 18, 24, 38, 104, 149]. Значительное внимание уделено зарубежным публикациям в связи с тем, что на русский язык переведены лишь немногие из них.
ОСНОВНЫЕ УСЛОВНЫЕ ОБОЗНАЧЕНИЯ
Решение контактных задач теории пластин и плоских негерцевских контактных задач методом граничных элементов тема диссертации и автореферата по ВАК РФ 01.02.04, кандидат физико-математических наук Малкин, Сергей Александрович
Оглавление диссертации кандидат физико-математических наук Малкин, Сергей Александрович
Глава 1. Изгиб пластин сложной формы.
§1.1. Исходные соотношения и гипотезы.
§ 1.2 Метод компенсирующих нагрузок при изгибе пластин.
§1.3 Граничноэлементное представление интегральных уравнений метода компенсирующих нагрузок.
§1.4 Аналитическое вычисление интегралов по элементам контура.
§1.5 Вычисление интегралов по области пластины.
§1.6 Тестовые задачи.
Глава 2. Контактные задачи взаимодействия пластин с жесткими телами.
§2.1 Методы решения контактных задач.
§2.2 Контактные задачи с известной областью контакта.
§2.3 Контактные задачи с неизвестной областью контакта.
§2.4 Аналитическое решение задачи контакта круглой пластины с жесткой плоскостью.
§2.5 Давление жесткого гладкого штампа на пластину.
Глава 3. Плоская задача теории упругости.
§3.1 Исходные соотношения.
§3.2 Фундаментальные решения для полуплоскости.
§3.3 Метод компенсирующих нагрузок в плоской задаче теории упругости.
§3.4 Граничноэлементная постановка задач плоской теории упругости.
§3.5 Аналитическое вычисление интегралов по элементам контура.
§3.6 Тестовая задача.
Глава 4. Плоские контактные задачи теории упругости.
§4.1 Условие контакта при взаимодействии упругого и жесткого тел.
§4.2 Контакт между двумя соприкасающимися цилиндрическими телами.
§4.3 Сжатие тел, ограниченных цилиндрическими поверхностями, радиусы которых почти равны (задача И.Я. Штаермана).
§4.4 Контакт балки с жестким основанием. Решение балочной теории.
§4.5 Контакта балки с жестким основанием. Решение плоской задачи теории упругости.
§4.6 Изгиб балки под действием жесткого криволинейного штампа.
§4.7 Расчет на прочность цилиндрических зубчатых передач.
Рекомендованный список диссертаций по специальности «Механика деформируемого твердого тела», 01.02.04 шифр ВАК
Контактные задачи взаимодействия мембраны сложной формы с жестким телом и жидкостью 2001 год, кандидат физико-математических наук Чумарина, Ольга Владимировна
Большие прогибы пластин и пологих оболочек со сложным контуром 1998 год, доктор физико-математических наук Грибов, Александр Павлович
Расчет ортотропных пластин и оболочек методом граничных элементов 2008 год, кандидат физико-математических наук Великанов, Петр Геннадьевич
Математическое моделирование процессов деформирования пластин и пологих оболочек методом граничных интегральных уравнений 2003 год, кандидат физико-математических наук Куканов, Николай Иванович
Статическое и динамическое контактное взаимодействие пластин и цилиндрических оболочек с жесткими телами 1983 год, кандидат физико-математических наук Кузнецов, Сергей Аркадьевич
Введение диссертации (часть автореферата) на тему «Решение контактных задач теории пластин и плоских негерцевских контактных задач методом граничных элементов»
Конструирование современных машин и механизмов неразрывно связано с проведением многовариантных прочностных расчетов. Высокие требования, предъявляемые к надежности конструкции, в настоящее время могут быть удовлетворены лишь при условии обеспечения процесса проектирования оперативной и достоверной информацией о ее напряженно-деформированном состоянии. Расчетные схемы исследуемых конструкций при этом должны быть максимально приближены к реальным объектам, учитывать сложность их конструктивных форм, структуры, характер нагружения и взаимодействия с окружающей средой, поведение материалов конструкции в экстремальных условиях и т. д.
В машиностроительных конструкциях передача усилий обычно осуществляется посредством контакта отдельных деталей. Однако при рассмотрении узлов, состоящих из системы взаимодействующих тел, явлениями в локальной зоне контакта зачастую пренебрегают. В этом случае, руководствуясь принципом Сен-Венана, проводят упрощение и схематизацию усилий, воспринимаемых отдельными деталями, и приходят к смешанной задаче теории упругости с заданными на границе силами и смещениями. Такие упрощения расчетной схемы приемлемы далеко не всегда. В большинстве реальных конструкций закон распределения истинных контактных давлений оказывает существенное влияние на НДС взаимодействующей пары, а иногда, как, например, во фланцевых соединениях с упругими прокладками, определяет работоспособность конструкции в целом. В таких случаях возникает необходимость решения контактных задач, где размеры и конфигурация площадок контакта, условия взаимодействия на них нелинейно зависят от приложенной нагрузки. Эти параметры являются искомыми и могут быть определены только в процессе решения задачи.
Теория контактного взаимодействия включает в себя различные классы задач [230]. Среди них выделяют статические и квазистатические, где не учитываются эффекты инерции, а также контактные задачи динамики, где рассматриваются различные режимы движения взаимодействующих тел, пульсирующее, ударное нагружение и т. п. В свою очередь эти задачи подразделяются на так называемые нормальные задачи без трения, где рассматриваются идеальные односторонние связи между телами, и задачи с трением. Для ряда случаев процесс трения аппроксимируется полным сцеплением. Различным проблемам контакта посвящена обширная литература, которая обобщена в подробных обзорах [90, 94, 101,170,203].
Исторически первыми, основополагающими работами в теории контактных задач явились исследования Герца, где впервые было получено распределение местных напряжений в районе контакта упругих тел. И хотя постановка задачи предусматривала ряд серьезных допущений, таких, как малость пятна контакта, отсутствие трения, однородность, изотропность и идеальная упругость материала, результаты исследований до сих пор не потеряли своей теоретической и практической ценности.
Значительный вклад в развитие аналитических методов решения контактных задач внесли фундаментальные труды отечественных ученых — Н.И. Мусхелишвили, И.Н. Векуа, Н.П. Векуа, С.Г. Михлина, J1.A. Галина, И.Я. Штаермана, Д.И. Шермана, B.J1. Рвачева, а также работы зарубежных математиков и механиков К. Каттанео, Н. Губера, Р.Д. Миндлина, А. Синьорини. Разработанные ими методы теории функций комплексной переменной и теории сингулярных интегральных уравнений оказались достаточно эффективными для решения смешанных задач упругости. Однако круг рассмотренных примеров при этом ограничивался в основном классическими смешанными задачами о внедрении жесткого индентора (штампа) в бесконечную или полубесконечную область, цилиндрическом изгибе пластин и стержней, осесимметричном контакте пластин.
В последние годы получили развитие методы решения контактных задач теории пластин и оболочек, связанные с работами Ю.П. Артюхина, И.А. Биргера, М.В.Блоха, Э.И. Григолюка, B.C. Гудрамовича, Б.Я. Кантора, С.Н. Карасева, М.М. Кира, С.А. Кузнецова, В.Н. Максименко, В.П. Ольшанского, Б.Л. Пелеха,
Г.Я.Попова, B.C. Саркисяна, М.А.Г. Сильвы, С.П. Тимошенко, В.М. Толкачева, М.М. Филоненко-Бородича, В.Ф. Чижова, Ф.Эссенбурга и др. Однако остался неразработанным вопрос о решении контактных задач для тел произвольной формы. Объясняется это тем, что для таких тел не удается получить в аналитическом виде функцию Грина, определяющую ядра интегральных уравнений, на основе которых построен интегральный подход решения контактных задач.
Впервые задача о стесненном изгибе тонкого первоначально прямого стержня вокруг круговой опоры или о распрямлении изогнутого стержня на плоской плите приложенными к концам стержня силами рассмотрена С.П. Тимошенко [255]. Решение построено на основе классической теории изгиба балок Бернулли—Эйлера. Был обнаружен скачок в поперечной силе на конце зоны контакта или сосредоточенная сила в составе реакции. После С.П. Тимошенко задачу стесненного изгиба стержня вокруг круговой опоры рассмотрел М.М. Филоненко-Бородич [194] на основе теории изгиба балок, учитывающей деформации поперечного сдвига без учета поперечного обжатия. В решении был устранен скачок в поперечной силе и, соответственно, исчезла сосредоточенная сила в реакции. Задача цилиндрического изгиба пластин и изгиба балок жесткими штампами в постановке М.М. Филоненко-Бородича позднее рассматривалась в работах [27, 145, 158, 178]. Изгиб пластин на упругом основании с помощью штампов изучен В.М. Александровым [11], Б.Л. Пелехом и Р.Д. Сысаком [157]. Эффект поперечного обжатия в рамках прикладных теорий пластин принят во внимание в работах Ю.П. Артюхина и С.Н. Карасева [103], М.В. Блоха [50], М.В. Блоха, Н.Г. Ващенко, А.А. Гинца [51], Э.И. Григолюка и В.М. Толкачева [91]. Поперечное обжатие в указанных работах учитывалось путем интегрирования по толщине пластины соотношения закона Гука для поперечной деформации. С.О. Саркисян [178] решал контактную задачу о цилиндрическом изгибе пластины с использованием теории С.А. Амбарцумяна [13] и теории П. Нагди [244], в которых учитывается как поперечный сдвиг, так и обжатие.
Задача цилиндрического изгиба пластин штампами на основе уравнений плоской теории упругости изучалась Э.И. Григолюком и В.М. Толкачевым [91], С.Н. Карасевым и Ю.Н. Артюхиным [102], С.О. Саркисяном [178], Киром и Сильвой [234]. В работе [234] рассмотрена задача цилиндрического изгиба бесконечной пластины периодической системой жестких штампов с круговой формой основания. Решение строится в тригонометрических рядах и сведено к парным уравнениям, которые решаются численно. Эта же задача описана в статье [91]. Решение строилось с помощью функции Грина и сведено к интегральному уравнению с периодическим логарифмическим ядром. Последнее путем обращения интеграла с логарифмическим ядром сведено к уравнению Фредгольма второго рода, которое решалось численно.
Первая работа, в которой рассмотрена осесимметричная контактная задача для круговой пластины, принадлежит К. Гиркману [224]. В ней с позиции теории пластин С. Жермен— Лагранжа—Кирхгофа предполагалось, что первоначально неизогнутая пластина покоится на абсолютно жестком плоском основании и прижимается к основанию равномерно распределенной нагрузкой (вес пластины, либо давление). Под действием прикладываемых к наружному контуру пластины равномерно распределенных по окружности изгибающих моментов некоторая кольцевая зона пластины, примыкающая к контуру, может оторваться от основания. В работе найдена зависимость между величиной зоны отрыва и моментом, а также напряжения в пластине. В этой же работе обсуждается аналогичная задача для балки, а также для круглой пластины, покоящейся первоначально на податливом основании Винклера. Задача Гиркмана рассматривалась Р. Хофманом [226] также на основе теории С. Жермен— Лагранжа—Кирхгофа. Кроме этого, в статье [226] рассмотрена пластина, которая первоначально опирается на круговую жесткую стенку достаточно мелкой цилиндрической полости с плоским дном. Под действием давления или собственного веса пластина прогибается и некоторой центральной зоной может войти в контакт с дном полости. Здесь в отличие от задачи Гиркмана на контуре пластины действуют только поперечные силы и отсутствует изгибающим момент.
Контактная задача изгиба круговой пластины жестким штампом с параболическим основанием на основе теории пластин с учетом поперечного сдвига без учета поперечного обжатия исследована JI. А. Розенбергом [175].
Ряд исследований по контактным задачам для круглых пластин выполнен Ф. Эссенбургом [215, 216, 217] на основе теории пластин, учитывающей деформации поперечного сдвига без поперечного обжатия. В статье [215] рассмотрена жестко защемленная по контуру пластина, изгибаемая жестким штампом с параболической формой основания. Приведены графики изменения величины зоны контакта в зависимости от нагрузки на штамп, от величины смещения штампа, графики распределения напряжений. Все это сравнивается с результатами, вытекающими из теории С. Жермен—Лагранжа—Кирхгофа. В статье [215] рассмотрена упомянутая задача Р. Хофмана [226] (опирание на днище полости). Решение строится по теории типа С. П. Тимошенко с небольшим видоизменением — учитываются нормальные поперечные напряжения при записи соотношений обобщенного закона Гука для изгибающих моментов. Прогиб не изменяется по толщине пластины и, таким образом, обжатием пренебрегается. Показано, что учет поперечного сдвига приводит к значительному изменению нормальных изгибных и касательных напряжений в пластине по сравнению с полученными по теории С. Жермен—Лагранжа— Кирхгофа. Что касается учета нормальных напряжений в формулах для моментов, то этот учет практически не изменяет ни картину, распределения напряжений, ни характер реакции. В работе [217] исследован контакт двух пластин, которые до деформации могут либо прилегать одна к другой, либо находиться на некотором расстоянии. Совместный изгиб двух пластин с позиций теории С. П. Тимошенко проанализирован также Ю. П. Артюхиным и С. Н. Карасевым [28].
Пластины на упругом основании, сжимаемые штампами, рассмотрены в работах В.М. Александрова [11], Б. Лаинга [237], Б.Л. Пелеха и Р.Д. Сысака [156], И. Снеддона [254] и др.
Наибольшие трудности при решении контактных задач вызывает проблема определения границы области контакта. Одномерные контактные задачи для тонкостенных элементов с неизвестной границей взаимодействия имеют решение [14, 90, 91, 121, 203, 168, 25]. Если контакт происходит по неизвестной площади (двумерная задача), то таких исследований немного. Поэтому особый интерес представляют задачи отыскания двумерных областей контакта. Основополагающей здесь является работа J1.A. Галина [84], в которой рассмотрена задача о давлении жесткого эллипсоида на защемленную круглую пластинку. Размеры области контакта удалось найти путем сведения задачи к обратной краевой задаче для одной аналитической функции. Автор показал, что давление, передаваемое штампом на пластинку, будет распределено по эллиптической границе области контакта, а внутри ее будет равно нулю.
Основная идея предложенного метода состоит в переходе от известного решения классической задачи о действии штампа на упругое полупространство к приближенному решению соответствующей задачи о внедрении того же штампа в упругое тело конечных размеров. Существенным преимуществом такого подхода является получение решения в зоне действия штампа в простой аналитической форме, удобной для инженерных расчетов.
Асимптотический метод и его модификации для решения различных смешанных задач был использован в работах И.И. Воровича [7, 77], В.М. Александрова [4, 5, 9], В.А. Бабешко [43, 44] и др.
Наряду с асимптотическими существует ряд методов сведения смешанной краевой задачи к бесконечным системам алгебраических уравнений. Например, в работах В.М. Александрова [8, 10], Г.Я. Попова [166,167], B.JI. Рвачева [172, 173] и др. широко используется метод ортогональных полиномов, с помощью которого производится разложение известной функции, входящей в правую часть интегрального уравнения. Регулярная часть ядра интегрального уравнения 1 рода также раскладывается в двойной ряд, после чего уравнение сводится к алгебраической системе. В работах Б.Л. Абрамяна [2], А.А. Баблояна [45,46] и др. предложены методы непосредственного сведения краевой задачи к бесконечной алгебраической системе, минуя интегральное уравнение.
Иногда интегральные уравнения смешанных задач удается привести к конечным алгебраическим системам. Это обычно достигается путем аппроксимации регулярной части их ядер вырожденными [6] либо применением метода коллокаций [76, 100], где контактное давление представляется определенным числом параметров, для определения которых используются условия связи, налагаемые на перемещения в конечном числе точек области контакта.
Широкое распространение получили методы, основанные на сведении смешанной краевой задачи к некоторым парным или тройным функциональным (интегральным) уравнениям (или рядам), которые в итоге преобразуются в интегральное уравнение Фредгольма II рода, решаемое одним из приближенных методов. Группа данных методов представлена в работах Ю.Н. Кузьмина и Я.С. Уфлянда [127,128], А. А. Баблояна [1], А. Ф. Улитко [192] и др.
Необходимо отметить, что при решении смешанных задач указанной группой методов снимается ряд упрощающих предположений классической теории. В частности, рассматриваются контактные задачи для неоднородных, анизотропных тел, в ряде случаев производится учет трения и микроструктуры контактирующих поверхностей. Существенно и то, что исследуемая область контактного взаимодействия для задач такого типа соизмерима с характерными размерами тел [75, 170, 182, 184].
Однако все указанные решения получены для частных, относительно простых областей, реологических свойств материала и условий контактирования. При этом решение каждой отдельной задачи сопряжено с большими, а порой и непреодолимыми трудностями математического характера, что требует высокой квалификации исполнителя. Поэтому в широкой инженерной практике распространение получила лишь малая часть аналитических методов, наиболее простых с вычислительной точки зрения.
Наряду с классическими постановками контактной задачи существует ее вариационная формулировка, впервые предложенная в работе А. Синьорини [253]. Для ее применения к рассматриваемым задачам строится функционал, достигающий минимума на решении исходной задачи и, кроме того, имеющий граничные условия в качестве необходимых условий экстремума.
На поверхности раздела контактирующих тел вводятся связи специального вида, способные передавать только односторонние сжимающие усилия в направлении общей нормали к контактирующим поверхностям. Взаимные перемещения соприкасающихся тел в том же направлении не могут быть произвольными и ограничиваются условиями непроникания контактирующих тел друг в друга.
В вариационной постановке решение контактной задачи без трения сводится к проблеме минимизации функционала полной энергии системы с линейными ограничениями в виде неравенств. С точки зрения методов оптимизации — это задача квадратичного программирования и для ее решения приемлемы известные процедуры градиентного спуска [72, 74], возможных направлений [126], множителей Лагранжа [73, 114, 164] и др.
Идея использования подходов квадратичного программирования для решения контактных задач впервые была предложена в работах В.М. Фридмана и B.C. Черниной [195, 196]. В дальнейшем вопросы применения квадратичного программирования изучались в работах [97, 110—112, 114, 246, 247]. Такой подход к решению контактных задач тесно связан с использованием современных численных методов, таких, как вариационно-разностный [99, 164] метод и МКЭ [124, 125, 177, 219, 221], которые базируются на эквивалентных вариационных формулировках задачи. Причем большинство авторов отдает предпочтение МКЭ благодаря его высокой универсальности и эффективности.
В настоящее время известен ряд подходов к решению контактной задачи методом конечных элементов. Наиболее прост с алгоритмической точки зрения прием, основанный на вычислении коэффициентов взаимного влияния точек контактирующих тел в нормальном и касательном направлениях. С помощью метода сил для составления равновесия каждого тела в отдельности находится распределение контактных напряжений. Полученные значения напряжений используются в качестве граничных условий для повторного вычисления по определению напряженного состояния контактирующей пары. Границы контактных площадок и участки проскальзывания находятся итерационным путем в процессе решения задачи. Такой подход использовался в работах [95, 250, 258]. Отметим, что наряду с относительной простотой такой метод не лишен недостатков, основным из которых является необходимость решения задачи на этапе определения коэффициентов податливости 2п раз, где п — число точек контакта.
Существует еще одна группа методов решения контактной задачи МКЭ, где условия взаимодействия между телами моделируются с помощью соотношений физически нелинейных задач механики твердого тела. Первыми работами, в которых механика контакта рассматривалась по аналогии с пластическим течением, явились исследования Р. Михайловского, 3. Мроза и В. Фридриксона. В работе [240] соотношения между силами и перемещениями в зоне контакта представлены в виде ассоциированного и неассоциированного законов скольжения. Несколько иной подход продемонстрирован в работах [221, 222], где использована аналогия между законами пластического течения и законами движения жестких или упругих блоков с сухим трением. Дальнейшее развитие этого направления представлено в работах А. Г. Кузьменко [124, 125], где проводится аналогия механики контактной среды с законами пластичности и ползучести. Достоинства такого подхода особенно ярко проявляются при решении упругопластических контактных задач.
Другой путь к решению контактных задач МКЭ открывается с использованием специальных стыковочных элементов, моделирующих диаграмму сила — смещения на поверхностях раздела взаимодействующих тел. Идея применения элементов особого типа принадлежит, очевидно, авторам работы [225], которые для моделирования трещин и швов горных пород применили разрывные контактные элементы. Дальнейшее усовершенствование контактных элементов проводилось в работах [52, 138, 155, 227, 252] и др.
Авторы работ [147, 212, 219, 242, 256] предлагают решение контактной задачи без использования каких-либо аналогий и стыковочных элементов. В отличие от предыдущего подхода, где контактные элементы объединяют взаимодействующие тела в одну систему, для работ данного направления характерно раздельное рассмотрение контактирующих тел. При этом общая система пополняется определенным количеством уравнений совместности, кратным числу контактирующих узлов. Для решения задачи обычно применяется пошаговый процесс нагружения [219, 223] с уточнением граничных условий на каждом шаге итерационным методом. Приращения нагрузки выбираются достаточно малыми [147] для сохранения линейной связи между перемещениями и деформациями в пределах каждого шага по нагрузке. Такой подход требует многократного решения краевой задачи, а также построения сложных итерационных алгоритмов корректировки граничных условий.
В ряде работ [176,177] для решения контактной задачи МКЭ предлагается использовать релаксационную процедуру. В этом случае континуальное тело предполагается состоящим из системы материальных точек, соединенных между собой упругими связями. Деформация в таком теле распространяется от ее источников равномерно во все стороны путем смещения материальных точек, что приводит к последовательному деформированию связей.
Наиболее современное состояние проблем, возникающих при конечно-элементной реализации приведено в книгах Wriggers Р. [257](2002) и Laursen Р. [238](2002). При решении контактной задачи методом конечных элементов основная сложность заключается в выполнении условий непроникновения, а также дополнительных кинематических условий в случае задачи с трением на общей неизвестной границе. Для выполнения условий контакта при конечно-элементной реализации получили распространения следующие методы: метод множителей Лагранжа, метод штрафа (penalty method), обобщенный метод множителей Лагранжа (Augmented Lagrange Method), прямое решение вариационной задачи ограничениями типа неравенств методами квадратичного программирования. Метод множителей Лагранжа основан на введении на неизвестной границе контакта дополнительных неизвестных, являющимися с механической точки зрения контактными усилиями. Основным неудобством является необходимость решения итерационной задачи с дополнительными неизвестными на контактной границе, что приводит к изменению глобальной матрицы жесткости и операциями с ней на глобальном уровне. Метод штрафа заключается в определении контактных усилий из дополнительно определяемых гипотез, включающих обычно параметр штрафа, при стремлении которого к бесконечности контактные условия выполняются асимптотически точно. Преимущество метода штрафа заключается в том, что дополнительные условия для контактных условий вводятся локально на элементе, что приводит к возможности построения, так называемых, контактных элементов локально. К недостаткам метода можно отнести сложность сходимости при больших значениях параметра штрафа. Объединяющим началом всех трех алгоритмов является алгоритм поиска зоны контакта. Таким стал, обоснованный с точки зрения решения экстремальных задач с ограничениями типа неравенств проекционный алгоритм, получивший название «алгоритм проекции ближайшей точки» (the closest point projection algorithm). Данный алгоритм позволяет построить контактные элементы, основанные на подходе, получившим название «мастер-слуга» («master-slave»).
В последнее время возросло количество публикаций, посвященных применению МГИУ к решению контактных задач. По мнению специалистов, использование МГИУ, обладающего высокой точностью результатов в зонах больших градиентов напряжений, простотой реализации и малым объемом входной информации, дает ряд преимуществ по сравнению с аналогичными схемами МКЭ. Наиболее применяемым для решения контактных задач является прямой вариант МГИУ, где решением являются неизвестные граничные значения переменных, выраженные в естественных физических величинах.
Постановки и подходы к решению контактных задач методом граничных интегральных уравнений во многом сходны со схемами МКЭ. В частности, в работе [209] развиваются идеи использования последовательных и параллельных блочных методов по аналогии с МКЭ для задач контакта нескольких тел. Решены задачи анализа напряжений в резьбовых соединениях с использованием постоянных, линейных и квадратичных граничных элементов. Внимания заслуживает исследование особенностей использования МГИУ для осесимметричных задач при наличии угловых точек на границе. Приведенные расчеты демонстрируют высокую эффективность предлагаемого подхода.
На основе вариационных неравенств и предложенных автором работы [248] полувариационных неравенств приводится постановка задач механики с односторонними ограничениями и ее решение непрямым МГИУ. В силу одностороннего характера взаимодействий вместо интегральных уравнений автором получены интегральные включения.
В работах [206,207] излагаются основные концепции, лежащие в основе формулировок и методов решения плоских контактных задач статической теории упругости. Описаны две методики решения плоских контактных задач, одна из которых применима при отсутствии сил трения, а другая — при их наличии. Рассматривается контакт двух тел, причем каждое из них независимо. Учет условий контакта позволяет связать две системы уравнений в одну. Для нахождения зоны контакта нагрузка прикладывается малыми приращениями, после каждого из которых зоны сцепления и проскальзывания определяются итерационным способом. В созданном программном обеспечении использовались простейшие кусочно-постоянные граничные элементы. Предложенный алгоритм демонстрировался на ряде конкретных задач. Однако рассмотрение контакта только двух тел и использование граничных элементов низкого порядка аппроксимации вводит существенные ограничения на класс и точность рассматриваемых прикладных задач, на возможность расчета НДС различных реальных конструкций.
В более поздних работах Т. Андерсона [208] описывается комплекс программ, реализованный на основе МГИУ для решения плоских контактных задач с учетом трения и без него. Решены задачи о давлении ролика на упругое основание, о контакте круглого диска с границей отверстия в бесконечной области, взаимодействии стальной заклепки с алюминиевым бесконечным листом.
Постановка задачи и возможности реализации МГИУ для решения смешанных и контактных двумерных задач рассмотрены в работе [211]. Особое внимание уделено перспективам применения МГИУ для решения контактных задач.
В работе [239] рассмотрены некоторые приложения МГИУ применительно к решению плоских задач теории упругости. Решена задача о контакте полукруга с полуплоскостью и контакте двух полукругов. Полученные результаты хорошо согласуются с решением по теории Герца.
Используя метод функций Грина и функций влияния, автор работы [259] строит систему МГИУ для решения плоских контактных задач теории упругости с учетом характерных вариантов граничных условий и условий взаимодействия в зоне контакта упругих тел (сцепление, проскальзывание, проскальзывание с трением). Для модельной задачи о контакте двух прямоугольников, в одном из которых имеется длинная неглубокая ступенчатая выемка, даются сравнения решений МГИУ, МКЭ и экспериментальных данных, подчеркивающие точность результатов, полученных МГИУ.
На примере контакта двух квадратных свободно опертых пластин, одна из которых нагружена поперечным давлением [210], рассмотрена задача контакта для изгиба пластин. Зона контакта определяется с помощью итерационной процедуры.
Методом граничных интегральных уравнений решен также ряд задач о внедрении штампов в упругие тела [104, 105, 162, 197]. В работах [104, 105] рассматриваются осесимметричные и плоские задачи о воздействии штампов на балочную плиту и о системе заглубленных штампов. Получены и реализованы системы граничных интегральных уравнений для задач такого класса. Решение сводится к реализации смешанной задачи теории упругости. С использованием методов последовательных приближений для решения граничного интегрального уравнения в работах [108, 167] решен ряд прикладных задач оценки прочности деталей прокатных станов. Подробно рассмотрены вопросы численной реализации для случая второй основной задачи теории упругости. Исследованы задачи о прессовой посадке составных цилиндров с учетом температурного воздействия, волочении проволоки из квадратного прута и т. д. Решение поставленных задач сводится к рассмотрению последовательности смешанных задач теории упругости.
В работе [115] задача о соприкосновении абсолютно жесткого гладкого штампа с изотропным упругим телом ставится как задача нелинейного программирования. Определение зоны контакта, контактного давления и НДС вытекает из минимизации соответствующего функционала. Решение поставленной задачи проводится методом потенциала.
Необходимо отметить, что применение численных методов, таких, как МГИУ и МКЭ, к решению контактных задач существенно расширило спектр их приложений. Благодаря индифферентности методов к описанию геометрии объектов и условий нагружения появилась возможность решения реальных, практически важных задач для резьбовых, фланцевых и замковых соединений различных типов, разнообразных узлов трения деталей машин, технологических посадок и других конструкций.
В большинстве публикаций, посвященных решению прикладных контактных задач, используется двумерная постановка краевой задачи, в которой НДС объектов определяется соотношениями осесимметричной либо плоской задачи теории упругости.
Необходимо отметить, что известные алгоритмы прикладных контактных задач не являются достаточно универсальными, поскольку ориентированы на решение задач определенного класса. Попытки построения более общих алгоритмов решения такого рода задач приводят, как правило, к наложению друг на друга ряда итерационных процедур. В этом случае вычислительная схема задачи становится чрезвычайно громоздкой, что отражается на сходимости процесса решения и затратах машинного времени. Поэтому поиск простых и эффективных методов решения контактных задач с учетом сложной геометрии, условий нагружения и характера деформирования по-прежнему остается актуальной задачей механики твердого деформируемого тела.
В последние годы для решения задач механики эффективно применяются методы граничных интегральных уравнений или методы потенциала. Они основаны на преобразовании исходной системы дифференциальных уравнений в систему граничных интегральных уравнений (ГИУ), из решения которой определяются некоторые определенные на границе функции плотности.
Искомое решение в области определяется граничными интегральными соотношениями с найденными функциями плотности. При численном решении дискретизация проводится не для области, а для ее границы.
Своими корнями метод ГИУ уходит в классический математический анализ. В XIX в. были развиты понятия о силах притяжения в ньютоновских гравитационных полях, получены функции Грина для некоторых частных конфигураций. В 1905г. вышла работа Фредгольма по исследованию интегральных уравнений [220].
До конца 50-х годов методы граничных интегральных уравнений интенсивно развивались математиками. Большой вклад в развитие этих методов был сделан Михлиным С.Г., Купрадзе В.Д., Мусхелишвили Н.И., Смирновым В.И. и др. [141, 142, 129, 130, 143, 144, 180, 181].
Купрадзе В.Д. введены векторные интегральные уравнения методов потенциала в задачах теории упругости [129, 130]. Он развил приближенные методы решения статических задач для однородных тел и динамических задач для кусочно-однородных тел. Сформулировал связь между перемещениями и напряжениями на границе среды, используя распределения поверхностной плотности источников.
В настоящее время теория линейных, а также некоторых классов нелинейных сингулярных уравнений хорошо разработана и изложена в известных монографиях Михлина С.Г., Гахова Ф.Д., Векуа Н.П., Пресдорфа 3., Чибриковой Л.И., Партона В.З., Перлина П.И., Купрадзе В.Д., Гегелиа Т.Г., Башелейшвили М.О., Барчуладзе Т.В. и др. [141,142, 85, 58, 169, 199, 153, 154, 130].
Из теории известно, что сингулярные интегральные уравнения в замкнутом виде решаются в очень редких случаях. Поэтому для приложений первостепенное значение приобретает разработка приближенных методов вычисления сингулярных интегралов и решения сингулярных интегральных уравнений.
В этом направлении выполнены фундаментальные работы Иванова В.В., Корнейчука А.А., Белоцерковского С.М., Лифанова И.К., Габдулхаева Б.Г., Бойкова И.В., Плещинского Н.Б. [98, 106, 48, 78, 53, 54, 55, 149, 160,161].
Методы решения граничных задач с помощью разложений по фундаментальным функциям разработаны в монографиях Купрадзе В.Д., Гегелиа Т.Г., Барчуладзе Т.Г., Башелейшвили М.О., Алексидзе М.А. [12, 130]. Идейно эти методы близки к методам ГИУ, где уравнения рассматриваются, как правило, на основной поверхности граничной задачи. Это приводит к интегральным уравнениям второго рода, но при этом ядро интегрального уравнения становится сингулярным. Решения граничных задач методом разложения по фундаментальным решениям приводит к интегральным уравнениям первого рода.
В работах Гавели С.П., Мельникова Ю.А. 79 строятся интегральные представления, ядрами которых служат не фундаментальные решения, а функции или матрицы Грина соответствующих уравнений или систем для областей, границы которых частично совпадают с рассматриваемыми. Идея построения таких неклассических потенциальных представлений принадлежит В.Д. Купрадзе и применялась С.П. Гавелей, Ю.А. Мельниковым и их учениками, где соответствующие функции и матрицы Грина строились приближенно.
Численной реализации методов потенциала в задачах теории упругости посвящены работы Партона В.З., Перлина П.И., Верюжского Ю.В., Угодчикова А.Г., Хуторянского Н.М., Бреббия К., Уокера С., Бенерджи П., Баттерфилда Р., Крауча С., Старфилда А., Круза Т., Риццо Ф., Теллеса Ж., Вроубела Л., Громадки Т., Лей Ч., Кузнецова С.В., Лившица И.М., Розенцвейга Л.Н. и др [153, 154, 71, 191, 49, 56, 57, 120, 214, 185, 122, 123, 92, 131].
Среди методов построения ГИУ можно выделить два основных направления. Прямой МГЭ, основаный на формуле Сомилиана, полученной из теоремы о взаимности работ, где неизвестные плотности в интегральном уравнении имеют реальный физический смысл. В теории оболочек такими величинами являются перемещения и усилия.
В непрямом МГЭ ядра интегральных уравнений представляют фундаментальное решение и его производные, распределенные на границе рассматриваемой области с некоторой плотностью. Функции плотности не являются решениями задачи на границе, но если они определены, то решение в области определяется вычислением граничных интегралов.
Непрямой МГЭ в задачах изгиба пластин известен как метод компенсирующих нагрузок. Функциям плотности придается смысл нагрузок, приложенных к бесконечной пластине и распределенных по границе области, или по некоторому контуру, внутри которого находится область. В задачах изгиба пластин первые работы в этом направлении выполнены Кореневым Б.Г. и дальнейшее развитие этот метод получил в работах Толкачева В.М., Артюхина Ю.П., Грибова А.П., Венцеля Э.С., Крамина Т.В., Крамина М.В и др. [188-190, 1623, 26, 29-33, 59-70,47, 88, 89].
Монография Верюжского Ю.В. [71] посвящена разработке методов исследования деформирования плоских и пространственных тел при статическом нагружении. Для основных видов граничных условий при дискретизации границы и аппроксимации плотностей выполнен переход от интегральных соотношений к их алгебраическим аналогам. Интегральные соотношения записаны на основе формулы Сомилиана. Предложена методика вычисления эластопотенциалов в виде замкнутых аналитических выражений для кусочно-линейных и сплайновых аппроксимаций функций плотностей. Рассмотрено решение задач изгиба пластин, плоской и пространственной задач теории упругости.
Существенным вкладом в развитие метода компенсирующих нагрузок для пластин сложной формы являются работы Толкачева В.М. 188, где особое внимание уделено теоретическому анализу. В частности, в связанной с контуром системе координат получены аналитические формулы для ядер интегральных уравнений, дан анализ поведения ядер в особых точках, определены предельные значения основных потенциалов, предложены способы устранения неинтегрируемых особенностей.
Работы Венцеля Э.С. и его соавторов 65 посвящены применению метода компенсирующих нагрузок к решению линейных задач теории упругости, пластин и пологих сферических оболочек. В некоторых из них компенсирующие плотности располагаются вне контура пластины, что приводит к интегральным уравнениям Фредгольма первого рода. Наиболее полно их исследования опубликованы в монографии [63].
Проблеме построения и анализа фундаментальных решений теории пластин и оболочек посвящены работы Лукасевича С., Ольшанского В.П., Шевченко В.П., Белоносова С.М., Артюхина Ю.П., Гурьянова И.Н. [132, 150, 151, 200, 198,47, 34, 93].
В постановке многих краевых задач предполагается, что решения считаются непрерывными. Однако при решении практических задач нагрузочные члены содержат различного рода особенности: например, сосредоточенные силы, сосредоточенные моменты, локальные температурные источники и т.д. Применение обобщенных функций решает эту проблему и расширяет возможности применения методов интегральных преобразований Фурье. В работах [198, 200] с использованием методов интегрального преобразования Фурье получены фундаментальные решения пластин и пологих оболочек, проанализирована их асимптотика при малых значениях аргумента и выделены особенности частных решений.
Работа Серазутдинова М.Н., Банцарева К.Н. [179] посвящена синтезу МГЭ и вариационного подхода при анализе пластин и оболочек. Отмечается, что сочетание достоинств этих методов позволяет снять ограничения, присущие каждому из них в отдельности.
Решению задач изгиба изотропных пластин сложной формы МГЭ посвящено большое число работ [87, 18, 22, 23, 65, 139, 146, 152] и др.
В этих работах рассмотрено применение прямого и непрямого МГЭ к решению задач изгиба пластин при различных силовых и температурных нагружениях; рассмотрены различные способы дискретизации границы, различные аппроксимации функций плотности; дан анализ ядер интегральных уравнений; разработаны методики вычисления сингулярных интегралов. Для различных граничных условий получены числовые результаты. Исследована практическая сходимость приближенного решения.
Пластины произвольного очертания при различных граничных условиях рассматриваются также в [249]. Изучается интегральная формулировка МГЭ, связанная- с интегральными представлениями для величин прогибов и изгибающих моментов. В [260] применяются прямолинейные граничные элементы с 9 степенями свободы. Сходимость исследуется на примере свободно опертой пластины. Проведено сопоставление с конечноэлементным решением для защемленной пластинки и прямоугольной пластины с квадратным отверстием.
Двумерным задачам теории упругости посвящены работы: [68], [86], [96], [109], [118], [119], [137], [148], [159].
В 1985 году применение МГЭ к решению плоской задаче теории упругости предложено У. Фишером и С.И. Богомоловым в статье [193].
Решение плоской задачи теории упругости МГЭ изложено в [120]. Особое место в ней уделено практическим задачам. Здесь иллюстрируются обобщения МГЭ и технические приемы для увеличения точности решения, построение вычислительных программ; проведено сравнение полученных результатов с аналитическим решением для тестовых задач.
Работа [68] посвящена исследованию плоской задачи теории упругости, однако только в случае разрывных и сосредоточенных нагрузок на границе области. Рассматривается вторая граничная задача в области с границей для системы уравнений Ляме в перемещения (плоская деформация). Основные теоретические и практические результаты применения МКН в граничных задачах обобщены в [69]. Рассмотрены и изучены вопросы, связанные с особенностями построения и исследования приближенных решений интегральных уравнений первого рода в качестве рабочего аппарата МКН. Это позволяет существенно расширить применение упомянутых интегральных уравнений. Разработана численная реализация метода, приведены некоторые численные результаты.
В работе [109] метод граничных интегральных уравнений используется для решения плоской задачи для тел с угловыми точками. Методы построения функций Грина и численной реализации метода граничных интегральных уравнений для уравнений равновесия плоской задачи теории упругости с использованием метода функции напряжений Эри рассматриваются в [261]. Показано применение численных методов решения граничных интегральных уравнений для задач теории упругости; анализа концентрации напряжений при растяжении прямоугольной пластинки с прямолинейной трещиной. Численные решения задач о плоском напряженном состоянии МГЭ представлены в [96], [243].
Прямой МГЭ в применении к плоским задачам теории упругости излагается в [228]; проводится сравнение результатов с данными, полученными МКЭ. Способ определения коэффициента интенсивности напряжений в плоской задаче предлагается в [159]. Применяется теорема взаимности к исходному решению краевой задачи теории упругости для области, ограниченной контуром с угловой точкой, и к вспомогательному решению, удовлетворяющему однородным краевым условиям. Последние имеют специальную, неинтегрируемую особенность в напряжениях. Выражение для коэффициента интенсивности напряжений получается в виде интеграла по контуру от произведения заданного на границе вектора напряжений на смещения вспомогательного решения.
В [245] проводится анализ двумерных задач теории упругости. Обсуждаются основы МГЭ, демонстрируются его преимущества при решении двумерных граничных задач. Вариант МГЭ для ПНС при использовании функции Эри представлен в [232]. Здесь, кроме этого, обсуждается алгоритм дискретизации. Граничные интегральные уравнения относительно производных и их применение в плоской задаче теории потенциала рассматриваются в [213]. В качестве неизвестной функции выступает тангенциальная производная потенциала
В [148] приведены плоские граничные интегральные уравнения, в которых в качестве неизвестных выступает вектор напряжения на границе и производные от перемещения на контуре. Такие неизвестные позволяют находить тангенциальные граничные напряжения без Численного дифференцирования перемещений на границе, что повышает точность. Приводятся результаты численного решения ряда модельных задач. Граничные интегральные уравнения плоской задачи в [140] получены при использовании потенциалов двойного слоя первого и второго рода. Эти уравнения, кроме Членов, возникших в соответствующих граничных интегральных уравнениях длк задачи с заданными граничными смещениями, содержат также операторы Вольтерра и конечномерные операторы. Представлены асимптотики решений граничных интегральных уравнений около угловых точек границы и указана их связь с соответствующими асимптотиками решения исходной краевой задачи. Конформным и неконформным граничным элементам в плоской задаче теории потенциала посвящена статья [119]. В [233] предлагается новый полигональный элемент (в виде ломаной линии) для решения двумерных задач эластостатики. Разработана техника аналитического интегрирования, позволяющая преодолеть трудности, связанные с сингулярностями. При использовании этого метода напряжения на границе могут быть получены без численного интегрирования. Для оценки точности предлагаемого элемента вычислены значения коэффициента концентрации напряжений для изгибаемой пластины с глубоким вырезом гиперболической формы. Отмечается более высокая эффективность этого элемента по сравнению с линейным элементом.
Развитие МГЭ для решения плоской задачи теории упругости дано в [86]. Исследуемая область может быть многосвязной с произвольным контуром при произвольном нагружении, включая концентрированные силы. На основе МГЭ разработана методика, алгоритм и программа, позволяющая разделить на подобласти с различными упругими и геометрическими характеристиками исследуемый объект. Затруднения вызывают угловые точки. Даны численные результаты, отражающие эффективность метода. Проблеме подразделения на подобласти посвящена работа [231]. Представление в теории потенциала и линейной теории упругости решения уравнений через решение и его производные на границе (при помощи фундаментального решения уравнений) используется для получения конечных соотношений для некоторых представителей решения и его производных в подобластях, на которые предварительно подразделена область определения решения. В работе даются оценки погрешности предлагаемого метода решения граничных задач, вводится регуляризующий параметр для невязок на границах подобластей. Приведены примеры двумерных задач, в которых область подразделена на две простейшие. МГЭ для подобластей рассматривается и в [235]. Представление Кельвина для перемещений статической линейной теории упругости путем подразделения области на подходящие части сводится к матричному уравнению для некоторых представителей перемещений и усилий в подобластях. Рассмотрены примеры плоской задачи теории упругости для областей, составленных из усеченных секторов и прямоугольников.
В [218] дается аналитический обзор применения МГЭ в теории упругости. Выделяется его преимущества: простота дискретизации, сокращение ненужной вычислительной информации, времени счета и машиной памяти, высокая точность определения напряжений и деформаций во внутренних точках, возможность рассмотрения бесконечных и полубесконечных тел. Из недостатков отмечаются: малое число разработанных коммерческих программных продуктов, необходимость высокого уровня подготовки пользователя, трудности решения задач с существенно вытянутыми областями, а также нелинейных задач. Проводится сравнение МКЭ и МКР.
В работе [116] также подчеркиваются преимущества МГЭ по сравнению с МКЭ. По мере увеличения размерности задач, совокупные расходы для схем МГЭ, связанные с ЭВМ, увеличиваются значительно менее резко, чем для схем МКЭ. Как только получены решения на границе, могут быть вычислены значения переменных, описывающих решение, в любых внутренних точках. Отмечается, что граничное интегральное уравнение, является формулировкой поставленной задачи, ведущей к точному ее решению, погрешности вследствие дискретизации и численных аппроксимаций возникают только на границах. Погрешности могут быть очень малыми, если процедура численного интегрирования сделана достаточно сложной. Кроме того, численное интегрирование всегда есть более устойчивый и точный процесс, чем численное дифференцирование.
Из приведенного обзора видно, что современные исследования в области решения контактных задач проводятся в основном в вариационной постановке с использованием метода конечных элементов. Интегральная постановка в решении контактных задач, с которой и получили начало исследования в этой области, ограничивается только задачами для тел с достаточно простой геометрией. Это объясняется трудностью построения функции Грина, которая выступает ядром в разрешающих интегральных уравнениях, для тел сложной формы. Использование метода граничных элементов позволяет построить функцию Грина в численно-аналитическом виде для тел произвольной формы. Отсюда становится актуальным развить МГЭ на решение контактных задач.
Решение большинства прикладных контактных задач, возникающих при прочностных расчетах элементов машиностроительных конструкций, выполнено по теории Герца, которая имеет некоторые ограничения на геометрию и способ взаимодействия тел. С развитием современных численных методов решения контактных задач становится актуальным исследовать применимость решений теории Герца.
Таким образом, перед автором диссертации были поставлены следующие цели:
– применить метод граничных интегральных уравнений к решению контактных задач для тел сложной геометрии путем построения численно-аналитической функции Грина;
– разработать и применить алгоритм на основе НМГЭ решения двумерных контактных задач для пластин произвольной формы с неизвестной областью контакта;
– принимая во внимание, что контактные нагрузки имеют распределение по границам тел, построить интегральные уравнения НМГЭ с фундаментальным решением для полуплоскости решения плоских задач теории упругости для тел сложной формы с замкнутым контуром;
– разработать и применить алгоритм на основе НМГЭ решения плоских контактных задач в отказе от положений теории Герца, сравнить полученные решения с решением по теории Герца.
Диссертация, в которой реализованы поставленные цели, состоит из введения, четырех глав, заключения и списка литературы из 261 наименования.