Контакт пластин, жесткие включения в которых выходят на границу Текст научной статьи по специальности «Математика»

ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

2011 Математика и механика № 3(15)

УДК 517.95

Т.А. Ротанова

КОНТАКТ ПЛАСТИН, ЖЕСТКИЕ ВКЛЮЧЕНИЯ В КОТОРЫХ ВЫХОДЯТ НА ГРАНИНУ1

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

Ключевые слова: вариационное неравенство, жесткое включение, пластина Кирхгофа - Лява, контактная задача, выход на границу.

Данная работа относится к изучению класса задач об одностороннем контакте упругих тел, в частности к контактным задачам для упругих пластин с жестким включением. Исследование задач о контакте упругих тел, содержащих жесткие включения, в настоящее время представляет огромный интерес в связи с активным изучением в последние годы композитных материалов.

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

Рассмотрим односторонний контакт двух упругих пластин под действием внешних сил (модель Кирхгофа - Лява [5]). Пластины расположены под углом а

друг к другу, где ае (0,п]. В естественном состоянии пластины контактируют

по прямолинейной незамкнутой кривой ненулевой меры у. Такая задача исследовалась в работах [6 - 8]. Особенностью данной работы является то, что каждая из пластин содержит жесткое включение, выходящее на область возможного контакта, и, таким образом, выполняется принцип виртуальных перемещений. Кроме того, предполагаем, что жесткое включение в верхней пластине выходит на её границу.

1 Работа выполнена при финансовой поддержке ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009 - 2013 годы (П597), гранта РФФИ № 10-01-00054 и гранта Президента РФ (МК-222.2010.1).

Рис. 1. Контакт пластин, жесткие включения выходят на границу

1. Постановка задачи

Пусть заданы ограниченные области О, О с К2 с гладкими границами Г и дв,

соответствующие срединным плоскостям контактирующих пластин. Пусть у с дО, упГ = 0 и у0 = (дО) \ О, тогда дО = у и у0. При этом нижняя пластина О деформируется в своей плоскости, а точки верхней (горизонтальной) пластины О допускают перемещение только в вертикальном направлении. Обозначим ц = (дъд2) - вектор внешней нормали к границе Г, V = ^,у2) - вектор нормали к у, расположенный в плоскости верхней пластины О, п = (п1,п2) - единичный вектор внутренней нормали к дО, расположенный в плоскости нижней пластины О. Кро-

заданные функции, описывающие действие внешних сил.

Разобьем область возможного контакта пластин на три части: у = у2 и у2 и у3. Каждая из пластин содержит жесткое включение (рис. 1). Для верхней пластины жесткое включение - это подобласть ю1 с гладкой границей 5ю1 = у2 и у3 и Е1, выходящая на внешнюю границу Q. При этом Q \ ю1 соответствует упругой части пластины. Для описания перемещения точек области ю1 введем пространство жестких перемещений:

Z(ro1) = {/11 (x) = a0 + a1 x1 + a2x2, ai = const, i = 0,1,2; x = (x1, x2) ею1}. (1)

Таким образом, перемещения точек жесткого включения представляют собой элементы пространства аффинных непрерывных функций. Будем предполагать, что на Г выполнено условие жесткого защемления пластин, тогда перемещения подобласти ю1 жесткого включения l0 = 0.

Жесткое включение в нижней пластине - подобласть ю2 с G с границей

9ю2 = у2 иЕ0, где Е0 является кривой класса С01. Таким образом, G\ ю2 соответствует упругой части пластины. С целью описать перемещение точек области ю2 введем пространство инфинитезимальных жестких перемещений:

R(®2) = {р = (Р1,Р2) |р(x) = Cx + D, x еи2},

ме того, wv

dw• Пусть также g=(gbg2), gt e L2(G), i=1,2; f e L2(Q)

где

где c, d1, d2 = const.

Будем рассматривать пространство Щ (G) х Н^ (Q), где

Я 1о (G) = {u е [Я\G )]2|u = 0 на у о},

H02(Q) = {w е Я2(Q) | w = — = 0 на Г}.

dq

Введем Кш ^ - множество допустимых перемещений:

