Метод нелинейных граничных интегральных уравнений для решения квазистатической контактной задачи о взаимодействии упругих тел при наличии кулонова трения

Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2016. T. 20, № 2. С. 306-327

ISSN: 2310-7081 (online), 1991-8615 (print) doi: http://dx.doi.org/10.14498/vsgtu1471 УДК 539.3

МЕТОД НЕЛИНЕЙНЫХ ГРАНИЧНЫХ ИНТЕГРАЛЬНЫХ УРАВНЕНИЙ ДЛЯ РЕШЕНИЯ КВАЗИСТАТИЧЕСКОЙ КОНТАКТНОЙ ЗАДАЧИ О ВЗАИМОДЕЙСТВИИ УПРУГИХ ТЕЛ ПРИ НАЛИЧИИ КУЛОНОВА ТРЕНИЯ

Ю. М. Стреляев

Запорожский национальный университет, Украина, 69600, Запорожье, ул. Жуковского, 66.

Аннотация

Рассмотрена пространственная квазистатическая контактная задача о взаимодействии двух линейно-упругих тел с учетом кулонова трения между ними. В граничных условиях задачи на каждом шаге дискретного процесса нагружения были приняты упрощения, основанные на модификации закона трения Кулона. Эта модификация состояла во введении в соотношения, выражающие закон трения Кулона, запаздывания контактных давлений, ограничивающих касательные контактные напряжения. В рассмотренной постановке задача бала сведена к последовательному решению серии однотипных систем нелинейных интегральных уравнений, описывающих взаимодействие тел на каждом из шагов нагружения. Для получения контактных напряжений на каждом шаге нагружения использован метод приближённого решения системы интегральных уравнений этого шага, который заключается в регуляризации этой системы уравнений, дискретизации регуляризированной системы и применении сходящегося итерационного процесса для получения решения дискретизированной системы. Предложенным методом получено численное решение контактной задачи о вдавливании упругого шара в упругое полупространство при возрастании и последующем убывании нормальной сжимающей силы.

Ключевые слова: упругое тело, контактная задача, трение Кулона, квазистатическая задача, интегральное уравнение, единственность решения, регуляризирующее уравнение, итерационный процесс.

Введение. Для решения контактных задач теории упругости с учетом трения часто используются численные методы, основанные на вариационной постановке задачи [1, 2] либо на её сведении к различным нелинейным операторным уравнениям [3-7]. Основная трудность реализации вариационных методов заключается в необходимости рассматривать сложные задачи нелинейного программирования. Иногда эту трудность удается устранить путем

© 2016 Самарский государственный технический университет. Образец для цитирования

Стреляев Ю. М. Метод нелинейных граничных интегральных уравнений для решения квазистатической контактной задачи о взаимодействии упругих тел при наличии кулонова трения // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2016. T. 20, № 2. С. 306327. doi: 10.14498/vsgtu1471. Сведения об авторе

Юрий Михайлович Стреляев (strelkiny@gmail.com), старший преподаватель, каф. математического анализа.

использования нелинейных операторных уравнений для моделирования контактного взаимодействия тел. Однако в большинстве случаев при использовании таких уравнений авторы либо не учитывают трение [3-5], либо его учет осуществляется при упрощенных граничных условиях, соответствующих полному проскальзыванию тел [6, 7]. Хотя метод нелинейных граничных интегральных уравнений, предложенный в работе [8], и позволяет рассматривать закон трения Кулона в сложном классическом виде, но при использовании этого метода трудно выразить условия единственности решения контактной задачи и доказать сходимость итерационных процессов, применяемых для получения этого решения [9]. В настоящей статье предложена модификация нелинейных граничных интегральных уравнений [8], позволяющая избавиться от указанных трудностей. Данная модификация основана на упрощении граничных условий квазистатической контактной задачи на каждом шаге процесса нагружения тел. Это упрощение позволяет свести контактную задачу к последовательному решению серии однотипных нелинейных операторных уравнений второго рода с квазисжимающим оператором в правой части. Такая особенность этих уравнений позволяет разработать эффективный метод их приближённого решения. Целью данной статьи является разработка такого метода и получение строгого математического обоснования основных его этапов.

1. Постановка контактной задачи. Предположим, что каждое из двух тел на некоторой части своей поверхности сцеплено с абсолютно твёрдым телом, которое будем называть жесткой опорой. Первоначально касаясь в одной точке, тела входят в контакт за счет того, что опора верхнего тела совершает заданное поступательное перемещение 7 (Ь), зависящее от времени Ь, в то время как опора нижнего тела остается неподвижной. Требуется найти поверхность соприкасания тел после установившего равновесия и распределение контактной нагрузки на этой поверхности. Относительно взаимодействующих тел принимаются следующие допущения:

1) оба тела являются линейно-упругими и изотропными;

2) упругие перемещения точек малы по сравнению с размерами поверхности контакта;

3) в каждой точке поверхности каждого из тел в зоне возможного контакта существует единственная касательная плоскость к этой поверхности;

4) операторы влияния поверхностных нагрузок на поверхностные упругие перемещения для каждого из тел являются известными и допускают интегральное представление.

Будем считать, что вектор-функция 7 (Ь) медленно изменяется с возрастанием Ь от 0 до фиксированного значения Т > 0, и пренебрегать инерционными и волновыми эффектами. Предположим, что процесс взаимодействия тел сопровождается поверхностным трением, подчиняющемся закону Кулона [10]. Будем ещё полагать, что в любой момент времени Ь > 0 поверхность соприкасания тел является плоской и содержится в известной ограниченной области О, принадлежащей общей для тел касательной плоскости П, проходящей через точку их начального касания. Граница областей сцепления и проскальзывания заранее не известна и подлежит определению в ходе решения. Введем прямоугольную систему координат Охху с началом в точке первоначального касания тел и осью х, направленной внутрь верхнего тела

Стреляев Ю. М.

перпендикулярно плоскости П.

Будем рассматривать процесс нагружения тел лишь в фиксированные моменты времени ¿о,¿1,¿2,...,¿г такие, что 0 = ¿0 < ¿1 < ¿2 < ••• < ¿г = Т. Будем также предполагать, что при каждом значении Ь = ^ для системы контактирующих тел устанавливается равновесие. Таким образом, процесс нагружения тел происходит пошагово и состоит из конечного числа I последовательных состояний равновесия. Пусть индекс г определяет номер шага нагружения и вектор-функция 7 (¿) принимает на каждом шаге нагружения значение 7» = 7 (¿¿) = (Дн, Д2г, Лэг) (рис. 1).

Граничные условия контактного взаимодействия тел на i-том шаге нагружения будем выражать с помощью двух вектор-функций

