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

Основы_акустики_Гринченко_Вовк

.pdf
Скачиваний:
56
Добавлен:
10.04.2023
Размер:
16.81 Mб
Скачать

Всесторонним растяжением (или всестороннее сжатие в зависи- мости от знака деформации) называют деформацию, при которой уд- линение одинаково вдоль всех трех осей, а деформации сдвига отсут- ствуют:

 

 

u11 = u22 = u33 ,

u12 = u13 = u23 = 0.

 

 

(6.88)

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

 

 

 

 

 

σ

= σ

22

= σ

33

=

λ + 2

μ u

, σ

= σ = σ

23

= 0.

(6.89)

11

 

 

 

3

kk

12

13

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отличны от нуля только нормальные напряжения. Величину

 

 

 

 

 

 

 

K = λ + 2

μ

 

 

 

(6.90)

 

 

 

 

 

 

 

3

 

 

 

 

 

называют модулем всестороннего сжатия или объемным модулем упругости.

Деформация чистого сдвига, скажем, в плоскости х1х2, есть де- формация, при которой отличны от нуля только компоненты u12 = u21 тензора деформаций. Из выражения (6.86) получим

 

 

 

 

 

u1

 

 

 

σ

= σ

21

= 2μu

= μ

+ u2

,

(6.91)

 

12

 

12

 

x2

x1

 

 

 

 

 

 

 

 

σ11 = σ22 = σ33 = σ13 = σ23 = 0.

 

(6.92)

Таким образом, коэффициент Ламе μ представляет собой модуль сдвига G, который определяется при чистом сдвиге из соотношения μ = σ12/(2u12) = G. Понятно, что постоянные Ламе λ и μ имеют размер- ность поверхностных сил.

Теперь рассмотрим однородное растяжение цилиндрического стрежня со свободной боковой поверхностью. Направим ось Ox1 вдоль оси стрежня. Единственной отличной от нуля компонентой на- пряжения будет σ11, поскольку на боковых стенках стрежня напря- жения равны нулю. В этом случае из уравнений (6.86) получим

σ11 = (λ + 2μ)u11 + λu22 + λu33,

(6.93)

σ22 = λu11 + (λ + 2μ)u22 + λu33 = 0,

(6.94)

σ33 = λu11 + λu22 + (λ + 2μ)u33 = 0.

(6.95)

Уравнения (6.94) и (6.95) выполняются при условии

331

u22 = u33 = −

λ

 

u11.

(6.96)

 

 

 

2(λ + μ)

Подставив (6.96) в уравнение (6.93), получим

 

σ

=

μ(3λ + 2μ)

u .

(6.97)

 

 

 

 

11

 

 

 

λ + μ

11

 

 

 

 

 

 

 

 

Величину отношения σ11/u11 обозначают

 

 

 

E =

μ(3λ + 2μ)

 

 

(6.98)

 

λ + μ

 

 

 

 

 

 

 

 

 

и называют модулем Юнга* для стрежня, а величину отношения

u22/u11 = u33/u11 обозначают

 

 

 

 

ν =

 

λ

 

 

(6.99)

 

 

 

 

 

 

2(λ + μ)

 

 

и называют коэффициентом Пуассона. Коэффициент Пуассона свя- зывает деформации в продольном и поперечном направлениях. Знак минусв (6.96) характеризует экспериментальный факт: при растя-

жении стрежня в продольном направлении происходит уменьшение его поперечного размера.

Определим физический смысл первого коэффициента Ламе λ. Снова рассмотрим растяжение стрежня вдоль оси Ох1, но при усло- вии, что боковая поверхность стрежня зажата и, следовательно, не- подвижна. В этом случае компоненты тензора деформаций равны ну- лю, кроме u11. Тогда вместо (6.93)-(6.95) следует записать

σ11 = (λ + 2μ)u11 ,

σ22 = λu11, σ33 = λu11.

Таким образом, коэффициент λ = σ22/u11 = σ33/u11 равен отношению поперечного нормального напряжения к деформации стрежня при растяжении стрежня с зажатыми боковыми гранями.

Поскольку в закон Гука для изотропного тела входят только две независимые характеристики вещества, то между модулями упруго- сти λ, μ, К, Е, ν, G может быть только линейная зависимость. Если основными являются постоянные Ламе λ и μ, то