кffli ,ffl2 = {(u, w) е Я;о (G) х Но2 (Q) | un sin а + w > 0 на у15

un > 0 на Y2 UT3; w |ш1 = 0 U |ш2 е ^(ю2)}-

Введем обозначения для изгибающего момента и перерезывающей силы верхней пластины:

d2 d д2

m(w) = qAw + (1 -q)—-, tv(w) = — (aw + (1 -q)—-), 5 ^з^) = (“v2,vi),

dv2 dv ds¿

где q - коэффициент Пуассона верхней пластины Для описания нижней пластины введем тензор модулей упругости B={bjkl}, i, j, k, l = 1, 2. Пусть имеет место положительная определенность коэффициентов bijkl е L (G):

bjki jJ5|2 V^j =| i, С0 > 0, (2)

а также их симметричность: b^ = bklij = bjikl , i, j, k, l = 1, 2. Также введем соответственно тензоры деформаций и напряжений e(u)={e,j(u)}, c(u)={cj(u)}, i, j = 1, 2,

стп = (c jUj,CT2¡n¡), sj.(u) =1 (Mi,j + мj,i), i, j = 1,2.

Полагаем, что с = с^пр, ст = en - cnn, ст = (cx1,cx2). Все величины с двумя нижними индексами предполагаются симметричными по этим индексам; по повторяющимся индексам проводится суммирование. Пусть функции w(y) и u(х) = (м1 (х), м2(х)), y = (y1,y2) е Q, х = (x1,x2) е G описывают перемещения точек верхней и нижней пластин соответственно.

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

функции u в G, w в Q такие, что выполняется

-div(Bs(u)) = g в G \ Й2, a2 w = f в Q\ ^ ; (3)

u = 0 на у0, w = wq = 0 на Г; (4)

u|®2 = Р0,гдеР0еR(ro2)> w|®1 =0; ()

[w] = [wv ] = 0 на у ; (6)

un sin а + w >0 на у1; un >0 на у2 иу3; (7)

[m(w)] = 0, [tv (w)] > 0, [tv (w)]n sin а =-cn, [tv (w)](un sin а+ w) = 0 на у1; (8)

cn < 0, cT = 0 un • cn = 0 на у3; (9)

- | an-y-J gy>-J [tv (w)](yn sina + ф) V(y^) e K^ ; (10)

I0 иуз ffl2 Yi

-J an-Po -J gPo = °- (11)

E0 ®2

Уравнения (3) суть уравнения равновесия, справедливые в области упругости нижней и верхней пластины соответственно. В области G \ ю2 выполняется закон Гука с= Be(u). Условия (4) обеспечивают жесткое защемление пластин. Соотношения (5) указывают на характер перемещений в жестких подобластях ю2 и ю1. Условие (4) - условие склейки на области возможного контакта. Неравенства в (7) есть условие непроникания пластин. Краевые условия (8) на упругой части контакта имеют вид системы равенств и неравенств. Соотношения (9) - условия

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

Отметим, что в областях жестких включений уравнения равновесия не выполняются, несмотря на то, что внешние силы приложены ко всем точкам пластин.

2. Вариационная формулировка

Вариационный подход позволяет исследовать вопросы существования и единственности решения. Для того чтобы его использовать, сформулируем задачу (3)-(11) как задачу минимизации функционала энергии пластин

E(u, w) = 2 J a(u)e(u) - J gu +1 an (w, w) - J fw

G G Q

на множестве допустимых перемещений Кш ш . Здесь aQ (w, v) - билинейная форма следующего вида:

aQ (w,v) = J (w,11v,11 + w, 22v,22 +^(w,11v,22 + w, 22 v,11) + 2(1 - ?)w,12v,12).

Q

Задача минимизации функционала inf E(u, w) эквивалентна следующему

(u,w)eK®1 ,®2

вариационному неравенству:

(u, w) e Kffl1,ffl2; (12)

J a(u)e(u - u)-J g(u - u) + aQXS1(w,w - w)-J f(w - w) > 0 V(u, w)e Kffl1,. (13)

g\®2 g q