где в = (ж, у) — произвольная точка области О, определяющая пару противолежащих точек (рис. 1) на поверхностях взаимодействующих тел (предполагается, что произвольная точка поверхности верхнего тела может вступить в контакт только с противолежащей точкой поверхности нижнего тела). Вектор-функция ^¿(в) определяет распределение удельной контактной нагрузки, передаваемой от нижнего тела к верхнему телу через точки области О на i-том шаге нагружения. Вектор-функция -¿(в) задает перемещения точек поверхности верхнего тела относительно противолежащих точек поверхности нижнего тела на i-том шаге нагружения. Зависимость компонент этих вектор-функций от времени выражается лишь в дискретном изменении их значений при переходе к следующему шагу нагружения. Отметим, что значения 1, 2, 3 индексов соответствуют направлениям координатных осей Оя, Ож, Оу. С учетом принятых предположений граничные условия контактного взаимодействия тел можно выразить следующей системой неравенств и уравнений:

где в € О, г = 1,1; положительная константа ^ — коэффициент трения, а -2г(в), ^зг(в) —компоненты вектора скорости относительного проскальзывания тел в точке в € О на i-том шаге нагружения.

Первое из соотношений (1) выражает знакопостоянство контактного давления, второе означает отсутствие взаимного проникания тел, третье означает отсутствие контактного давления за пределами площадки контакта. Последние три соотношения системы (1) выражают закон трения Кулона.

Будем считать, что компоненты р^, вектор-функций р^(в), ¿¿(в) принадлежат гильбертову пространству ¿2(0), а сами вектор-функции р^(в), -¿(в) являются элементами гильбертова пространства ¿2(0) [9]. Связь между вектор-функциями р^(в) и -¿(в) выражается соотношениями

Рг(в) = (р1г(в),Р2г(в),Р3г(в)) , -(в) = (-1г(в),-2г(в),-3г(в)) ,

3

"^¿(в) = ^ Лу (р,г)8 - /ь(в) = ^ (й, Ук = 1, 3,

у=1

Рис. 1. Схема контактного взаимодействия упругих тел [Figure 1. The scheme of contact interaction of elastic bodies]

где Akj: ¿2(0) ^ L2 (0) —операторы влияния поверхностных нагрузок на поверхностные упругие перемещения; fi(s) = (fu(s), f2i(s), f3i(s)) £ L;](0) — известная вектор-функция, для которой fii(s) = —50(s) — Aii, f2i(s) = — A2i, f3i(s) = — A3i, а функция 50(s) £ L2(0) задает первоначальный зазор между телами, определяемый уравнениями их поверхностей. Выразим компоненты V2i(s), v3i(s) вектора скорости относительного проскальзывания тел в точке s £ 0 на i-том шаге нагружения с помощью конечных разностей:

. ,, Vki(s) — Vki-l(s) 1 г

Vki(s) =-h-= hFk(Pi, Aki)s

hi hi

i = 1,l, k = 2,3,

где Ъ = и - Аы(в) = Аы - Рк (рг-1, Аы-\)3, а Рк (рг, Аы)3 определяется соотношениями (2). Тогда систему (1) можно записать в виде

Pii(s)^0_;_Fi(Pi, Aii)s > 0; Pii(s)Fi (pi, Au)s = 0;

\Jp^i(s)+ p3i(s) < rni(s);

F22(pi, A2i)s + F32(ppi, A3i)sP2i(s) + №ii(s)F2(pi, A2i)s = 0; F22('Pi, A2i)s + F2(pi, A3i)sP3i(s) + VPii(s)F3(pi, A3i)s = 0;

где в £ О, г = 1,1.

Таким образом, рассматриваемая контактная задача на г-том шаге нагружения состоит в определении неизвестных функций Ркг(в) £ Ь2(О), к = 1, 3, которые почти всюду на О удовлетворяют соотношениям (4). Отметим, что из соотношений (3) следует, что для нахождения решения системы (4) на г-том шаге нагружения должно быть известно решение этой системы на (г - 1)-м шаге. В качестве начальных условий на нулевом шаге (в момент начального касания тел) примем тождественное равенство нулю всех компонент вектор-функций 70 и р0(в).

Для упрощения постановки контактной задачи рассмотрим следующую модификацию соотношений, выражающих закон трения Кулона в системе (4).

Вначале запишем эту систему в эквивалентном виде:

рн(в) ^ 0; ¿1 (р, Ди)в ^ 0; рн(в)^\ (р, Ди)в = 0;

^р^ + р^(в) < ^ (Р1»-1(в) + ^(в));

/--

^(р, Л+ ^(р, Л^^¿(в) + ^ (рН-1(в) + ^(в)) ^(р, Л2^ = 0; У 7

\/¿2 (р^, Л2») + ¿Кр^, Л3г)sP3г(в) + ^ (р1^-1(в) + ^¿(в)) ¿з(р^, Дз^ = 0,

где в € О, г = 1,1; ^¿(в) = рн(в) — р1^-1(в). Затем будем предполагать, что при достаточно большом числе шагов нагружения I разности рн(в) — р^-^в) малы по сравнению с величинами контактных давлений рн(в) на каждом i-том шаге нагружения. Наконец, учитывая последнее предположение, отбросим в системе (5) все приращения ^¿(в). Тогда система (5) примет следующий модифицированный вид:

рн(в) ^ 0; ¿1 (р, Дн)5 ^ 0; рн(в)^\ (р, Ди)в = 0;

Л/рМ + р3»(в) < №;-1(в);

/-

^Р22

у/¿Кр», Л2г)s + ¿Кр», Л3г)sP3г(в) + ДОН-1 (в)^^, Л3;)« = 0, где в € О, г = 1,1.

При использовании указанной модификации закона трения Кулона рассматриваемая контактная задача на i-том шаге нагружения сводится к определению неизвестных функций рн(в), р2г(в), р3г(в) € ¿2(О), которые почти всюду на О удовлетворяют соотношениям (6). Принятую модификацию закона трения Кулона, выраженную последними тремя соотношениями системы (6), можно истолковать как наличие запаздывания в выражениях дон(в), входящих в последние три соотношения системы (6) и ограничивающих модуль удельной касательной контактной нагрузки. При этом очевидно, что на i-том шаге нагружения ограничивающее контактное давление рн-^в) является известным, что существенно упрощает граничные условия контактной задачи. Если система (6) решена для i-того шага нагружения, то конфигурацию области контакта на этом шаге можно определить следующим образом:

о0° = {в € О | ¿1 (р, Дн)в = 0},

ОС = {в € о01) | ^2(рг, Л+ ¿32(рг, Л^ = 0}, оП) = {в € о01) | ^2(рг,Л+ ¿32(рг,Л^ > 0},

где о0^ —область контакта, О С —зона сцепления, О^) —зона проскальзывания. Отметим, что при изменении индекса г от 1 до I соотношения (6) представляют собой серию систем, подлежащих последовательному решению. Окончательным решением рассматриваемой квазистатической задачи условимся считать функции ри(в), к = 1, 3 и

области О01), О^, О}?, найденные

на последнем шаге нагружения I.

2. Операторное уравнение контактной задачи. Для неизвестных вектор-функций р1(в),р2(в),... ,р[(в) € Ь3(и) серия систем соотношений (6) эквивалентна серии нелинейных операторных уравнений [9]

Рг

г_1

в которых Е — произвольная положительная константа, ¡г € Ь2(и) и определяются равенствами

¡г = (-бо(в) - Аи, -Аф), -Аф)).