* Юнг (Young) Томас (1773—1829)английский физик.

332

 

 

μ(3λ + 2μ)

2

 

 

 

 

 

λ

 

 

 

E =

 

 

 

, K = λ + 3

 

μ ,

 

ν =

 

 

 

;

 

(6.100)

λ + μ

 

2(λ + μ)

 

если основной парой есть Е и ν, то

 

 

 

 

 

 

 

 

 

 

 

λ =

 

νE

,

μ = G =

 

E

 

,

K =

E

;

(6.101)

(1+ ν)(12ν)

2(1+ ν)

3(12ν)

если основной парой есть К и G, то

 

 

 

 

 

 

 

 

 

 

 

 

λ = K 3G ,

E =

9KG

 

,

ν =

3K 2G ;

 

(6.102)

 

3K +G

 

 

 

 

2

 

 

 

 

 

 

6K + 2G

 

 

если основной парой будет G, ν, то

 

 

 

 

 

 

 

 

 

 

 

 

λ =

2νG

,

E = 2G (1+ ν),

K =

2G (1+ ν)

.

 

(6.103)

 

 

 

 

 

 

12ν

 

 

 

 

 

 

 

 

3(12ν)

 

 

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

E < 3K , 2μ < E < 3μ , K < λ + 2μ < 3K , 0 < ν < 0,5.

Коэффициент Пуассона, близкий к нулю, имеет пробка: при растяже- нии и сжатии поперечные размеры тела из пробки практически не изменяются. Предельный случай ν 0,5 соответствует μ 0, т.е. пе- реходу от твердого тела к жидкости.

В таблице приведены характерные величины для ряда упругих

материалов (о величинах ce и ct

будем говорить позднее).

 

 

Материал

ρ 10–3,

Е 10–10,

 

ν

G 10–10,

ce , м

сt , м

 

 

кг/м3

Н/м2

 

 

Н/м2

с

с

 

Алюминий

2,7

7,0

 

0,35

2,6

6420

3110

 

Вольфрам

18,7

36,0

 

0,29

13,3

5230

2860

 

Железо

7,87

21,0

 

0,28

8,2

5850

3230

 

Золото

19,3

8,0

 

0,42

2,9

3140

1170

 

Кварц

2,2

7,4

 

0,17

3,2

5970

3780

 

Латунь

8,5

9,0

 

0,33

3,0

4240

2140

 

Медь

8,9

12,0

 

0,33

4,6

4560

2250

 

Олово

7,18

5,4

 

0,33

2,0

3320

1670

 

Свинец

11,34

1,6

 

0,43

0,6

2120

740

 

Сталь

7,86

2,2

 

0,29

8,5

5890

3210

 

Стекло

2,5

7,2

 

0,25

2,9

5800

3350

 

 

 

 

 

 

 

 

 

333

6.6. Общие уравнения распространения волн в твердых телах

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

Рис. 6.11. Мысленная элементарная частичка среды твердого тела

Рассмотрим компоненту движения вдоль одной из трех осей коор- динат, например оси Ох2 (рис. 6.11). Сила, которая действует по второ- му закону Ньютона, приложена к граням элементарного параллелепи- педа из окружающей среды. Она выражается через механические напряжения на гранях параллелепипеда. На рис. 6.11 показаны про- екции напряжения на ось Ох2. Запишем проекцию силы на ось Ох2:

dF2 = (σ12′ − σ12 )dx2dx3 + (σ′22 − σ22 )dx1dx3 + (σ32′ − σ32 )dx1dx2 . По-

скольку σ12(x1 + dx1,x2,x3 ) , σ′22(x1,x2 + dx2,x3 ), σ′32(x1,x2,x3 + dx3 ), то,

применяя формулу Лагранжа о приращении функции при малом из- менении аргумента, получаем

 

 

∂σ

∂σ

∂σ

32

 

 

 

∂σ

 

 

dF

=

12 +

22

+

 

dx dx dx

=

i2

dx dx dx

. (6.104)

 

 

 

xi

2

 

x1

x2

x3

1 2 3

 

1 2 3

 

Под действием этой силы частичка будет двигаться вдоль оси Ох2 с ускорением 2u2/t 2, где u2 проекция на ось Ox2 вектора смещения u = (u1, u2, u3) частицы. Масса частицы равна ρdx1dx2dx3, где ρ плотность среды. Итак, уравнение движения частицы вдоль оси Ox2