Покажем, что выполнены все условия теоремы Вейерштрасса, тогда задача минимизации, или, что то же, вариационное неравенство (12), (13), имеет решение. Исходное пространство рефлексивно.

Множество Kra ш выпукло и по теореме вложения замкнуто, следовательно,

является слабо замкнутым.

Из выпуклости и непрерывности функционала E(u,w) вытекает его слабая по-лунепрерывность снизу. Последовательно применяя неравенство (2), первое неравенство Корна [9] и оценку функций в пространстве hQ (Q) [10], получаем

(ст(и), є(и))е + ап (м, м) > с (є(и), є(и))е + с2 (м „■, м„-) >

> с3

4,(0) +с4

м н 2

1 "я02( □)

Здесь и далее используемые константы положительны. В силу неравенства Коши имеем оценку

(g, и)о +(/, м)о ^ с5

4(°) +Сб

V N 2 1 '»¿(О)

Полученные оценки обеспечивают коэрцитивность функционала Е(и,м>). Таким образом, условия теоремы Вейерштрасса выполнены.

Решение задачи (12), (13) единственно. Чтобы убедиться в этом, предположим, что (и1^1), (и2,^2) - решения задачи (12), (13), и подставим в (13) (иь^1) в качестве решения (и,^), а в качестве пробной функции берем (и2,^2) и наоборот. Путем сложения полученных неравенств извлекаем (иь^1)= (и2,^2).

Ниже покажем эквивалентность вариационной и дифференциальной постановок задачи, предполагая, что функции (и,^) являются достаточно гладкими.

Пусть (и,^) - решение вариационного неравенства (12), (13). Чтобы получить уравнения равновесия пластин (3), необходимо подставить в (13) в качестве пробных функций (и+ф,^+у) и (и-ф, ^-у), где (ф,у) е [С° (О \ Ю2)]2 х С° (О\ Ю1) и

ф, у продолжаются нулем в ю2 и Ю1 соответственно. Уравнения равновесия будут выполнены в смысле распределений. Краевые условия (4) - (7) вытекают из определения множества допустимых перемещений ш .

Для того чтобы вывести краевые условия (8) - (11), нам потребуется формула Грина [10], справедливая для произвольной области О с границей Г класса С1,1, к которой V - внешняя нормаль:

| фд2м = аа(м,ф) + |Ґ4(м)-ф-|т(м)~фу Уфє Н2(0). (14)

□ гг

Предполагаем, что кривую у можно продолжить до кривой уиЕ2 так, что область О разбивается на две упругие подобласти О1 и О2 и жесткое включение ю1, как указано на рис. 2. Применяя формулу Грина к □1 иП2 и учитывая склеивание функций и их производных на Е2, получим следующее выражение:

/ фЛ м aQs^—l (м,ф)-| [^ (м)]ф + | [ш(м)]фу, (15)

□\—1 71 71

где [tv(w)] = tv(w)+ - tv(w)-, [m(w)] = m(w)+ - m(w)-, а значения tv(w)±, m(w)±

соответствуют положительному и отрицательному (по отношению к нормали v) берегам.

Также нам потребуется формула Грина, которая справедлива для области G\ю2 с внутренней нормалью п к границе:

| ст(и)е(у) = - | divCT(u)y- | стп-у. (16)

G\ffl2 G\®2 d(G\Ш2)

В (15), (16) пробные функции выбираются из множества допустимых перемещений.

Вывод условий (8). Чтобы получить условия [m(w)] = 0, [tv (w)] > 0 на y1, подставим в вариационное неравенство (12), (13) тестовую функцию в виде (uw+ф), где феНО(^), ф> 0 на уь Ф = 0 на raí:

ап\ш1(^,Ф) > j f Ф-0\®1

Далее применим к левой части неравенства формулу (15) и, с учетом уравнения равновесия, имеем