Линейный ограниченный оператор влияния А: ^ Ь2(и), входящий

в уравнения (7), задается соотношениями

Х = (Х1,Х2,Хз), у = (У1,У2,Уз) € ЬЦО); 3

у = А(х), ук = ^ Аку (X]) У к = .

3 = 1

(8

Входящий в правую часть соотношений (7) нелинейный непрерывный операнд : £>3(и) ^ Ь^ (

тор С^д : Ь3(0) ^ Ь3(и) определяется равенствами

X = (Х1 ,Х2,Хз), у = (У1,У2,Уз) € Ь^(иу,

V = (x),

VI (в) = Н( Х1 (в)), (9) У2 (в) = д(х2(в),хз(в),рд(в)),

Уз (в) = д(Хз(в),Х2(в),ц,д(в)), в € и,

в которых д(в) € Ь2(и) —заданная, почти всюду на и неотрицательная функция, а функции Н(х) и д(х, У, г) определяются соотношениями

, , . | Х, если Х ^ 0; п(х) =

0, если Х < 0;

д(Х,У,г) = I хг

если д/х2 + У2 ^ г;

\/х2 + У2

, если д/х2 + У2 > г.

Для получения оператора С№1: Ь2(и) ^ Ь2(и), который фактически используется в выражении правой части уравнения (7), нужно в соотношениях (9), задающих оператор С^д: Ь3(и) ^ Ь2(и), положить д(в) = р1 г-1(в) почти всюду на и.

Важной особенностью [9] каждого из уравнений (7) является то, что множество его решений в Ь:3(и) не изменяется при изменении входящего в это уравнение положительного параметра Е.

Оператор С^д: Ь3(и) ^ Ь2(и), определяемый соотношениями (9), является метрическим проектором [11] гильбертова пространства Ь2(и) на выпуклое замкнутое множество

и = {(х1,Х2,хз) € Ь3(и) | Х1 (в) ^ 0, ^х2(в) + х3(в) ^ ^д(в), в € и}

Стреляев Ю. М.

этого пространства. Этот факт следует из очевидных свойств функций Л,(ж) и ^ (ж, у, г), определяемых соотношениями (10). Действительно, функция Л,(ж) является метрическим проектором одномерного евклидова пространства М1 на промежуток [0; С М1, а отображение ^у: М2 ^ М2 вида

ж = (ж1,Ж2), у = (У1,У2) € М2; у = ^У (ж);

У1 = 9(Ж1,Ж2 , г);

У2 = ?(Ж2,Ж1, г); г ^ 0,

осуществляет метрическую проекцию двумерного евклидова пространства М2 на замкнутый круг

|(Ж1,Ж2) € М2 | ^ж2 + Ж2 ^ г|

этого пространства. Отметим, что входящая в (7) известная неотрицательная функция р1 ¿-1 € ¿2(О) определяется из решения рг-1 € ^(О) этого уравнения на (г — 1)-м шаге нагружения, а известный элемент / € ^(О) также зависит от решения р^^. Таким образом, решение р1 € ^(О) уравнения (7) при фиксированном значении индекса г можно получить после нахождения решений р1,р2,..., р!-1 € ¿2(О) этого уравнения на предшествующих шагах нагружения (при получении р1 предполагается, что в момент начального касания тел вектор-функция контактной нагрузки и жесткие сближения тел тождественно равны нулю на О). Окончательным решением квазистатической контактной задачи, описываемой серией уравнений (7), будем считать решение р^ € ^(О), полученное на последнем ¿-том шаге нагружения. Для многих важных для практики случаев входящие в соотношения (8) линейные ограниченные операторы А^-: ^(О) ^ ¿2 (О) допускают интегральное представление. В этих случаях каждое из операторных уравнений (7) при фиксированном значении индекса г представляет собой систему трех нелинейных интегральных уравнений относительно трех неизвестных функций р1 ¿(в), р2¿(в), р3¿(в), задающих распределение удельной контактной нагрузки, передаваемой от одного тела к другому.

3. Единственность решения контактной задачи. Для доказательства единственности решения контактной задачи, описываемой серией операторных уравнений (7), достаточно доказать единственность решения каждого из уравнений этой серии. Фиксируя индекс г в (7) и обозначая р = € ¿2(О), д = р1 ¿-1 € ¿2(О), Ь = / € ¿3(О), получим более простой вид каждого операторного уравнения из совокупности (7):

р =

В дальнейшем будем полагать, что линейный ограниченный оператор влияния А: (О) ^ ¿3(О), входящий в уравнение (11), является самосопряженным и положительным:

(А(ж), у) = (ж, А(у)) Уж,у € ¿3(О),

где (x,y) —скалярное произведение элементов x, y пространства ¿¡(П), $L3(n) —нулевой элемент пространства ¿2(П). Самосопряженность оператора A(p) обусловлена свойством взаимности упругой среды, а его положительность гарантирована положительностью энергии упругих деформаций системы взаимодействующих тел, подверженных ненулевому нагружению.

Сформулируем и докажем теорему о единственности решения уравнения (11) при условии, что оператор влияния A удовлетворяет соотношениям (12).

Теорема 1. Пусть линейный ограниченный оператор A : ¿¡(П) ^ ¿¡(П) удовлетворяет условиям (12). Тогда для любого элемента b G ¿¡(П) любой заданной и почти всюду неотрицательной на П функции g G ¿2(П) и любого действительного числа ц > 0 уравнение (11) не может иметь более одного решения в пространстве ¿2(П).

Доказательство. Предположим противное. Пусть p и p — два различных решения уравнения (11) в пространстве ¿2(П). Для оценивания нормы разности элементов p и p воспользуемся нерастягиваемостью на этом пространстве метрического проектора : ¿2(П) ^ ¿¡(П) [11]:

IIP - P II = (P - E(A(p) - b)) - G^ (p - E(A (p) - b)) || <

< IIp - p - EA (p - p)||.

Возведя обе части полученного неравенства в квадрат и перейдя в нём от норм элементов к скалярным произведениям, получим оценку

IIp - p II2 < IIp - p II2 - 2E(p - p , A (p - p)) + E2 IIA (p - p)II2 ,

означающую справедливость неравенства

E

(p - p , A (p - p)) < - IIA (p - p)I2 .

Очевидно, что полученное неравенство будет выполняться для любого положительного значения E. Переходя в этом неравенстве к пределу при E ^ +0, получим соотношение

(p - p , A (p - p)) ^ 0,

которое противоречит второму из условий в (12). Полученное противоречие означает, что доказываемое утверждение верно.

Доказанная теорема означает, что при выполнении её условий рассматриваемая квазистатическая контактная задача, описываемая серией уравнений (7), не может иметь более одного решения в пространстве ¿¡(П).

4. Регуляризация операторного уравнения контактной задачи. Разработка численного алгоритма решения уравнения (11) осложняется тем, что во многих случаях входящий в него линейный оператор влияния A : ¿¡(П) ^ ¿2(П) является вполне непрерывным [9,11]. Это обстоятельство приводит к тому, что задача нахождения неизвестного элемента p G ¿2(П) по известным элементам b G ¿¡(П), g G Uo из решения уравнения (11) является некорректно поставленной в смысле А. Н. Тихонова [12] (здесь и ниже символом Uo обозначен конус неотрицательных функций пространства ¿2(П) [13]). Поэтому

вместо уравнения

р = О^д

в котором £ > 0 есть некоторое фиксированное число — параметр регуляризации. Сформулируем и докажем две теоремы, из которых следует, что уравнение (13) является регуляризированным аналогом уравнения (11).

Теорема 2. Пусть линейный ограниченный оператор А: Ь^(и) ^ Ь^(и) удовлетворяет условиям (12). Тогда для любого фиксированного положительного числа £ задача нахождения неизвестного элемента р € Ь3(и) по известным элементам Ь € Ь33,(и), д € и0 из решения уравнения (13) является корректно поставленной в смысле А. Н. Тихонова.

Доказательство. Для доказательства теоремы достаточно показать, что для любых элементов Ь € Ь33,(и), д € и0 уравнение (13) однозначно разрешимо в пространстве Ь2(и) и его решение р € Ь2(и) непрерывно зависит от элементов Ь € Ь'2(и), д € и0 [12].

Пусть £ — фиксированное положительное число и Ь € Ь3(и), д € Щ. Выберем положительное число Е в уравнении (13) так, чтобы выполнялось неравенство

Е<-^т",

£ + II »и V 7

где ||A||* —норма оператора A: L^3(Q) ^ [14].

Докажем, что из неравенства (14) следует оценка

II — E

где I — тожественный оператор, действующий в пространстве L2(Q). Из неравенства (14) и условий (12) следует, что линейный ограниченный оператор (I — E (eI + A)) : L2,(Q) ^ L2,(Q) является самосопряженным и положительным. Найдя его норму с помощью известного соотношения [14], получим для этой нормы очевидную оценку:

Ц1 — E (eI + А)Ц^ = sup (x — E(ex + A(x)),x) = l|x|| = 1

= 1 — Ee — E inf (A(x), x) < 1 — Ee = q< 1,

llxll = 1

из которой следует неравенство (15). Из этого неравенства и нерастягиваемости метрического проектора G^g: L2,(Q) ^ L2,(Q) следует сжимаемость оператора F: L3(Q) ^ L3(Q), стоящего в правой части уравнения (13). Действительно, для любых элементов x,y £ L3(Q) можно получить оценки

UF (x) — F (у)Ц =

= \\Gm (x — E(ex + A(x) — b)) — G^g (y — E(ey + A(y) — b)) || < ^ \\x — y — E(e (x — y) + A (x — y)) \\ ^

где 0 < q < 1. Таким образом, однозначная разрешимость уравнения (13) в пространстве L2(O) для любых элементов b £ Lg (О), g £ Uo следует из принципа сжимающих отображений [15].

Пусть элементы p',p'' £ L2(O) являются решениями уравнения (13), полученными при b = b' £ L2(O), g = g' £ U0 и b = b" £ L|(n), g = g'' £ U0 соответственно. Тогда, учитывая нерастягиваемость оператора Gw: L2(O) ^ L2(O) и соотношения (10), (15), можно получить неравенства

IIP' - Р"! =

= ||GV (p' - E (ep' + A(p') - b')) - G^ (p'' - E (ep'' + A(p'') - b'')) || = = (p' - E (ep' + A(p') - b')) - (p'' - E (ep'' + A(p'') - b'')) + (p" - E (ep'' + A(p'') - b'')) - G^ (p'' - E (ep'' + A(p'') - b'')) || <

< ||G^g/ (p' - E (ep' + A(p') - b')) - G^g (p'' - E (ep'' + A(p'') - b'')) || + + ||G^ (p'' - E (ep'' + A(p'') - b'')) - G^ (p'' - E (ep'' + A(p'') - b'')) || <

< || (p' - p'') - E (e(p' - p'') + A(p' - p'')) + E(b' - b'') || + ^||g' - g"Hi <

из которых для нормы ||p' - p''|| следует оценка

e Ee

Здесь и ниже символами ||  || и |  Hi обозначены нормы элементов пространств L2(O) и L2(O) соответственно. Эта оценка, полученная при фиксированном e > 0 для любых элементов b', b'' £ L2(O) и g',g'' £ U0, означает, что решение уравнения (13) в пространстве L2(O) непрерывно зависит от элементов b £ Lg(O), g £ Uo.

Теорема 3. Пусть линейный вполне непрерывный оператор A: Lg(O) ^ L2(O) удовлетворяет условиям (12) и для некоторых фиксированных элементов b £ L2(O), g £ U0 уравнение (11) имеет решение p* £ L2(O). Тогда

£^+0

где p£ £ Lg(O) является решением уравнения (13), полученным при заданном положительном значении параметра регуляризации e.

Доказательство. Из очевидных равенств p* = (p* - E (A(p*) - b)), p£ = (pe - E(ep£ + A (pe) - b)) вытекает оценка

Hp* - p£H2 < Hp* - p£ - E (A(p* - p£) - ep£)H2 ,

справедливая для любого числа E > 0. Переходя в этой оценке от норм элементов к скалярным произведениям, получим неравенство

(p* - p£,A (p* - p£) - ep£) ^ E ¡A (p* - p£) - ep£112.

Переходя в этом неравенстве к пределу при Е ^ +0, получим условие

Используя положительность оператора А и неравенство Коши—Буняковского [15], получим из условия (16) оценки

0 < (р* - р£,А (р* - р£)) < (р* - р£,ер£) < е 1Ы (||р* || - ||р£И), влекущие за собой неравенство

Ир£И < Ир*И.

Очевидно, что неравенства (16) и (17) выполняются для любого числа е > 0.

Предположим теперь, что {еп} — произвольная последовательность положительных чисел, для которой

lim en = 0.

n—>оо

Введя обозначение pSn = pn, покажем, что

lim \\pn - p*\\ = 0.

Это и будет означать справедливость доказываемой теоремы.

Для доказательства сильной сходимости последовательности {pn} элементов L2,(Q) к элементу p* £ L^(^) достаточно показать, что из любой подпоследовательности {pnk } последовательности {pn} можно извлечь другую подпоследовательность {pnk }, которая сильно сходится к элементу p*. Покажем это. j

Пусть {pnk } — произвольная подпоследовательность последовательности {pn}. Поскольку подпоследовательность {pnk } является ограниченной в силу неравенства (17), существует такая её подпоследовательность {pnk.} и такой

элемент p £ L^(^), что {pnk,} слабо сходится в L^fä) к p (см. [11, теорема 1.7]). Из оценки (16) при е = enk, вытекает неравенство

{p* - pnk. , A(p* - pnk. )) < enk. (p*, pnk. ) - e'nk. Ы

'nk . \\ ,

справедливое для всех значений индекса . Переходя в этом неравенстве к пределу при nkj ^ ж, получим неравенство

(р* - р, А (р* - р)) ^ 0,

которое в силу второго из условий (12) означает, что р* = р. Следовательно, {рПк.} слабо сходится к элементу р*. Но поскольку ещё из неравенства (17) следуют оценки

урпк. И < ур*И УЩз , подпоследовательность {рпк.} сильно сходится к р* (см. [11, теорема 1.4]).  316

Из доказанных теорем 2 и 3 следует, что уравнение (13) является регуля-ризирующим аналогом уравнения (11).

5. Дискретизация регуляризированного уравнения. Наличие слагаемого ер в уравнении (13) позволяет построить эффективный устойчивый численный алгоритм приближенного решения этого уравнения. Такой алгоритм основывается на следующей аппроксимационной теореме.

Теорема 4. Пусть b € g е U0, а линейный ограниченный опе-

ратор A: ^(П) ^ L2(^) удовлетворяет условиям (12). Пусть последовательность {An} линейных ограниченных операторов, действующих в а также последовательность {bn} элементов и последовательность

{gn} элементов U0 удовлетворяют следующим условиям:

lim ||A - An И* = 0,

га^-те

га^-те v 7

lim ||g - gnli = 0.

га^те

Тогда для любого числа е > 0 существует натуральное число п0, зависящее только от е, такое, что для всех натуральных п ^ п0 уравнение

Р =

является однозначно разрешимым в ^(П) и справедливо равенство

га^те

где pn, р£ е ¿2(П) являются решениями уравнений (19) и (13) соответственно.

Доказательство. Выберем для заданного числа е > 0 номер п0 и число E > 0 в уравнении (19) так, чтобы выполнялись следующие неравенства:

е 1

2 е +

Используя эти неравенства, можно получить оценку

F

Действительно,

X =1

2

г е

< 1 - Ее + Е- = 1 - Е- < 1. 2 2

Из неравенства (20) следует сжимаемость оператора Ега: (О) ^ сто-

ящего в правой части уравнения (19) (см. доказательство теоремы 2), которая

влечет за собой однозначную разрешимость этого уравнения в пространстве Ь2(О) для всех натуральных п ^ п0.

Докажем теперь сильную сходимость в пространстве Ь^(О) последовательности {рп} решений уравнений (19) к решениюре € Ь^(О) уравнения (13). Для этого рассмотрим последовательность {ап} элементов Е^(О), определяемую равенством

ап = Ап (рп) - А (рп)+ Ь - Ьп, п ^ по.

Учитывая условия (18) и ограниченность в Ь^(О) последовательности {рп}, имеющую место в силу ограниченности в Ь^(О) сходящейся последовательности {Ьп}, а также неравенства (20), получим для ап оценки

1КН < \\Ап - А\1 \\рп\\ + \\Ь - Ьп\\ < \\Ап - А\1  в + \\ь - Ьп\\, где в — некоторое положительное число. Из последних оценок следует, что

Иш \\ап\\ = 0.

п—^^о

Поскольку из теоремы 2 следует непрерывная зависимость решения уравнения (13) от элементов Ь € Ь2(О), д € Щ, из очевидных соотношений

Ре = С^д (ре - Е(ере + А (р£) - Ь)) , Рп = 0^дп (рп - Е(ерп + А (рп) - Ь + ап)) ,

а также из последнего равенства (18) и доказанного соотношения (21) следует,

что Иш \\рп - ре\\ =0.  п—<х

Из теорем 1-4 следует, что приближенное решение рассматриваемой контактной задачи на каждом шаге нагружения можно получить, решив уравнение (19) при достаточно малом значении е > 0 и достаточной близости оператора Ап и элементов Ьп, дп к оператору А и элементам Ь, д соответственно.

Построим такие удовлетворяющие условиям теоремы 4, аппроксимирующие последовательности {Ап}, {Ьп} и {дп}, для которых уравнение (19) можно решить численно с любой степенью точности. Для этого зададим область О в виде открытого квадрата, ограниченного прямыми, параллельными осям некоторой декартовой системы координат, введенной на общей для взаимодействующих тел касательной плоскости. Далее для каждого натурального п разобьем область О на п2 непересекающихся квадратных областей ш\,ш2,... ,шп2 равной площади, ориентированных подобно квадрату О.

Будем полагать, что в выражениях (8) линейные ограниченные операторы А.ц : Ь2(О) ^ Ь2(О) допускают интегральное представление

Ац(х)ь = Кц(Ь, 8)х(8)(1,8, Ь € О, 1,2 = 1,3 Jn

и являются вполне непрерывными. Зададим линейный вполне непрерывный оператор Ап: Ь3,(О) ^ Ь2(О) следующим образом:

ж = (Ж1,Ж2,Жз) , y = (У1,У2,Уз) е

3

y = ^^ yi = Y^ а

j=1

(Xj)t = ^ K^s^-(s)ds, t е П, i,j = 173.

Входящие в соотношения ( ) функции Kj°(í,s) определяются равенствами

L j s) mes(wfc) mes(wm) I \ I K j

I ( I Kj(t,s)ds)dt,

J u>k \ Jum J

если t е , s е wm; k, m = 1, n2; i, j = 1,3,

в которых mes(wfc) есть площадь квадрата . Элементы bn € ¿3(0) и gn € U0 зададим в виде b„ = (P„(bi),P„(b2), РЦЬз)), g„ = Pn(g), где (ЬьЬ2,Ьз) = b и оператор Pn: ¿2(0) ^ ¿2(0) определяется соотношениями

Pn(x)t =-1—т i x(s)ds, t € wfc, k = 1,n2.

mes(wfc ) ,/ Wfc

Очевидно, что в случае полной непрерывности линейных интегральных операторов Aj: ¿2(0) ^ ¿2(0) построенные аппроксимирующие последовательности {An}, {bn} и {gn} удовлетворяют условиям (18). Кроме этого, эти последовательности заданы так, что оператор : L;](0) ^ L;](0), стоящий в правой части уравнения (19), отображает пространство ¿3(0) в конечномерное подпространство = Pn (L2(0)) х Pn (L2(0)) x Pn (L2(0)) пространства ¿3(0). Следовательно, каждое решение € L;](0) уравнения (19) содержится в и может быть найдено из решения следующей системы 3n2 скалярных уравнений с 3n2 неизвестными:

