Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Механика сплошной среды. Т.4 Основы механики твёрдых сред

.pdf
Скачиваний:
25
Добавлен:
19.11.2023
Размер:
41.7 Mб
Скачать

§ 5.7. Плоские волны в пластических средах

561

решения, этот случай реализуется, когда обращается в нуль определитель

системы (5.7.32):

 

det ( ) = —p d X s + с1dt2dX = 0.

(5.7.33)

Из (5.7.33) находим уравнения двух семейств характеристик:

 

dX = ±adt,

(5.7.34)

где

(5.7.35)

a2(k) ее c\{k)/p

— скорость звука в пластической среде.

Подставляя (5.7.34) в первое и второе уравнения системы (5.7.28), полу­

чаем условия на характеристиках:

 

 

 

 

 

о

dt

. dT

77

 

dt

. dv

 

pdv =

dl —— =

± — ,

dk =

dv—— =

± — ,

 

или

dX

a

 

 

dX

a

 

pa dv =F dT = 0,

a dk =F dv = 0.

 

(5.7.36)

 

 

Введем функции

 

 

 

 

 

 

 

 

k

 

 

 

k

 

 

Л

л iJ

,„(

„ i J \

m J

(5.7.37)

= ci(k')dk',

cp(k)

=

a(fc')

dkf,

 

l

 

 

 

l

 

 

для которых d^ = ci dk, dip = adk, тогда уравнения (5.7.36) можно проинте­

грировать:

y?(fc) =F с = const.

(5.7.38)

£(fc) —Т = const,

Таким образом, имеем два семейства характеристик и по два условия на каждой из них:

dX — a(k) dt = 0,

 

< v — ip(k) =

рР+ = const,

(5.7.39)

Т — £{k) =

^

= const;

 

{ dX + a(fc) dt =

0,

 

c +

y?(fc) =

 

= const,

(5.7.40)

T -

f (fc) =

^

= const,

 

Выберем точку At на фазовой плоскости (X, t). Для нее всегда су­ ществуют две характеристики, причем для одной из них касательная яв­ ляется положительной (dt/dX = 1/a > 0), а для другой — отрицательной (dt/dX = - 1 /a < 0).

Рассмотрим такую точку At, для которой обе характеристики пересека­

ются с осью ОХ

(рис. 5.7.2) в некоторых точках Р и Q. Поскольку

точки

Р и Q принадлежат области задания начальных условий (5.7.30), то

в них

v = 0, k = I и Т

= 0. Тогда из условий (5.7.39) и (5.7.40) следует, что вдоль

§ 5.7. Плоские волны в пластических средах

 

563

образом, для рассматриваемой точки М вдоль характеристик М Р и M Q из

(5.7.39) и (5.7.40) имеем следующие соотношения:

 

 

МР:

v=<p(k) + <p°+,

Т = № ) + $_,

(5.7.44а)

MQ: v = -<p(k),

T = £(k).

(5.7.446)

Из уравнений (5.7.44а, б) находим соотношения в точке М:

 

 

 

v = (р+/2 = const,

= 0,

 

 

к = к0 = const,

Т = £(к°) = Т° = const.

(5.7.45)

Если выбрать произвольную точку М ' на «+»-характеристике

ЛЛР и

построить ее «—»-характеристику

M 'Q ', то для точки М ' получим

те же

соотношения (5.7.44а, б) и (5.7.45), причем с теми же константами

 

и fc°,

поскольку константа

одна и та же вдоль всей РЛЛ. Следовательно, значе­

ния функций v, к и Т, хотя уже и отличаются от значений в области покоя, но остаются постоянными вдоль всей «+»-характеристики. Это означает, что скорость звука а(к) = а(к°) = const также постоянна и характеристика М Р является прямой линией, уравнение которой

t

X

= С, С = const.

 

а(к°)

Поскольку к0 и Т° — значения функций k(X,t) и T(X,t) вдоль характе­ ристики М Р , то для различных значений X и t, удовлетворяющих уравнению t — (Х/а{к°)) = (7, выполняются соотношения

к0

 

X

= const,

 

= к(С) = k ( t -------Q-)

 

 

 

а(ки)'

 

 

г° = t(k°) = m e )) = е

X

= const.

(5.7.46)

 

 

 

а(/г)

 

Если применить соотношения (5.7.46) для точки Р, где X = 0 и кроме соотношений (5.7.45) имеет место граничное условие (5.7.29а), то из (5.7.46)

и (5.7.29а) получим

 

Р : m t ) ) = ~Pe(t).

(5.7.47)

Из этого соотношения находим функцию k(t)\

 

k(t) = c4~Pe(t)),