имеет вид ∂σi2 = ρ 2u2 . Рассматривая аналогично движение частицы

xi t2

вдоль других осей координат, получаем систему из трех уравнений:

334

∂σik

= ρ

2uk

,

k = 1,2,3.

(6.105)

 

xi

 

t2

 

 

Это и есть уравнения движения (для ускорения использована частная производная 2/t 2 вместо d 2/dt 2, что обусловлено пренебрежением ускорения переноса (см. п. 4.1.2)).

Уравнение (6.105) связывают две функции координат и времени: тензор напряжений σik (x1,x2,x3,t ) и вектор смещения u(x1,x2,x3,t). Для того чтобы число уравнений соответствовало числу неизвестных функций, необходимо иметь еще одно уравнение, которое связывает σik и u. Таким уравнением является обобщенный закон Гука

(см. (6.86)):

 

 

ue

 

ui

 

uk

 

 

σ

= λδ

+ μ

+

.

(6.106)

 

xk

 

ik

ik xe

 

 

xi

 

Подстановка (6.106) в (6.105) позволяет получить уравнение относи- тельно проекций вектора u:

 

 

 

 

ue

 

 

 

 

 

 

ui

 

 

2

 

 

 

 

2

 

 

 

 

λ

 

 

 

 

δ

+ μ

 

 

+ μ

uk

 

= ρ

uk

,

k = 1,2,3,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t2

 

 

 

xi

xe

 

ik

 

 

xi

xk

 

 

xi xi

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ue

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

λ

 

 

 

+ μ

 

 

ui + μ

uk

 

= ρ

uk

,

 

k = 1,2,3.

(6.107)

 

 

 

 

 

 

 

xi xi

 

 

 

 

 

 

xk

 

xe

 

xk

xi

 

 

 

t2

 

 

 

Поскольку ue

divu , 2uk

u

(провели суммирование по индексу i),

 

 

 

xe

 

 

 

 

 

 

xi2

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

то уравнение (6.107) можно переписать в виде

 

 

 

 

 

 

 

 

 

 

 

(λ + μ)

 

divu + μΔu

= ρ 2uk ,

k = 1,2,3.

(6.108)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xk

 

 

 

 

 

 

k

 

t2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Умножим каждое из трех уравнений (6.108) на соответствующий ба- зисный вектор ej, j = 1,2,3, и просуммируем. Поскольку e1 x1 + e2 x2 +e3 x3 grad , то искомое уравнение запишем в

виде

(λ + μ)graddivu + μΔu = ρ

2u

.

(6.109)

t

2

 

 

 

335

При решении конкретных задач часто используют другую форму урав- нения движения. Учитывая тождество u = graddivu — rotrotu, пере- писываем (6.109) так:

(λ + 2μ)graddivu − μrotrotu = ρ

2u

.

(6.110)

t

2

 

 

 

Известно [8, 52], что любое векторное поле можно представить в виде суммы потенциального (от латинского слова potentia сила) вектора, ротор которого равен нулю, и соленоидального (от греческого слова swlhu (солен) — трубка) вектора, дивергенция которого равна нулю. Если записать в таком виде вектор смещения u, то можно по- лучить отдельно уравнение потенциальной и соленоидальной частей смещения. Действительно, пусть u = ue + ut, считая rotue = 0 и divut = 0; подставляя в (6.110), получаем

(λ + 2μ)graddivu

− μrotrotu

= ρ

2u

2u

(6.111)

e + ρ

t .

e

t

 

t2

t2

 

 

 

 

 

Поскольку разложение вектора на потенциальную и соленоидальную части является единственным, то согласно (6.111) имеем

(λ + 2μ)graddivu

= ρ

2u

(6.112)

e ,

 

 

e

 

t2

 

 

 

 

 

 

−μrotrotu

= ρ 2ut .

(6.113)

 

t

 

t2

 

 

 

 

 

Применяя тождество rotrot = graddiv –

к векторам ue и ut, находим

graddivue = ue, rotrotut = – ut, (6.114)

 

 

 

кроме того, введем обозначения

 

 

 

 

ce =

λ + 2μ

,

ct =

μ.

(6.115)

 

ρ

 

 

ρ

 