жзк-2 = ^жзк-2 - E^^ Ö3fc-2jж^- - Ь'зк-^^j ; жзк-1 = ^жзй-1 - E^^ азк-Уж^' - bjjfc-^ ,

жз& - ¿ ^ж - ^, rn);

жзк = Л жзк - E í ^ азку ж^- - Ь'з^,

жзк-1 - ^ ¿ азк-1, ж - Ь'зк-^, rn);

k = 1, n2.

С т р е л я е в Ю. М.

Решение рп = (р1п,р2п,р3п) € Нп уравнения (19) выражается через решение х1,х2,... ,х3п2 системы (24) соотношениями х3к-2 = р1п(8), х3к-1 = р2п(8), хзк = рзп(8) при 8 € Шк Ук = 1,п2. Неотрицательные постоянные дк в этой системе являются известными, а ^ — коэффициент трения. Значения параметров Ь\, ац, дк в (24) с учетом (22), (23) определяются следующими соотношениями:

Ьзк-2 = шек) /,Шб8 ЬЬзк-1 = к(8>Л8-

Ь3к = ШЛ8' дк-ШЪЗ I1 Ф)"8'

а3к-й 3]-е = Ше8(Шк) ! ! К3-й 3-е(Ь,8)6^6Ь + е8зк-<1 3]-в,

