О косом ударе жестким телом, имеющим плоскую границу, по нелинейному упругому полупространству Текст научной статьи по специальности «Физика»

13. Дударев В.В. Об уточненной модели изгиб-ных колебаний предварительно напряженной балки // Современные проблемы механики сплошной среды:

Тр. XII Междунар. конф. 2008. Т. 2. С. 56-59. 14. Филиппов А.П. Колебания деформируемых систем. М.: Машиностроение, 1970. 734 с.

УДК 539.3

О КОСОМ УДАРЕ ЖЕСТКИМ ТЕЛОМ, ИМЕЮЩИМ ПЛОСКУЮ ГРАНИЦУ, ПО НЕЛИНЕЙНОМУ УПРУГОМУ ПОЛУПРОСТРАНСТВУ

О.В. Дудко, Д.А. Потянихин

Институт автоматики и процессов управления ДВО РАН, Владивосток,

лаборатория механики деформируемого твердого тела E-mail: dudko@iacp.dvo.ru, potyanikhin@mail.ru

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

Ключевые слова: теория упругости, ударные волны, контактное взаимодействие.

About Oblique Impact by Perfectly Rigid Body with Plane Boundary on the Nonlinear Elastic Half-Space

O.V. Dudko, D.A. Potyanikhin

Institute of Automation and Control Processes FEB RAS, Vladivostok,

Laboratory of a Deformable Solid Mechanics E-mail: dudko@iacp.dvo.ru, potyanikhin@mail.ru

In this paper the impact interaction of perfectly rigid body and nonlinear elastic solid, which have plane boundaries, are investigated. Suppose that the moving rigid body has constant velocity, resulting in self-similar formulation of the problem. Possible variants of wave combinations, arising from such interaction, are discussed. The existence condition for evolutionary shock waves and the thermodynamic discontinuities compatibility condition serve as criterions for choosing the wave pattern. The scheme for the solution of a self-similar boundary-value problem includes checking for the mentioned criterion during computations.

Key words: theory of elasticity, shock wave, shock interaction.

ВВЕДЕНИЕ

Плоские автомодельные задачи нелинейной динамической теории упругости были рассмотрены в работах [1-5]. Автомодельное динамическое деформирование упругопластических твердых тел достаточно подробно изучено в [6-8]. Постановка задачи, представленной в настоящей работе, близка к [1]. Но в отличие от статьи [1], где рассматривается ударное взаимодействие двух нелинейно-упругих тел, здесь одно из тел полагается абсолютно твердым.

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

1. ОПРЕДЕЛЯЮЩИЕ СООТНОШЕНИЯ

Система модельных уравнений динамического деформирования нелинейной упругой среды в прямоугольной декартовой системе координат в переменных Эйлера представляется соотношениями:

— и г + , "г — уг + 1] , 2аг] — иг] п ик п ик ] ,

Р Ш (О 2 ) (1'1}

О] — —т;-(Ок] - 2ак]), О]— Р"г,

Р о да к]

Ж = -/? + ^ + К/! 12 + */? + П/3 + . . . , (1.3)

— = 1 - 2/1 + 2/2 - 2/2 - 4/3 + 4/1 /2 - 8/3,

Ро V 1 1 2 3 1 1 2 3 3 (1.2)

/1 = акк, /2 = ац, «ы, /3 = ауг •

