книги из ГПНТБ / Кристенсен, Р. Введение в теорию вязкоупругости
.pdf§ 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р |
Ѵ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 0ф |
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 = Ä ( sü) = А. |
(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ц |
З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) |
и |
|
|
|
|
Я — Яо — (Яі/sü) -(- |
(Яг/а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 и х = =УюУ- Величины скачков даются формулами