к, 2 = 1,п2, 6,е = 0,2.

Здесь Ь = (Ь1,Ь2,Ь3), 5ц —символ Кронекера, а положительное число е — параметр регуляризации. Предложенный способ построения аппроксимирующих последовательностей {Ап}, {Ьп}, {дп} гарантирует, что решение рп € Ь2(О) операторного уравнения (19) (существование и единственность которого в Ь2 (О) при достаточно большом п вытекает из теоремы 4) непременно должно содержаться в конечномерном подпространстве Нп пространства Ь2(О). Это означает, что при выполнении всех условий теоремы 4 (а также при полной непрерывности линейных интегральных операторов Ац : Ь2(О) ^ Ь2(О)) система уравнений (24) является совместной, если только значение п достаточно велико. Таким образом, искомое решение рп уравнения (19) при достаточно большом п может быть восстановлено по найденным значениям неизвестных х1,х2,... ,х-3п2, удовлетворяющих системе уравнений (24).

Для многих контактных задач теории упругости используется аппроксимация взаимодействующих тел упругими полупространствами [10]. В этом случае для последовательности {Ап}, определяемой формулами (22), (23), можно вместо формул (23) использовать следующие соотношения [9]:

Кйе(^п, 8Ц), если Ь € ш., 8 € Шц, г = ]]

