CC BY
37
УДК 517.958
ОБ ОБРАТНОЙ граничнои задаче деформирования ПЛАСТИЧЕСКОГО СЛОЯ ПРИ ЕГО СЖАТИИ СО СКОЛЬЖЕНИЕМ
© В.Л. Дильман, А.Н. Дияб
Ключевые слова: пластический слой; плоская деформация; осадка; напряженно-деформированное состояние; гипотеза плоских сечений; линии тока; обратная граничная задача; система уравнений гиперболического типа.
Исследуется напряженно-деформированное состояние пластического слоя в процессе его сжатия при наличии скольжения между плоскими контактными поверхностями в условиях плоской деформации. Полученные аналитические зависимости позволяют находить все механические и геометрические параметры процесса либо в виде явных приближенных формул, либо на основе итерационных вычислительных схем.
Сжатие пластического слоя впервые изучалось в работе [1] и затем многими авторами. Обычно напряжения на контактной поверхности между слоем и сжимающими поверхностями неизвестны. Для определения сжимающего усилия необходимо решить обратную граничную задачу нахождения нормальных напряжений на контактной поверхности. Система уравнений этой задачи квазилинейна. При плоской деформации она имеет гиперболический тип. Рассматриваются различные математические модели процесса сжатия. Обычно они основываются, во-первых, на технологических особенностях: величине трения при скольжении контактных поверхностей, наличии деформации контактных поверхностей, вовлечении сжимающих тел в пластическое деформирование. Во-вторых, вид модели зависит от введения ограничений на классы функций, в которых ищется решение [2]. Последнее позволяет упростить математическую модель и доопределить обратную граничную задачу. Экспериментально обоснованным ограничением деформационного типа является известная гипотеза плоских сечений [2-5]
уу _ W (у). (1)
Здесь уу - скорость перемещения точки слоя в его поперечном направлении. Упомянутые ограничения обычно распространяются на часть слоя (которая для толстых слоев может и отсутствовать), не включающую окрестности свободных поверхностей. Около свободных поверхностей напряженно-деформированное состояние слоя определяется на основе решения задачи Коши для нелинейной системы уравнений в частных производных гиперболического типа методом характеристик [2, 6]. В работах [2-6] рассмотрен процесс поперечного растяжения слоя. Сжатие со скольжением контактных поверхностей относительно друг друга в рамках таких моделей ранее не рассматривалось. Целью данной работы является изучение напряженно-деформированного состояния пластического слоя в процессе его сжатия при наличии скольжения между плоскими контактными поверхностями. Когда слой достаточно тонкий или становится таким в процессе сжатия, применяется гипотеза (1). Полученные, явные и неявные, аналитические зависимости позволяют находить все механические и геометрические параметры процесса либо в виде явных приближенных формул, либо на основе итерационных вычислительных схем.
Напряженно-деформированное состояние пластической среды при плоской деформации в безразмерных переменных задается системой уравнений:
дах + дтху _ о дау + дтху _ о (2)
дх ду ' ду дх '
(ах - ау)2 + 4тХу = 4, (3)
дих диу о (4)
дх ду '
°х - у у _ дих/дх - диу/ду (_)
2тху дУх/ду + диу / дх'
Здесь ах,ау и тху - напряжения, ух,уу - скорости перемещений. Функции из уравнений (2) - (5) определены на сечении слоя толщиной 2к, к € (0; 1]; к _ Н/1 , Н и I - толщина (высота) и ширина (по контактной поверхности) пластического слоя с двумя осями симметрии (осями координат). Обозначим
Г (у) _ 0, 5Ж'' (у)/Ш '(у). (6)
Из (1), (3), (5) и (6) следует:
т _ У (у)х ; _ 2
ху + у 2(у)х2' х у + у 2(у)х2'
В работах [2, глава 3; 3] показано, что функцию У можно приближенно найти как решение уравнения У" + 4УУ' _ 0 при начальных условиях У(0) _ 0, У'(0) _ —а2/2 :
У(у)_ -(а/2^(ау). (8)
Параметрами математической модели являются: 1) толщина слоя Н, которая рассматривается как условное время; 2) длина слоя I; относительная толщина слоя к _ Н/1; 3) скорость сжатия (движения контактной поверхности к срединной плоскости слоя) ^ _ _ (Н)|; 4) параметр а, который характеризует возникающий в процессе деформирования прогиб поперечных (по слою) координатных линий; 5) т _ тахтху(х, Н) - наибольшее значение касательных напряжений; оно достигается на контактной поверхности всюду на некотором промежутке [х*; I], причем на этом участке касательные напряжения постоянны [2, гл. 3]; 6) ближайшая к поперечной оси симметрии слоя точка х* , в которой касательные напряжения достигают наибольшего значения т; в работе [2, с. 83] показано, что
х*/I _ 1 - 2к/(сов и + 8т и) _ 1 - 2к/\/1 + т, (9)
где и - угол поворота характеристики от свободной до контактной поверхности; формула (9) имеет место, когда вычисляемая по ней величина х* > 0 (например, если к < 0, 5); 7) коэффициент трения у; условием скольжения предполагается равенство
т _ ц\ау(х,Н)|, х € [х*; I]. (10)
Параметры а, т и х* являются внутренними и вычисляются через внешние параметры. Известно [2, с. 70], что при х € [х*; I]
т _ |тху(х,Н)| _ 8ш(2и), |уу(х,Н)| _ 1 + 2и + сов(2и). (11)
В начальный момент нагружения поле характеристик является простым однородным, всюду в слое | _ 2 , тху _ 0, и _ 0. Затем трение контактных поверхностей приводит к возникновению ненулевых касательных напряжений в слое и, как следствие, возникновению так называемого контактного упрочнения. В поле характеристик появляется участок центрированного поля (то есть сектор) с вершиной в особой точке - точке выхода контактной границы на свободную. Угол этого сектора и увеличивается, пока не достигнет
своего наибольшего значения и* , определяемого коэффициентом трения у. В дальнейшем, в процессе деформирования величина и* остается неизменной. Из (10) и (11) следует трансцендетное уравнение для вычисления и* :
sin(2u) = ^(1 + 2ш + cos(2u)). (12)
Из (11) следует, что наибольшие (по всем точкам слоя) значения касательных и нормальных напряжений остаются инвариантными в процессе деформирования. Для возникновения внутри слоя зоны, к которой можно применить гипотезу (1), необходимо, чтобы ввычисля-емая по формуле (9) величина x* > 0, что равносильно условию
к < sinш + cosш. (13)
С момента, когда начинает выполняться условие (13), начинается второй этап деформирования слоя. Так как, в силу (7), (8) и (11),
sin(2u*) = - QX*tg(Qh) =, (14)
л/4 + a2 (x*)2 tg2(ah)
из (14) следует трансцендентное уравнение для нахождения параметра a:
ax* tg(ah) = 2tg(2u). (15)
Запишем это уравнение в другой форме. Пусть Ф - функция, обратная к функции y = = x tg x, x € [0; п/2) . Тогда уравнение (15) равносильно уравнению
a =1ф(2Щ2и)
h у x* '
Здесь x* вычисляется по формуле (9), и находится решением уравнения (12). Из (6) и (8) следует уравнение для вычисления функции W :
0, 5W''(y)/W'(y) = -(a/2) tg(ay). (16)
Учитывая очевидное условие vy |y=0 = 0, получим зависимость величины скорости деформации в каждой точке от ее ординаты слоя в направлении внешнего усилия:
vy = — vh sin(ay)/sin(ah). (17)
Из гипотезы (1), условия несжимаемости (5) и следующего из симметрии условия |x=o = 0 можно получить скорость деформации точек слоя в направлении, ортогональном внешнему усилию:
vx = —vhax cos(ay)/sin(ah). (18)
Линии тока являются интегральными кривыми уравнения dy/dx = vy/vx . Подставляя сюда правые части уравнений (17) и (18) и интегрируя, получим уравнение линии тока, проходящей через точку (x0; y0) : x sin(ay) = x0 sin(ay0) .
Решением обратной граничной задачи на участке контактной границы при x € [0; x*] является функция для вычисления напряжений ay :
2a2y2 1
ay =--2 — 2--+ - in | cos(ax)| + 2 + c.
cos2(ax)y4 + a2y2 tg2(ax)(y4 + a2y2 tg2(ax) +2) 2
Для ее получения следует проинтегрировать уравнения (2), используя формулы (7) и (8). Постоянную интегрирования с можно найти, приравняв полученную функцию второму выражению в (11) при х _ х* .
Чтобы получить результаты в приближенной аналитической форме, следует найти явные аналитические зависимости для параметров и и а , которые в неявной форме заданы уравнениями (12) и (15). Например, приближенное решение уравнения (15) можно получить, используя следующие утверждения.
Теорема1. Функцию Ф можно представить в виде Ф(х) _ л/хф(-х), где функция ф аналитична на всей числовой прямой, причем
I / \ Л 1 4 2 16 3 256 4 ф(х)_1 + 3 х +45 х2 + 945 х3 + 127575 ^ + -
Теорема2. Функцию Ф можно представить в виде Ф(х) _ р(1/х), где функция р аналитична в интервале (- 1;1), причем
р(г)_ |(1 -г + г2 - (1 - П2 )г3 +
ЛИТЕРАТУРА
1. Прандтль Л. Примеры применения теоремы Г. Генки к равновесию пластических тел // Теория пластичности / под ред. Ю.Н. Работнова. М.: Издательство иностр. литературы, 1948. С. 102-113.
2. Дильман В.Л., Ерошкина Т.В. Математическое моделирование критических состояний мягких прослоек в неоднородных соединениях. Челябинск: Издат. центр ЮУрГУ, 2011. 276 с.
3. Дильман В.Л., Остсемин А.А. Напряженное состояние и статическая прочность пластичной прослойки при плоской деформации // Проблемы машиностроения и надежности машин. 2005. № 4. С. 38-48.
4. Дильман В.Л. О напряженно-деформированном состоянии при растяжении пластического слоя с двумя осями симметрии // Изв. РАН. MTT. 2001. № 6. С. 115-124.
5. Дильман В.Л., Носачева А.И. Математическое моделирование критических состояний пластического слоя // Вестник Тамбовского университета. Серия Естественные и технические науки. Тамбов, 2013. Т. 18. Вып. 5. С. 2502-2504.
6. Дильман В.Л. Напряженное состояние и прочность неоднородной пла-стической полосы с дефектом в более прочной части // Изв. РАН. МТТ. 2010. № 2. С. 89-102.
Поступила в редакцию 5 мая 2015 г.
Dilman V.L., Dheyab A.N. ON THE INVERSE BOUNDARY VALUE PROBLEM OF THE PLASTIC LAYER DEFORMATION UNDER ITS COMPRESSION WITH SLIDE
The stress-strain state of a plastic layer in the process of compression with slide between flat contact surfaces in a plane strain is under consideration. The obtained analytical dependences allow to find out all the mechanical and geometrical parameters of the process, either in the form of explicit approximate formulas or on the basis of iterative computational schemes.
Key words: plastic layer; plane strain; sediment; stress-strain state; the hypothesis of plane sections; streamlines; inverse boundary problem; the system of equations hyperbolic type.
Дильман Валерий Лейзерович, Южно-Уральский государственный университет, г. Челябинск, Российская Федерация, доктор физико-математических наук, доцент, заведующий кафедрой прикладной математики, e-mail: dilman49@mail.ru
Dilman Valery Lazerovich, South-Ural State University, Chelyabinsk, the Russian Federation, Doctor of Physics and Mathematics, Associate Professor, the Head of the Applied Mathematics Department, e-mail: dilman49@mail.ru
Дияб Аус Нидал, Южно-Уральский государственный университет, г. Челябинск, Российская Федерация, аспирант кафедры прикладной математики, e-mail: Aws.nth@gmail.com
Dheyab Aws Nidhal, South-Ural State University, Chelyabinsk, the Russian Federation, Post-graduate Student of the Applied Mathematics Department, e-mail: Aws.nth@gmail.com
УДК 517.958
КЛАССИФИКАЦИЯ КРИТИЧЕСКИХ СОСТОЯНИИ НЕОДНОРОДНОГО ПЛАСТИЧЕСКОГО СЛОЯ НА ОСНОВЕ ИССЛЕДОВАНИЯ ОДНОГО ФУНКЦИОНАЛЬНОГО УРАВНЕНИЯ
© В.Л. Дильман, Д.А. Трунова
Ключевые слова: неоднородный пластический слой; напряженное состояние; гипотеза разделения переменных; системы нелинейных уравнений в частных производных; функциональное уравнение.
Исследуется краевая задача для системы нелинейных уравнений в частных производных, моделирующая напряженное состояние неоднородного пластического слоя. Слой находится под растягивающей нагрузкой в условиях плоской деформации. В предположении разделения переменных для касательных напряжений задача сведена к некоторым нелинейным системам обыкновенных дифференциальных уравнений. Сведение основано на полученной в работе полной классификации решений некоторого чисто функционального уравнения.
Рассматривается математическая модель напряженного состояния неоднородного пластического слоя. Пластически деформируемые слои возникают при нагружении сварных соединений ( сварные швы, зоны термического влияния, диффузионные прослойки) и при осадке заготовок жесткими матрицами. В работе предполагается что слой имеет прямоугольную форму, расположен между жесткими участками соединения и находится под сжимающей или растягивающей нагрузкой в условиях плоской деформации. Цель работы — разработка вычислительной схемы нахождения напряженного состояния слоя в критический момент нагружения, а также получение в аналитической форме явных зависимостей напряжений от координат для некоторых характерных частных случаев.
Расположим оси координат по осям симметрии слоя [-1; 1] х [—к; к] , здесь 0 < к < 1 — толщина слоя. Математическая модель содержит систему нелинейных уравнений «пластического равновесия» гиперболического типа [1—5]:
ж + ^ = 0.^ + ^ = °' е. — -у)2+<4, = 4Я2 м
Условие на свободной границе Г в форме Сен-Венана:
^ -.¿У = 0. (2)
Неоднородность слоя определим функцией
Я (х,у) = и (ж)У (у), (3)
дифференцируемой по каждой переменной.