В соотношениях (1.1), (1.2) иг, г>г, — компоненты векторов перемещений, скоростей и ускорений точек среды; агу, сту — компоненты тензора деформаций Альманси и тензора напряжений Эйлера - Коши; ро, р — плотность материала в начальном и текущем состоянии соответственно; <5гу — символ Кронекера. Латинский индекс после запятой обозначает частное дифференцирование по пространственной координате (/,у = д//джу), точкой обозначена частная производная по времени / = д//д^. Здесь и далее принято правило суммирования по повторяющимся индексам; латинские индексы принимают значения 1, 2, 3.

Система уравнений (1.1)—(1.2) становится замкнутой, если задан упругий потенциал Ж. В адиабатическом приближении функция Ж для нелинейной упругой среды зависит только от деформаций агу и принята в виде

А

2

где А, д — параметры Ламе; к, х, п — упругие модули третьего порядка.

Решение системы (1.1)—(1.3) может не являться непрерывным на некоторой движущейся поверхности £(£). В этом случае на £(£) необходимо потребовать выполнения условий совместности разрывов. Если при непрерывных перемещениях иг на £(£) терпят разрыв их первые производные, то в среде распространяется со скоростью С поверхность сильных разрывов — ударная волна, на которой параметры напряженно-деформированного состояния и движения точек среды связаны динамическими и кинематическими условиями совместности первого порядка [9]:

[стгз ] П3 = Р+ (л+ Ну - С) [V] , [иг] Ну = -С [и*, у ] Ну при [иг] = 0. (1.4)

Если же на £(£) непрерывны не только перемещения, но и их первые производные, а скачком изменяются вторые производные, то на такой поверхности необходимо потребовать выполнения условий совместности второго порядка [9]:

[сту ,у ] = Р+ И , [и] Ну = С2 [игу ] Нун^ при [иг] = 0, [иг,у ] = 0. (1.5)

В этом случае С — скорость распространения поверхности слабых разрывов.

В условиях (1.4), (1.5) квадратными скобками обозначен скачок функции на поверхности £(£): [т] = т+ - т-, т+ — значение функции перед поверхностью, т- — сразу за ней; Ну — компоненты единичной нормали к поверхности волны (нуНу = 1).

Анализ условий совместности (1.4) показывает, что в нелинейной упругой среде с потенциалом (1.3) в условиях плоского деформированного состояния возможно существование двух типов плоских ударных волн — квазипродольных, на которых преобладает изменение объемных деформаций, и квазипоперечных, вызывающих преимущественно сдвиговые деформации [10].

2. ПОСТАНОВКА ЗАДАЧИ

Рассмотрим косой удар абсолютно твердым телом, движущимся с постоянной скоростью, по покоящемуся нелинейно-упругому телу (рис. 1). Пусть упругое тело занимает полупространство х1 > 0. Его границу обозначим через . Второе тело с плоской границей Ь2, двигаясь как жесткое тело, сталкивается с первым так, что появляется общая граница двух тел ОЬ. Движение второго тела определяется постоянной скоростью его точек vo = {г>10 ,г>20}. На участке ОЬ выполняется закон сухого трения.

Рис. 1. Волновая картина: а - с двумя ударными волнами; б - с ударной и простой волнами

Пусть S — скорость движения точки O — начала подвижной системы координат — вдоль оси Ox2. Величина S выражается через заданные параметры задачи v10, v20 и угол столкновения ф соотношением S = v2o + vio • ctg ф. Будем считать эту скорость настолько большой, чтобы возникающие вследствие ударного столкновения волновые фронты не могли отделиться от точки O.

Ввод автомодельной переменной в = x1 /(St — x2) и представление компонент вектора перемещений в виде u1 = (St — x2)F(в), u2 = (St — x2)Ф(в), u3 = 0, где F(в), Ф(в) — первая и вторая компоненты безразмерного вектора перемещений в координатном пространстве в, позволяют перейти в системе определяющих соотношений (1.1)-(1.3) от дифференциальных уравнений в частных производных к однородной системе обыкновенных дифференциальных уравнений относительно неизвестных функций F(в) и Ф(в):

f A • F''(в)+ B • Ф''(в) = 0, C • F''(в) + D • Ф''(в) = 0,

где A, B, C, D — функции от в, F, Ф, F', Ф' и упругих модулей среды. Условие существования нетривиального решения системы (2.1)

(2.1)

AD — BC = 0

(2.2)

может выполняться либо при некотором значении 0 — 0* (что соответствует ударной волне), либо в целом интервале 0 Е [0 + ,0-] (0+ и 0- определяют положение двух слабых волн разрывов ускорений, которые являются передним и задним фронтами простой волны Римана). В областях тривиального

решения между волновыми фронтами ненулевые параметры движения точек среды и напряженно-деформированного состояния остаются постоянными и могут быть выражены через произвольные константы интегрирования а, Ь, е, / системы (2.1):

^(в) = ав - Ь,

Ф(в) = ев - /, -£Ь

VI =

=

и1,1 = а, и,2 = Ь,

-£(/ - а/ + Ье)

(1 - а)(1 - /) - еЬ'

(1 - а)(1 - /)- еЬ' а11 = а - ^ (а2 + е2) , а12 = 2 (Ь + е - аЬ - е/),

Ро

2

= (1 - а)(1 - /) - еЬ.

и2,1 = е, И2,2 = /, =0, и>2 = 0,

«22 = / - 1 (Ь2 + /2) ,

(2.3)

Компоненты тензора ст^- (а, Ь, е, /) здесь не приводим из-за их очевидной громоздкости.

Условия совместности разрывов (1.4), (1.5), записанные с использованием введенных автомодельных параметров процесса ударного деформирования, допускают существование в рамках поставленной задачи плоских ударных волн (квазипродольных и квазипоперечных) и простых волн (передними и задними фронтами которых являются плоские волны разрывов ускорений), а также областей постоянства параметров напряженно-деформированного состояния, каждой из которых соответствуют свои значения констант интегрирования а, Ь, е, /.

Таким образом, волновая картина, возникающая в области упругого тела, с математической точки зрения может состоять из различных комбинаций плоских ударных волн и простых волн Римана. Их характер определяется параметрами задачи г>10, г>20, ^ и коэффициентом трения к на границе ОЬ и может быть конкретизирован проверкой термодинамического условия совместности разрывов и условия эволюционности ударных волн.

Анализируя соотношения модели (1.1)—(1.3) совместно с условиями (1.4), (1.5), можно показать, что для нелинейной среды (1.3) скорость распространения объемных деформаций всегда оказывается больше, чем скорость распространения дополнительного сдвига. Значит, передним фронтом всегда оказывается продольная ударная волна Е1, положение которой определяется значением = tg в1 (см. рис. 1). Следом за ней может распространяться либо квазипоперечная ударная волна Е2, соответствующая значению в2 = tgв2 (см. рис. 1, а), либо центрированная простая волна в е [в+,в-] (см. рис. 1, б).

Решение краевой задачи в случае возникновения двух ударных волновых фронтов Е1 и Е2 заключается в определении соответствующих значений в1, в2 и констант интегрирования а, Ь, е, / в областях 1 и 2 (см. рис. 1, а) из системы алгебраических уравнений, составленных из условий совместности разрывов и непрерывности перемещений (1.4):

а1в1 - Ь1 = 0, е1в1 - /1 =0,

а102 - Ь1 = а202 - Ь2, е102 - /1 = е2 02 - /2,

стЦ + ст52} в1 + М2 в!^ =0,