Подставляя (6.114), (6.115) в (6.112), (6.113) получаем волновые урав- нения для ue и ut :

 

 

 

1

2u

 

 

 

 

1

2u

 

u

e

=

 

 

 

e ,

u

t

=

 

 

t .

(6.116)

 

 

 

 

 

 

 

ce2 t2

 

 

ct2 t2

 

 

 

 

 

 

 

 

Для гармонических волн, т.е. при временной зависимости exp(–iωt), уравнения (6.116) превращаются у уравнения Гельмгольца:

u

+k2u

e

= 0,

u +k2u

t

= 0,

(6.117)

e

e

 

t t

 

 

336

где ke = ω/ce, kt = ω/ct волновые числа соответственно для ue и ut.

6.7. Продольные и поперечные волны

В пятом разделе мы исследовали решение волнового урав- нения в виде плоской продольной волны в идеальной сжимаемой жидкости. Рассмотрим теперь, что представляют собой плоские вол- ны ue и ut в твердом теле, как решение уравнений (6.116). Очевидно постоянные ce и ct (см. (6.115)) в волновых уравнениях (6.116) имеют смысл скорости распространения волн. Рассмотрим плоские гармо- нические волны, которые распространяются вдоль оси Ох1 декарто- вой системы координат:

u

= A(e) exp[iω(t x

c

e

)],

(6.118)

e

1

 

 

 

u

= A(t ) exp[iω(t x

c

t

)],

(6.119)

t

1

 

 

 

 

где A(e) = e1A1(e) + e2A2(e) + e3A3(e) и A(t ) = e1A1(t) + e2A2(t ) + e3A3(t ) векторы

амплитуд.

Определим характер колебаний частиц среды в этих волнах, опи- раясь на основные свойства векторов ue (rotue = 0) и ut (divut = 0). Ис- пользуя (6.118), вычисляем rotue и записываем уравнение

 

e

 

A(e) + e

 

A(e)

iω

exp[iω(t x

c

e

)] = 0.

 

 

 

 

 

2

3

3

2 сe

1

 

 

Отсюда имеем A3(e) = 0, A2(e) = 0. Итак, вектор ue, направлен вдоль оси

Ox1, т.е. движение частиц среды происходит в направлении распро- странения волны. Таким образом, волна ue(x1, x2, x3, t) является про- дольной волной.

Используя (6.119), вычисляем divut и записываем уравнение

A(t ) iω exp[iω(t x ct )] = 0.

1 ct

Отсюда имеем A1(t) = 0. Итак, вектор ut ориентирован перпендикуляр-

но к направлению распространения волны (оси Ох1). Таким образом, волна ut(x1, x2, x3, t) является поперечной волной. Поскольку в попе- речной волне частицы среды двигаются в плоскости, перпендикуляр- ной к направлению распространения волны, то возникает вопрос о характере движения частиц в этой плоскости. Это свойство называ-

ют поляризацией поперечной волны.

Пусть ось Ох3 декартовой системы координат совпадает с направ- лением волнового вектора поперечной волны. Тогда произвольное смещение частиц u в такой волне можно представить в виде суммы

337

двух взаимно ортогональных векторов u1, u2, направленных вдоль осей соответственно Ox1 и Ox2. При этом возможны ситуации, кото- рые показаны на рис. 6.12.

Рис. 6.12. Формы поляризации поперечной волны

Если u1 = 0, то имеем волну, в которой смещение лежит в плоско- сти х2Ох3 и изменяется по закону (рис. 6.12, а):

u2 = A2e2 cos(ωt kt x3), u1 = u3 = 0.

(6.120)

Выражение (6.120) определяет линейно поляризованную волну. Аналогично определяется линейно поляризованная волна, смеще-

ние частиц которой лежит в плоскости х1Ох3 (рис. 6.12, б) и задается выражением

u1 = A1e1 cos(ωt ktx3), u2 = u3 = 0.

(6.121)

Произвольная линейная комбинация выражений (6.120) и (6.121) также определяет волну, которая поляризована в плоскости х1Ох2. Ес- ли взять суперпозицию движений (6.120) и (6.121) с произвольным сдвигом фаз α, то получим движение в виде

u = A2e2 cos(ωt kt x3 ) + A1e1 cos(ωt kt x3 + α),

(6.122)

где со временем конец вектора смещения описывает

эллипс