К^мМ 1 С „ ,п ,,

шв8

г / Кае(^п,8)й8, если Ь,8 € ш.

где г, 2 = 1,п2, 6,е = 1,3. В соотношениях (26) 8п — центр квадрата ш., а ядра Кц (Ь,8) интегральных операторов Ац определяются в соответствии с формулами Буссинеска—Черрути [10].

Таким образом, в случае аппроксимации взаимодействующих тел упругими полупространствами коэффициенты ац в системе (24) можно вычислить без использования соотношений (25) по следующим простым формулам:

а3к-2 3к-1 = а3к-2 3к = а3к-1 3к-2 = а3к-1 3к = а3к 3к-2 = а3к 3к-1 = ° азк-4 3ц-е = Шв8(шц)К3-а 3-е^, 8~п),

а3к-а 33к-а = е + К-3Ч 3Ч (8^, 8)6,8,

Ук,з = 1,п2 (к = з), У 6, е = 0,2. 320

Отметим, что при аппроксимации взаимодействующих тел упругими полупространствами линейные интегральные операторы А. : ¿2(0) ^ ¿2(0) являются вполне непрерывными [15], а порождаемый ими вполне непрерывный оператор А: ¿2(0) ^ ¿2(0) вида (8) удовлетворяет всем условиям теорем 1-4.

6. Итерационный процесс для решения дискретного аналога регуляризи-рованного уравнения. Из всех доказанных выше утверждений следует, что нахождение приближенного решения операторного уравнения (11), соответствующего определенному шагу нагружения рассматриваемой квазистатической контактной задачи, сводится к решению системы нелинейных скалярных уравнений (24) при достаточно малом е > 0 и достаточно большом натуральном п. Для приближенного решения системы (24) можно предложить следующий итерационный процесс:

(х(10),х20),...,х20„)2) € м3"2;

> ^2

(ш+1) _ , / ГШ) Е х3к-2 _ 'Ч х3к-2 Е

3га2

Е<

¿=1

3га2

(т+1) ( (т) 3^ (т) ,/ \

Х3/с-1 _ 9(х3*-1 - Е( 2^-1 . ) - Ь3к-1 ) ,

( 3"2 \ \

Е^^. -ь3И, №); \?=1 / /

„

^ - . -/ / (т) 7-1 / (т) т./

Ж^ _ ^ - Е^^к.X - Оч

3га2

Е<

¿=1

3к

3к- 1

/3 \ \

к _ 1, п2, т _ 0,1,2,...

Сходимость итерационного процесса (28) к решению системы (24) будет иметь место независимо от выбора начального приближения (ж!0),ж20),..., ) в 3п2-мерном евклидовом пространстве М3"2 при условии, что матрица Ао _

/ ч 3га2,3га2

_ (а.)г=1 ¿=1 является симметрической и положительно определенной, а входящая в (28) положительная константа Е удовлетворяет неравенству