(2.4)

ст21) + ст22} в1+м2 в^ = 0,

ст<;> - *1?) + V1) - ■ (2) ■ 12

'ст21) - ст<?>) + - ■ (2)

М 2= р£2 = ,

Л + 2д

5?)

■22)

£2

Л + 2^'

~(к) V

(к) г_

£

¿(к) = = (1 - ак)(1 - /к) - Ькек, Ро

(к)

С1 =

'Л +

Ро

1

2

1

2

2

ст

2

с

1

В соотношениях (2.4) и далее верхним индексом в скобках при О] и обозначена соответствующая область постоянства параметров напряженно-деформированного состояния.

Для замыкания системы (2.4) необходимо учесть отсутствие разрыва нормальной к полуплоскости контакта ОЬ составляющей скоростей точек ударяемого и ударяющего тел:

у(2) - УЮ) + (у22) - У2о) п^ — 0. (2.5)

Разрыв касательной составляющей этих же скоростей запрещаем лишь в случае выполнения нера-

■ ■ (ОЬ)(ОЬ) ■ ■

венства

(2) о • •

г]

(2) О • •

г]

(ОЬ) (ОЬ)

П п]

у(2) - Ухо) т1ОЬ) + (у22) - Т2(ОЬ) — 0. (2.6)

В противном случае условие (2.6) должно быть заменено выражением, следующим из закона трения Кулона:

П(ОЬ) п]ОЬ). (2.7)

(2) О • •

г]

