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

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

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

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

187

Отраженная P -волна будет простой, если отношение

Lp jMp является действительным. Для

этого в соответ­

ствии с (4.123) требуется, чтобы

 

 

(К* + 2р*) было действительным числом и

р*/(А* + 2 p *)^ s in 2 щ.

 

(4.124)

Если условие (4.124) не выполняется,

то

отраженная

P -волна не будет простой. Этот случай мы рассмотрим позже. Записывая для р* и X* действительные и мнимые

части, из требования, чтобы отношение в (4.124)

было

действительным, можно заключить, что

 

Т/Дсоо) /А/(<йо) = i i " ( ( l>o) / [ i '((O q) ,

(4.125)

где р' и X' — действительные, а р" и X" — мнимые ча­

сти р* и Я* соответственно. Соотношения (4.124)

вместе

с (4.125) приводят к требованию

 

sin2a i ^ p /(o»o)/[^'(coo) + 2 р /((о0)]-

(4.126)

Неравенство (4.126) определяет критический угол паде­ ния. Если а,- меньше критического угла и если удовлет­ воряется (4.125), то отраженная P -волна будет затухать в направлении распространения. Если хотя бы одно из этих условий не выполняется, то отраженная P -волна не является простой. Этот случай мы сейчас и рассмотрим.

Если условие (4.125) не удовлетворяется, то следует применить более общий подход. Комбинируя последнее

из соотношений (4.119) с условием

из (4.99),

имеем

 

Lp = рсооДЯ f 2р ) — (р<»о/р*) sin2

(4.127)

Угол отражения P -волны определяется формулой

tgcsP = MplLp,

(4.128)

где Мр и L'p— действительные части Мр

и Lp . Отно­

шение мнимых частей M"p 'iL'p определяет направление затухания. Для вычислений по формуле (4.128) Lp бе­ рется из (4.127), а Мр — из (4.122).

Более общую информацию можно получить из пол­ ного решения задачи с граничными условиями (4.115) и

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

(4.116). С использованием уже полученных результатов

эти два уравнения принимают вид

 

 

 

(М2 - Lie) Us + 2MLP Up = -

(М2 -

L b ) USi

(4.129)

и

 

 

 

2Ls t MUs + { M * - L l t)U p =

2Ls tMUsr

(4.130)

Из этих формул можно определить величины Us

и Uр ,

подставив Lp из (4.127), M = MSi = М р

из (4.122) и L Si

также из (4.122), но с заменой sin а* на cos а*. Они будут комплексными постоянными, которые определяют ампли­ туды и фазовые углы отраженных волн.

Аналогичные результаты можно получить и в случае, когда падающая волна является не рассмотренной выше S -волной, а P -волной. Можно показать, что отраженная P -волна является простой и имеет угол отражения, рав­ ный углу падения. Для того чтобы отраженная S -волна была простой, необходимо, чтобы К"/к' = ц"І\і' и чтобы угол падения был меньше некоторого критического зна­ чения, которое можно найти. Когда Х"/к'ф\і"/\і', отра­ женная S -волна не будет простой, т. е. будет затухать в направлении, отличном от направления распростране­ ния.

Более общую задачу отражения и преломления гар­ монических волн на поверхности раздела между двумя вязкоупругими средами рассматривал Локкет [4.17]. Задачи этого вида исследовали также Купер и Рейсс [4.8, 4.9]. В этих работах рассмотрена также проблема возникновения поверхностных волн. Относительно по­ верхностных волн типа Рэлея в вязкоупругой среде опу­ бликованы экспериментальные результаты с соответст­ вующим теоретическим анализом [4.22].

§ 4.5. Движущиеся нагрузки на вязкоупругом полупространстве

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

§ 4.5. Движущиеся нагрузки на полупространстве

189

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

Пусть в декартовых координатах граница полупрост­ ранства определяется уравнением х2 = 0. Нагрузка, а следовательно, и решение задачи не зависят от х3. Гра­ ничные условия имеют вид

<т22 (хи х2, t) = — ph (хг +

vt) при х2= 0

и

(4.131)

од (Ч, Ч , 0 = 0

при х2 = 0.

Таким образом, давление в виде ступенчатой функции перемещается по поверхности полупространства со ско­ ростью и в направлении отрицательной оси х\.

Обозначим через щ компоненты скорости. Вязкоупру­ гие соотношения между напряжениями и деформациями (2.2) в сочетании с уравнениями движения (2.6) дают

f\i(t — T)v.'j j (xl,x)dx +

+J \Mt—x)-\-\k{t—x))v.'{.(x.-\ T)dT=pdvi(xi,t)ldt, (4.132)

—- oo

где 7=1, 2.

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

