CC BY
27
2440
Фундаментальные проблемы теоретической и прикладной механики Вестник Нижегородского университета им. Н.И. Лобачевского, 2011, № 4 (5), с. 2440-2442
УДК 539.3
К ВЫБОРУ ФИЗИЧЕСКИ ОБОСНОВАННЫХ РЕШЕНИЙ АВТОМОДЕЛЬНЫХ КРАЕВЫХ ЗАДАЧ НЕЛИНЕЙНОЙ ДИНАМИЧЕСКОЙ
ТЕОРИИ УПРУГОСТИ
© 2011 г. Д-А. Потянихин, О.В. Дудко
Институт автоматики и процессов управления ДВО РАН, Владивосток
potyanikhin@mail.ru
Поступила в редакцию 24.08.2011
Уравнения нелинейной теории упругости относятся к классу гиперболических систем, выражающих законы сохранения. Известно, что существуют такие гиперболические системы уравнений, для которых построение решений автомодельных задач с использованием непрерывных функций и ударных волн оказывается неоднозначным [1-3]. Для выделения единственного, физически обоснованного решения необходимо учитывать дополнительные ограничения. В настоящем исследовании выбор единственного решения из числа математически возможных связан с термодинамическим условием совместности сильных разрывов, выражающим неубывание энтропии при необратимых процессах на ударных волнах, а также с условием эволю-ционности ударных волн.
Ключевые слова: нелинейная упругость, ударная волна, автомодельная волна Римана, контактное взаимодействие.
В твердых телах скорость распространения тепла существенно меньше скорости распространения ударных возмущений. Поэтому в задачах динамики деформирования упругих тел процессом теплопередачи часто пренебрегают, полагая коэффициенты теплопроводности равными нулю. Однако исключение из полной системы соотношений теории термоупругости уравнений, связанных с балансом температуры, приводит к неединственности решения даже простейших краевых задач, допускающих автомодельную постановку. Иногда удается удовлетворить начальным и граничным условиям задачи, предполагая появление в среде различных комбинаций волновых фронтов, вызванных ударным нагружением. С математической точки зрения все такие решения равноправны, ни одному из них нельзя отдать предпочтение.
Предлагается способ, позволяющий определять единственное решение плоских автомодельных задач динамики деформирования нелинейной упругой среды при помощи дополнительных критериев — условия эволюционности и условия неубывания энтропии в необратимом процессе на ударной волне. Этот способ продемонстрирован на решении плоской автомодельной задачи об ударе абсолютно твердым телом по упругому полупространству [4].
Динамическое деформирование среды нелинейно-упругого тела представляется в перемен-
ных Эйлера соотношениями:
V, = йг + } , Wt = V, + ]■,] ,
2аг] = иг,] + и], - ик,гик,] , г],] = Р™г =
р дЖ
= -
р0 да
гк
<^к] - 2ак])
(1)
Ро
4
8
=, 11 - 2/, + 212 - 212 - -13 + 4/, 12 - -1
3
[1
Ч±2
3
[3 :
Л
Ж = -/2 = ц/2 + К/1 /2 + Х/13 +П/3 + •••,
/1 = акк, /2 = а гкакг, /3 = агкак]а ]г ■
В соотношениях (1) иг, wi — компоненты
векторов перемещений, скоростей и ускорений
точек среды; аг], О,, — компоненты тензора дефорУ У
маций Альманси и тензора напряжений Эйлера-Коши; р0, р — плотность материала в начальном и текущем состоянии соответственно; 5] — символ Кронекера. При адиабатическом деформировании изотропной среды упругий потенциал Ж зависит от инвариантов тензора деформаций / / /3; X, Д — параметры Ламе; к, %, п — упругие модули третьего порядка.
На волновых фронтах скачки параметров напряженно-деформированного состояния и движения точек среды связаны динамическими и кинематическими условиями совместности [5]. Условие неубывания энтропии на поверхности сильных разрывов приводит к термодинамическому
Р
условию совместности [6], аналогом которого в газовой динамике является теорема Цемплена:
а+[V, ]п]-р+ (У+П] -с/т |>0. (2)
j j
Пусть упругое тело I занимает полупространство х1 > 0. Его границу обозначим через Ь1. Аб-солютно твердое тело II с плоской границей Ь2, двигаясь с постоянной скоростью V) = {Vl0 ,V20 } , сталкивается с упругим телом так, что появляется общая граница двух тел ОЬ (рис. 1, 2). На участке ОЬ выполняется закон сухого трения.
Рис. 1
Рис. 2
Скорость движения точки О (начала подвижной системы координат) определяется соотношением S = ^0 + v20ctgф. Полагаем эту скорость достаточно большой, чтобы возникающие при столкновении волновые фронты не могли отделиться от точки О. Для этого необходимо считать скорость соударения v0 достаточно большой, а угол ф — малым. Ввод автомодельной переменной 0 = = х1/( — х2) и представление компонент вектора
перемещений в виде и1 = ^ — Х2)Д0), и2 = (5У — — Х2)Ф(0), и3 = 0 позволяют перейти в системе определяющих соотношений (1) к однородной системе обыкновенных дифференциальных уравнений относительно неизвестных функций F(0) и Ф(0). Анализ решений этой системы показывает, что волновая картина может состоять из двух ком -
бинаций плоских ударных фронтов и простых волн Римана (см. рис. 1, 2), причем из-за отсутствия предварительных деформаций передним фронтом всегда оказывается продольная ударная волна S1 (О1 = tgP1). Следом за ней может распространяться либо квазипоперечная ударная волна Х2, отвечающая значению О2 = tg^2 (см. рис. 1), либо центрированная волна 0є [0+,0- ] (см. рис. 2).
Получив решение задачи при всех возможных постановках, проверяем выполнение термодинамического условия совместности (2) на ударных волнах. Если соотношение (2) не выполняется при некоторой постановке, то такую постановку исключаем из числа возможных. Еще одним ограничением на существование ударных волн является условие эволюционности. Если в результате сопоставления двух решений оказывается, что фронт ударной волны Ъ2 занимает промежуточное положение внутри веера простой волны 0 є [0+,0- ], то считаем ударную волну £2 неэволюционной и выбираем решение с волной Римана.
Серия вычислительных экспериментов позволила изучить зависимость напряженно-деформированного состояния от параметров соударения v10, v20, ф и коэффициента трения к на границе OL. Оказалось, что при различных наборах параметров реализуются обе волновые картины (см. рис. 1, 2). При этом каждая волновая картина может соответствовать как жесткой сцепке двух тел, так и проскальзыванию на границе OL.
Список литературы
1. Буренин А. А., Лапыгин В .В. Об отражении пло ской продольной ударной волны постоянной интенсивности от плоской жесткой границы нелинейной упругой среды // ПМТФ. 1985. №5. С. 125-129.
2. Агапов И.Е., Буренин А. А., Резунов А.В. О соударении двух нелинейно-упругих тел с плоскими границами // Прикладные задачи механики деформируемых сред. Владивосток: ДВО АН СССР 1991. С. 206-215.
3. Kulikovskii A.G., Chugainova A.P., Sveshniko-va E.I. Nonuniqueness of solution to nonlinear of the elasticity theory // Journal of Engineering Mathematics. 2006. Vol. 55, No 1-4. P. 97-110.
4. Дудко О.В., Потянихин Д. А. О косом ударе жес -тким телом, имеющим плоскую границу, по нелинейному упругому полупространству // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2009. Т. 9. Вып. 4, ч. 2. С. 32-40.
5. Быковцев Г.И., Ивлев Д.Д. Теория пластичности. Владивосток: Дальнаука, 1998. 523 с.
6. Буренин А. А., Чернышов А. Д. Ударные волны в изотропном упругом пространстве // ПММ. 1978. Т. 42. Вып. 4. С. 711-717.
ON THE CHOICE OF PHYSICALLY FEASIBLE SOLUTIONS OF SELF-SIMILAR BOUNDARY PROBLEMS OF NONLINEAR DYNAMICAL ELASTICITY THEORY
D.A. Potyanikhin, O. V Dudko
Equations of nonlinear elasticity theory belong to hyperbolic systems resulting from conservation laws. It is well known that for certain hyperbolic systems of equations the solutions of self-similar boundary problems composed of continuous functions and shock waves are non-unique. In order to choose unique physically feasible solution it is necessary to take into account additional contingencies. In this paper the choice of unique solution from among all the probable solutions is associated with the thermodynamical compatibility condition and the existence condition of evolutionary shock waves. The former condition follows from the second law of thermodynamics for shock waves.
Keywords: nonlinear elasticity, shock wave, self-similar Riemann wave, contact interaction.