п(ОЬ)т(ОЬ) — к

(2) О • •

г]

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

О+ [Уг] П] - Р+ (у+ П] - С) (^ + > 0, (2.8)

аналогом которого в газовой динамике является теорема Цемплена о термодинамической невозможности существования ударных волн расширения.

Проверку эволюционности ударной волны 02 проводим, определяя положение переднего фронта 0+ простой волны из условия (2.2) с подстановкой в него предварительных (предшествующих ударному фронту 02) постоянных параметров напряженно-деформированного состояния ах, , вх, Д, вычисленных из системы (2.4). Нарушение эволюционности ударной волны означает, что ее фронт 02 оказывается внутри веера простой волны: 02 < 0+. Такая ударная волна должна быть заменена на простую волну Римана 0 е [0+,0-] (см. рис. 1, б).

Решение в области простой волны строится путем интегрирования по автомодельному параметру 0 системы обыкновенных дифференциальных уравнений, состоящей из условия (2.2) и одного из равенств (2.1):

'л • ^''(0) + В • Ф''(0) —0,

У 7 У 7 ' (2.9)

А£ - ВС — 0.

Положение заднего фронта 0- простой волны заранее неизвестно, его приходится уточнять непосредственно во время счета, согласовывая получаемые значения с граничными условиями задачи.

3. РЕЗУЛЬТАТЫ ЧИСЛЕННЫХ ЭКСПЕРИМЕНТОВ

В процессе расчетов изменялись параметры задачи уХ0, у20, к, а также значения упругих постоянных А, д, к, х, п. Там, где это специально не оговаривается, механические свойства среды задаются безразмерными константами: А — 0.4, Д — 0.3, к — -10.0, х — -10.0, п — -10.0 (все упругие модули нормированы по (А + 2д)).

В каждом случае решение проводилось сначала в постановке с двумя ударными волнами (см. рис. 1, а), а затем уточнялось путем проверки условия (2.8) и эволюционности ударного фронта Е2.

Серия вычислительных экспериментов показала, что и при проскальзывании, и при сцепке возможно появление как ударной, так и простой волны.

На рис. 2, а представлены диаграммы распределения параметров напряженно-деформированного состояния, полученные при ^ — 1°, уХ0 — 0.3 сх, у20 — 0.2 сх, к — 0.3 и соответствующие скольжению

]

тел. При этом возмущения в среду распространяются посредством продольной ударной волны Е1 и волны Римана

в е [в+,в- ].

Изменение коэффициента сухого трения приводит к изменению характера контакта тел на границе ОЬ. Так, например, при ^ = 1°, г>10 = 0.3 с1, г>20 = 0.2 с1, к = 0.5 происходит жесткая сцепка взаимодействующих тел. Волновая картина при этом остается прежней. В обоих случаях продольная ударная волна вызывает сжатие среды, простая волна — расширение.

Выбирая другие параметры задачи, можно получить волновую картину с двумя ударными волнами. При ^ = 2°, г>10 = 0.2 с1, г>20 =0.1 с1, к = 0.2 на границе соударяющихся тел происходит проскальзывание, а при ^ = 3°, г>10 = 0.2 с1, г>20 = 0.1 с1, к = 0.2 — жесткая сцепка (рис. 2, б). Ударные волны Е1 и Е2 вызывают сжатие среды.

Изменение напряжений, плотности и компонент вектора скорости в зависимости от коэффициента трения к показано на рис. 3, 4. Для иллюстрации влияния физической нелинейности для каждого параметра приведено по три графика. Графики с индексом 1 построены при оговоренных значениях упругих модулей третьего порядка к, х, п; с индексом 2 — когда х, п уменьшены на порядок, и с индексом 3 — для случая геометрической нелинейности, когда К = х = П = 0. Исходными параметрами экспериментов выбраны ^ = 1°, г>10 = 0.1 с1, г>20 = 0.1 с1. Для каждого фиксированного набора параметров г>10, г>20, ^ и упругих модулей существует критическое значение коэффициента трения к*, в котором графики функций имеют излом. При к < к* происходит скольжение граничных плоскостей соударяющихся тел. Значение к = к* соответствует наступлению жесткого контакта. Для остальных параметров задачи г>10, г>20, ^ можно аналогично указать критические значения, при которых меняется характер взаимодействия тел. Механика