ется, чтобы выражение для щ зависело от Х\ и (

только

в комбинации Х\ + v t, откуда

 

ѵі = fi(x 1 + vt, х2).

(4.133)

Поскольку Ѵі имеет вид (4.133), естественно произвести замену переменных

190

Гл. 4.

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

 

 

X =

Xi -\-vt,

y = X2.

(4.134)

Система координат x, у движется вместе с нагрузкой, причем точка х = 0 определяет положение разрыва дав­ ления. С учетом замены координат (4.134) уравнения движения (4.132) принимают вид

Іх-

д

дѵх (l,у)

, дѵу (l,у)

d l =

+

ді

ді

ду

 

 

рѵ*дѴх (х'у)- (4.135)

дх

И

X

 

д

дѵх (l,i/) 1

дѵу(1,у)

I

л>

+ (J.

"V M )l_ ду

[

dl

 

ду

J d l

 

 

 

 

=

рп2 дѴу (х’^ - t

(4.136)

 

 

 

 

 

 

дх

 

где ѵх и Ѵу — компоненты вектора

скорости, отнесенные

к системе координат х, у.

 

 

 

 

 

Соответствующие граничные

условия

теперь

будут

выглядеть так:

 

 

 

 

 

 

 

оѵу =

—ph (x)

при

у =

0

(4.137)

 

 

 

 

 

 

 

и

О х у =

0

при

у =

0.

 

 

 

В соответствии с подходом классической теории уп­ ругости уравнения движения являются несвязанными. Используя представление Гельмгольца в двух измере­ ниях, получаем

§ 4.5. Движущиеся нагрузки на полупространстве

191

ѵ

_

dq> (х,у)

дф (х,у)

 

 

х

дх

ду

 

и

 

 

 

(4.138)

X )

_

д ф ( * , у )

&ty(x,y)

 

у

 

ду

дх

 

Подставив (4.138) в уравнения движения (4.135) и (4.136), получим

 

 

 

X

 

 

 

 

 

 

О

 

 

 

+ 2(х

Ѵ2Ф (l,y ) d l

д

 

 

ду

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.139)

рѵ*

д<р

 

 

 

+

 

ду

дх )

 

 

 

дх

 

 

 

 

 

 

 

 

 

X

V24>(I.y) rfS

+

Ѵ2ФЫ

d %— ~

J

 

 

 

 

 

 

(4.140)

где

у 2 =

(д2/<?|2) +

(d 2f d y 2) .

Нижний

предел — оо

в (4.135) и (4.136) заменен на 0, так как предполагает­

ся, что разрыв давления движется

достаточно быстро

и ѵх— ѵу= 0

при xsg:0.

Кроме того, при переходе к фор­

мам (4.139)

и

(4.140)

использовались соотношения вида

- 5

 

 

 

 

і Г > ' <*•««.

U

 

 

І т

П

 

 

 

ü

 

 

а также зависимость ф(0, у) = 0 .

х,

 

Складывая

производную по

определяемую из

(4.139), с производной по у, определяемой из (4.140), находим,что

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

рѵ2

2p,

Ѵ2Ф (b y ) dl. (4.141)

дх

 

 

Подобным же образом из (4.139) и (4.140)

следует, что

р°2

^

1 fl (“ 7^") ѵ2 ^ ^

^4’142^

 

 

6

 

Уравнения (4.141)

и

(4.142) удобно решать с использо­

ванием умноженного на s преобразования Лапласа; при таком преобразовании имеем

f(s) = s \ f(x )e sxdx.

(4.143)

о

 

Умноженные на s преобразования Лапласа по перемен­

ной X от уравнений (4.141)

и

(4.142) имеют вид

 

(А, + 2р)

=

(рп2— А — 2р) s2<p

(4.144)

и

 

 

 

 

Р _р!_ =

(рп2

— |Г) s Y

(4.145)

В этих соотношениях умноженное на s преобразование Лапласа функции р(х/ѵ) дается формулой

% (х/ѵ)} = [I (sn) = р,

( 4 . 1 4 6 )

где p(s) — умноженное на s преобразование Лапласа по времени функции р (0 - Аналогичное (4.146) соотно­ шение существует и для А:

S? [ A ( * / ü) 1 = Ä ( ) = А.

(4 .1 4 7 1

При выводе уравнений (4.144) и (4.145) используется за­

висимость ф — ф= 0 при

х ^ О ,

совместная с условием

Ѵі = 0 при х^іО.

 

 

Преобразование Лапласа граничных условий (4.137)

дает

 

 

<Ууѵ($, У) =

— Р

при у = 0

и

 

(4.148)

§ 4.5. Движущиеся нагрузки на полупространстве

193

Oxu(s, у) = 0

при

у = О.

 

Напряжения, необходимые для подстановки в граничные условия, берутся из (2.2), (4.138) и соотношений между деформацией и перемещением (или между скоростями деформаций и скоростями). Преобразование этих напря­ жений имеет вид

З2 ф

3 ф

 

 

 

д2 ф

 

sv

~diF

ду )

 

 

 

Щ/*

(4.149)

и

 

 

 

 

 

 

 

аху

J L

2s

dq>

,

З2ф -S2l|)

 

 

 

SV

 

ду

 

Зу2

 

 

Решение уравнений (4.144) и (4.145), удовлетворяю­

щее граничным

условиям

(4.148),

определяется форму­

лами

 

 

 

 

 

 

 

 

Ф

рѵ

(t2 _ ,) 2 + 4 ^

esV fy

(4.150)

[iS

Ts

 

 

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

ф =

рѵ

Г ______ 2у/_________

 

(4.151)

 

 

1)2+

4У/ vs

 

 

(Ts

L

 

 

где

 

 

 

 

 

 

 

 

 

Y /=

[рц2/(Х +

2 ,Г ) - 1 ] 1/2

 

(4.152)

и

 

 

 

 

 

 

 

 

 

Ys = (ро2/ф— 1)1/2,

 

(4.153)

а две другие постоянные, входящие в решение уравне­ ний (4.144) и (4.145), принимаются равными нулю в си­ лу ограниченности решения при г/-ѵоо. Преобразованное решение для скорости, отвечающее (4.150) и (4.151), согласно (4.138), имеет вид

' (у2 —

1) е syfV — 2yf ys e'

Ejl

(4.154)

V

(Vs— !)2+ 4VfVs

 

13—851

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

и

Vf( Ys

e-Wfy + 2yf e ~ sV

(4.155)

 

И- L

- lY+ 4yys

Легко получается также преобразование решения для напряжений.

Применяя теорему о начальном значении, соответст­ вующую умноженному на s преобразованию Лапласа, имеем

Ііш р. (s) — р,(0) = р()

 

 

 

 

 

(4.156)

 

lirn^(s) = Х(0) =

Я0.

 

Примем следующие обозначения:

 

 

Cso =

Р ) '7 ’.

Cfo

K V

I ■ Spo)^]1/2

(4.157)

и

 

 

 

 

 

Yso =

[ ( « )

- i r /2,

Y/o = [ ( y2/c /o)

l ] l/2.

 

 

 

 

 

(4.158)

Теперь для достаточно больших s параметры реше­ ния приобретают простую форму, которая приведена ни­

же. Например, для ѵх имеем

ѵхÄi Сг e~ syf«y -f- Сa

 

) s >

1,

(4.159)

где С{, yfo и Yso не зависят от

s.

Используя

теорему

о сдвиге и теорему о начальном

значении,

находим, что

обращение двух членов в (4.159)

дается

выражениями

C x h ( x Y/oУ) I*—ѵ/0і/=0

и

C2h(x YsOУ) ljc-Ys04f = о-

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

х = уІ0 У и X = Yso У

§ 4.5. Движущиеся нагрузки на полупространстве

195

Эта ситуация изображена на рис. 4.3. Припоминая, что система координат х, у движется вместе с нагрузкой, мы замечаем, что эти две линии, х = у^у и х = увоУ, явля­ ются волновыми фронтами, первый из которых связан

Р и с . 4.3. Движущиеся нагрузки на границе полупространства.

с продольными возмущениями, а второй — с поперечны­ ми. Теперь мы видим, что если скорость не столь велика, чтобы радикалы в (4.158) были действительными, волны могут распространяться и впереди движущегося раз­ рыва давления. Этот тип задачи поддается анализу с значительно большими трудностями, чем только что рассмотренная задача о быстродвижущейся (сверхсей­ смической) нагрузке.

Чтобы сделать результаты более наглядными, разло­ жим функции релаксации в ряды Тейлора в окрестно­ стях точки ^=0; получим

ц ( 0

=

Цо — Ult +

( М 2/ 2 ) — ...

и

 

 

(4.160)

А ( 0

=

Ао h t +

(h t2/2) — ... .

13*

196

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

Умноженное на s преобразование Лапласа по времени этих разложений дает

p(s) = р0 — (ці/s) +

W s 2) — ...

и

(4.161)

Ms) — Яо— (U/s) +

(Яг/s2) — ••• •

Отсюда р и Я, согласно (4.146) и (4.147), теперь опреде­ ляются зависимостями

|1 = Ро —

(pi/sü) +

 

(б2 ls2V2) — ...

(4.162)

и

 

 

 

Я — Яо — (Яі/) -(-

(Яг/а2^2) — ... .

 

Подставляя (4.162)

в (4.153)

и используя биномиаль­

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

 

 

 

 

( 7so+ 0 ( M-i /Fq)

(4.163)

s Y s = 5 Yso

+

2Yso

 

 

 

где уso берется из (4.158), а члены, опущенные в (4.163), имеют более высокий порядок по 1/s. Используя (4.163), видим,что

lime sV = lime ^ W x p

(— ( ^so +

) | ^

(4.164)

S^oo

^ I

2VsoMo

J

Обозначим через [i\]s разрыв в vx при переходе через фронт волны сдвига x = y soy. Из (4.154), (4.163), (4.164) и из теоремы о начальном значении и теоремы сдвига следует, что

lvx]s _

 

c x p j ^ '0 + 1)

1X J .

(4.165)

0

Mo { (т?о — l ) 2 + 4Т/о Vso}

1 2АюИо

I

 

Аналогичным образом находятся разрывы ѵх при пере­ ходе через x = yj0y и ѵу при переходе через х = уs0p и х = =УюУ- Величины скачков даются формулами