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

книги из ГПНТБ / Кристенсен, Р. Введение в теорию вязкоупругости

.pdf
Скачиваний:
38
Добавлен:
22.10.2023
Размер:
10.86 Mб
Скачать

§ 4.3. Гармонические термовязкоупругие волны

177

Биквадратное уравнение (4.70) определяет ц в виде

г] = ± [— Ы2а + (112а)(Ь2— 4ас)Чі\,

(4.72)

где а, b и с находятся по формулам (4.71). Параметр т] является комплексным, и его действительную и мнимую части можно непосредственно связать с фазовой скоро­ стью и затуханием в зависимости от расстояния, пробе­ гаемого волнами. Можно считать т] заданной действи­ тельной величиной. Тогда со следует искать в комплекс­ ной форме. Эта процедура соответствует ситуации, рассмотренной в § 2.11 на примере задачи о свободных колебаниях, когда было обнаружено затухание волн со временем.

Для удобства интерпретации найдем теперь для г) асимптотическое решение. Для простоты допустим далее, что комплексные модули ф* и т* являются действитель­ ными константами, причем

Ф * ( ц о ) = ф, Т0т*(і(р) = т.

(4.73)

Обозначим через г\м и г\Т асимптотические разложения значений гр соответствующие механическим и тепловым

.волнам. Сначала рассмотрим механические волны и по­ ложим

Ч ]м = ± [(P®2/G3О-®))7, + Cj8 + с2е2 -1----- ],

(4.74)

где

 

е = ф2/ т |G* (tö)|,

(4.75)

а Сі подлежат определению. Зависимость (4.74) выра­ жает разложение по безразмерному параметру связи е, который для большинства материалов может считаться малым. Если им пренебречь, то связи между механиче­ скими и тепловыми явлениями не будет.

Возьмем вторые и четвертые степени цм согласно вы­ ражению (4.74) и подставим эти члены в (4.69), исполь­ зуя (4.73). Далее соберем коэффициенты при одинако­ вых степенях е и положим их равными нулю, чтобы удов­ летворить (4.69). Член нулевого порядка обращается в нуль тождественно. Приравнивание нулю коэффициен­ та при е дает

12-851

178

 

 

Гл.

4.

Распространение волн

 

 

 

 

 

 

_