О Х2

Р/Ро

V /е1

V, /с.

1.1 1.05 1

0.3 0.2 0.1 0

0.15 0.1 0.05 0

ап/(Х+2|д) 0 -0.5 -1

а12/(Х+2|) 0 -0.1 -0.2 -0.3

а22/(Х+2|) 0 -0.2 -0.4 -0.6

Ь 2 5* Х1

-

-

N -

-

- У -

-

б

Рис. 2. Диаграмма распределения параметров напряженно-деформированного состояния и движения: а - постановка с ударной волной и волной Римана; б - постановка с двумя ударными волнами

а

0+2Ц)

-0.15

-0.2

-0.25

-0.3

12 -0.01

-0.02

-0.03

-0.04

-0.05

-0.06

0.1

-0.04 -0.06 -0.08 -0.1 -0.12 -0.14 -0.16

0.1

Р1"/Рп

э<1> 1.11

1.1 1.09 1.08 1.07 1.06 1.05 1.04

0.1

11111

-!-•—

3 2

1

0.2 0.3 0.4 0.5

0.2 0.3 0.4 0.5

0.2 0.3 0.4 0.5

—•—

......7

1

\ 1

2 2 3

.... / __—•-

11111

к

0+2Д)

-0.12 -0.14 -0.16 -0.18 -0.2 -0.22 -0.24 -0.26 -0.28

0.1

^/(Х+гц)

-0.05

-0.1

-0.15

-0.2

Р(2)/РП

1.09 1.08 1.07 1.06 1.05 1.04

0.1

0.2 0.3

0.4

/

2. у..

3

1

0.2 0.3 0.4 0.5

Рис. 3. Зависимость компонент тензора напряжений и относительного изменения плотности от коэффициента

трения при различных упругих модулях

к

к

к

к

к

1

2

3

0.12 0.118 0.116 0.114 0.112 0.11 0.108 0.106 0.104

2

1

0.06 0.055 0.05 0.045 0.04 0.035 0.03 0.025

%—

2

/ ----

0.35

0.4

0.45

0.5

0.55

0.35

0.4 0.45 0.5 0.55

С1

0.1008 0.1007 0.1006 0.1005 0.1004 0.1003 0.1002 0.1001 0.1

3.....

- -1 у

0.1 0.09 0.08 0.07 0.06

1

^ / 3

0.35

0.4

0.45

0.5

0.55

0.35

0.4

0.45

0.5

0.55

Рис. 4. Зависимость компонент вектора скорости от коэффициента трения к при различных упругих модулях

Выясним влияние параметров задачи на скорости возникающих волновых фронтов. Зафиксируем ^ = 1°, г>ю =0.3 сь г>2о = 0.2 С1 и начнем изменять коэффициент трения к. При малых значениях к возникают две ударные волны, тела проскальзывают (рис. 5, а). При увеличении к, на отрезке [к**, к*], вторым отраженным фронтом становится волна Римана, тела по-прежнему проскальзывают. Дальнейшее увеличение к приводит к жесткой сцепке с ударной и простой волнами.

Если теперь при неизменных скоростях г>10 = 0.3 с1, г>20 = 0.2 с1 увеличить угол соударения ^ до значения 3°, то и проскальзывание, и сцепка происходят с двумя ударными волнами (рис. 5, б).

4 3.5

3

2.5

2