Е < (

3га2

\1^3га2^ . ¿=1

-1

Этот факт следует из принципа сжимающих отображений.

7. Численные результаты. Предложенным методом было получено численное решение пространственной квазистатической контактной задачи о взаимодействии упругого шара радиуса 0.3 м с упругим полупространством. Процесс нагружения тел осуществлялся в соответствии со следующим законом

изменения жёстких смещений Ai¿, Аэ^ шара, выраженных в метрах:

-0.000002  i, если 1 ^ i ^ 50;

Ali = ^

-0.0002 + 0.000002  i, если 51 ^ i ^ 65;

A2i = A3i = 0, если 1 ^ i ^ 65.

График зависимости абсолютной величины нормального жесткого смещения шара |Ali| от номера шага нагружения i представлен на рис. 2. В соотношениях (29) величины A2¿ и A3Í есть касательные составляющие жесткого смещения шара в направлении осей Ox и Oy соответственно.

Рис. 2. Кривая нагружения [Figure 2. A loading curve]

В расчетах использовались следующие модули Юнга Е1, Е2 и коэффициенты Пуассона VI, ^ шара и полупространства: Е1 _ 3  105 МПа, Е2 _ _ 105 МПа, v1 _ v2 _ 0.2. Коэффициент трения ^ _ 0.12375. Значение безразмерного параметра в, определяемого согласно [10] формулой

_ 1 Е1 (1 - 2У2) (1 + У2) - Е2 (1 - 2У!) (1 + VI) в 2 Е1 (1 - V2) + Е2 (1 - V2) ,

составляет для рассматриваемой контактной задачи 0.1875. Отношение ^/в _ _ 0.66. При дискретизации зона контакта разбивалась на 41 х 41 _ 1681 квадратных граничных элемента одинаковой площади (п _ 41, шв8(^г) _ _ 7.84  10-8 м2). Для вычисления элементов а^. матрицы податливости использовались соотношения (27), в которых функции (з,^') задавались при помощи формул Буссинеска—Черрути, а значение е определялось соотношением

т—31 7/„га

е = 10-3 •/ Kn(sn,s)ds. Jwk

С помощью итерационного процесса (28) в настоящей работе были получены два численных решения контактной задачи, которые соответствуют различным этапам процесса нагружения. Первое решение было получено на

тридцать пятом шаге нагружения, что соответствует монотонному нормальному внедрению шара в полупространство на 0.00007 м (отрезок ОА, рис. 2). Второе решение, полученное на последнем, шестьдесят пятом, шаге нагру-жения (ломаная ОВС, рис. 2), соответствует первоначальному вдавливанию шара в полупространство на 0.0001 м (отрезок ОВ, рис. 2) и последующей разгрузке (отрезок ВС, рис. 2). Отметим, что одноименные значения относительных жёстких смещений тел на тридцать пятом шаге совпали с их значениями на шестьдесят пятом шаге процесса нагружения тел.

Полученные из первого решения значения касательных напряжений тхх, действующих в точках поверхности шара, соприкасающихся с полупространством по оси Ох, сопоставлялись с результатами решения соответствующей осесимметричной контактной задачи, приведенными в работе [10] для случая ^/¡3 = 0.66. В таблице показана зависимость безразмерной величины тгх/(вртах) от безразмерного параметра х/а.

Первый столбец таблицы соответствует решению [10], второй — численному решению, полученному на тридцать пятом шаге процесса нагружения (29), третий — численному решению, полученному на шестьдесят пятом шаге этого процесса. Здесь символ ртах обозначает максимальное значение контактного давления, а символ а — радиус площадки контакта. Данные, приведенные в первом и втором столбцах таблицы, показывают хорошее соответствие решения [10] и численного решения, полученного при монотонном на-гружении. Максимальная относительная погрешность в этом случае составляет приблизительно 6 % по отношению к наибольшему значению величины Тгх/(вртах), найденному согласно [10], что свидетельствует о корректности предложенного метода решения контактной задачи. Результаты, приведенные в третьем столбце таблицы, означают, что при разгрузке наблюдается уменьшение достигнутых значений касательных напряжений и смена знака этих значений вблизи границы пятна контакта.

Зависимость безразмерных касательных напряжений от безразмерного расстояния [The dependence of the dimensionless shear stresses on the dimensionless distance]

x/a Solution from [10] Solution by the Eqs. (28), (29) at the step 35 Solution by the Eqs. (28), (29) at the step 65

0.00 0.00 0.00 0.00

0.06 0.10 0.07 0.07

0.12 0.18 0.17 0.17

0.18 0.25 0.22 0.23

0.24 0.30 0.29 0.30

0.30 0.35 0.35 0.35

0.36 0.39 0.39 0.40

0.42 0.42 0.42 0.43

0.48 0.45 0.45 0.46

0.54 0.47 0.47 0.48

0.60 0.48 0.47 0.50

0.66 0.48 0.46 0.48

0.72 0.45 0.42 0.44

0.78 0.40 0.37 0.40

0.84 0.33 0.30 0.32

0.90 0.23 0.22 0.12

0.96 0.03 0.05 -0.14

Стреляев Ю. М.

На рис. 3 показаны конфигурации пятен контакта и зон сцепления, соответствующие полученным численным решениям. Здесь внешний замкнутый контур есть граница пятна контакта, заштрихованная область — зона сцепления, незаштрихованная область — зона проскальзывания. Рис. 3, а соответствует решению, полученному на тридцать пятом шаге, рис. 3, Ь — решению, полученному на шестьдесят пятом шаге процесса нагружения (29).

-0.005 0 х, m -0.005 0 х, m

a b

Рис. 3. Границы областей контакта и зон сцепления шара и полупространства при монотонном

loading]

Анализ представленных на рис. 3 результатов показал, что качественные характеристики контакта тел невозможно установить по конечным значениям жёстких смещений этих тел, не учитывая, каким образом эти конечные значения были достигнуты в процессе нагружения. Полученные данные также показывают, что после прохождения этапа увеличения нагрузки до максимального значения (отрезок AB, рис. 2) и последующей разгрузки (отрезок BC, рис. 2) граница центрального участка зоны сцепления несколько сместилась к центру и одновременно вблизи границы пятна контакта образовалась вторая кольцевая зона сцепления. Отметим, что подобную конфигурацию зон сцепления на контактном пятне получил J. R. Turner [16] при исследовании случая разгрузки в задаче о взаимодействии цилиндрического штампа, имеющего плоскую подошву, с упругим полупространством.

Выводы. Разработан метод получения приближённого решения трёхмерной квазистатической контактной задачи о взаимодействии двух линейно-упругих тел с учетом кулонова трения между ними при неизвестной заранее площадке контакта и неизвестной границе раздела зон проскальзывания и сцепления на этой площадке. Анализ полученных численных результатов для задачи о вдавливании упругого шара в упругое полупространство свидетельствует о корректности предложенного метода и его применимости к решению контактных задач рассматриваемого класса.

ORCID

Юрий Михайлович Стреляев: http://orcid.org/0000-0002-4400-7824

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Кравчук А. С. Вариационный метод в контактных задачах. Состояние проблемы, направления развития// ПИИ, 2009. Т. 73, №3. С. 492-502.

2. Reina S. A., Dini D., Hills D. A., Lida Y. A quadratic programming formulation for the solution of layered elastic contact problems: Example applications and experimental val-idation// European Journal of Mechanics - A/Solids, 2011. vol.30, no. 3. pp. 236-247. doi: 10.1016/j.euromechsol.2010.12.003.

3. Галанов Б. А. Метод граничных уравнений типа Гаммерштейна для контактных задач теории упругости в случае неизвестных областей контакта// ПИИ, 1985. Т. 49, №5. С. 827-835.

4. Галанов Б. А. О приближенном решении некоторых задач упругого контакта двух тел // Изв. АН СССР, ИТТ, 1981. №5. С. 61-67.

5. Александров В. М., Kalker I. I., Пожарский Д. А. Пространственная контактная задача для двухслойного упругого основания с заранее неизвестной областью контакта // Изв. РАН, ИТТ, 1999. №4. С. 51-55.

6. Чебаков М. И. Трехмерная контактная задача для слоя с учетом сил трения в области контакта// Изв. РАН, ИТТ, 2002. №6. С. 59-68.

7. Александров В. М., Пожарский Д. А. Трехмерные контактные задачи при учете трения и нелинейной шероховатости // ПИИ, 2004. Т. 68, №3. С. 516-527.

8. Александров А. И., Стреляев Ю. М. Метод нелинейных граничных интегральных уравнений для контактных задач теории упругости // Восточно-Европейский журнал передовых технологий, 2014. Т. 3, №7(69). С. 36-40.

9. Александров А. И. Метод решения пространственной контактной задачи о взаимодействии двух упругих тел при наличии трения между ними // Иатематичт методи i фiзико-механiчнi поля, 2013. Т. 56, №3. С. 29-42.

10. Johnson K. L. Contact Mechanics. Cambridge: Cambridge University Press, 19S5, xii+452 pp. doi: 10.1017/CBO9781139171731.

11. Александров А. И. Неподвижные точки непрерывных операторов в гильбертовом пространстве. Запорожье: Запорож. гос. ун-т, 2002. 77 с.

12. Тихонов А. Н., Арсенин В. Я. Иетоды решения некорректных задач. М.: Наука, 1986. 288 с.

13. Красносельский М. А., Вайникко Г. М., Забрейко П. П., Рутицкий Я. Б., Стеценко В. Я. Приближенное решение операторных уравнений. М.: Наука, 1969. 456 с.

14. Ахиезер Н. И., Глазман И. М. Теория линейных операторов в гильбертовом пространстве. М.: Наука, 1966. 543 с.

15. Канторович Л. В., Акилов Г. П. Функциональный анализ. М.: Наука, 1984. 752 с.

16. Turner J. R. The frictional unloading problem on linear elastic half-space // IMA J. Appl. Math., 1979. vol.24, no. 4. pp. 439-469. doi: 10.1093/imamat/24.4.439.

Поступила в редакцию 24/I/2016; в окончательном варианте — 27/IV/2016; принята в печать — 27/V/2016.

Vestn. Samar. Gos. Techn. Un-ta. Ser. Fiz.-mat. nauki

[J. Samara State Tech. Univ., Ser. Phys. & Math. Sci.], 2016, vol. 20, no. 2, pp. 306-327 ISSN: 2310-7081 (online), 1991-8615 (print) doi: http://dx.doi.org/10.14498/vsgtu1471

MSC: 74M15, 74M10

A NONLINEAR BOUNDARY INTEGRAL EQUATIONS METHOD FOR THE SOLVING OF QUASISTATIC ELASTIC CONTACT PROBLEM WITH COULOMB FRICTION

Yu. M. Streliaiev

Zaporizhzhya National University,

66, Zhukovskogo st., Zaporizhzhya, 69600, Ukraine.

Abstract

Three-dimensional quasistatic contact problem of two linearly elastic bodies' interaction with Coulomb friction taken into account is considered. The boundary conditions of the problem have been simplified by the modification of the Coulomb's law of friction. This modification is based on the introducing of a delay in normal contact tractions that bound tangent contact tractions in the Coulomb's law of friction expressions. At this statement the problem is reduced to a sequence of similar systems of nonlinear integral equations describing bodies' interaction at each step of loading. A method for an approximate solution of the integral equations system corresponded to each step of loading is applied. This method consists of system regular-ization, discretization of regularized system and iterative process application for solving the discretized system. A numerical solution of a contact problem of an elastic sphere with an elastic half-space interaction under increasing and subsequently decreasing normal compressive force has been obtained.

Keywords: elastic body, contact problem, Coulomb friction, quasistatic problem, integral equation, uniqueness of solution, regularizing equation, iterative process.

ORCID

Yurii M. Streliaiev: http://orcid.org/0000-0002-4400-7824

REFERENCES

1. Kravchuk A. S. The variational method in contact problems. The present state of the problem and trends in its development, J. Appl. Math. Mech., 2009, vol. 73, no. 3, pp. 351-357. doi: 10.1016/j.jappmathmech.2009.07.004.

2. Reina S. A., Dini D., Hills D. A., Lida Y. A quadratic programming formulation for the solution of layered elastic contact problems: Example applications and experimental

© 2016 Samara State Technical University. Please cite this article in press as:

Streliaiev Yu. M. A nonlinear boundary integral equations method for the solving of quasistatic elastic contact problem with Coulomb friction, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. & Math. Sci.], 2016, vol. 20, no. 2, pp. 306-327. doi: 10.14498/vsgtu1471. (In Russian) Author Details:

