Полный текст

Введение Контактные задачи и их аналитические решения вызывают большой интерес в современной мировой науке и очень важны для решения многих проблем в машиностроении [1-4]. Математический аппарат контактных задач теории упругости был заложен в [5]. Аналитические методы решения задач со смешанными краевыми условиями были рассмотрены в [6]. Исследовались плоские и пространственные контактные задачи для однородных [7, 8] и анизотропных [9] тел. Развиты методы решения плоских и осесимметричных контактных задач для неоднородных тел [10, 11]. Ряд контактных задач для функционально-градиентных тел и покрытий решен в [12-17]. Рассматривалось влияние различных факторов (динамика, микрорельеф, адгезия) на поверхности контакта [18-21]. Одним из основных приближенных методов в неоднородной теории упругости является метод возмущений [22]. Для случая переменного коэффициента Пуассона точные фундаментальные решения были получены для полупространства и слоя (непрерывная зависимость коэффициент Пуассона от глубины) [23-25], для цилиндрических тел (зависимость от радиальной координаты) [26, 27] и для плоского клина (зависимость от угловой координаты) [28, 29]. Изучались контактные задачи для однородного и составного плоского и пространственного клина [5, 7, 8, 29, 30]. Для плоского клина с постоянным коэффициентом Пуассона исследовались случаи зависимости модуля упругости [31-33] или сдвига [10] от угловой координаты. Были получены приближенные решения контактных задач без трения для плоского клина с переменным коэффициентом Пуассона, когда область контакта выходит на угловую точку, основанные на специальных аппроксимациях символов ядер интегральных уравнений [28, 29]. Однако эти решения эффективны лишь для малых изменений коэффициента Пуассона, когда погрешность аппроксимации невысока. В настоящей работе рассмотрен случай, когда область контакта не выходит на угловую точку клина, для решения применяется регулярный асимптотический метод. Предлагаемый подход эффективен для областей контакта, относительно удаленных от угла клина. Для учета трения и шероховатости используется метод коллокаций. 1. Контактные задачи без трения Рассмотрим плоские контактные задачи о взаимодействии бесконечного линейно-упругого клина угла раствора a с жестким штампом по области j = , < a £ r £ b (рис. 1). Для решения задачи удобно ввести полярные координаты r, j с полюсом в вершине клина. Принимается, что вне области контакта грань j = не нагружена, трение не учитывается. Краевые условия задачи формулируются следующим образом. Грань клина j = a жестко заделана (задача А) или свободна от напряжений (задача Б). Материал клина является неоднородным: коэффициент Пуассона n = n(j), модуль сдвига постоянен: G = const. Без ограничения общности будем считать, что штамп имеет плоское основание и вдавливается без перекоса на величину d под действием силы P. В задаче А при заданных величинах a, b, a, d, G и известной функции n(j) требуется определить контактное давление sj(r,0) = -q(r) в области контакта. Затем может быть определена сила P и точка ее приложения. В задаче Б невозможно установить связь между приложенной к штампу силой P и осадкой d [7, 28, 29], поэтому вместо d заданной считается величина P. Этот известный парадокс связан с плоской постановкой задачи для клина или полуплоскости (частный случай клина) [7] и устраняется при рассмотрении пространственной контактной задачи (типа Герца) для клина с одной свободной гранью и полупространства (частный случай пространственного клина) [29]. j = 0 j = a P а б a b j = a j = P n(j), G a b n(j), G Рис. 1. Контакт штампа и клина при жесткой заделке (а) и отсутствии напряжений (б) Fig. 1. Contact of the punch and the wedge for a rigid fixation (а) and without stress (b) Для сведения контактных задач к интегральным уравнениям относительно q(r) рассмотрим вспомогательные задачи о действии заданной нагрузки на грань описанного выше неоднородного клина при разных граничных условиях на другой его грани (жесткая заделка или отсутствие напряжений). Точные решения вспомогательных краевых задач находятся на основе общего представления типа Фрайбергера [34] решения уравнений упругого равновесия в полярной системе координат при переменном коэффициенте Пуассона [28, 29]. Задачи сводятся к векторному уравнению Лапласа и скалярному уравнению Пуассона, правая часть которого зависит от функции n(j). Решения этих уравнений находятся при помощи интегрального преобразования Меллина [35]. Нормальное упругое перемещение поверхности клина под приложенной нагрузкой представляется в форме интеграла Меллина (нагрузка также находится под знаком интеграла). Приравнивая найденное упругое перемещение осадке штампа, получаем интегральное уравнение (ИУ) (1) где для задач А (m = 1) и Б (m = 2) , (2) При n(j) = const ИУ (1), (2) совпадают с известными ИУ для однородного клина [5, 7]. Предположим, что функция h(j) разлагается в ряд Фурье, и далее без ограничения общности оставим в этом разложении два члена: (3) В отличие от задачи А, в задаче Б символ ядра ИУ имеет двукратный полюс при s = 0. Это и будет препятствовать возможности найти связь между P и d в задаче Б. Переходя в ядре ИУ (1) к интегрированию по вещественной оси, используя теорию вычетов для задачи Б, после введения безразмерных величин (4) придем к ИУ (m = 1 и 2 для задач А и Б соответственно) (5) , , . Безразмерный параметр l вида (4) характеризует относительную удаленность области контакта от вершины клина. Для решения ИУ (5) применим регулярный асимптотический метод [5-8], эффективный при больших значениях параметра l (вдали от угловой точки). Как известно [6], структура асимптотики существенно зависит от поведения символа ядра на бесконечности. Для однородного материала (n1 = 0) функция экспоненциально убывает при u®¥ [7]. Здесь же имеем (6) Поведение (6) приводит к появлению в асимптотическом решении логарифмических членов, которые отсутствуют для однородного клина [6,7]. Именно на основании (6) можно показать, что для ядра ИУ (5) имеет место степенно-логарифмическое разложение по малому параметру [6, 36]: (7) . Для задачи Б разложение (7) содержит бесконечную постоянную ввиду расходимости интеграла в нуле. Регулярный асимптотический метод для обеих задач начинается с почленного дифференцирования ИУ (5) по x. После дифференцирования ИУ (5) бесконечная постоянная в ядре для задачи Б исчезает (интеграл в ядре становится сходящимся), а главный логарифмический член в (7) для обеих задач переходит в известное сингулярное ядро Коши. Асимптотическое решение разыскивается в форме степенно-логарифмического разложения по 1/l с неизвестными функциями - коэффициентами этого разложения. Группируются и приравниваются члены при одинаковых степенях . Для нахождения каждого последующего члена асимптотики возникает сингулярное ИУ с ядром Коши, правая часть которого зависит от предыдущих, уже найденных членов. Поскольку производная правой части ИУ (5) равна нулю, для первого члена асимптотики имеем однородное сингулярное ИУ с ядром Коши, нетривиальное решение которого зависит от одной произвольной постоянной [6]. Эта произвольная постоянная окончательно находится из интегрального условия равновесия штампа (8) Таким образом, контактное давление оказывается прямо пропорциональным приложенной к штампу безразмерной силе . Затем в задаче А для определения связи между силой и осадкой штампа полученное решение следует подставить в исходное ИУ (5). В задаче Б этого сделать нельзя из-за расходимости первого интеграла (7). В результате приближенное асимптотическое решение контактных задач получим в виде ( ) (9) , , . В задаче А интегральная характеристика контактного давления (8) (безразмерная вдавливающая сила) вычисляется по формуле . (10) Погрешность решения (9), (10) при l ³ 5 и не слишком малых углах a не превосходит 5 %. Это решение можно рекомендовать, когда область контакта (штамп) относительна удалена от вершины клина. p/2 3p/4 p 5p/4 a a30 2 1 p/2 3p/4 p 5p/4 a p/2 3p/4 p 5p/4 a a31 0,4 0,2 -0,2 P0 1,1 1,0 0,9 0,8 б в а Рис. 2. Графики величин (7) - а, б и (10) - в, P0 = P*/f, в задаче А: сплошные линии при n1 = 0,1, n = 1; пунктир при n1 = 0,1, n = 3; точки при n1 = -0,1, n = 1; точка-пунктир при n1 = -0,1, n = 3 Fig. 2. Plots of parameters (7) - (a, b) and (10) - (c, P0 = P*/f), in problem A: solid lines for n1 = 0,1, n = 1; dashed for n1 = 0,1, n = 3; dotted for n1 = -0,1, n = 1; dashed-dotted for n1 = -0,1, n = 3 p/2 3p/4 p 5p/4 a а -a31 0,8 0,4 -0,2 б j0 0,33 0,32 0,31 0,30 p/2 3p/4 p 5p/4 a Рис. 3. Графики величин (7) - а и (11) - б в задаче Б: сплошные линии при n1 = 0,1, n = 1; пунктир при n1 = 0,1, n = 3; точки при n1 = -0,1, n = 1; точка-пунктир при n1 = -0,1, n = 3 Fig. 3. Plots of parameters (7) - (a) and (11) - (b) in problem B: solid lines for n1 = 0,1, n = 1; dashed for n1 = 0,1, n = 3; dotted for n1 = -0,1, n = 1; dashed-dotted for n1 = -0,1, n = 3 Для задачи А на рис. 2 показано поведение постоянных a30 и a31 (7) и отнесенной к осадке силы (10) в зависимости от угла a при разных n1 и n (n0 = 0,3; в (10) взято l = 5). На основе решения (9) для давления при x = 0 получим (11) Для задачи Б на рис. 3 приведены графики зависимостей постоянной a31 (7) и отнесенного к силе давления (11) от угла a при разных n1 и n (n0 = 0,3; в (11) взято l = 5). 2. Учет трения и шероховатости Рассмотрим контактную задачу Б, которая отличается от задачи А тем, что в зоне контакта будем учитывать силы трения и шероховатость поверхности клина. Пусть силы трения подчиняются закону Кулона: где m - коэффициент трения. Шероховатость поверхности в контактных задачах будем моделировать покрытием винклеровского типа в предположении, что коэффициент постели линейно зависит от радиальной координаты. Такое покрытие вносит дополнительный вклад в нормальное перемещение в области контакта (снова используем размерные обозначения) (12) Условие контакта имеет вид (13) где - упругое нормальное перемещение точек поверхности под действием контактного давления и сил кулоновского трения, которое определяется с использованием фундаментальных решений вспомогательных краевых задач; d - осадка плоского штампа. Для изменяющегося по угловой координате коэффициента Пуассона снова принимаем закон (3). Добавляя к обозначениям (4) величину безразмерное ИУ задачи В на основании формул (12) и (13) можно записать в форме (14) Здесь ядра имеют вид (15) В формулах (15) использованы обозначения, принятые в формулах (5). В отличие от уравнения (5), ИУ (14) является уравнением второго рода. Несмотря на логарифмическую особенность ядра, его решение уже не обладает корневыми особенностями и принимает конечные значения на краях области контакта при Символы (15) характеризуются следующим поведением в бесконечности: (16) Для численного решения ИУ (14) используем метод коллокаций на основе квадратурной формулы Гаусса. Метод эффективен при не слишком малых углах a и значениях параметра l ³ 1. Для улучшения сходимости интегралов в ядре (14) на основе поведения (16) выделим из них главные части: (17) Здесь использованы значения интегралов [6, 36] (18) В правых частях формул (17) подынтегральные функции убывают на бесконечности как . На рис. 4 для задачи Б приведены графики интегральной характеристики контактных давлений в зависимости n1 при разных уровнях шероховатости (коэффициентов постели основания Винклера) для случая m = 0,2 и и различных значениях l. -0,2 0 0,2 n1 а P0 2,4 2,2 1,8 б -0,2 0 0,2 n1 P0 2,0 1,8 1,4 Рис. 4. Графики силы, отнесенной к осадке (P0 = P*/f), в задаче Б: a - при l = 1, m = 0,2; б - при l = 2, m = 0,2: сплошные линии при A* = 0,01; пунктир при A* = 0,001 Fig. 4. Plots of the force related to the settlement (P0 = P*/f) in problem C, a - for l = 1, m = 0,2; b - for l = 2, m = 0,2: solid lines for A* = 0,01; the dashed line for A* = 0,001 Как показывают расчеты, в диапазоне -0,2 £ m £ 0,2 отличие соответствующих значений P0 от приведенных на рис. 4 не превосходит 1 %. Зависимость P0 от n1 близка к линейной. Заключение Неоднородность клина по угловой координате меняет характер поведения символа ядра ИУ контактных задач на бесконечности, что приводит к появлению логарифмических членов в регулярном асимптотическом решении (9), (10). Такие члены отсутствуют в решениях для однородного клина [7]. Ранее аналогичное изменение поведения символа было отмечено в контактных задачах для неоднородной по глубине полосы (экспоненциальное убывание на бесконечности заменяется на степенное) [10]. Коэффициент Пуассона n(j) при законе изменения (3) возрастает при n1 > и убывает при n1 < от значения (от грани контакта j = 0). В задаче А при возрастающем коэффициенте Пуассона по угловой координате (при этом модуль продольной упругости также возрастает) контактные давления и сила, отнесенные к осадке штампа f, больше, чем при убывающем. С ростом l давления и сила убывают ввиду отдаления области контакта от грани с жесткой заделкой. При росте n (частоты осцилляций изменения коэффициента Пуассона) в задаче А сила возрастает при фиксированном n1 > и убывает при фиксированном n1 < . В задаче Б для неоднородного плоского клина с одной свободной от напряжений гранью невозможно установить связь между силой, приложенной к штампу, и осадкой штампа. Это находится в соответствии с известными контактными задачами для однородного плоского клина и для упругой полуплоскости (частный случай клина) [7, 37]. В случае задачи Б интегральная характеристика давлений слабо зависит от коэффициента трения, который влияет в основном на асимметрию (перераспределение) давлений по области контакта. При приближении области контакта к вершине клина (с уменьшением l) значения интегральной характеристики увеличиваются. При возрастании параметра шероховатости контактные давления снижаются. Значения P0 монотонно возрастают с ростом n1, что объясняется повышением коэффициента Пуассона вблизи грани клина в контакте. Для всех трех задач достоверность выведенных ИУ подтверждается их совпадением в частных случаях с известными для однородного материала клина [5-7, 29]. Достоверность полученных асимптотических решений (9)-(11) для задач А и Б подтверждается их совпадением в частных случаях с известными решениями ИУ для однородного материала клина [5-7]. Численные решения задачи Б, полученные по методу коллокаций, верифицировались путем увеличения числа узлов коллокаций.

Об авторах

Д. А Пожарский

Донской государственный технический университет

Е. Д Пожарская

Донской государственный технический университет

Список литературы

Статистика

Просмотры

Аннотация - 305

PDF (Russian) - 222