Full Text

В современном мире контактные задачи возникают постоянно, например, при касании пальцем экрана смартфона [1], пальпационной томографии [2]. Задачи механики анизотропных материалов, контактные задачи важны для развития современной промышленности и занимают важное место в современной мировой науке [3-9]. Аналитические методы и подходы к решению контактных задач были развиты в исследованиях[10-14]. Контактные задачи для функционально-градиентных материалов исследовались в работах [15-18]. Проблемы множественного, дискретного и периодического контакта рассматривались в работах [19-21]. Учет дополнительных факторов при контактном взаимодействии (температура, адгезия) проводился в статьях [22-25]. Изучались контактные задачи для трансверсально-изотропного полупространства как с параллельными [26] или перпендикулярными [27; 28] границе плоскостями изотропии, так и с произвольно ориентированными плоскостями изотропии [29]. Отметим, что для случая параллельных границе плоскостей ядро интегрального уравнения контактной задачи отличается только множителем от соответствующего уравнения для изотропного случая, тогда как для случая перпендикулярных границе плоскостей ядро усложняется и теряет симметрию, осесимметричная формулировка становится невозможной [27]. В работе [28] ядро интегрального уравнения «нетрадиционной» [27] контактной задачи для трансверсально-изотропного полупространства получено в форме, свободной от квадратур, на основе теории обобщенных функций. Исследовались как взаимодействие жестких штампов на транстропном полупространстве [30], так и контакт двух упругих транстропных тел [31]. Асимптотические решения контактных задач для изотропного слоя приведены в работе [11]. Задача о действии кругового штампа на транстропный слой (плоскости изотропии параллельны его граням) на основе оригинального метода последовательных приближений рассматривается в исследовании [32]. Интегральные уравнения контактных задач для транстропного тела связаны с интегральными уравнениями задач о трещинах в подобных телах [33]. В настоящей статье фундаментальное решение для транстропного полупространства в форме, свободной от квадратур [28], используется для выделения главной части ядра интегрального уравнения трехмерной контактной задачи для транстропного слоя, что позволяет эффективно рассчитать ядро. В связи с тем что эффективность регулярного асимптотического метода решения контактных задач [11], как правило, падает с ростом анизотропии, для плоских задач предлагается метод специальной аппроксимации символа, приводящий к замкнутому решению, эффективность которого увеличивается при возрастании анизотропии. 1. Пространственные контактные задачи Рассмотрим трансверсально-изотропный упругий слой {|x| < ¥, |y|<¥, 0 £ z £ h}, в верхнюю грань которого без перекоса вдавливается жесткий штамп под действием силы P. Форма основания штампа описывается функцией f (x, y), осадка штампов равна d. Нижняя а b Рис. 1. Схема контакта в задачах А (а) и Б (b) Fig. 1. Scheme of contact in problems A (a) and B (b) грань слоя лежит без трения на недеформируемом основании (скользящая заделка). Плоскости изотропии z = const параллельны граням слоя (задача А, рис. 1, а). Рассмотрим аналогичную задачу Б для слоя {0 £ x £ h, |y| < ¥, |z| < ¥}, которая отличается от задачи А направлением осей координат и тем, что плоскости изотропии z = const перпендикулярны граням слоя (рис. 1, b). Для обеих задач закон Гука имеет вид [4] (1) В частном случае изотропного тела в формулах (1) следует положить (G - модуль сдвига, n - коэффициент Пуассона) (2) Граничные условия задачи А имеют вид (W - область контакта) (3) , а для задачи Б , (4) . При заданных параметрах упругости Amn, известных величинах h, d, функции f (x,y) (задача А) или f (y,z) (задача Б) требуется определить область контакта W и контактное давление sz (x, y, 0) = -q (x,y) ((x,y) ÎW, задача А) или sx (0, y, z) = -q (y,z) ((y,z) ÎW, задача Б), затем при использовании интегрального условия равновесия штампов может быть найдена сила P. Для сведения задач с граничными условиями (3) и (4) к интегральным уравнениям следует рассмотреть вспомогательные задачи типа Буссинеска о действии заданной сосредоточенной силы на грани слоя при разных ориентациях плоскостей изотропии. Решения таких задач находятся при помощи двойных преобразований Фурье. Получив вспомогательные решения, интегрируем выражения для нормальных перемещений по области контакта и удовлетворяем первому граничному условию (1.3) или (1.4). Для задачи А придем к интегральному уравнению ((x,y) ÎW) , (5) Для задачи Б аналогичное уравнение имеет вид ((y,z) ÎW) , , , (6) Здесь gl (l=1,2) - корни характеристического уравнения (возможно, комплексно сопряженные, Re gl> ) [4; 27] (7) Ввиду положительной определенности матрицы упругих постоянных q> [5]. При переходе к изотропному материалу в формулах (5) следует учесть соотношения (2) и пределы , (8) что приводит к известному интегральному уравнению контактной задачи для изотропного слоя [11]. Соответствующий предельный переход возможен и из формул (6). При h®¥ функции-символы ядер (5) и (6) переходят в известные для транстропного полупространства [27; 28; 32] ( переходит в 1, в исчезают котангенсы). Для символа в задаче А выполняется асимптотика . (9) Символ в задаче Б обладает важными асимптотическими свойствами: поведение в бесконечности по одному из аргументов связано с поведением в нуле по другому аргументу. Именно, если один из аргументов фиксирован, (10) Для решения задач А и Б применим метод нелинейных граничных интегральных уравнений (ГИУ) [14; 34; 35], в соответствии с которым уравнения (5), (6) по неизвестной области W сводятся к нелинейным уравнениям по заданному прямоугольнику S со сторонами 2b и 2c (b ³ c), содержащему W. При численном решении нелинейных уравнений используется модифицированный метод Ньютона, область W определяется по точкам разбиения S, в которых решение положительно. В ходе метода требуется рассчитать ядра (5), (6) в узлах решетки. Для эффективного проведения расчетов ниже выделим из ядер главные члены. Введем безразмерные обозначения , (11) и т.д., звездочки далее опускаем. Параметр l характеризует относительную толщину слоя. Используя известную технику [14], асимптотику (9) и обозначения (11), преобразуем ядро (5) для задачи А к форме (12) где J0 (u) - функция Бесселя [36]. Для задачи Б, используя ядро соответствующей контактной задачи для полупространства в форме свободной от квадратур [28], представим ядро в виде , (13) Важно, что в формулах (12) и (13) функции в квадратных скобках экспоненциально убывают на бесконечности. Для расчетов выберем штампы в форме круговых параболоидов (задача А) и (задача Б). Возьмем следующие транстропные материалы (параметры упругости Aij см. [4, с. 22-23]): древесина (ель Дугласа, ), гнейс влажный ( ). На рис. 2 и 3 показаны зависимости безразмерной интегральной характеристики контактных а b Рис. 2. Графики зависимостей P0 (l) для древесины (а) и гнейса влажного (b) для задач А (сплошные линии) и Б (пунктир) при d=A=B=10-3 Fig. 2. Plots of relation P0 (l) for wood (a) and gneiss rock water-saturated (b) for problems A (solid lines) and B (dashed) for d=A=B=10-3 а b Рис. 3. Графики зависимостей P0 (d) для древесины (а) и гнейса влажного (b) для задач А (сплошные линии), Б (пунктир) при l=1, A=B=10-3 Fig. 3. Plots of relation P0 (d) for wood (a) and gneiss rock water-saturated (b) for problems A (solid lines) and B (dashed) for d=A=B=10-3 давлений от относительной толщины слоя l и осадки штампа d соответственно (e = 1). Как видим, вдавливающая сила может быть больше в одной из задач в зависимости от материала слоя. Интегральная характеристика уменьшается с ростом l в связи с усилением влияния скользящей заделки нижней грани слоя (рис. 2) и возрастает с ростом осадки (рис. 3). 2. Осесимметричная контактная задача Частный случай задачи А, когда выполнены условия осевой симметрии, назовем задачей В. Используя обозначения (11) при e = 1, вводя полярные координаты [11], преобразуем интегральное уравнение задачи В к безразмерной форме , (14) где a - радиус области контакта, - контактное давление, - форма основания штампа. Вводя новые функции и по формулам , сведем уравнение (14) к парным уравнениям (15) Применяя для решения парных уравнений (15) известный метод [10], приходим к интегральному уравнению Фредгольма второго рода , (16) для численного решения которого применим метод механических квадратур. При условии ограниченности давления на краю области контакта неизвестный радиус a можно определить из условия [10] . (17) Зная решение уравнения (16), безразмерную интегральную характеристику контактных давлений можно найти по формуле [10] (18) Расчеты интегральной характеристики по формулам (16)-(18) показывают хорошее совпадение с численными результатами для задачи А, полученными при помощи метода нелинейных ГИУ (отличие около 1 %). На рис. 4 показаны зависимости радиуса области контакта a, найденного из условия (17), от осадки d для кругового параболоида (). При l = ¥ (полупространство) упругие характеристики не влияют на зависимость a (d), поскольку ядро интегрального уравнения (16) тождественно обращается в нуль. В случае слоя (l = 1) область контакта расширяется, а значения a (d) для древесины превышают соответствующие значения для изотропного материала. Рис. 4. Графики зависимостей a (d), задача В, при A=10-3: l=¥ - сплошная линия; l=1 - пунктир для изотропного материала, точки для древесины Fig. 4. Plots of relation a (d) in problem C for A=10-3: solid line for l=¥; dashed for isotropic material (l=1), dotted for wood (l=1) 3. Плоские контактные задачи Пусть трансверсально-изотропные слои, изображенные на рис. 1, а и 1, b, взаимодействует с жесткими полосовыми штампами соответственно по заданным областям (задача Г) и (задача Д). Формы оснований штампов в задача Г и Д описываются соответственно функциями f (x) и f (z). Интегральные уравнения контактных задач Г и Д относительно соответственно контактных давлений q (x) и q (z) можно вывести из уравнений (5) и (6), используя теорию обобщенных функций, а для задачи Д еще и асимптотику (10) функции при . После введения безразмерных обозначений l=h/a, Г) , (19) Д) , (звездочки далее опускаем) интегральные уравнения обеих плоских контактных задач можно записать в форме (n= 1 для задачи Г, n=2 для задачи Д) (20) Здесь функции-символы ядер определяются по формулам (5) и (10). Замечание. Если для слоя, показанного на рис. 1, b, полосовой штамп направить вдоль оси z, то в каждом плоском сечении, перпендикулярном оси штампа, будет изотропия упругих свойств. В этом случае на основании асимптотики (10) функции при интегральное уравнение плоской контактной задачи также будет иметь вид (19), где следует в качестве символа взять функцию L (u) из формул (8), которая соответствует контактной задаче для изотропного слоя [11]. В результате интегральное уравнение будет отличаться от известного [11] только множителем, учитывающим анизотропию (он исключается введением обозначений типа (20)), для решения такой задачи можно рекомендовать асимптотические методы [11]. Регулярный асимптотический метод [11], основанный на разложении решения в ряд по степеням малого геометрического параметра 1/l, может применяться при l> 1 в контактной задаче для изотропного слоя, при в задаче Г или в задаче Д, когда сходятся соответствующие разложения ядер (20). С ростом анизотропии (действительные корни ) область применимости метода в задачах Г и Д существенно снижается. Заметим, что в этом случае членом с одним из котангенсов в функциях-символах Ln (u) ядер (20) можно пренебречь ( в задаче Г и в задаче Д), что позволяет получить решение уравнения (20) в замкнутой форме. Именно при учете асимптотического поведения символов Ln (u) в нуле и бесконечности аппроксимируем их при £ u < ¥ выражением ; Г) Д) . (21) Интегральное уравнение (20), (21) сводится к сингулярному интегральному уравнению с ядром Коши, имеющему точное решение [14]. Для штампов с плоским основанием (f (x) =f=const) приближенные решения контактных задач Г и Д имеют вид (22) где K (e) - полный эллиптический интеграл. Относительная погрешность решения (22) не превышает погрешности аппроксимации (21) и уменьшается с ростом анизотропии. Важно, что решение (22) справедливо для любой относительной толщины слоя l. Корни характеристического уравнения (7) и погрешность аппроксимации (21) Roots of the characteristic equation (7) and error of the approximation (21) Материал g1 g2 g1/g2 e, % Изотропный 1,000 1,000 1,00 9,5 Гнейс влажныйi 1,621 0,5476 2,96 7,1 Композит (60 % волокон) 4,724 0,6888 6,86 4,1 Ламинат T650/950-1 5,413 0,6458 8,38 3,5 Древесина (ель Дугласа) 3,713 0,3504 10,6 2,8 В таблицу включено несколько известных трансверсально-изотропных материалов ([4, табл. 1.4]), а также изотропный материал, указана относительная погрешность e аппроксимации (21), практически одинаковая для задач Г и Д. Как видно из данных таблицы, относительная погрешность e монотонно снижается с ростом параметра g1 / g2 (с увеличением анизотропии). Заключение Показана связь интегральных уравнений пространственных контактных задач для трансверсально-изотропного слоя с уравнениями соответствующих плоских задач и задачи при осевой симметрии. Можно сделать вывод о существенном влиянии ориентации плоскостей изотропии на характеристики контакта. Для случая перпендикулярных граням слоя плоскостей изотропии фундаментальное решение соответствующей задачи для транстропного полупространства в форме, свободной от квадратур, позволяет эффективно регуляризовать ядро интегрального уравнения контактной задачи для слоя. При уменьшении l (относительной толщины слоя) возрастает отличие контактных характеристик от соответствующих величин для случая трансверсально-изотропного полупространства. В случае осевой симметрии радиус заранее неизвестной круговой области контакта для слоя оказывается больше, чем в соответствующей задаче для полупространства, а для транстропного слоя область контакта может быть больше, чем для изотропного слоя. В то время как эффективность регулярного асимптотического метода для рассмотренных плоских контактных задач снижается при увеличении анизотропии, эффективность метода, основанного на специальной аппроксимации ядра интегрального уравнения, наоборот, возрастает с ростом анизотропии.

About the authors

D. A Pozharskii

Don State Technical University

N. B Zolotov

Don State Technical University; Southern Federal University

References

Statistics

Views

Abstract - 357

PDF (Russian) - 141