CC BY
10
УДК 539.3
Ю. П. АРТЮХИН, С. А. МАЛКИН
КОНТАКТ ПЛАСТИН С ЖЕСТКИМ ОСНОВАНИЕМ
Решается контактная задача об изгибе пластин, ограниченном гладким жестким плоским основанием. Решение проводится в интегральной постановке аналитически и численно непрямым методом граничных элементов (НМГЭ). Аналитически задача решается для круглой пластины с учетом и без учета обжатия. В случае неучета обжатия задача становится математически некорректной и дает скачкообразное распределение контактных напряжений. Учет обжатия приводит к уравнению Фредгольма ¡¡рода, которое решается сведением к краевой задаче. В этом случае контактные напряжения принимают гладкий вид. Для решения задачи для пластин произвольных очертаний предложен численный алгоритм на основе НШЭ. Сравнение аналитического и численного решений для круглых пластин показало хорошую согласованность результатов как по областям контакта, так и по контактным напряжениям.
Контактные задачи теории пластин могут быть сформулированы в виде эллиптических вариационных неравенств с односторонним ограничением внутри области [1]. В настоящей работе предлагается алгоритм решения контактных задач в интегральной постановке [2].
Рассмотрим задачу об изгибе круглой пластины радиуса и под равномерно распределенной нагрузкой q, ограниченном гладким жестким плоским основанием, отстоящим от пластины на расстоянии а. Требуется определить напряженно-деформированное состояние пластины, контактные напряжения а(х, у) и область контакта S.
Решим задачу без учета обжатия пластины. Общее решение уравнения осе-симметричного изгиба пластины
О имеет вид
¿5 Ггггф^Л С/ (1)
Неизвестными здесь будут четыре константы общего решения и радиус контакта Ь. Для их определения используем следующие условия. Два граничных условия жесткого защемления краев пластины ю(я}=0, н/(а)= 0. (2)
На границе контакта известен прогиб и угол поворота пластины, изгибающий момент равен нулю, тогда
= (£■}=<). Ш
Решение системы уравнений (2), (3) дает радиус контакта Ь. Контактные 29 Вестник УлГТУ 1/2001
напряжения внутри области контакта будут равномерно распределены с интенсивностью q, а распределение напряжений по линии границы области определяется из условия статического равновесия пластины и будет иметь вид погонного усилия. Это решение является математически некорректным, поэтому решение задачи будем проводить с учетом обжатия пластины в интегральной постановке. Функция Грина защемленной круглой пластины для кольцевой нагрузки, приложенной по радиусу £„ имеет вид [3]
^Ь1-?)- а»* И * Щ 1\+
где Щг -£) - единичная функция Хевисайда. В такой постановке условие контакта запишется следующим образом [4]: Ь
©О(!■)+ = -Ш ■+ Опч{г\ | Г !<*. (5)
о
Здесь ^ ^ частное решение уравнения осесимметричного изгиба пла-
стин, 6 = Ю; к - коэффициент постели, G(r,4) = DG0(г,£,). Учет обжатия пластины моделировался основанием Винклера, коэффициент к равен [4]
2ЕЩ '
4 ' где п - толщина пластины; Ъ - модуль упругости; V - коэффициент Пуассона. Интегральное уравнение (5) можно свести к краевой задаче на интервале (0,Ь) [4]:
од
где Б Однородное уравнение (6) можно раз-
бить на два: - Л (9)
Уравнение (8) имеет решение
Где Ьп(г) и Ки(г) - модифицированные функции Бесселя мнимого аргумента первого и второго рода; Л, - постоянные интегрирования. Решением уравнения (9) является
(-1?"* рг^+Г,^ I.функция /щ/*) имеет особен.
ность при г 0, и поэтому постоянные при ней приравниваем нулю. Введем функции Томпсона
м*}=/4Д(-1Г4
Видно, что они удовлетворяют следующим соотношениям: Тогда регулярное решение уравнения (6) имеет вид
Таким образом, неизвестными являются две константы выражения контактного напряжения и радиус контакта Ь. Рассмотрим условие (7) и выпишем его члены, учитывая, что при А = 0 оно выполняется тождественно:
«М-, - р++
■НИМФ
= 1.
16*
- ЪЩ + ^-ЬефьУсАЬегШ я -0125^=^+0.2+С, р Ц16),
4 Ш = №25Ь{Ь2 -а2}; * (С2Ы'(р&)+ С.бйгОД),
4 П
Из приведенных выражений видно, что условие (7) представляет собой полином по г, в котором присутствуют только нулевая и вторая степени. Так как это условие должно выполняться при любых гс^,Ь), то приравниваем нулю коэффициенты при обеих: степенях г. В результате получаем два условия для определения неизвестных, в которые постоянные С, и С2 входят линейно. Эти условия позволяют выразить постоянные через Ь. Для отыскания области контакта необходимо поставить на ее границе условие равенства ну-
лю контактных напряжений в силу плавного сопряжения форм поверхностей пластины и основания
ст(Ь)= 0. (11)
Это уравнение является трансцендентным относительно Ь. Все приведенные вычисления проводились в пакете программ Mathematica 4. Уравнение (11) имеет настолько громоздкий вид, что встроенные процедуры пакета не позволяют его решить. Для его решения была создана процедура, реализующая метод половинного деления, со следующим кодом:
flb]j>0.на-While[MBS|N|gn], NMstns-cpa, u -N№1]
Iffu«t™[|, i-fi-tsyj]; Щц S=ft+i)ßJ; * N[t]j]; Z).
Входные параметры процедуры bisect имеют следующий смысл: f- заголовок функции; а,Ь - интервал функции, в которых она имеет разные знаки; eps - погрешность решения.
Для случаев неучета и учета обжатия была решена задача для пластины радиуса а = 0,5м, толщиной h = 0,05л*, контактирующей с жестким основанием, отстоящим на расстояние а = 0,01л*. Свойства материала пластины: модуль упругости Е = 2-Ю"Па, коэффициент Пуассона v = 0,3 • На рис. 1,2 показаны прогиб пластины и распределение изгибающих моментов для задачи без учета обжатия при нагрузке q = 20 • 107 Па.
На рис. 3 представлена зависимость радиуса контакта и нагрузки, приложенной к пластине: 1 - без учета обжатия, 2 - с учетом обжатия.
9ДЯ ■ _ ' *------ ------
If
Ряс. 3. Зависимость ркшгуса Контакта я лагруз*н
На рис. 4 показано развитие контактных напряжений для этой задами с учетом обжатия при изменении нагрузки от 2,4-107 Па до 20*1078 Па.
Рис.4. Развитие контактных напряжений
При рассмотрении задачи для пластин произвольной формы аналитическое решение получить не удается. В этом случае решение проводится численно непрямым методом граничных элементов [5]. Прогиб пластины w(x,y) записывается в интегральном виде
О I
Г г
Здесь Q, Г - область и граница пластины; pfe,г|) и /и(Д,л) - неизвестные компенсирующие нагрузка и момент, распределенные по границе Г; п(£,т\) -нормаль к границе Г. Функция G(x,y£,r[)= —г2 1п г является фундаментальным решением изгиба пластины, где г - расстояние между точками (*,>') и т]), которое дает значение прогиба бесконечной пластины в точке (л-,;;) от действия единичной силы в точке (4,л)- Таким образом, имеем три неизвестные подынтегральные функции стЛу), т(лл) и неизвестную область S, по которой
ведется интегрирование. Для определения этих неизвестных используем следующие условия. Два граничных условия закрепления пластины, которые в случае жесткого защемления имеют вид
Условие контакта с учетом обжатия пластины
Условие равенства нулю контактных напряжений на границе области контакта
0. (х^)е Г,
Первые три условия имеют интегральный вид, а последнее - конечный. Подставив выражение для уу(х^у) в приведенные выше условия, получим систему интегральных уравнений
Граница Г разбивается на N постоянных граничных элементов, интегралы по граничным элементам вычисляются аналитически. Постоянный элемент предполагает постоянные значения компенсирующих нагрузки и момента в своих пределах. При дискретизации интеграла по 5 область контакта аппроксимируется трехузловыми треугольными элементами. В пределах такого элемента контактное напряжение берется в виде линейной функции. Искомая граница области контакта разбивается на ^ узлов, координаты которых определяются в полярной системе координат с центром в некоторой точке, лежащей в области контакта. В качестве этой точки рекомендуется брать точку максимального прогиба рассматриваемой пластины. Для каждой точки разбиения границы Г, угловая координата полагается постоянной, это позволяет отыскивать только радиусы точек рк,к = 1 ,...,№. Получаемая дискретная система уравнений является нелинейной относительно рк и решается методом Ньютона. Для пластин простых форм (круг, прямоугольник) решение сходится при задании начального приближения области контакта в виде круга. Для более сложных форм пластин такое приближение дает сходимость только в случае малых областей контакта. При исследовании задач с развитыми областями контакта требуется использовать пошаговое нагружение, беря за начальное приближение области контакта решение с предыдущего шага. 34 Вестник УлГТУ 1/2001
В таблице 1 представлено сравнение результатов аналитического и численного методов, где ст0 - контактные напряжения в центре пластины; спш -максимальные напряжения в области контакта. При решении МГЭ граница пластины разбивалась на 120 граничных элементов, граница области контак-та - на 36 узлов, область контакта аппроксимировалась 1332 треугольными элементами.
Таблица 1
Сравнение аналитического и численного решений
¿у,1(Г Па АНДлигичиооч рсинк Чиелкин« годен МГЭ
6, Л/ а0,1О*/Та Ь,м *Па
2.4 0 ,С]«9 ЗДТКЗ 0,01649 1,22Ш
J 4.75655 Р.0ЯМ7 4,79245 4,79245
6 О,!05Н 0,311522 0.323981 Э,2Т4Т$
С,357!04 (ШЙО ■>,156941 З.ЙР!
оилга 0.660395 4,01106 0,15865 №¿44668 л,««*
,5 оимга 5.1П&85 о.кмм 5, N532
Приведенная таблица демонстрирует хорошее совпадение результатов численного решения с аналитическим. Предложенный численный алгоритм позволяет решать контактные задачи для пластин произвольных очертаний.
СПИСОК ЛИТЕРАТУРЫ
1.Кравчук А. С. Вариационный метод исследования контактного взаимодействия и его реализация на ЭВМ // Расчеты на прочность. 25. - М.: Машиностроение, 1984. - С. 33-50.
2.Григолюк Э. И., Толкачев В. М. Контактные задачи теории пластин и оболочек. - М.: Машиностроение, 1980.- 411 с.
3.Тимошенко С. П., Войновский-Кригер С. Пластинки и оболочки. - М.: Физма-тгиз, 1963.-636 с.
4.Артюхин Ю. П. Контактные задачи для упругих элементов конструкций. - Наб. Челны: КамПИ, 1988. - 78 с.
5.Бенерджи П., Баттерфилд Р. Метод граничных элементов в прикладных науках: Пер. с англ. - М.: Мир, 1984. - 494 с.
Артюхин Юрий Павлович, доктор физико-математических наук, профессор кафедры «Теоретическая механика» Казанского государственного университета, окончил физико-математический факультет Казанского государственного университета. Имеет монографии и статьи в области теории и методов расчета пластин и оболочек.
Малкин Сергей Александрович, аспирант кафедры «Теоретическая механика» Казанского государственного университета, окончил механко-математический факультет Казанского государственного университета. Имеет статьи по теории пластин и оболочек.
