CC BY
35
94 Вестник СамГУ — Естественнонаучная серия. 2006. №6/1(46).
МЕХАНИКА
УДК 539.3
ЛУЧЕВЫЕ РАЗЛОЖЕНИЯ В ИЗУЧЕНИИ ЗАКОНОМЕРНОСТЕЙ РАСПРОСТРАНЕНИЯ НЕПЛОСКИХ УДАРНЫХ ВОЛН1
© 2006 Е.А.Герасименко, В.Е.Рагозина2
На основе обобщения теории геометрических и кинематических условий совместности разрывов на случай задания движущейся поверхности разрывов в криволинейной системе координат излагается способ построения приближенных решений краевых задач ударного деформирования лучевым разложением в прифронтовой области. Обсуждаются нелинейные эффекты в закономерностях распространения одномерных цилиндрических ударных волн в зависимости от характера предварительных деформаций, построено приближенное решение с цилиндрической продольной расходящейся ударной волной.
Введение
Приближенный метод построения решений краевых задач для линейных и квазилинейных систем уравнений в окрестностях поверхностей разрывов, называемый лучевым методом, широко используется в теории пластического течения и в динамике деформирования. Его развитие непосредственным образом связано с исследованиями ученых Самарского государственного университета [1, 2]. Не случайно наиболее полный и наиболее квалифицированный обзор работ, обобщающий этот метод [3], был посвящен светлой памяти Г.И. Быковцева, долгое время возглавлявшего кафедру механики деформируемого твердого тела Самарского госуниверситета. Именно профессором Г.И. Быковцевым вместе с учениками была разработана теория рекуррентных условий совместности разрывов функций и их производных на движущихся поверхностях разрывов, обобщающая представления Т.Томаса [4] для геометрических и кинематических условий совместно-
1 Представлена доктором физико-математических наук, профессором Ю.Н. Радаевым.
2Герасименко Екатерина Андреевна ([email protected]), Рагозина Виктория Евгеньевна ([email protected]), Институт автоматики и процессов управления ДВО РАН, 690041, Россия, г.Владивосток, ул.Радио, 5.
сти разрывов. С использованием таких рекуррентных ограничений, накладываемых на разрывы геометрией и кинематикой движущейся поверхности разрывов, были построены приближенные решения целого ряда краевых задач идеальной пластичности [1], гиперболической теории теплообмена [2], динамики линейной [5] и нелинейной [6] упругой среды. Использовались эти зависимости и для изучения условий существования ударных волн [7, 8]. Позднее на данной основе был предложен способ построения приближенных решений краевых задач ударного деформирования [9, 10]. Оказалось, что с помощью построенных таким образом прифронтовых разложений решений возможно построение алгоритмов отслеживания положений волновых фронтов в процессе численного решения краевых задач ударного деформирования [11, 12].
Часто решение краевой задачи вследствие ее симметрии существенно упрощается при использовании соответствующей криволинейной системы координат. Однако теория, построенная Г.И. Быковцевым [13], изначально подразумевает прямоугольную декартову систему координат. Отдельные условия совместности в криволинейных координатах записаны в [14]. Ранее [15] нами были получены рекуррентные соотношения, аналогичные [13] в произвольной системе координат, что потребовало новых корректных определений для 6-производных по времени от функций, определенных как в пространстве, так и на движущейся поверхности. Настоящей статьей иллюстрируем использование этих соотношений при изучении закономерностей распространения одномерных ударных возмущений при наличии цилиндрической симметрии.
1. Условия совместности разрывов в криволинейных координатах
Для описания положения частиц сплошной среды и их движения в евклидовом пространстве Б3 введем две пространственные стационарные системы координат: произвольную криволинейную систему координат хг (г = = 1,2,3) и декартову систему координат Всюду в дальнейшем, если не оговорено иное, латинские индексы принимают значения 1,2,3. Ко- и кон-травариантный метрический тензор (§-у и а также символы Кристоф-феля 1-го и 2-го рода (Г,д и Гд) вычисляются по стандартным формулам:
_ дт^дт^ ¡к _ ;
8и~ дх^дхГ 88к1~°Р
г»- _ г - I + _ ^
]к~8 2 + ы дхт
В (1.1) и в дальнейшем принято суммирование по повторяющимся индексам. Рассмотрим в Б3 гладкую поверхность Е(г), положение которой с течением времени задается уравнениями х1 = х1 [у1, у2, ^, и каждая ее точка
движется в направлении своей единичной нормали с сохранением постоян-«12
ных значений у , у , так что
X (у«, г) = Оу1
(1.2)
В (1.2) и далее греческие индексы принимают значения 1,2; О — скорость поверхности Е в направлении нормали.
Во многих задачах механики и математической физики часто необходимо проследить изменение со временем различных по типу тензорных полей в данной точке Е(г). В качестве характеристики такого изменения используют 6-производную по времени. Исходно определение 6-производной записывалось для декартовой пространственной системы координат. Однако есть много случаев движения Е с ненулевой гауссовой кривизной, для которых особенности решения удобнее изучать в пространственных криволинейных координатах. Один из простейших примеров — движение одномерных цилиндрических или сферических волн. Для решения таких задач необходимо уточнение операции 6-дифференцирования, но в согласовании со следующими требованиями: 6-производная должна определять тензорный объект; различные виды 6-дифференцирования не должны противоречить друг другу; для предельного случая перехода от х1 к zi 6-производные должны вычисляться по известным классическим формулам [13]. В зависимости от типа тензорного поля 6-производная должна определяться по-разному.
Рассмотрим, к примеру, как вычисляется 6-производная для поля тензора, заданного всюду в рассматриваемой области и, в частности, на Е, т.е. тензора вида А](хк (уа, г), ^. Для этого проведем дифференцирование инвариантного объекта А]е3 е;, где е[ и е3 — ко- и контравариантные базисные векторы системы х'. Для такого объекта получим:
уа=еоп81
йА'. _]_
йг
уа=еопв1
7 л г ае] е]е; + А • —
3 Л
е;+
а1 3 йе' 3 Л
уа=еопв1
(дА]
~дГ
ЗА1.
_I
д хк
уа=еопв1
] тк т к]
Оук\ е] е;,
(1.3)
где при выводе учтено, что е; = е; х , е; = е; х и
де1_гк _ _Г,- к
().\ 1 ~ <)х 1 ~ '
Из (1.3) следует, что величина в фигурных скобках может быть названа компонентами 6-производной. Она складывается из мгновенного изменения в данной точке пространства и производной вдоль выбранной лучевой траектории.
'2(уа ]*1]*2"-] г '
ные только на Е, но имеющие пространственные индексы (например, вектор
Другой часто встречаемый случай — тензоры вида А]?1 '^"Ч (уа, г), задан-
+
нормали V1 (уа, г)). На основании соображений, аналогичных предыдущим, можно получить:
bAlll2."'h. дЛ'У2.-h. ,
ЛД-Jr _ JU2-Jr Jr'l Amh—is , , r's ihi2...m
{n
A"U2--is + . . . + p's An
gt dt \ mk jl j2-jr mk jl j2 ...jr
—Ym A'1'2'"'s — — rm A'1 i2'"is \Qyk
jlk mj2'.. jr jrk jl j2'-m\ '
(1.4)
Отметим, что (1.3) может быть сведена к (1.4), поскольку Aj^X (ya, t), t) = = Al. (ya, t), поэтому в качестве основной формулы будем использовать (1.4). Ее интерпретация: (1.4) определяет производную по времени вдоль лучевой траектории.
Для тензоров, которые заданы только на Е и имеют только поверхностные индексы, к примеру первая aaß (yY, t) и вторая baß (yY, t) квадратичные формы поверхности, в качестве аксиомы примем
~öT " ~дГ'
И последнее. Могут быть тензоры, имеющие пространственные и поверхностные индексы, к примеру, касательные векторы ^а- Обобщением (1.4, 1.5) будет для них формула
gA'X'-'s ax.'.am
jl-jr ßl-ßn
gt
dA'l-.'s ai ■■•am
jl-jr ßl-ßn
dt
is lii—lai".am
+
{Г1
il Ali2... is al...am
Ik jl".jr ßi-ßn
+
+r„A
- Г,. A'lr's:
il...is al...am
lk jl".jr ßl -ßn jlk lj2'''jr ßl-ß.
-----Г' A'l""
il...is al...am
jrk jl -l ßl-ßn
Gvk
(1.6)
n
Перечислим некоторые важные свойства операции, определенной в (1.6), не останавливаясь на их доказательствах, найти которые можно в [15].
Во-первых, 6-производная является тензором того же порядка и типа, что и исходный. Во-вторых, для нее выполняется правило Лейбница дифференцирования произведения. В-третьих, в рассматриваемом случае евклидова пространства вычисление 6-производной перестановочно с операцией тензорного дифференцирования по уа [15]. В-четвертых, вычисление 6-производной перестановочно по отношению к операции свертывания по латинским или греческим индексам.
Для геометрических характеристик Е на основании (1.6) получим
= (Оч%,
5 1дх1
"б?
6а,
дув
ав
6г
= - 2ОЪ,
ав'
6г
~бГ
дув
= 2ОЬав,
6Ь,
ав
(1.7)
6г
~5Г
§с«р
~бГ
6г
= 3Осав + аа ав°О
уОг
= —Осав + О,ав'
(ЪХвО,ух + 2ОЬав
суа) + аа (ЪхаО:ОХ + 2ОЪуас1а) ■
= аа^
,аа + ЪуаО,ва) '
В формулах (1.7) сав — компоненты третьей квадратичной формы Е.
Рассмотрим теперь функцию /(х\ ^, которая соответствует компонентам некоторого тензорного поля, непрерывного и непрерывно дифференцируемого требуемое число раз в рассматриваемой области Б3, кроме поверхности Е. На Е функция / или ее производные определенного порядка могут иметь разрывы первого рода. Область, в сторону которой движется Е, и функции в этой области в дальнейшем обозначаются знаком " +", остальную часть Б3 и функции в ней обозначим знаком " —". На Е / не определена, но определены ее предельные значения /+ и /—, взятые по разные стороны от поверхности Е. Тогда, вычисляя ковариантную производную и частную производную /±, получим
г* = + а^х1 г-х-
ду 1 ¿± 6/±
V '3 ',а'
/± =
6г
—/ру]'
(1.8)
ду ']
.к х] = Я;кх ■
Здесь и далее индексом после запятой обозначено ковариантное дифференцирование, если не оговорено иное. Из (1.8) можно получить величины, определенные на Е, — разрывы функции / и ее производных:
ил =
д/
ду
ы
[/] = /+ — Г
V] + аав [/],в х,а,
и
д/
ду
О,
(1.9)
Соотношения (1.9) носят название геометрических и кинематических условий совместности 1-го порядка.
Для исследования поведения решений в динамических задачах в областях, примыкающих к поверхностям, где имеют разрыв искомая функция или ее производные, применяются различные приближенные методы, в частности, лучевой метод [3, 10]. В нем необходимо знать рекуррентные соотношения на разрывы производных от поля / произвольного порядка. Такие соотношения были получены в [13] для случая декартовой системы координат. В [15] было показано, что в криволинейной пространственной системе координат их обобщением будут соотношения:
Ш = ] V+^{[¿'Д. -
- 2 2 ей -1)! [ *г2"'*Т)+
¡=3 р=0
+Ьав^ 2 2 С-<5 - 1)' [№I <2-ЧР ¡=2 р=0
[ЯЦ = - с [Щ + № - -
- 2 2 С'-2-1<5 - 1)' [¿Т"],а кТ2-ркр+1)+
5=3 р=0
+<к - 1)ЬавО,в 2 С'-2-1<5 - 1)! /Т4]а
5=2 р=0
А = ] ...vffl
(1.10)
причем в формулах для и V свертка с компонентами V^... V™ прово-
дится <к-1)-раз, — обозначена к-я ковариантная производная. В (1.10)
К и К2 — обозначения, принятые для главных кривизн Е. Отметим, что область применения полученных формул не ограничивается тематикой настоящей статьи. Они имеют универсальный характер в евклидовом пространстве и могут применяться при решении динамических задач со слабыми волнами или же для задач, включающих стационарные поверхности разрывов и т.д.
2. Модельные соотношения упругой среды с малой нелинейностью
Движение нелинейной изотропной упругой среды в криволинейной пространственной системе координат Эйлера х1 задано уравнениями:
V1 = и' + и1]У],
aij = р (V + j),
1 I k \ i P dW /„ ; . ; \
ij = ■= {Ui,j + Hi - uk4u)), a) = — — (6< - 2a\),
a i P dw
a-ij = т IЩ,] + Ujti - ukju ■), a -=---
2V j j po dak
/ 4 8 \1/2
P = Poll - 2/i + 2/2 - 2/2 - -/3 + 4hh ~ з/зj , (2-1)
/1 = ai, I2 = ajaj, /3 = а)а[ ak,
,•/• - dul - ^ + F' и - дщ г* и
u = -, u • = -: + 1 u , u; ; = -: — 1 • -ui.
dt •■> dxJ Jk 4 dxJ l]
Здесь ul и vl — компоненты векторов перемещения и скорости; a,j — компоненты тензора деформаций Альманси; aj — смешанные компоненты тензора напряжений Эйлера-Коши; р и po — плотность среды в текущем и свободном состоянии; W — функция упругого потенциала среды. Система (2.1) определяет движение среды всюду, где входящие в нее функции непрерывны и дифференцируемы требуемое число раз. Этой системе соответствует адиабатическое приближение, означающее постоянство энтропии в области, где выполняется (2.1). Для конкретизации типа нелинейно-упругой среды зададим функцию W ее разложением в ряд Тейлора в окрестности свободного состояния. Если считать среду изотропной, то это разложение можно провести по инвариантам тензора деформаций:
X
w=-lj + \!h + Ihh + ml\ + nh + III + цф2 + xhh + -x/t + ■ ■ ■ ■ (2-2)
В (2.2) X, ц — упругие модули второго порядка; l, m, n — модули третьего порядка; п, к, х — модули четветого порядка; многоточием обозначены слагаемые более высоких порядков малости по деформациям. Ограничимся в дальнейшем представлением (2.2), тем самым считая, что среда имеет малую нелинейность. Объединение (2.1) и (2.2) приводит к замкнутой системе уравнений нелинейно-упругой изотропной среды. Для некоторых из рассматриваемых в дальнейшем задач, в которых нас будут интересовать прежде всего сдвиговые деформационные процессы, принимается модель несжимаемой изотропной нелинейно-упругой среды (р = po). Для несжимаемой среды W = W (Ii, /2), а в представлении степенным рядом
W = (a - ц)/1 + a/2 + bl2 - к/1 /2 - 6/3 + c/f + d/f + k/2/2 + ••• , (2.3)
в которой учтено, что в этом случае /1 ^ 0, /2 ^ 0, так что все упругие постоянные в (2.3) ц, a, b, к, 6, c, d, k положительны.
В основе построения модели упругого тела лежат законы сохранения массы, импульса и энергии, следствиями которых в дифференциальной форме являются уравнения (2.1). Предположение о разрывах функций, входящих в (2.1), имеющих место на движущихся в пространстве поверхностях, позволяет определить понятие обобщенного решения. Для такого типа решения малая по толщине область быстрого изменения характеристик задачи заменяется на поверхность разрывов. В дальнейшем будем рассматривать поверхности разрывов градиента перемещений — ударные или сильные волны. На ударных волнах не выполняются уравнения (2.1), но справедливы динамические условия совместности, следующие из интегральных законов сохранения:
[Р(- а)] = о,
[с'] V у = р+ ' у - ,
' [V;] Vу = р+ 'у - в)
[V,-] [V]
(2.4)
2
- [е]
к!
V у.
Первое уравнение в (2.4) — следствие закона сохранения массы, второе — закона сохранения импульса, третье — закона сохранения энергии, причем в нем е — удельная плотность внутренней энергии, к — компоненты вектора теплового потока. Кроме динамических условий совместности на Е разрывы функций связаны геометрическими и кинематическими условиями совместности (1.9), а разрывы их производных к-го порядка условиями (1.10).
3. Возможные скорости и типы одномерных цилиндрических ударных волн
Рассмотрим в цилиндрической системе координат х1 = г, х2 = ф, х3 = = г случай, когда физические проекции вектора перемещений иг, Иф, щ зависят только от г, t, причем это выполнено как для предварительных деформаций, так и для послеударного процесса. В соответствии с типом задачи вектор нормали возможных ударных волн имеет компоненты vг = = ±1, Vф = vz = 0. В дальнейшем рассмотрим случай расходящихся волн, т.е. vг = 1. В качестве основной характеристики разрыва на Е выберем вектор
[и У] V' = XV ' + уц',
[и ' V ■'Vц'ц,- = 1, V 'ц, = 0,
(3.1)
т =
который выделяет нормальную составляющую разрыва т и ортогональное дополнение уц', вычисляемое проектированием на касательную к Е плос-
д х'
кость. Проекции уи,' на касательные векторы —- обозначим та, тогда для
дуа
нашей задачи
т = [иг,г] , т1 =
ди
дг
7 К']'
т2 =
ди3
дг
= [иг,г] ,
диг диф диг
Ыг г ~ , ^ ? " \у
дг дг дг
В дальнейшем будем работать с более удобными физическими проекциями:
тг = т, тф = гт1, тг = т2. (3.2)
Вычислив согласно (2.1) тензоры деформаций а,у и напряжений с,у и подставляя их наряду с (3.1) и (3.2) в динамическое условие совместности, следующее из закона сохранения импульса (второе из уравнений (2.4)), приходим к следующей системе уравнений:
Fтг + Ктф + Ьт1 + Мтгтф + Итгт1 = 0,
Ртг + Qтф + Ятг + S тг тф + Ттфтг = 0, (3.3)
Утг + Wтф + Хтг + 7тгтг + Zтфтz = 0.
Коэффициенты, входящие в (3.3), являются функциями предварительных деформаций иг/г, иф/г, иг/г, иг,г, иф,г, иг,г, §1 = иг/С1, §2 = иф/С1, §з = иг/С2, #4 = иф/С2, компонент волнового вектора разрывов тг, тф, тг, скорости волны С, а также упругих модулей среды; С\ = -у/Л + 2(х/ро и Сг = л/Ц/ро соответствуют скоростям продольной и поперечной волн в линейном приближении.
Для большей наглядности формул знак "+" у всех величин, соответствующих деформациям перед Е, опущен. При записи (3.3) предполагалось, что
иг иф «V
каждая из величин ипг, а -¡г. -у- является малой как до 2., так и
после, а также то, что тг, тф, тг — малые величины одного порядка малости с деформациями. Из (3.3) получаем три возможные скорости ударных волн. Предполагая в ~ С1 +..., где многоточием обозначены все невыписан-ные нелинейные слагаемые, из последних двух уравнений в (3.3) получим следующие соотношения для компонент разрыва:
( иф N
Тф ~ I Ах — + АгМф,г Iтг, тг ~ А2и^гхг,
Зц-/-§И (3-4)
А\ = \ + Ао, А2 = -2-А0, А О = —Л-
Л + ц
откуда следует, что на такой волне поперечные составляющие разрыва имеют второй порядок малости по отношению к продольной. Из (3.4) следует, что на этой волне поперечные составляющие изменяются неравноправно:
если приближенно считать аГ1 ~ и агф ~ — ^Иф,г--— то из (3.4) полу-
иф
чим Тф ~ 2А2а«ьТт--~ 2А2аггтг- Определим плоскость поляризации
г
волны положением векторов V и ц. Для ц из (3.4) следует:
Цф Тф 7* 7*
А2и7,
иф иф
(3.5)
г
Если перед ударной волной присутствуют только скручивающие деформации (ы1,г = 0, Ыф,г Ф или же только антиплоская деформация
Ы1гГ Ф 0, иф,г = иф = 0^, то в этих случаях из (3.5) следует неизменность типа деформации за волной. Такое же заключение невозможно сделать в общем случае произвольных предварительных деформаций: для них может поменяться направленность предварительного сдвига, причем изменение это зависит от текущего положения волны (от г) и от констант материала, что отличает рассматриваемый случай от плоских ударных волн, где вектор ц полностью определялся только предварительными сдвиговыми деформациями [10]. В соответствии с принятой в литературе терминологией, назовем эту волну квазипродольной. С учетом (3.4) из (3.3) для ее скорости получим:
С ~ Сх \ 1 + а\Хг + а2игг + + а^--1- а$х +
г
Ыг 2
+(ц,иг<гхг + а-1$\Хг + «8—хг + адигг + «шИг,п?1 + г,
Ыг 2 Ыф 2 (3.6)
+«ц—+ а\2Чфг + «1з—+ а\4и1Г+
2 Ыг Ыф / иг \2 / Ыф\2
+«15,? 1 + «16—+ «17-82 + «18 — + «19 -
г г г г
Коэффициенты щ здесь являются функциями упругих модулей. Скорость О1 в линейном приближении зависит от поля предварительных объемных деформаций и от тг. Зависимость О1 от сдвиговых деформаций и их разрывов является эффектом более высокого порядка малости. Если предварительные сдвиговые деформации отсутствуют, то на основании (3.4) волна становится чисто продольной.
Теперь рассмотрим волны, для которых продольная составляющая разрыва мала относительно поперечных. Эти волны в дальнейшем будем называть квазипоперечными. Считая для них О ~ С2 + ..., из (3.3) можно получить две различные формулы для скорости. Обозначим О = О11 или О = О111 в соответствии с номером уравнения в (3.3), откуда была вычислена скорость. Из первого уравнения (3.3) для тг получим:
т, = - АоуТф + • • • , ¥ = [и\ г + му . (3.7)
Из (3.7) следует, что для квазипоперечных волн для цилиндрических задач условие ¥ = 0 не влечет за собой тг = 0, что выполнялось для плоских одномерных волн [10]. Отметим, что в (3.4), (3.7) везде присутствует константа А0, определяемая упругими модулями X, ц, I, п, которая в рассматриваемом приближении отвечает за взаимодействие различных процессов
(3.8)
деформации. С учетом (3.7) из остальных уравнений (3.3) получим: 11 I Ыг 2
С « с211 + Ьщ^г + Ъг%ъ + Ъъ— + ¿4Тф+
Ыф 2 2
Н + + Ь%и1<гх1 + 69и,.г + Ь10иг^3 +
Ыг 2 Ыф 2
+Ьп — иг<г + ¿12"ф,г + ¿13 —"ф,г + ¿14"г,г +
Ыф , (Ыг \2 (Ыф \2
+¿15—£4+¿16 — +¿17 — +• г V г /
Формула для О111 аналогична (3.8), но следует помнить, что входящие в них функции различны. Коэффициенты в формуле для О111 обозначим dí•. Как и прежде Ь, di являются алгебраической суммой упругих модулей.
На каждой из волн Е11, Е111, которым соответствуют найденные скорости О11, О111, кроме (3.7) есть еще одно условие, связывающее между собой тг, Тф, тг. Эти условия определяются из (3.3) после исключения из этой системы О11 или О111. Так, для значения скорости О11 выполняется:
( Ыф ] ( Ыг (Ыг\2
ТфМг,г|/г1Тф + й2Мф,г + й3— |+тг|/г4—■ + й5 1 — 1 +
Ыг (Ыф \2 Ыф 2 2 ("Х ОЛ
+Й6 — ит,т + ¿7 I — I + ¿8~"ф,г + ¿9"ф,г + ¿10"г,г+
+ЙцМг,гТг} + ТфТг |й12Мф,г + + • •' = 0.
В (3.9) = Ь (к, ц, I, т, п, к, х, п, , г = 1.....13. За исключением частно-
3П
го случая, когда Й4 = 4 - — обращается в ноль, следствием (3.9) будет
невозможность на Е11 совместно присутствия Тф и тг одного порядка малости. Поскольку экспериментальные данные дают отрицательные значения модулей I, т, п, то можно считать, что это выполняется для любых нелинейно-упругих сжимаемых материалов. Для них наличие предварительных объемных деформаций, созданных за волной Е1, приводит к условию
Мф
и + ' ' ' • (З.Ю)
Й1Тф + /г2Мф,г + —
/г4-т + •••
Из (3.10) можно сделать вывод, что баланс между Тф и тг на Е11 также определяется не только предварительными деформациями, но и видом среды (ее модулями) и расстоянием до нагружаемой поверхности. Также из полученного условия следует, что на Е11 основной составляющей разрыва оказывается Тф, а тг и Тг имеют относительно нее второй порядок малости,
Ыф Ыг
если считать, что щ1Г, иг<т, —, — — величины одного порядка малости.
По аналогии с (3.9) для скорости О также есть связь между Тф и тг:
Ыг\2 Ыг Ыф / Ыф\2
Н\е\— + — + еъ — иг<г + е4 — 1иг + е5\ —
1 г \ г / г г \ г /
2 2 Ыф
+ + е1и1Г + Тф + едИф.ДфГ +
( Ыф Ыф ^ (3.11)
+Т£ + +£цИф,,-Тг + е\2-и .г + £13—тЛ +
+^14 ТфТгЫг,г + ••• = 0,
в1 = е(- (Л, ц, I, т, п, к, х, п, , г = 1,..., 14
причем в1 = -Й4, что определяет на Ъ111 ситуацию, обратную рассмотренной выше:
Ыф , Ыф
Тф ~--—-тг + • • • . (3.12)
е\ г + • • •
Т.е. на Ъ111 основная составляющая разрыва — тг, а остальные — малые более высокого порядка.
До сих пор был оставлен открытым вопрос о том, какая из квазипоперечных волн движется быстрее. С целью ответа на него рассмотрим два наиболее типичных случая предварительных деформаций в среде. Предположим, во-первых, что до момента ударного воздействия в среде есть объемные и антиплоские деформации, т.е. до прохождения Ъ1 ы+ ф 0, ы+ ф 0, Ыф = 0. Из (3.5) в этом случае получим, что за Ъ1 поле деформаций не меняет своего типа. Провести точный сравнительный анализ О11 и О111 невозможно, поскольку предварительные деформации, входящие в формулы для них, в (3.8), вообще говоря, различные. Тем не менее, известно, что на каждой из этих волн Тг меняется очень мало. Для послеударных времен, когда Ъ11 не отходит далеко от Ъ111, это позволяет приближенно считать ы+|ъи ~ ы+|ъШ и ы+г|ъ77 ~ ы+г|ъи7. Тогда, отбрасывая в формулах скоростей квадратичные слагаемые, получим
Си-Сш «(й3-<*з)-, (3.13)
г
откуда можно считать О11 < О111, если только иг > 0 (т.е. для ударных волн сжатия). Полученная оценка будет справедлива не только для рассматриваемого частного случая предварительных деформаций. В соответствии с ней квазипоперечная волна, на которой в основном меняется антиплоская деформация, движется быстрее, чем волна, которая изменяет прежде всего кручение среды. Пусть теперь для нашего случая за Ъ1 идет Ъ111; учитывая, что ыф = 0, из (3.12) получим
Тф!Ъш = 0, тгУп ф 0, Тг\ъш ^ тг. (3.14)
Следом за Ъ111 движется Ъ11, на которой выполняется условие (3.10).
Ы
Остановимся на вопросе о возможности существования нейтральных волн для нашей задачи. Напомним [10], что в теории плоских одномерных волн нейтральной названа волна, на которой Т = 0 и предварительный сдвиг меняет свое направление, не меняя при этом интенсивности. Для таких волн скорость распространения зависит только от предварительных деформаций среды. Нейтральной волне предшествует волна, изменяющая интенсивность предварительного сдвига, не меняя направления. В нашей задаче Ъ111 действительно меняет только интенсивность сдвига. Проверим, может ли волна Ъ11 быть нейтральной в смысле приведенного выше определения. Если это так, то кроме (3.10) необходимо считать ¥ = |ыфг + ы2г| =
= 0. Из условия ы+ = ы+ = 0 и из (3.7) это эквивалентно требованию
•> ф, г ъ11 ф ъ"
тг = 0. Но совместное выполнение (3.10) и ¥ = 0 приводит к тому, что Тф|ъ„ = тг\ъп = 0, что невозможно. Таким образом, для приведенной задачи невозможно появление нейтральной волны.
Рассмотрим еще один вариант предварительного деформирования, считая теперь ы+ ф 0, ыф ф 0, ы+ = 0. Здесь на Ъ1 остается тот же вариант деформаций, что и перед ней. Но теперь у нас нет квазипоперечной волны, меняющей только направление сдвига: действительно, за Ъ1 идет Ъ111, на которой из (3.12):
Ыф
--—г-+ (3-15)
+ •••
Поэтому за Ъ111 появляются все виды деформаций. Тогда на следующей волне Ъ11 имеем условие (3.10). Ни одна из этих волн не может быть нейтральной, поскольку в нашем случае условие тг = 0 означает, что ¥ = Тф = 0 и тем самым тг = 0. Если же ослабить требования и считать, что на нейтральной волне наряду с поворотом предварительного сдвига может присутствовать слабое изменение продольного скачка тг, то все равно получаем, рассматривая остальные требования, очень жесткие условия на предварительные деформации, скачки и модули среды, связанные, кроме того, с текущим положением волны. Эти условия вряд ли можно выполнить при наличии произвольного ударного воздействия. Проведенное исследование позволяет говорить о том, что для одномерных цилиндрических ударных волн понятие нейтральной волны (ударной волны круговой поляризации) не должно иметь место и в общем случае.
4. Лучевое разложение решений одномерных цилиндрических задач с ударными волнами
Лучевым методом называют способ построения краевых задач динамики сплошных сред путем разложения искомого решения в ряд по типу ряда Тейлора за подвижной поверхностью разрывов. Такое разложение может
проводиться по временной или лучевой координате. Последняя представляет собой расстояние, отсчитываемое от поверхности разрывов до данной точки пространства по соответствующему этой точке лучу. В случае слабых волн можно получить и последовательно интегрировать обыкновенные дифференциальные уравнения на коэффициенты лучевого ряда. Для задач с ударными волнами прямое интегрирование невозможно, но модификация метода, предложенная А.А.Бурениным [9, 10], позволяет и в этом случае получить решение, когда послеударное время мало. Таким образом были получены решения ряда одномерных краевых задач с плоскими ударными волнами. Однако такой метод решения можно с успехом перенести на задачи с криволинейными граничными поверхностями, в том числе многомерные задачи, когда необходимо определять геометрию лучевых координат [11]. В настоящей статье в качестве иллюстрации применения лучевого метода и полученных ранее условий совместности разрывов рассмотрим более простую одномерную задачу о нормальном ударе по внутренней границе цилиндрической полости в пространстве, заполненном сплошной средой.
В результате нормального нагружения на внутренней поверхности бесконечного кругового цилиндра радиуса го его граница L(t) начинает двигаться по закону
иг L(t) = g(t) , t ^ 0 ,
» 1 (4.1)
g(0 = > , -rgnt", gi > о, L{t) = r0 + g(t). ¿—t n!
n=1
В (4.1) условие gi > 0 означает, что воздействие на L(t) сразу создает расходящуюся продольную ударную волну. Функция g(t), определяющая закон движения, представлена рядом Маклорена. Количество членов этого ряда можно ограничить в зависимости от требуемой точности решения. Например, если в искомом решении приемлемо оставить слагаемые до второго порядка по времени, то для g(t) можно принять квадратичную зависимость
at2
g(t)=v0t+—. (4.2)
С целью упрощения дальнейших выкладок предположим отсутствие в среде предварительных деформаций. Тогда на фронте ударной волны Ъ(г) получим два краевых условия:
иг = 0
t
. |ад , I(t) = го + f
тг |вд = -и
Искомое решение для поля перемещений представим рядом
ы(г, г) = -
то 1
п=1
(г - гЕ)п,
^е(г) =
г г
(Щ,
щу
г=Е
Кп =
дпЫ
дгп
(4.3)
г=Е
В нашем случае лучевая координата совпадает с радиальной г, под ы(г, г) будем понимать единственную ненулевую компоненту Ыг вектора перемещений, относительно нее предполагаем сколь угодно высокую гладкость в окрестности Ь(г).
Из геометрических и кинематических условий совместности (1.9) в рассматриваемом случае следует
[ЫМ*)\ = -О
- _г~3—(— -
|".г V I = С 1 - Яц+Х
б2
бг \ бг
Кг+Л + О
2
Кг
б г2
-2
Ы
+ Кг+2
(4.4)
Закон сохранения импульса на поверхности разрывов (второе из соотношений 2.4) сводится к единственному уравнению, из которого скорость О продольной ударной волны можно выразить через ее интенсивность в виде
' К1
О = С1
к=1
р1 =-7 +
3(1 + т + п) 2(Х + 2[л) '
(4.5)
1 1 4 2(Х + 2ц)
3(1 + т + п) 7 2(Х + 2ц) 4'
На Е(г) внутренние координаты удобно выбрать так, что у1 = ф, у2 = г. Уравнение движения выпишем до квадратичных слагаемых включительно:
Ы Ы,г Ы ЫЫ,г Ы2 Ы,г
и гг 1 + а\и г + а2— н—:---Н--7т + аз— + 014 — +
г г г2 г2 г3 г
1--+ 2мм г +
а! = -7 + 6-
,1 + т + п к + 2ц
а2 =
С1
-4к - 2ц + 2(1 + 3т)
(4.6)
5к + 7ц - 41 - 6т - 3п
аз =-г—^-, а4 =
к + 2ц
к + 2ц
-6к - 9ц + 41 + 6т + 3п к + 2ц
Предположение о гладкости функции ы(г, г) в исследуемой области позволяет продифференцировать уравнение (4.6) до произвольного к-го порядка и записать полученные уравнения, начиная с к = 0, в разрывах с учетом
2
условий совместности (4.4). В [10] показано, что в результате получается бесконечная система нелинейных дифференциальных уравнений, связывающая неизвестные функции кп с их 6-производными. В отличие от случая слабых волн, в задачах с ударными волнами в каждом из уравнений к-го шага присутствует неизвестная функция к + 1-го шага, поэтому непосредственное интегрирование становится невозможным. Однако представив неизвестные функции Кп рядами по 6-производным в окрестности момента времени г = 0 и включив эти разложения в (4.3), коэффициенты лучевого ряда можно определить из граничных условий с требуемой степенью точности [10]
^ 1 ЪкКпо к §кХп0 ОкК
о = 2^—,
, „ к! Огк ' Огк Огк
к=0
, Кпо = Кп(0). (4.7)
г=0
Подстановка (4.7) в (4.5) позволяет представить скорость волны О в виде степенного ряда по времени
Ог
2^
1 О"К|0 /ОК|0\
+
г2 +..
(4.8)
Г « р1С-1К10 + р2С-2к20, Г1 « р1 + 2р2С-1К10, Г2 « р2.
Последнее представление скорости О позволяет определить время, через которое ударная волна будет занимать положение г = г (г^)
^ = С-1 (1 + р1С-1К10 + |32С-2х20 + ...)-1 (г - Г0) -
X (1 + р1С-1К10 + р2С-2х20 + ...)-3 (г - Г0)2 ) + ....
Входящие в (4.8)—(4.9) неизвестные величины должны быть определены в ходе решения задачи. На первом шаге метода, записывая уравнение движения (4.6) в разрывах, получим:
ы ~ сх ' ' ' 1 ' ;
Г0 + / О(1Ш Г0 + / О&Ш
00
где А, В, D — безразмерные коэффициенты, сложным образом зависящие от упругих модулей среды и от искомой функции К1. Уравнения типа (4.10) обычно называют уравнениями затухания интенсивности разрыва. В отличие от плоской задачи в (4.10) входит величина г-1, определяющая изменение кривизны волнового фронта. В пределах квадратичного по времени
анализа необходимо определить только величины хю, 10 и Я2о- Подставок
ляя (4.3), (4.7)-(4.9) в граничное условие (4.2), получим
| у0 1 аг2 1 бхю Г^г2 |
——— + --
С1(1 +Г) 2 С1(1 +Г) 2 6г [С1 (1 +Г)]
бкш уо (л уо 1 2
3
,1--^-и2- (4.11)
Ы С1(1+Г)\ С 1(1+Г)/ у 7
к2о Л уо \2 д , аг2
— 1--г = у0г + —.
2 \ Сх (1 + Г)/ 2
В первом приближении отсюда следует
К20 « -уо. (4.12)
Положив в (4.10) г = 0 получаем алгебраическое соотношение, связывающее
хю, Х2о, ——• Это позволяет исключить из (4.11) одну из неизвестных ог
бкш „
величин, например —— и наптп выражение для кю- £> силу громоздкости ог
. бкш
последнего здесь его приводить не будем, но —-— должно быть определено
ог
из (4.10) с учетом найденного Х20.
Приведем решение для поля перемещений до квадратичных слагаемых включительно
И = Vо 11 + (АоС^хго - —vo + — | ь + ... I (г - г-£) -
I \ Г0 г°/ ) (4.13)
12 ~2 ("20 + ...)(?- ьг + ■■■ ■
Индекс "0" соответствует величинам, вычисленным при г = 0. Разложение
и(г, г) в ряд (4.3) можно продолжить и дальше в зависимости от принятого
граничного условия. Например, если во всех разложениях учесть еще сла-
3 г г г б2К10 6X20
гаемые с г, то необходимо было бы определить также хзо, -т— и —-—.
ог2 ог
Дифференцируя уравнение движения (4.6) и записывая его в разрывах,
6к2
можно получить связь между —— и хз, а дифференцируя (4.10), выра-
ог
6К2 62К1
зить —-— через ——. Рассматривая полученные зависимости в привязке к ог ог2
начальному моменту времени, из (4.11) можно определить кзо в зависимо-
„ 6К20 62кш
сти от #1, #2) 8з, а значит, будут найдены —— и -—.
ог ог2
Аналогичную задачу можно поставить для расходящихся сферических продольных волн. Для них в уравнение затухания (4.10) будет входить средняя кривизна поверхности, что приведет к более быстрому затуханию интенсивности, чем для цилиндрических волн. Для среды с внутренней геометрической связью, которую накладывает условие сохранения объема,
авторами рассматривались задачи антиплоского деформирования в нелинейно-упругой несжимаемой среде, а также задача о скручивающем воздействии на границе цилиндрической полости. В каждом из рассмотренных случаев основные этапы решения те же, что и в приведенном примере, а решение может быть выписано с любой степенью точности.
Полученные формулы для перемещений и волнового фронта необходимы для выделения поверхностей разрыва в численных расчетах [11]. Здесь вся область от нагружаемой границы до прифронтовой в качестве решения использует какой-либо конечно-разностный метод, а неизвестные константы, входящие в лучевое разложение, определяются исходя из значений поля перемещений в нескольких узлах сетки, близких к прифронтовой области. Уточнение констант позволяет определить и текущее положение волнового фронта. Такая реализация метода расчетов может быть альтернативой методам сквозного счета, особенно в случаях, когда необходимо знание о точном положении поверхности ударной волны, без размытия последнего.
Заключение
Приближенные решения многих краевых задач ударного деформирования в прифронтовой области строятся методом лучевых разложений. Основные результаты здесь получены для плоских волновых картин. Поэтому закономерно возникает задача обобщения полученных для плоских волн результатов на случай движения волн произвольной геометрии. В настоящей статье эта задача представлена на примере одномерных цилиндрических ударных волн. Применение лучевого метода для волн произвольной геометрии потребовало обобщения понятия дельта-производной тензорных полей, а также рекуррентных соотношений, связывающих разрывы производных произвольного порядка по времени и пространственной координате, на случай криволинейной пространственной системы координат. Кроме того, для таких волновых процессов определены возможные типы и скорости ударных волн. Показано, что на скорости распространения ударных волн влияют как интенсивность разрыва, так и характер предварительных деформаций в среде. По своему типу волны разделяются на квазипродольную волну и две квазипоперечных: на каждой из волн присутствует продольная и поперечная составляющая разрыва, причем одна из них является доминирующей, а вторая имеет относительно нее более высокий порядок малости.
Работа выполнена при финансовой поддержке программы "Государственная поддержка ведущих научных школ" (НШ-890.2003.1).
Литература
[1] Быковцев, Г.И. Особые линии и поверхности в пространственных течениях идеальных жесткопластических сред / Г.И. Быковцев, И.А.Власова // Мех. деформ. тв. т. (динамика сплошной среды). -Новосибирск, 1979. - Вып. 41. - С. 31-43.
[2] Шаталов, А.Г. Разрывные решения в связанной задаче термоупругости / А.Г. Шаталов // Механика деформ. сред. Куйбышевский ун-т. -1979. - Вып. 6. - С. 85-90.
[3] Rossikhin, Yu.A. Ray method for solving dynamic problems connected with propagation of wave surfaces of strong and weak discontinuities / Yu.A. Rossikhin, M.V. Shitikova // Appl. mech. rev. - 1995. - V. 48. -No. 1. - P. 1-39.
[4] Томас, Т. Пластическое течение и разрушение в твердых телах / Т.Томас. - М.: Мир, 1964. - 308 с.
[5] Rossikhin, Yu.A. On construction of uniformly fit ray decompositions for solving dynamical problems of linear viscoelasticity (Engl transl) / Yu.A. Rossikhin, M.V. Shitikova // Soviet Appl. Mech. - 1991. - V. 27. -No. 1. - P. 77-82.
[6] Бестужева, А.П. К исследованию нестационарных поверхностных волн в нелинейно-упругих средах / А.П. Бестужева, Г.И. Быковцев,
B.Н.Дурова // Прикл. механика. - 1981. - Т. 17. - No. 12. - С. 27-33.
[7] Заварзина, А.А. Об ударных волнах в деформированной упругой среде / А.А. Заварзина, Г.Ф.Филатов // Нелинейные волны деформации: матер. Международного симпозиума. Таллинн, - 1978. - Т. 2. -
C. 70-73.
[8] Гринфельд, М.А. Лучевой метод вычисления интенсивности волновых фронтов в нелинейно-упругом материале / М.А. Гринфельд // ПММ. -1978. - Т. 42. - Вып. 5. - С. 883-898
[9] Буренин, А.А. Лучевой метод решения одномерных задач нелинейной динамической теории упругости с плоскими поверхностями разрывов /
A.А.Буренин, Ю.А. Россихин // Прикладные задачи механики деформируемых сред. Владивосток: ДВО АН СССР, 1991.
[10] Буренин, А.А. Об одной возможности построения приближенных решений нестационарных задач динамики упругих сред при ударных воздействиях / А.А. Буренин // Дальневосточный мат. сборник. 1999. Вып. 8. C. 49-72.
[11] Буренин, А.А. Выделение поверхностей разрывов лучевым методом в задачах динамики упругих сред / А.А. Буренин, П.В. Зиновьев,
B.Е.Рагозина // Сб. докладов международной научной конференции "Фундаментальные и прикладные вопросы механики". - Хабаровск: ХГ-ТУ, 2003. - С. 64-66
[12] Буренин, А.А. К проблеме выделения поверхностей разрывов в численных методах динамики деформируемых сред / А.А. Буренин, П.В.Зиновьев // Проблемы механики. Сборник статей к 90-летию
A.Ю.Ишлинского. - М.: Физматлит, 2003. - С. 146-155.
[13] Быковцев, Г.И. Теория пластичности / Г.И. Быковцев, Д.Д. Ивлев. -Владивосток: Дальнаука, 1998. - 528 с.
[14] Гринфельд, М.А. Методы механики сплошных сред в теории фазовых превращений / М.А. Гринфельд. - М.: Наука, 1990.
[15] Герасименко, Е.А. Геометрические и кинематические ограничения на разрывы функций на движущихся поверхностях / Е.А. Герасименко,
B.Е.Рагозина // Дальневосточный мат. сборник. - 2004. - Т. 5. - №1. -
C. 100-109.
Поступила в редакцию 19/У/2006;
в окончательном варианте—19/V/2006.
RAY SERIES IN STUDY OF PROPAGATION OF NON-PLANAR SHOCK WAVES3
© 2006 E.A. Gerasimenko, V.Y. Ragozina4
In this paper a way of forming approximate solutions of shock deformation boundary problems is discussed. The analysis is based on the theory of geometric and kinematic discontinuities compatibility equations when the moving discontinuity surface is determined in a curvilinear coordinate system. We also consider nonlinear effects in regularities of propagating one-dimensional cylindrical shock waves depending on type of initial deformation. Approximate solution with cylindrical longitudinal diverging wave is also obtained.
Paper received 19/V/2006. Paper accepted 19/V/2006.
3Communicated by Dr. Sci. (Phys. & Math.), Prof. Y.N. Radayev.
4Gerasimenko Ekaterina Andreevna ([email protected]), Ragozina Victoriya Yevgenjevna ([email protected]), Institute of Automatics and Control Processes Far Eastern Branch Russian Academy of Science, Vladivostok, 690041, Russia.
