ПРИКЛАДНАЯ МАТЕМАТИКА И МЕХАНИКА
Том 72. Вып. 4, 2008
УДК 539.3
© 2008 г. С. М. Айзикович, Л. И. Кренев, И. С. Трубчик
ДЕФОРМИРОВАНИЕ ПОЛУПРОСТРАНСТВА С ГРАДИЕНТНЫМ УПРУГИМ ПОКРЫТИЕМ
ПРИ ДЕЙСТВИИ ПРОИЗВОЛЬНОЙ ОСЕСИММЕТРИЧНОЙ НАГРУЗКИ
Предлагается аналитико-численный метод решения граничной задачи Неймана для упругого полупространства с градиентным упругим покрытием. Дается постановка задачи и излагается процесс построения фундаментального решения (функции Грина). Метод позволяет найти решение задачи для достаточно широкого класса законов неоднородности среды и аналитически исследовать эффекты, связанные с неоднородностью. Описывается процедура расчета полей смещений, напряжений и деформации. Особое внимание уделяется анализу механических характеристик в зоне перехода от покрытия к упругой подложке.
Неоднородные покрытия с плавно изменяющимися по глубине упругими свойствами, так называемые градиентные (или функционально-градиентные, ФГП) покрытия, дают возможность избежать концентрации напряжений в зоне сопряжения покрытия с упругой подложкой, поэтому применение таких покрытий - одно из перспективных направлений в трибологии. Однако развитые ранее математические модели слоистых материалов недостаточно точно описывают ФГП [13], наблюдаются как количественные, так и качественные отличия в поведении функционально-градиентных материалов по сравнению с однородными или слоистыми.
Среди неоднородных покрытий с плавно изменяющимися по глубине упругими свойствами особый интерес представляют ФГП, у которых градиент изменения упругих свойств меняет знак по одной из координат. Так, например, это покрытия, жесткость которых на поверхности максимальна и уменьшается с глубиной до некоторого значения, а затем, наоборот, увеличивается по мере приближения к подложке (покрытия с немонотонными законами изменения по глубине упругих свойств). В связи со сложностями математического моделирования сред, свойства которых меняются непрерывно хотя бы по одной из координат, в большинстве известных публикаций исследовались, в большинстве случаев, слоистые материалы (или кусочно-однородные).
Решение контактных задач для функционально-градиентного полупространства с произвольным законом изменения неоднородности по глубине было построено [4^9] двусторонним асимптотическим методом. Было доказано [10], что построенное приближенное аналитическое решение - асимптотически точное решение исходного уравнения как для малых, так и для больших значений характерного геометрического параметра задачи.
Внимание будет сосредоточено на анализе полей упругих смещений и напряжений, в частности, на исследовании влияния изменения упругих свойств материала на внутренние напряжения в поверхностных слоях материала со слоистым или функционально-градиентным покрытием при индентировании. Используемые методы обеспечивают устойчивые численные результаты для разных видов слоистых и функционально-градиентных материалов.
В зависимости от степени отличия упругих свойств покрытия от упругих свойств основания или от градиента упругих свойств в функционально-градиентном материале напряжения внутри материала имеют разный характер. Например, было показано [11, 12], что для некоторых материалов и условий работы самые высокие растягивающие напряжения появляются в зоне контакта покрытия с подложкой, а не на поверхности покрытия.
Ниже проводится анализ механических характеристик функционально-градиентных материалов в зоне перехода от покрытия к упругой подложке, в том числе и для задач контактного взаимодействия.
г
H
Фиг. 1
1. Постановка граничной задачи для непрерывно-неоднородного по глубине полупространства при заданных на его поверхности усилиях. Рассмотрим упругое полупространство, упругие характеристики которого непрерывно меняются с глубиной в пределах прилегающего к поверхности слоя толщиной Н, а затем стабилизируются и остаются постоянными. С полупространством свяжем цилиндрическую систему координат (г, ф, г). Будем полагать, что прилагаемые к полупространству усилия вызваны или воздействием штампа, имеющего форму выпуклого тела вращения, или распределенной внутри круга радиуса а силой р(г) и поэтому симметричны относительно оси координат (фиг. 1). Таким образом, смещения, деформации и напряжения не зависят от угловой координаты ф. Полагаем также, что коэффициенты Ламе М = М(г) и Л = Л(г), связанные с модулем Юнга Е и коэффициентом Пуассона V соотношениями
М Е Л =_EV_. Е М ( 2М +3 Л) V = Л
2(1+ V)' (1+ v)( 1 - 2v)' М + Л ' 2 (М + Л)
изменяются по глубине в приповерхностном слое толщиной Н, а при г < -Н остаются постоянными:
М( г) = М (-Н) = М5, Л( г) = Л(-Н) = Л5, г < - Н
М( г) = Мс( г), Л( г) = Лс (г), -Н < г < 0
(1.1)
Мс (-Н) = М5, Лс (- Н) = Л5
шшЛ( г) > Л * > 0, тахЛ( г)<Л * , ттМ( г )> М* > 0, тахМ( г )< М*
Минимум и максимум берутся при г 6 0). Индекс 5 соответствует подстилающему однородному полупространству, а индекс с - неоднородному слою, Л*, Л*, М*, М* - постоянные.
В этом случае уравнения равновесия имеют вид
Э , N дтгг ^ Э , , д°г ^ Э , 2 , 2Этфг „ ,-,ч
эГ(г°г) +г аГ - °ф = 0, ЭГ(гтгг) +г ЭТ = 0, эГ(^ +г "эТ = 0 (1.2)
Первые два уравнения системы (1.2) описывают осесимметричное напряженное состояние, возникающее, например, под действием нормальной к поверхности нагрузки, а
последнее - уравнение равновесия полупространства, скручиваемого касательным усилием.
Используя закон Гука, связь упругих напряжений со смещениями можно записать в виде
и , п „,,и .„ ~.дц> .„ „ ди и дw
сг = + Л0, аф = 2М- + Л0, о7 = 2М (г) -=--+ Л0, 0 = -— + - + --—
г дг ф г г дг дг г дг
(1.3)
Тгг = М( д- + д-], ТгФ = М(ди- и], V = мди
*( ди + д^]
( д г + д г ),
Используя формулы связи компонент тензоров деформации и напряжений в осесим-метричной задаче, получаем систему уравнений Ламе
М(г)(у2и - и] + (М(г) + Л(г))|0 + М'(г+ | 1 = 0
М(г)У+ (М(г) + Л^))^ + 2М'(г) ^ + Л'(г)0 = 0 (1.4)
ог ог
2
= (г I-1+ К М' (г) = ^М^, Л' (г) = гдг( дг) дг йг йг
В силу изложенного выше поверхностная нагрузка прилагается в ограниченной области, а именно внутри круга радиуса а. Рассмотрим случай произвольной нормальной к поверхности г = 0 нагрузки
-р(г), 0 < г < а, 0, а < г
г = 0: ог( г, 0) = ' тгг(г, 0) = 0, 0 < г (1.5)
На границе неоднородного слоя с однородным полупространством, при г = -Я, должны выполняться условия сопряжения по смещениям и напряжениям:
< (г, -Я) = 0*( г, -Я), ТСГг(г, -Я) = ТСГг(г, -Я) (1 б)
г = -Я: (1.6)
иС( г, - Я) = и* (г, -Я), wc (г, -Я) = г, - Я)
При (г, -г) ^ ^ смещения, деформации и напряжения исчезают.
Таким образом, сформулирована задача Неймана для неоднородного полупространства: найти смещения, деформации и напряжения внутри полупространства, удовлетворяющие уравнениям (1.4), заданному распределению усилий на границе (1.5), (1.6) и условию затухания на бесконечности.
2. Построение фундаментального решения для неоднородного по глубине полупространства. Будем разыскивать решение для смещений и и ^ в виде интегралов Ганкеля
и (г, г) = -1 и (у, г) J1 (уг )у йу, ^ (г, г) = | Ж (у, г) J0 (уг )у йу (2.1)
0 0
Подставим данные выражения в систему (1.4) и приравняем нулю подынтегральные выражения. Получим систему обыкновенных дифференциальных уравнений, которую, используя векторное представление для трансформант
т
х = (х1, х2, х3, х4), х1 = и, х2 = и', х3 = Ж, х4 = Ж'
запишем в матричном виде, при этом явно выделив части, соответствующие поверхностному неоднородному слою и однородному основанию:
dx
dz
= Acxc, -H < z < 0
Ac =
0 1 0 0
2 К +1 YK-1 -n -Yn 2 -Y- к -
0 0 0 1
3-к Yn-_ к +1 Y IIIII2IIIIIII Y 1С-- 1 2 К - 1 Y К + 1 -n
dx_
dz
A
= A x , -ro < z < - H
Y
0
2K -I- 1
1-C- 1 0
0
Y
0 0 0
2 К - 1
0
-Y
К + 1
Y
К-1 1
0
К + 1
к = 1 + 2М/(М + Л), п = М'/М = Е'/Е Граничные условия (1.5), (1.6) принимают вид
(7-к)(к-1)
x4(Y, 0) - Y(К+Кx1(Y, 0)
2 E( 0 )(к + 1)
P(Y)
(2.2)
(2.3)
(2.4)
Yx3(Y, 0) + x2(Y, 0) = 0, xck(y, 0) = 4(y, -H), k = 1, ..., 4; P(y) = ¡p(p)J0(YP)Pdp
Общее решение системы (2.3) для однородного полупространства г < -Н, М' = Л' = 0, М > 0, Л > 0 записывается следующим образом:
3 'Уг 3 'Уг
х1(у, г) = (й 1 + уг^) е , Х2(У> г) = (¿1 + (1+ Уг)й2)Уе
з V г з V г
*э(7> г) = (й 1 + (- к + уг)¿2)е , Х4(у, г) = (йх + (1 - к + уг)¿2)уе'
где ¿1, й2 - произвольные функции параметра у.
Применяя метод модулирующих функций [13], ищем решение системы (2.2) в виде
хс(У, г) = й 1 (у)а:(у, г)етг + ¿2(7)^(У, г)етг (2.6)
Векторы а1(у, г), а2(у, г) определяются из решения задач Коши
(2.5)
da- c
-_■ = Aca; - ya¡, -H < z < 0,
dz
1, 2
(2.7)
22 z = -H: ax(Y, -H) = (1, y, 1, Y), a2(Y, -H) = (-yH,y - Y H, - к - yH, y - ky - Y H)
a
0
Постоянные d1, d2 определяются из условий (2.4), где ak¡ (у, г) (к = 1, 2, 3, 4) - к-я компонента вектора а;(у, г) (г = 1, 2). Имеем
¿1(7) = Р *(у)М2(у)А-1, d2(y) = -Р*(у) (у)А-1 А = М2 (у) N 1(у) - М^у) N 2(7), Р* (7) = -Р (7) / (2М (0) + Л( 0))
(7) = а^у, 0) - К+К 7 а^у, 0), Мк (у) = а2 (у, 0) + у а3 (у, 0)
Окончательно получаем следующие выражения для компонент вектора хс(у, г) при -Н < г < 0:
хк(у, г) = Р*(у)А-1Ак(у, г)в1г, к =1, ..., 4
к к (2.8)
Ак(у, г) = М2 (у)а1 (у, г) - М^у)а2 (у, г) Введем обозначения
1кт = 1кт(Г, г) = |хк(У, г)]т(чг)у¿у, /кт = ОГ г) = |хк(7> г)(7Г)72¿7
0 0
Тогда, используя формулы (1.3), (2.1), распределение смещений, деформаций и напряжений для неоднородного по глубине полупространства при действии на его поверхности произвольной осесимметричной нагрузки можно представить в виде
и (г, г) = -1П, V (г, г) = 131)
еДг, г) = -1 ^ГI п( г, г), ег( г, г) = 140, еф(г, г) = -Г I lo, е^ = -1 (/21 +1 ^)
, , 2М ( г ,, ЛТ1 1Т °г( г,г) =1-^ (у/40-(1-у) 110 + Г ^ (2.9)
°г(Г, г) = пМЬ((1 + 2V)Г111+ у(/40(Г, г) -110))
°Ф( Г, г) = 12М^(( 1- у) /40 - У110 ), тгг( Г, г) = -М (/21 +
3. Численный анализ решения для некоторых характерных распределений модуля
Юнга. Значения компонент полей смещений, деформаций и напряжений строятся с помощью численного интегрирования по формулам (2.9) в заданных точках (г, г) неоднородного по глубине полупространства.
Для дальнейшего прочтения статьи необходимо приобрести полный текст. Статьи высылаются в формате PDF на указанную при оплате почту. Время доставки составляет менее 10 минут. Стоимость одной статьи — 150 рублей.