(—f-m/2fe) [I g; [/(pG‘ )’^]

 

 

(4.76)

 

 

 

 

 

1 +

(imGg/wftp)

 

 

 

 

 

 

 

 

 

 

 

 

Если

используется

представление

Gg (г®) = G3 (со) +

+iGg

(to),

то зависимость

(4.76)

можно

разложить на

действительную и мнимую части,

что дает

 

 

 

 

 

 

 

Сі = с '+

гс",

 

 

 

 

(4.77)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m [Gg|

GgV ^l—mG3/ (akp)y\ -G ln mG'3/(a k p )

(4.78)

 

 

 

2kp<

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m |бз| <$И ( 1 — mGg/(cüftp))— G,(v mG3 /(wßp)

 

 

2kp 1

 

 

 

 

 

 

 

 

(4.79)

 

 

 

 

 

 

 

 

 

 

 

здесь

 

 

 

 

 

 

 

 

 

 

 

 

 

Щ11 I

1

mG3 1

■ о'? —

 

 

 

 

 

 

Д =

(ükp

 

 

 

 

 

 

 

 

 

 

 

 

(ükp

 

 

 

 

 

 

 

 

 

 

 

 

f l

mG3

\

 

mG3

 

 

 

 

 

+

 

GJ11-

 

 

 

 

 

G‘v |l

coép

/

 

 

 

 

 

 

 

\

 

 

ö)£p

При этом выполняется условие

 

 

 

 

 

 

 

[G3 (со) +

iG” (со)]1/2= Gg11 (со) + (G'v (со).

(4.80)

Используя

(4.77) — (4.80),

можно разложить выражение

(4.74)

на действительную и мнимую части

 

 

■м

(Ор^г Gg11

+ Cj е + ■

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

і

■ op■ '■ ОТ'

Cj е + '

 

; (4.81)

 

 

 

(oSn)a +

(oSv)a

 

многоточия

заменяют

члены

более

высокого

порядка

по е.

 

 

мы видим, что и(х,

t) можно записать так:

Из (4.66)

u (x t t) = « е х р ( т ] 2 ^ ) е х р [ і т і і ( х + ® / r i i ) ( ] ,

( 4 . 8 2 )

§ 4.3. Гармонические термовязкоупругие волны

179

где г]м = г)і—it]2. Ситуация, представленная формулой (4.82), характеризует волну, движущуюся в отрицатель­ ном направлении х с фазовой скоростью о = со/г|і и зату­ ханием т]2* Используя действительную и мнимую части (4.81), находим, что фазовая скорость и затухание дают­ ся соответственно формулами

и

 

 

 

 

 

 

где

 

 

 

 

 

 

А, =

[CS»-(m/fflÄp)(G;G“4-G;<?8V)]* +

 

 

 

 

+ [G3V -

(m/(okp)(—G' G‘n + G" G'v)]2;

здесь

опущены

члены более

высокого

порядка

по е

и G* (іи) =Л *(іи )+2р .*(ісо).

При отсутствии связи,

ког­

да е = 0, первые

члены в (4.83)

и (4.84)

дают фазовую

скорость и затухание продольных изотермических вязкоупругих волн. Следующие члены (4.83) и (4.84) дают по­ правки первого порядка к характеристикам изотермиче­ ского распространения волн при наличии термомехани­ ческой связи. В противоположность связанной термоупругости эти поправочные члены могут как увели­ чивать, так и уменьшать изотермическую фазовую ско­ рость в зависимости от величины отношения G3/G3 .

Для тепловых волн параметр т)г принимается в форме

г]т= ± [(— m m jk)и + dx8 + d2е2 + •••], (4.85)

где е также берется из зависимости (4.75). Следуя про­ цедуре, подобной только что описанной, можно получить фазовую скорость и затухание для тепловых волн. Они оказываются равными

2k \7.

— 8

°зI[(°з — Gl)(ö>m/kp) + со2]

f

V — со

Ä

+

(-COOT/

 

 

 

( 4 . 8 6 )

12*

180

 

Гл.

4.

Распространение волн

 

И

 

 

 

 

 

_

(^HL\k _i

/

 

ыт 1G3 1[ ( ° з + ° з ) («а/п//~) - со2]

_____

П2 “

[ 2k ) +

8 [ 2k

)

АрА2

' ' ’

 

 

 

 

 

(4.87)

где

А2 = [[Gj — G])(com/&p) + w2]2 + [(Gg + G’3)((am/kp) — со2]2.

Вторые члены в (4.86) и (4.87) дают вызванные связно­ стью поправки первого порядка для фазовой скорости затухания несвязанных тепловых волн.

Этим завершается анализ распространения гармони­ ческих волн в связанной теории термовязкоупругости. Хотя с точки зрения механических характеристик рас­ сматриваемые волны являются волнами расширения, мы видим, что распространение гармонических волн сдвига, т. е. волн, для которых Ыі,і 0, не будет сопровождаться связанными термическими эффектами. Действительно, повторяя анализ, проведенный здесь для волн расшире­ ния, можно показать, что фазовая скорость и затухание плоской гармонической волны сдвига даются первым членом в (4.83) и (4.84) с той разницей, что G\ (Ію) за­

меняется на

р*(ію). Таким образом, фазовая скорость

и затухание

гармонических волн сдвига

определяются

формулами

 

 

 

 

V =

1[I* |/|ХШ Р^2

(4.88)

 

Т]2 =

юр^рІѴ/| ц*|,

(4.89)

 

[р* (Ію)]'/г = рШ (ю) + ірП (ю).

(4.90)

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

§ 4.4. Отражение гармонических волн

181

В только что проведенном анализе продольных волн использовалась линейная теория, описываемая уравне­ ниями § 3.1, и приходилось пренебрегать эффектом дис­ сипации энергии и зависимостью вязкоупругих механи­ ческих свойств от температуры, за исключением зависи­

мости от

фиксированной

исходной

температуры Т0.

В задаче

рассмотренного

здесь типа

можно ожидать,

что роль этих эффектов будет значительной, и в силу этого предшествующие рассуждения являются в высо­ кой степени идеализированными. Исследования, в кото­ рых учитывались эти нелинейные эффекты, связанные с энергией диссипации и зависимостью механических свойств от температуры, провели Петров и Грач [4.19], Волосевик и Грач [4.26] и Хантер [4.13].

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

§ 4.4. Отражение гармонических волн

В задачах, обсуждавшихся до сих пор в этой главе, не рассматривалось непосредственно влияние граничных условий на распространение волн. Такого рода эффекты учитывались лишь косвенно в примерах § 4.2. Исследуем теперь общую задачу об отражении гармонических волн от свободной поверхности вязкоупругого полупростран­ ства. В дальнейшем изложении мы следуем работе Лок-

кета [4.17].

Преобразование Фурье уравнений движения таково:

Р* Uk,jj + [Я* (гео) + р* (ІО))] и P j k = — PO)2 uk. (4.91)

Эти уравнения движения можно разделить, как это бы­ ло сделано в (4.30) и (4.31), поочередно, полагая

182

Гл. 4. Распространение волн

Щ,і{хи

t ) = 0 и rot и(Хі, 0 = 0 , где u — вектор перемеще­

ния. Получающиеся в результате уравнения движения имеют вид

р* ик,ц = — pco2uk,

(4.92)

(X* + 2р*) tikjj

рсоaMfe.

(4.93)

Соотношения (4.92) и (4.93)

описывают

распростране­

ние соответственно поперечных и продольных волн. Эти

волны иногда называют

соответственно S -волнами и

P -волнами. Равенства

(4.92)

и

(4.93)

можно

записать

в виде

 

 

 

 

 

 

 

[cp fu k.jj

= — а

uk,

[ c s f Uk.tj = — со2uk,

(4.94)

где

 

 

 

 

 

 

 

('cp)2 =

р

+

2р*),

(cs)2 =

р 7 р -

(4.95)

Для гармонических волн вектор перемещения и за­

дается формулой

 

 

 

 

 

 

 

u =

U ехр [і (Lx -f- Му -f- <а00 ],

(4.96)

где U — вектор комплексной амплитуды, L и М — комп­ лексные параметры, а ©= соо — частное значение интере­ сующей нас частоты. При записи (4.96) предполагается, что движение волны в декартовой системе координат X, у, z не зависит от координаты г . Преобразование Фурье уравнения (4.96) дает

ц = иеі(і*+л1у) [2я6 (со — со0)J ,

(4.97)

где 6( ) — дельта-функция. Подставляя (4.97) в (4.94), получаем

(c;)2(La+ M

2) = o)2,

(4.98)

где с*0 можно отождествить с

или с*р при (о==(оо. В дей­

ствительности уравнение (4.98) можно получить и не прибегая к преобразованию Фурье, поскольку в данном случае рассматриваются только гармонические волны, однако здесь оно применено для удобства.

Параметры L и М в общем случае являются комп­ лексными. Однако перед тем, как рассматривать самый общий случай, рассмотрим сначала частный случай,

§ 4.4. Отраокение гармонических волн

183

в котором отношение L/М является действительным числом.

Случай действительного отношения ЫМ

Когда отношение комплексных чисел L/M действи­ тельно, решение (4.98) дает

L — (о)0/со) cos а

 

и

М — (со0/со) sin а.

(4.99)

Отношение &oIcq записывается в виде

 

 

 

 

щ!со =■

а (1 + іф),

 

 

(4.100)

который определяет а

и

ф

через

to0

и

с*0. Из

(4.100)

и (4.99) следует, что

 

 

 

 

 

 

 

L = /(1 + іф)

и

М =

/п( 1

+

іф),

(4.101)

где

 

 

 

 

 

 

 

 

/ =

а cos а

 

и

m ~ a sin а.

(4.102)

С учетом (4.101)

уравнение

(4.96)

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

u — U ехр[—ф(/л: + т у)]ехр[і(1х +

Ѵш/ +

и0і)].

(4.103)

Из (4.103) мы видим, что в этом частном случае волна затухает в зависимости от пройденного расстояния и на­ правление распространения совпадает с направлением затухания. Волна, обладающая этим свойством, по тер­ минологии работы [4.17] называется простой волной. Из (4.102) мы видим, что tg a = m//, где угол а, отсчиты­ ваемый от оси X , определяет направление распростра­ нения. Параметр ф определяет характеристику затуха­ ния. Рассмотрим теперь общий случай, когда отношение L/M не является действительным.

Случай комплексного отношения L/M

В этом случае отношение комплексных постоянных L и М является комплексным. При таких обстоятельст­ вах затухание происходит не в направлении распрост­

ранения

волны, и формула, соответствующая

(4.103),

имеет вид

 

u =

U exp (ßx + уу) exp [i (lx + тг/ + сй0і )],

(4.104)

184

Гл. 4. Распространение волн

где теперь

ß/v Ф. Чш.

Этот тип волны является более общим, чем гармониче­ ские волны, рассмотренные в § 4.3, где направления рас­ пространения и затухания считались совпадающими.

В дальнейшем будет удобнее использовать не вектор­ ные, а компонентные обозначения. Таким образом, пре­ образование компонент перемещения и и ѵ в направле­ ниях X и у принимается в виде

 

[и, ѵ) =

2яб (со —шо) (U, V) еҢ1х+Му),

(4.105)

где U и V — комплексные постоянные.

Два уравнения

движения (4.91) и (4.105) принимают вид

 

 

р* (L2 +

М2) U +

(А,* +

р*) L (LU +

MV) =

рсо2 U

(4.106)

и

 

 

 

 

 

 

 

 

р* (L2 +

М2) V +

(А* +

р*) М (LU +

МѴ) =

рсо2 V.

(4.107)

Для 5-волн, комбинируя

(4.98)

и

(4.106), находим

 

Ѵ =

— (L /М) U,

5-волны.

 

(4.108)

Для P -волн, комбинируя (4.98) и (4.107), после не­

которых упрощающих преобразований получаем

 

 

V = (M /L)U ,

Р-волны.

 

(4.109)

Теперь, используя (4.105), (4.108) и (4.109), можно пе­ реписать компоненты перемещения в виде

(“ . ѵ ) I s -волны =

( M t f . - Щ

И

(4.110)

(«, Дя.волны =

(LU, MU) е ^ х+ мУ+*J).

Решения (4.110) далее будут конкретизированы для сле­ дующего случая.

Отражение падающей 5-волны от свободной границы

Пусть плоскость х = 0 определяет свободную от на­ пряжений границу полупространства. Найдем характе­ ристики отражения падающей 5-волны от этой границы, если угол падения равен сц.

§ 4.4. Отражение гармонических волн

185

Обозначим через us i, vS(

компоненты вектора пере­

мещения падающей S -волны, а через us , vs , ир

и

ѵр

компоненты перемещения для отраженных S- и Р-волн.

Тогда в соответствии с (4.110)

получим

 

 

{ustt üSi! = 2л8((о щ ) (M s iU s i, L s i U s i )

 

 

 

Xexp [i(L siX + MSiy)\,

 

(4.111)

(«s, os) =

2яб(ю — co0) {Ms Us, —Ls Us ) exp [i{Ls x +

Msy)\

 

 

 

 

(4.112)

и

 

 

 

 

[up, vp) =

2я8(ю — co0) [Lp Up, Mp Up) exp [i[LPx +

MPy)\.

 

 

 

 

(4.113)

Вектор полного перемещения получается суперпозицией компонент всех трех типов (4.111) — (4.113):

U = U si Т - U s + Up И V — V st + V s + V p .

Уравнения движения удовлетворяются в силу условия (4.98), для выполнения которого в данном случае тре­ буется, чтобы

Lsi + М%і= Р І +

М% =

pcoo/fx »

(4 114)

Lp +

Mp =

р (Oo/(L +

2[а ).

 

 

 

Условия отсутствия внешних усилий на границе х = 0

записываются следующим образом:

 

 

 

вхх{х,

у,

0

=

0

при

X

=

0

(4.115)

и

 

 

 

 

 

 

 

 

 

Оху (X,

г/,

0

=

0

при

X

=

0.

(4.116)

Условие (4.116), переписанное для преобразований пере­ мещений, имеет вид

ди(х, у, а )/д у + дѵ(х, у, а )/д х = 0 при х = 0 . (4.117)

Подставляя выражения (4.111) — (4.113) в (4.117), по­ лучаем уравнение, содержащее экспоненциальную зави­

186 Гл. 4. Распространение волн

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

MSi = Ms = Мр = М,

(4.118)

где М подлежит определению. Выполнение (4.118) необ­ ходимо также для того, чтобы удовлетворялось гранич­ ное условие (4.115). Подставляя (4.118) в (4.114), на­

ходим

/4 . = L| = p o ^ * - M 2

и

(4.119)

Ь% = рсо2/(Г +

2

Поскольку L Si относится

к падающей поперечной

волне, а L s — к отраженной поперечной волне, они дол­

жны.иметь противоположные знаки, в силу чего, исполь­ зуя первую из зависимостей (4.119), получаем Ls = —LSj. Учитывая это соотношение и (4.118), мы видим, что

Ls/Ms = - L s JM sr

(4.120)

Однако для простой падающей поперечной волны отно­ шение LS .IMS[ является действительным согласно опре­

делению, данному в связи с (4.103). Из (4.120) видно, что отраженная поперечная волна также является про­ стой и что эта волна испытывает затухание в направле­ нии распространения. Из отношения двух выражений

(4.99) и из (4.120) находим, что

a s == Ms/Ls — MS{/LS[— tg a ;

(4.121)

и что для поперечных волн угол отражения равен углу падения.

Находим теперь характеристики отраженной Р-волны. Поскольку падающая волна является простой, уравне­ ния (4.99) и (4.118) дают

Мр = (p1/2« 0sina.)/(p*)‘/*.

(4.122)

Это равенство вместе с последним уравнением (4.119) приводит к зависимости

1. (4.123)

(X* + 2 (і*) sin2 а,-