(5.7.48)

полагая, что существует обратная функция к = <^_1(Т°).

Подставляя (5.7.48) в (5.7.46) и (5.7.446), определяем значения к и Т для точки М с координатами X и t:

к° = к ( х , t) = щ - J) = Г 1(-Ре(* - J ) ) .

т° = Г(Х, t) = z(kit - J ) ) = - p e(i - J ) ,

(5.7.49)

564 Глава 5. Пластические среды

v(X,t) = -<p(k(X,t)) = -<p

В силу произвольности точки М , полученные соотношения (5.7.49) пред­ ставляют собой искомое решение в возмущенной области для таких значений t, при которых аргумент функций (5.7.49) неотрицателен, т. е. t > Х/а.

Для

моментов времени t

^

Х /а

возмущение не доходит до рассмат­

риваемой

точки М и вместо

(5.7.49)

получаем решение,

соответствующее

состоянию покоя:

 

 

 

 

 

к = 1,

Т =

0,

v = 0, t < Х/а.

(5.7.50)

Построенное решение (5.7.49), (5.7.50) имеет место только до тех значений времени t, при которых головная волна не доходит до тыльной поверхности пластины X = h®, т. е. при t < h°{/aQ.

5 .7 .3 . Сравнит ельный анализ реш ения для разны х м оделей А £

Проанализируем полученное решение (5.7.49). При активном нагружении,

(п) ^

когда dpe/dt ^ 0, имеем С\\ = 0. Тогда для t > Х /а из (5.7.49) следует, что da\\/dt = dT/dt ^ 0 и нагружение пластины происходит в области сжатия. Следовательно, соотношения (5.7.20), (5.7.21) и (5.7.25) в упругой области

принимают вид

 

 

_

_

h+ = 0,

h- = 0,

h+ =

h- = 0,

^ — 1

^

П. _1 97Щ 72(п—ш —1)

CQ = 1,

с\

= (1\ +

212)кх

 

,гг—I I I —1

(5.7.51)

 

если Т и — Т 22 > —л/3сг5.

Н

(п)

Подставляя в формулу (5.7.18) выражение (5.7.19) для Т 22, а вместо Т и его выражения через к, согласно (5.7.12), (5.7.38) и (5.7.49):

(п)

Т и

находим предельное значение к = ks < 1, при котором появляются пластиче­ ские свойства при сжатии:

2(n —ill; —1

(Д —//Д/i —LLi)Ks

l\

-f Zt2

(5.7.52)

 

Здесь, как и ранее, г/ = 1\/{2{1\ + ^)) — коэффициент Пуассона.

566

Глава 5. Пластические среды

не возрастает (для AjV остается постоянной, а для Ау убывает). Для всех моделей Ап функция а{к) в рассматриваемом случае имеет разрыв при к = = ks, поскольку было сделано допущение о линейном упрочнении среды (см. формулу (5.7.20)).

 

 

5 .7 .4 .

П лоски е волны в м оделях AjV и А у

В силу

того, что

£(fc) является отрицательной при k < 1 и монотонно

убывает

в

диапазоне

значений к от 1 до 0 для всех моделей Ап, функ­

ция k(t)

= <^_1(—Pe(t))

монотонно убывающая (здесь к(0) = <^_1(—_ре(0)) = 1,

а также учтено, что по предположению pf(t) ^ 0). Следовательно, на лицевой поверхности пластины X = 0 функция k(t) монотонно убывает, а функция а(к) возрастает для моделей А\, А^ и убывает для моделей AjV, Ay.

Тогда на фазовой плоскости (t,X) в области возмущения для «+»-ха- рактеристик (как было показано выше, они являются прямыми) тангенс их угла наклона к оси ОХ с ростом t в моделях AjV, Ау увеличивается, а в моделях Aj, AJJ — уменьшается (рис. 5.7.6). Уменьшение же тангенса угла наклона характеристик t = Х/а{к) означает, что характеристики с разными значениями к могут пересекаться на головной волне. В результате возникает неоднозначность решения, что недопустимо в рамках сделанных нами допу­ щений об отсутствии для исследуемого решения скачков функций к, Т и v (допустимы только скачки их первых производных). Таким образом, для моделей Aj и AjXполученное решение неприменимо и оно будет построено иначе (см. п. 5.7.5).

Для моделей AjV и Ау характеристики не пересекаются и полученное решение (5.7.49) действительно имеет место. Графический способ построения решения (5.7.49) по заданным значениям функции pe(t) показан на рис. 5.7.6.

Рассмотрим частный случай нагружения, когда нагрузка ре имеет скачко­

образный характер:

 

Pe=P°e h(t),

(5.7.54)

где h(t) — функция Хевисайда. Тогда Т и к на лицевой поверхности пластины вначале изменяются от 0 и 1 соответственно до конечных значений Т° и fc° < 1, а затем остаются постоянными для всех t > 0. Следовательно, на фазовой плоскости (£, X) образуется угол, ограниченный характеристиками X = aQt, X = a(kP)t и покрытый характеристиками X = a{k)t, к0 < к < 1 («веер» характеристик). Волны, соответствующие этим характеристикам, на­ зывают волнами Римана (центрированными волнами) по аналогии с волна­ ми в газовой динамике (см. т. 3, п. 1.3.16). Эти волны характеризуются тем, что движутся без изменения амплитуды значений Т и v (рис. 5.7.7).

568

Глава 5. Пластические среды

Однако удовлетворить граничному условию (5.7.296) удается, только если допустить существование скачка функций, в области за которым снова спра­ ведливо тривиальное решение

Т = 0,

k = 1, v = 0.

(5.7.556)

Иначе говоря, для моделей А\ и

реализуется

решение в виде ударной

волны.

 

 

Функция X = X]j{t), разделяющая на фазовой плоскости (t,X) два ре­ шения (5.7.55а) и (5.7.556), представляет собой уравнение фронта ударной волны. Для нахождения этой функции, а также для расчета значений к и v (соотношения на характеристиках здесь уже не выполняются) используют соотношения (т. 2, (4.4.17)) на поверхности сильного разрыва в материальном описании. Для рассматриваемой задачи эти соотношения сводятся к следую­ щим:

Mv — p = 0,

(5.7.56а)

< М ( к2- l) + fw = 09

(5.7.566)

+ [е]) - p v = 0.

(5.7.56B)

Здесь [v] = v, [Р] = Рп = сгп, [F] = F\ \ — 1 = к — 1, так как по одну сторону

поверхности разрыва, согласно (5.7.556), среда является невозмущенной.

о

Массовая скорость М определяется по формулам (т. 2, (4.4.4) и (4.2.45)):

М = Dp,

D = с = dx^/dt = dXp/dt.

(5.7.57)

Положим, что на ударной волне скачок температуры равен нулю: [в] = 0, тогда, используя результат упр. 10 к § 5.1, для скачка внутренней энергии [е] имеем

 

 

[е\ = е - е 0 = Г jf ( C e) +

| h ( C e)2 .

(5.7.58)

 

 

 

 

 

2Р

 

Р

 

 

 

 

Здесь учтено, что в невозмущенной области, согласно (5.7.556), е = ео-

Из

(5.7.6)-(5.7.8)

получаем

следующие

выражения для

инвариантов

(п)

(п)

 

 

 

 

 

 

 

 

 

 

/!(Се) И/!(С2):

 

 

 

 

 

 

 

 

 

 

 

(п)

(п)

 

( п )

( п )

 

( и ) „

 

( п )

 

 

 

h { C e) = C \ {

+ 2 с 22 =

с е

+ С Рп = С и ,

 

 

 

 

 

11

 

 

(п)9

( п )

( п )

 

(п)

(п).

2

1

(п)

 

 

(С е

)2 + 2(6*22, 2 = (С п

 

 

, 1

(nV

(5.7.59)

 

1\{С I)

- с*п у

+ № 4 0

 

 

11

 

 

 

 

 

 

 

 

 

Подставляя (5.7.59) в (5.7.58), находим выражение для скачка внутренней

(п)(п)

энергии через деформации С*ц и С \ х\

 

 

§ 5.7. Плоские волны в пластических средах

569

[е] —

( п )

( п )

60

1

 

 

 

С\\

+ 2h { { C п - Cipn )2

+ х(С 'п)2)) —

 

 

 

 

z

 

 

 

 

 

=

1

(п)

(п)

(п)

(п)

 

 

4((*1 +2l2) C 2n - 4 1 2С и С ри + 312(С ри )2). (5.7.60)

 

 

 

2Р

 

 

 

 

При

первоначальном

пластическом

сжатии,

(п) *

= 0, подставляя

когда С \\

 

 

 

(п)

 

 

 

(п)

в (5.7.20) формулу (5.7.10) для Г ц, получаем следующее уравнение для С рп :

( п )

1 /

 

(п)

 

(n)

( n )

 

г -

 

 

 

 

щц

w

+

(5.7.61)

с ри = J-((Z i

+2l2) C n

- 2 1 2C PU - h C u

ask),

 

Hk

 

 

 

 

 

 

 

 

 

 

из которого находим пластическую деформацию

 

 

 

 

 

 

 

 

 

( п )

 

 

 

 

 

 

 

%

_

V3 crsk + ЪС п

Ь =

1

(5.7.62)

 

11

3(6 + Щк/2)

 

2(1 -

^)

 

 

 

 

 

 

 

 

(п)

и (5.7.4) для

( п )

Подставляя в (5.7.60) формулу (5.7.62) для

С рп

С п, полу­

чаем скачок внутренней энергии [е] как функцию от к:

 

 

 

 

 

[е] =

[е](к).

 

 

 

(5.7.63)

 

 

 

 

 

 

 

 

 

 

( п )

 

Если нагружение происходит только в упругой области, где С р = 0, то

 

M W = J

1+Я1 Щ

’‘~Ш -

Щ

 

(5.7.64)

 

 

 

2р(п -

III)2

 

 

 

 

 

С учетом выражения

(5.7.63) три

соотношения

(5.7.56а-в)

позволяют

 

 

 

 

 

 

 

О

через р. Для этого выразим

определить три неизвестные функции: к, v и М

v из (5.7.566):

 

 

v = М( 1 -

 

 

 

 

 

 

 

 

к)

 

 

 

(5.7.65)

 

 

 

 

 

Р

 

 

 

 

 

и, подставив (5.7.65) в (5.7.56а),

найдем

 

 

 

 

 

 

 

 

 

°

/

РР

 

 

 

(5.7.66)

 

 

 

М =

1

- к

 

 

 

 

 

 

 

 

 

 

 

 

Перед корнем выбран знак «+», исходя из физического смысла решения: движение ударной волны должно происходить в положительном направлении оси ОХ пластины.

Подставляя (5.7.66) в (5.7.65), получаем

^ = ^ р ( 1 - к ) / р ,

(5.7.67)