J [tv (^)]ф-J [m(w)^ > 0.

Y1 Y1

В силу отсутствия ограничений на фу, получаем требуемое.

Для вывода следующего условия на у1 подставим в (12), (13) пробные функции (и±у,^±ф), такие, что (у,ф)еUffl2 и y|ffl2 = 0, кроме того, предполагаем, что

ф + уп sin а = 0 на у. Имеем

J ст(и)е(Ж) + апх-(^,ф) = J gy + J f^

G\Ю2 G\Ю2 Ю1

Воспользовавшись формулами Грина (15) и (16), а также уравнениями равновесия (3), получим следующее соотношение:

- J стп - у + J [tv (^)]ф = 0,

Y1 Y1

откуда, с учетом выбора пробных функций, извлекаем

- J (стп + [tv (w)]n sinа)у = 0.

Y1

Таким образом, на у1 выполнено условие [tv(w)] п sin а = -сп.

Теперь покажем выполнение условия [tv(w)] (ип sin а+w) = 0 на у1. Выведем это равенство локально. Рассмотрим произвольную точку х0 еу1 и предположим, что в некоторой окрестности B выполняется ип sin а+w > 0, тогда можно взять

тестовые функции (и ± v, w±£), где (v, |) е [C“ (G)]2 х C“ (Q) - произвольная

функция, такая, что supp(v, £) с B и в Bx<¡ выполнено условие

(и±у)п sin а + (w±£) > 0. Подставляя тестовые функции в (13) и применяя формулы (15) и (16), получаем

- J стп - v + J [tv (w)]| = 0.

Y1 Y1

Воспользовавшись условием [t(w)] n sin а = -en, находим

j [tv (w)] (vn sin a + |) =0,

Yi

откуда [tv(w)] = 0 на y1. Предполагая [tv(w)] > 0 на y1, имеем выполнение условия un sin a + w = 0, так как иначе, следуя предыдущим рассуждениям, придем к противоречию. Таким образом, условие (8) получено.

Вывод условия (9) на у3. Подставим в вариационное неравенство (13) тестовые функции (U, w) = (u + у, w) е K , у=0 на Y1 ^Y2. Применяя к полученному неравенству формулу Грина (16), находим, что -(an, у) > 0 . Последнее означает, что

СТт = 0 an ^ 0 на Y3 •

Теперь покажем выполнение условия uncn = 0 на у3. Предположим, что в окрестности Bнекоторой точки х0 е y3 выполняется un>0, тогда можно взять

тестовую функцию (u±v,w), где v е [С0°° (G)]2 - произвольная функция, такая, что suppv с Bx0 и (u±v)n > 0 в B Получим, что (an, v)Y3 = 0, откуда сп = 0 на у3.

В случае, когда сп > 0, имеем un = 0, так как иначе придем к противоречию.

Вывод условий (10) - (11). Подставим в вариационное неравенство пробные функции (u+yw+ф), (у,ф) е Кш ш и применим формулы (15) и (16). С учетом

выполнения уравнений равновесия, получим

- J an-y + J [tv (w)ft-J [m( w)]фv - J gy > 0.

Y1UY3 u2ö Y1 Y1 ®2

Воспользовавшись условиями (8) на y1 и (9) на у3, получим (10). Подставляя в (10) пробные функции (у,ф) = (0,0) и (у,ф) = 2(u,w), имеем (11).

Обратно, из (3) - (11) можно получить вариационное неравенство (12), (13).

3. Предельный случай

Выполним предельный переход, устремив жесткость пластины Q. к бесконечности. Для этого вместо уравнения

д2 w = f в Q \ Й1

рассмотрим семейство уравнений, характеризующихся параметром а:

1 2

—д w = f в Q \ ю1 , а> 0, а

и выполним предельный переход при а ^ 0 , который соответствует случаю, когда верхняя пластина переходит от упругого состояния к жесткому.

Для каждого а > 0 вариационная задача выглядит следующим образом:

(ua, wa) е Kffl1,ffl2; (17)

Ja(ua)е(й-ua)-Jg(U-ua) + aQ(wa,w- wa)-J f(w - wa)>0 V(U,п)еK^, • (18)

G G Q

Подставляя в качестве пробных функций (0,0), 2(ua,wa), извлекаем из (18):

I ст(и“ )е(и“) + aQ (wa, wa) = | gua +| (19)

о о □

Из (19) при 0<а<ао следует

II И“ Н21 + II Н22 < С7,

11 "я!0(О) " "И о2(П) 7’

Выбирая для сходящейся подпоследовательности прежнее обозначение (и“^а), имеем при а — 0 следующие сходимости:

w слабо в И0 (^),

иа —— и слабо в И1 (О).

Уо

Из (19) также следует оценка:

—1| wа ||22 < с8,

а 11 "ио2(а\®1) 8’

wa —— 0 сильно в HO (Q\ roj).

что означает сильную сходимость

V

Предельная функция w = 0 во всей области В силу выполнения для любого a > 0 условия непроникания, для предельной функции и имеем неравенство

un > 0 на у.

Введем обозначение

K = {u £ H1 (G) | un > 0 на у; и |Ш2 е Д(га2)}.

Подставив в (18) тестовую функцию (и, 0), где u е K, найдем

I CT(ua )е(й) -1 g(u - ua) +1 fwa> I CT(ua )e(ua) + — a“x_(wa, wa).

G \ roo G Q G \ So

Осуществляя предельный переход к нижнему пределу в последнем неравенстве, получим следующее вариационное неравенство:

u е K : I a(u)e(u - u) -| g(u - u) > 0 Vu е K. (20)

G\ro2 G

Задача (20) соответствует контакту нижней пластины G, содержащей жесткое включение ю2, с неподвижным препятствием.

Мы получили, что решение задачи (17), (18) при a — 0 сходится к решению задачи (20).

ЛИТЕРАТУРА

1. Алексеев Г.В., Хлуднев А.М. Трещина в упругом теле, выходящая на границу под нулевым углом // Вестник НГУ. 2009. Т. 9. № 2. С. 15-29.

2. Лойгеринг Г., Хлуднев А.М. О равновесии упругих тел, содержащих жесткие включения // ДАН. 2010. Т. 430. № 1. С. 1-4.

3. Хлуднев А.М. Задача о трещине на границе жесткого включения в упругой пластине // Изв. РАН. 2009. № 5. С. 98-110.

4. Khludnev A.M., Novotny A.A., Sokolowski J., Zochowski A. Shape and topology sensitivity analysis for cracks in elastic bodies on boundaries of rigid inclusions // J. Mechan. Phys. Sol. 2009. V. 57. No. 10. P. 1718-1732.

5. Вольмир А.С. Нелинейная динамика пластин и оболочек. М.: Наука, 1972.

6. Неустроева Н.В. Жесткое включение в контактной задаче для упругих пластин // СибЖИМ. 2009. Т.12. № 4. С. 92-105.

7. Неустроева Н.В. Односторонний контакт упругих пластин с жестким включением // Вестник НГУ. 2009. Т. 9. № 4. C. 51-64.

8. Хлуднев А.М. Об одностороннем контакте двух пластин, расположенных под углом друг к другу // Журнал ПМТФ. 2009. Т. 49. № 4. С. 553-567.

9. Темам Р. Математические задачи теории пластичности. М.: Наука, 1991.

10. Хлуднев А.М. Задачи теории упругости в негладких областях. М.: Физматлит, 2010.

Статья поступила 18.01.2011 г.

Rotanova T.A. CONTACT PROBLEM FOR PLATES WITH RIGID INCLUSIONS INTERSECTING THE BOUNDARY. This paper deals with the unilateral contact problem for two elastic plates located at an angle to each other. Both plates contain rigid inclusions intersecting the contact area. The rigid inclusion in the top plate intersects also its external boundary. The lower plate is deformed in its plane with the top plate being vertically deformed. Using a variational method the solvability of the problem is established. Assuming that the solution is sufficiently smooth, the differential statement being equivalent to the variational formulation is justified. We analyze the limit case corresponding to the unbounded increase of the bending rigidity of the top plate.

Keywords: variational inequality, rigid inclusion, Kirchhoff-Love plate, contact problem, boundary intersection.

ROTANOVA Tatiana Alexandrovna

(Lavrentyev Institute of Hydrodynamics, Siberian Branch of Russian Academy of Sciences) E-mail: t.stekina@gmail.com