Yurii M. Streliaiev (strelkiny@gmail.com), Senior Lecturer, Dept. of Mathematical Analysis.

validation, European Journal of Mechanics - A/Solids, 2011, vol.30, no. 3, pp. 236-247. doi: 10.1016/j.euromechsol.2010.12.003.

3. Galanov B. A. The method of boundary equations of the Hammerstein-type for contact problems of the theory of elasticity when the regions of contact are not known, J. Appl. Math. Mech., 1985, vol.49, no. 5, pp. 634-640. doi: 10.1016/0021-8928(85)90084-X.

4. Galanov B. A. Approximate solution to some problems of elastic contact of two bodies, Mech. Solids, 1981, vol. 16, no. 5, pp. 61-67.

5. Alexandrov V. M., Kalker J. J., Pozharskii D. A. Three-dimensional contact problem for a two-layered elastic base with an unknown contact area, Mech. Solids, 1999, vol. 34, no. 4, pp. 41-45.

6. Chebakov M. I. 3D contact problem with friction for a layer, Mech. Solids, 2002, vol. 37, no. 6, pp. 48-56.

7. Aleksandrov V. M., Pozharskii D. A. Three-dimensional contact problems taking friction and non-linear roughness into account, J. Appl. Math. Mech., 2004, vol. 68, no. 3, pp. 463472. doi: 10.1016/S0021-8928(04)00061-9.

8. Aleksandrov A. I., Streliaiev Yu. M. Nonlinear boundary integral equations method for contact problems of the elasticity theory, Eastern-European Journal of Enterprise Technologies, 2014, vol. 3, no. 7(69), pp. 36-40 (In Russian).

9. Aleksandrov A. I. Method of solution of three-dimensional contact problem on interaction of two elastic bodies in the presence of friction between them, Matematichni metodi i fiziko-mekhanichni polia, 2013, vol. 56, no. 3, pp. 29-42 (In Russian).

10. Johnson K. L. Contact Mechanics. Cambridge, Cambridge University Press, 1985, xii+452 pp. doi: 10.1017/CBO9781139171731.

11. Aleksandrov A. I. Nepodvizhnye tochki nepreryvnykh operatorov v gil'bertovom prostranstve [Fixed points of continuous operators in the Hilbert space]. Zaporozhye, Zaporozh. State Univ., 2002, 77 pp.

12. Tikhonov A. N., Arsenin V. Yu. Solutions of ill-posed problems, Scripta Series in Mathematics. New York, John Wiley & Sons, 1977, xiii+258 pp.

13. Krasnosel'skij M. A., Vainikko G. M., Zabreiko P. P. Rutitskij Ya. B., Stetsenko V. Ya. Approximate solution of operator equations, Wolters-Noordhoff Series of Monographs and Textbooks on Pure and Applied Mathematics. Groningen, Wolters-Noordhoff Publ., 1972, xii+484 pp.

14. Akhiezer N. I., Glazman I. M. Theory of linear operators in Hilbert space, vol. I, II, Monographs and Studies in Mathematics, vol. 9, 10. Edinburgh, Scottish Academic Press, 1981.

15. Kantorovich L. V., Akilov G. P. Functional analysis. Oxford, Pergamon Press, 1982, xiv+589 pp.

16. Turner J. R. The frictional unloading problem on linear elastic half-space, IMA J. Appl. Math., 1979, vol.24, no. 4, pp. 439-469. doi: 10.1093/imamat/24.4.439.

Received 24/I/2016;

received in revised form 27/IV/2016;

accepted 27/V/2016.