(рис. 6.12, в); имеем поперечную волну с эллиптической поляризаци- ей. В частном случае А2 = А1, α = ±π/2 получим поперечную волну с круговой поляризацией (рис. 6.12, г).

Согласно (6.118) и (6.119) продольные и поперечные волны рас- пространяются с фазовыми скоростями соответственно ce и ct. Из формул (6.115) видно, что скорость продольных волн всегда больше

338

скорости поперечных. В конце параграфа 6.5 в таблице указаны величины скоростей ce и ct для некоторых веществ.

Целесообразно привести ряд формул, которые определяют ce и ct через различные упругие постоянные:

 

λ + 2μ

 

E

 

(1 v)

 

 

 

 

3K + 4G

 

2G

1

− ν

)

 

 

ce =

=

 

 

=

 

=

 

(

 

,

(6.123)

ρ

ρ (1 +v)(1

2v)

 

3ρ

ρ(1 2ν)

 

 

 

 

 

 

 

 

 

 

 

ct =

μ

=

E 1

=

G .

 

 

 

 

 

 

(6.124)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ

2ρ1+v

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ

 

 

 

 

 

 

 

Отсюда между скоростями продольной и поперечной волн существует такая связь:

 

ce

=

2

 

1v

.

(6.125)

 

 

 

ct

 

12v

 

Это отношение для разных металлов изменяется довольно существен- но, например, для цинка (ν 0,25) и свинца (ν 0,44) отношение ско- ростей, вычисленное по (6.125), определяется числами 1,73 и 3,3.

В конце параграфа сделаем важный вывод. Как видим, уравнение движения (6.111) можно заменить двумя волновыми уравнениями (6.116). Решение первого представляет собой продольную волну, а решение второго поперечную. Таким образом, произвольная пло- ская волна в упругой изотропной среде может быть представлена в виде суперпозиции продольной и поперечной волн, каждая из кото- рых распространяется со своей фазовой скоростью. Определить фа- зовую скорость для волнового движения в виде суперпозиции этих волн невозможно. В сущности, здесь имеем два волновые движения, которые не взаимодействуют друг с другом.

6.8. Граничные условия для твердых тел

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

339

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

Рассмотрим важные модели граничных условий на границе твер- дого тела с другими телами или с вакуумом; разнообразие здесь большее, чем для жидкостей. Будем обозначать ось, которая совпада- ет с нормалью к границе, индексом n, а две взаимно перпендикуляр- ные оси в плоскости границы индексом α, α = 1,2.

1. Свободная граница. Компоненты тензора напряжений: равны

нулю σnn = σnα = 0.

2. Граница с абсолютно жестким телом при наличии склеива-

ния”. Все компоненты смещения точек границы равны нулю:

un = uα = 0.

3. Граница с абсолютно жестким телом при наличии смазки”.

Нормальное смещение точек границы и касательные напряжения на границе равны нулю: un = 0, σnα = 0.

4. Граница с жидкостью. Нормальное напряжение на границе рав- но давлению в жидкости, взятому с противоположным знаком; каса- тельные напряжения равны нулю; нормальные скорости частиц твердого тела и жидкости на границе равны между собой:

σnn = –p, σnα = 0, un = υn .

t

5. Граница с другим твердым телом при наличии склеивания”.

Все компоненты смещения обоих тел на границе и одноименные ком- поненты тензора напряжения попарно равны: ui(I) = ui(II), i = 1,2,3;

σnn(I) = σnn(II), σn(I)α = σn(II)α .

6. Граница с другим твердым телом при наличии смазки”. По-

парно нормальные смещения и нормальные напряжения обоих тел на границе равны; касательные напряжения в обоих телах равны нулю:

un(I) = un(II), σnn(I) = σnn(II), σn(I)α = σn(II)α = 0.

6.9. Отражение плоской гармонической волны от свободной границы полупространства

Рассматривая волновое движение в однородном бесконеч- ном пространстве, мы поняли очень важную вещь: любое волновое движение в твердом теле всегда можно представить в виде суперпо- зиции продольных и поперечных волн. Наличие двух типов волновых движений в твердом теле значительно усложняет решение конкрет- ных волновых задач по сравнению с подобной ситуацией для акусти- ческих волн в жидкой среде. Ярким и очень важным примером явля- ется задача падения плоской упругой волны на плоскую границу раз-

340