(7-, / С1

^ г " 1......

........ о2/ И

и2с1

3.53 -2.52 -

&1/с1

О 0.05 0.1 £** 0.2 0.25

0.3

0.4

0.05 0.1 0.15 0.2 0.25

б

0.3 0.35

Рис. 5. Зависимости скоростей волн от коэффициента трения к

Отметим, что значения к**, при которых изменяется волновая картина, соответствуют точкам максимума графиков зависимости относительного изменения плотности р(2) /р0 от коэффициента трения (см. рис. 3).

Решение задачи удается получить только при небольших углах соударения. Углы в1 и в2 (или в+ и в—) быстро увеличиваются с увеличением и задача перестает быть автомодельной. Так, например, уже при ^ = 4°, г>10 = 0.3 с1, г>20 = 0.2 с1 угол в1 близок к п/2.

к

4

0

к

а

ЗАКЛЮЧЕНИЕ

Основной постановочной проблемой в автомодельных задачах ударного деформирования остается неоднозначность в определении возникающей волновой картины. Иногда удается удовлетворить начальным и граничным условиям задачи, предполагая появление в среде различных комбинаций ударных и простых волн [1-4], то есть построить более одного решения задачи. Выбор реализуемого варианта распространения деформаций связан с законами термодинамики и с условием эволюционно-сти сильных разрывов. Согласно принятому адиабатическому приближению для нелинейной упругой среды, энтропия оказывается кусочно-постоянной функцией пространственных координат и времени. Пренебрегая теплопроводностью среды при формировании модельных соотношений, считаем единственным необратимым процессом распространение ударной волны. Именно на поверхности разрывов деформаций энтропия может изменяться скачкообразно. Термодинамическое условие совместности разрывов (2.8) связывает производство энтропии с параметрами среды. Его следствием является аналог [10] теоремы Цемплена о термодинамической невозможности скачков разрежения. Более жесткие ограничения на существование поверхностей сильных разрывов накладывают условия эволюцион-ности разрывов [5]. В работе проверка критериев производилась в процессе численых расчетов на примере серии вычислительных экспериментов в задаче о косом ударе абсолютно твердым телом по упругому полупространству. Проведенные вычисления также, как и в [11], не выявили различия между двумя критериями. Два ограничения на существование поверхностей разрывов деформаций оказались равноправными.

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

Работа выполнена при финансовой поддержке РФФИ (проект 08-01-00001-а). Библиографический список

1. Агапов И.Е., Буренин А.А., Резунов А.В. О соударении двух нелинейно-упругих тел с плоскими границами // Прикладные задачи механики деформируемых сред: Сб. науч. тр. Владивосток: ДВО АН СССР, 1991. С. 206-215.

2. Буренин А.А., Лапыгин В.В., Чернышов А.Д. К решению плоских автомодельных задач нелинейной динамической теории упругости // Нелинейные волны деформаций: Материалы междунар. симпозиума. Таллин, 1978. Т. 2. С. 25-28.

3. Буренин А.А., Лапыгин В.В. Об отражении плоской продольной ударной волны постоянной интенсивности от плоской жесткой границы нелинейной упругой среды // ПМТФ. 1985. №5. С. 125-129.

4. Куликовский А.Г., Свешникова Е.И. Нелинейные волны в упругих средах. М.: Московский лицей, 1998. 412 с.

5. Чугайнова А.П. Автомодельная задача о действии бегущей нагрузки на границу нелинейно-упругого слабоанизотропного полупространства // ПММ. 1993. Т. 57, вып. 3. С. 102-109.

6. Баскаков В.А., Быковцев Г.И. Об отражении пло-скополяризованной волны от свободной поверхности в упрочняющейся упругопластической среде // ПММ. 1971. Т. 35, вып. 3, № 1. С. 71-72.

7. Быковцев Г.И., Колокольчиков А.В., Сыгуров А.В. Автомодельные решения уравнений динамики идеального упругопластического тела при условии пластичности Треска // ПМТФ. 1984. №6. С. 148-156.

8. Быковцев А.Г. О преломлении ударных волн чистого сдвига в упругопластическое полупространство // ПММ. 1989. Т. 53, вып. 2. С. 309-318.

9. Быковцев Г.И., Ивлев Д.Д. Теория пластичности. Владивосток: Дальнаука, 1998. 523 с.

10. Буренин А.А., Чернышов А.Д. Ударные волны в изотропном упругом пространстве // ПММ. 1978. Т. 42, вып. 4. С. 711-717.

11. Дудко О.В., Потянихин Д.А. Автомодельная задача нелинейной динамической теории упругости о взаимодействии продольной ударной волны с жесткой преградой // Вычисл. мех. сплош. сред. 2008. Т. 1, № 2. С. 27-37.