Научная статья на тему 'Рассеяние наклонно падающей плоской звуковой волны упругим цилиндром с неоднородным покрытием'

Рассеяние наклонно падающей плоской звуковой волны упругим цилиндром с неоднородным покрытием Текст научной статьи по специальности «Физика»

CC BY
400
147
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
РАССЕЯНИЕ / ЗВУКОВЫЕ ВОЛНЫ / УПРУГИЙ ЦИЛИНДР / НЕОДНОРОДНОЕ УПРУГОЕ ПОКРЫТИЕ

Аннотация научной статьи по физике, автор научной работы — Толоконников Лев Алексеевич

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Рассеяние наклонно падающей плоской звуковой волны упругим цилиндром с неоднородным покрытием»

Известия Тульского государственного университета Естественные науки. 2013. Вып. 2. Ч.2. С. 265-274 Механика

УДК 539.3:534.26

Рассеяние наклонно падающей плоской звуковой волны упругим цилиндром

*

с неоднородным покрытием *

Л. А. Толоконников

Аннотация. Рассматривается задача о рассеянии звука упругим цилиндром с радиально-неоднородным упругим покрытием при наклонном падении плоской волны. Получено аналитическое выражение, описывающее рассеянное акустическое поле.

Ключевые слова: рассеяние, звуковые волны, упругий цилиндр, неоднородное упругое покрытие.

Существуют различные подходы к изменению звукоотражающих характеристик тел в определенных направлениях. Изменение характеристик рассеяния звука упругих тел можно осуществить с помощью специальных покрытий. Представляет интерес исследовать звукоотражающие свойства тел с покрытиями в виде непрерывно неоднородного упругого слоя. Такой слой легко реализовать с помощью системы тонких однородных упругих слоев с различными значениями механических параметров (плотности и упругих постоянных).

Рассеяние звуковых волн на однородных изотропных упругих цилиндрах рассматривалось в ряде работ. Например, случай нормального падения плоской волны исследовался в [1], а случай наклонного падения — в [2]. В работах [3, 4] решена задача о рассеянии плоских звуковых волн на неоднородном упругом цилиндре. Исследованию рассеяния звуковых волн трансверсально-изотропным неоднородным цилиндрическим слоем посвящена работа [5]. Дифракция плоской звуковой волны на неоднородном анизотропном полом цилиндре в общем случае анизотропии рассмотрена в [6]. Решение задачи о рассеянии плоской звуковой волны неоднородным упругим полым цилиндром в вязкой жидкости получено в [7]. В работе [8] изучена дифракция плоской звуковой волны на неоднородном термоупругом цилиндрическом слое, граничащем с невязкими теплопроводными жидкостями. Рассеяние звуковых волн абслютно жестким цилиндром с неоднородным упругим покрытием рассмотрено в [9].

* Работа выполнена при финансовой поддержке РФФИ (проект № 13-01-97514).

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

1. Постановка задачи. Рассмотрим бесконечный однородный упругий

цилиндр радиуса т2, материал которого характеризуется плотностью р2 и упругими постоянными А 2 и ^2. Цилиндр имеет покрытие в виде неоднородного изотропного упругого слоя, внешний радиус которого равен Т\. Цилиндрическая система координат r,ip,z выбрана таким образом, что координатная ось z является осью вращения цилиндра. Полагаем, что

модули упругости А и ц материала неоднородного цилиндрического слоя

описываются дифференцируемыми функциями цилиндрической радиальной координаты т, а плотность р — непрерывной функцией координаты т. Окружающая цилиндр жидкость является идеальной. Ее плотность и скорость звука соответственно равны pi и с.

Пусть из внешнего пространства на цилиндр произвольным образом падает плоская звуковая волна, потенциал скоростей которой равен

Фо = Ao exp{i[(k ■ r) - ut)]},

где Ao — амплитуда волны; k — волновой вектор падающей волны; r —

радиус-вектор; и — круговая частота. В дальнейшем временной множитель в-гш1 будем опускать.

В цилиндрической системе координат падающая волна запишется в виде

Ф0 = A0 exp{ik[r sin 90 cos (^ — <^0) + z cos 00]},

где 00 и ^>0 — полярный и азимутальный углы падения волны; к = и/с — волновое число во внешней области.

Определим отраженную от цилиндра волну, а также найдем поля смещений в упругом цилиндре и неоднородном слое.

2. Определение волновых полей. Распространение малых возмущений в идеальной жидкости в случае установившихся колебаний описывается уравнением Гельмгольца [10]

ДФ(1) + к2Ф(1) = 0,

где Ф(1) — потенциал скоростей полного акустического поля, Ф(1) = Ф0 + Ф«, Ф5 — потенциал скоростей рассеянной волны. При этом скорость частиц v и акустическое давление p в жидкости определяются по формулам: v = gradФ(1), p = ip1 иФ(1).

Потенциал скоростей падающей плоской волны представим в виде [11]

ГО

Ф0(т,^)= А0в^ injn(er)ein(v-v0),

П= — ГО

где Jn(x) — цилиндрическая функция Бесселя порядка n; а = kcos00; в = = к sin 00.

С учетом условий излучения на бесконечности функцию Ф8 будем искать в виде

ГО

Ъ8(т,<р,г) = ега ]Т А„Ип(вт)ет^о), (1)

п=-го

где Ип(х) — цилиндрическая функция Ханкеля первого рода порядка п.

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

Представляя вектор смещения и(2) частиц упругого изотропного однородного цилиндра в виде

и(2) = gradФ(2) + го1Ф,

где Ф(2) и Ф — скалярный и векторный потенциалы смещения, из уравнения Ламе получаем два волновых уравнения [12], которые для установившегося режима движения переходят в волновые уравнения Гельмгольца

дф(2) + к2ф(2) = о,

ДФ + к2 Ф = 0, (2)

где к = ш/в[ — волновое число продольных упругих волн; кТ = ш/ет — волновое число поперечных упругих волн; с = л/(А2 + 2^2)/р2 и ст = л/^2/р2

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

Векторное уравнение (2) в цилиндрической системе координат в общем случае не распадается на три независимых скалярных уравнения относительно проекций вектора Ф, а представляет собой систему трех уравнений, решение которой сопряжено со значительными математическими трудностями.

Представим вектор Ф в виде

Ф = го1(Ьех) + кт Мёг,

где Ь и М — скалярные функции пространственных координат г, р, г; ~ёх — единичный вектор оси г.

Тогда векторное уравнение (2) заменится двумя скалярными уравнениями Гельмгольца относительно функций Ь и М

ДЬ + к2 Ь = 0,

ДМ + к2 М = 0.

С учетом условия ограниченности функции Ф(2), Ь и М будем искать в виде

го

Ф(2)(г, р, г) = е*°“ ^2 БпЫЬг'У^-^, (3)

п=-го

го

Ь(г,р,г)= е™ ^ CnJn(k2г)ein(^-^0), (4)

п=-го

М(г, р, г) = е™ £ 7„(к2г)е

Іп()—) о)

(5)

где

кі = \/кг2 — а2; к2 = д/к2 — а2.

Компоненты вектора смещения и(2), записанные через функции Ф(2), Ь и М в цилиндрической системе координат, имеют вид

,.®

и.

(2) =

дФ(2)

дг

+

д2Ь

дгдг

кт дМ + —■

г др ’

д Ф(2) дг

(2) = 1 д Ф(2) др

дЬ Ч

1

+ -

д 2Ь

1 д_ г дг

гдг /

г дрдг

1 д2Ь

г2 др2

к

дМ

дг

дФ(2)

дг

д2Ь ,2г

+ д^2 + к^'

Соотношения между компонентами тензора напряжений стг(2) и вектора смещения и(2) в однородном изотропном упругом цилиндре записываются следующим образом [13]:

ст® = А2^уи(2) + 2^2

ди®

дг

ст® = А2ё1уи(2) +2^^1 +

(2)

иг

г др

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

ст® = А2ё1уи(2) + 2^2

ди®

дг

ст,

(2)

Г)

= ^2

1 ди® г др

+

ди

(2)

и

(2)

дг

ст

(2)

= ^2

'ди®

дг

+

диГ2)

дг

ст

(2)

= ^2

ди

(2)

дг

+ -

1 ди®

г др

где <Ь-и('2) = ^ + „И) + я .

дг г \ др / дг

Используя приведенные выше формулы и и учитывая, что divu(2) = = ДФ(2) = —кг2Ф(2) , выразим компоненты тензора напряжений ст®, ст® и ст^ через функции Ф(2), Ь и М. Получим

ди

(2)

ди®

ст® = — А2к2Ф(2) + 2^2 ^

д2Ф(2) д3Ь кт д2М

----------1----------+ —---------

дг2 дг2 дг г дгдр

кт дМ

г2 др

ст

(2) =

Г)

2^2 / д2Ф(2) г \ дгдр

1 дФ(2) д3 Ь

+

—^2кт

г др дгдрдг

/ д2М 1 дМ 1 д 2М Ч

дг2 г дг г2 др2

1 д2Ь Ч

г дрдг у

п= — гх

г

)

)

г

(2) _ (2 9 2 ф(2) 2 д3Ь к2 дЬ кт д2М \

СТгх д2 \ дгдг + дгдг2 + т дг + г дрдг у ’

Уравнения движения неоднородного изотропного упругого цилиндрического слоя в случае установившихся колебаний в цилиндрической системе координат имеют вид [13]

д@тт , 1 д^гш , дстгх , СТ22 2 / \

-------1---------------------1-----— _ -и2р(г)пг;

дг г др дг

д@Тф ! 1 дст 22 | дС(рг ! 2

ГГ + ГГ + ^ + _ ^ Р(' )п2;

дг г др дг г

^ ^ ^ _ -^2р(г)и^; (6)

г др дг г

дСТг2 . 1 дСТщх . дСТхх . 1

•'гг . 1 ди 2х хх 1 2 / \

+-------ГГ2- + -----+ - СТГх _ -^ р(г)и2

дг г др дг г

где пг, п2, пх — компоненты вектора смещения и частиц неоднородного слоя; — компоненты тензора напряжений в неоднородном слое.

Соотношения между компонентами тензора напряжений и вектора смещения и в неоднородном упругом цилиндрическом слое аналогичны соотношениям для однородного упругого цилиндра. Только упругие

постоянные Л2 и д2 следует заменить на функции А — А(г) и д — д(г).

Используя упомянутые соотношения, запишем уравнения (6) через компоненты вектора смещения и. Получим

/л ^ ч д2пг Л, „ , А + 2д \ дпг д д2пг . д2пг д2пг

(А + 2д)^г-2т + А' + 2д' +----------^ + А^-^ + д^1 +

дг2 \ г у дг г2 др2 дгдг дг2

А дпг А + д д2п2 А дп2 1 / ' А + 3д\ дп2

г дг г дгдр г дрдг г \ г у др

ч д2пх ,, дпх /А' А + 2д 2 і

+(А + д) ———+ А —-----------+ I----------2----+ ш р) пг _ 0

дгдг дг \ г г2

А + д д2пг 1 / . А + 3д\ дпг д2п2 А + 2д д2п2 / . д \ дп2

+ Г д + : ""я77 + + 12 + д + г “к;г +

г дгдр г г др дг2 г2 др2 г дг

д2п2 А + д д2пх / д' д 2 \ „

д^ТХ +------- 7ПТ + - ^ + ^Р К _ 0, (7)

дг2 г дрдг \ г г2 у

д2пг / ' А + д \ дпг А + д д2п2 д2пх д д2пх

дгдг + \д + г у дг + г дрдг +д дг2 + г2 др2

^д' +

(А + 2д)+ (д' + д) ^7 + ^Рпх _ 0

^1 / , и\ии'х1 ,2

г.

где штрих означает дифференцирование по г.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Так как неоднородность материала покрытия проявляется лишь в радиальном направлении, то зависимость компонент вектора смещения и от координаты г, согласно закону Снеллиуса, имеет вид егаг. Компоненты

вектора смещения u в неоднородном упругом слое являются периодическими функциями координаты р с периодом 2п. Поэтому функции ur(r, р, z), u2(r, р, z) и uz(r, р, z) будем искать в виде

го

Ur(r,p,z) = eiaz ^ um(r)ein(^o),

n=-ro

го

u2(r,p,z) = eiaz ]T U2„(r)ein(^-^0), (8)

n=-ro

ro

uz(r,p,z) = eiaz ]T u3„(r)em(^0).

n=-ro

Подставляя выражения (8) в уравнения (7), получим следующую систему линейных обыкновенных дифференциальных уравнений второго порядка относительно неизвестных функций u1n(r), u2n(r) и u3n(r) для каждого n:

AnUn + BnUn + CnUn = 0, (9)

где Un = (uin,u2n,u3n)T; An = (anij )3x3,Bn = (bnij )3x3,C<n = (Cnij )3x3

матрицы третьего порядка;

«nil = (A + 2^)r2, fln.22 = fln33 = дг2, an? =0 (i = j),

bnii = (A; + 2д')r2 + (A + 2д)г + iaA, bni2 = in(A + д)г, bni3 = ia(A + д)г2,

bn21 = in(A + д)г, bn22 = ^V2 + ДГ, bn23 = 0, bn3i = ia(A + д)г2, bn32 = 0, bn33 = д'г2 + дг,

cnii = A'r — (1 — iar)A — (n2 + 2 + а2г2)д + w2pr2,

Cni2 = in (A'r — A — 3д), Cni3 = iaA',

Cn2i = in (д'г + A + 3д),

Cn22 = —д'г — n2A — (2n2 + a2r2 + 1)д + w2pr2, c,n23 = —na(A + д)г,

Cn3i = ia^'r2 + (A + д)г], Cn32 = —na(A + д)г,

cn33 = —(n2 + 2а2г2)д — a2Ar2 + w2pr2.

Коэффициенты An, Bn, Cn, Dn разложений (1), (3)-(5) и функции uin(r),u2n(r),u3n(r) из разложений (8) подлежат определению из граничных условий.

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

и тангенциальные напряжения. Имеем г = Г1 : — шиг = -иг, стгг = —р, = 0, = 0;

(2) (2) (2) г = Г2 : и = иг ', % = ,

= (2) = (2) = (2)

(7"гг 0Уг ; .

Из условия равенства нормальных скоростей при г = Г1 находим коэффициенты Ап, выраженные через величины и1П(г1):

, Аоггав^П(вг1) + г^Шп(г1) (10)

П вНП(вг1) ’ ( )

где штрихи означают дифференцирование по аргументу.

Из условий непрерывности составляющих вектора смещения при Г = Г2 находим коэффициенты Вп, Сга и Дга, выраженные через величины и1П(г2), и2п(Г2) и иэ„(Г2):

q-S

Bn = [binuin(r2) + b2n U2n(—2) + b3nU3n(—2 )]q,-Cn = [cin Uin(—2) + C2„U2„(r2) + C3nU3n(—2)] q-1, (11)

D„ = [dinuin(f2) + d2nU2„(r2) + d3nU3n(—2)] q—,

где

bin = fcrk2(k? - a2)—|Jn(k2—2) Jn(k2—2); b2n = mkr(k^ - a2)—2Jn(k2—2);

Ьз„ = iakr[n2Jn(k2—2) - (k2r2)2Jn2(k2^2)]; cin = —iafcr k2 —Jn(ki—2) J^k2—2); C2n = ankr r2Jn(kir2) Jn(k2 —2)

C3n = [1 - iaJn(kir2)b3n]/[(k2 - a2) Jn(k2^2)];

din = mfc^fc,2 - a2)—2Jn(ki—2) Jfe^); d2n = (k2 - a2)r2 Jn(k2r2)[ki(k2 - a2) J^ki—2) Jn(k2—2) + a2k2 Jn(ki—2) J^(k2—2)]; d3n = an(k2 - a2)—2 Jn(k2—2)[k2 Jn(ki—2) Jn(k2^2) - ki Jn(ki—2) Jn(k2—2)]; qn = krk2r|jn(k2r2)[ki(k2 - a2)jn(kir2)Jn(k2—2) + a% Jn(ki—2) J^(k2—2)] -

-n2k3 Jn(kir2)Jn(k2r2).

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

—2 Anun + KUn) = (ini

r 1 Л Л \ri (12)

-9 Anun + —nUn) = 0|

—2 n^n

' / Г=Г2

где элементы матриц ' = (enj )3х3, —n = (/mj )3х3, Gn = (/nj )3xi

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

Л w2piHn(e—i) in—

enii — + 0,r. (n \ ; eni2 — ; eni3 — iaA;

—i eHn(e—i) —

en2i = —^ ; en22 = - — ; en3i = iaW en23 = en32 = en33 = 0;

.nii = — (7inbin + Y2n cin + Y3n din)qn ;

I in— / , ,\-1

/ni2 =-------(7in02n + Y2nC2n + 73n«2n)qn ;

,/ni3 = iaA - (Yinb3n + Y2nc3n + Y3nd3

n qn- i ;

,/n2i — — (Y4nbin + Y5ncin + Y6ndin)qn ;

,/n22 = - - (74nb2n + 75nc2n + 76nd2n)qn i;

»/*n23 = - (74nb3n + 75nc3n + 76nd3n)qn i; ,/n3i = ia^ - (Y7nbin + Y8ncin + Y9ndin)qn i; /n32 = - (Y7nb2n + Y8nc2n + Y9nd2n

n -i

/'n33 = - (77nb3n + 78nc3n + 79nd3^ " i'

3n n

Здесь

_ _ 2inAowpi _ _ _ _

/ni = пв—i Hn(e—i); gn2 = /n3 = 0'

Yin = 2д2^П(ki— 2) - —2k2Jn(ki—); Y2n = 2i^2ak|J^(k2—2);

2i 2i

Y3n = -0 ^2nkr [k2—2Jn(k2—2) - Jn(k2—2)]; Y4n = “o №n[ki—2 Jn(ki— 2) - Jn(ki—2)]; —2 —2

2

Y5n = -“o^2na[k2—2Jn(k2—2) - Jn(k2 —2)];

—22

Y6n = Л^2kr[-(k2— 2)2Jn(k2—2) + k2—2Jn(k2—2) - n°Jn(k2—2)];

—22

Y7n = 2i^2akiJn(ki—2); Y8n = ^2k2Jn(k2—2)(2a2 + k°);

Y9n = -—^2kr naJn(k2—2)'

—2

Краевая задача (9), (12) может быть решена различными методами, например, методами, приведенными в [14].

После решения краевой задачи (9), (12) по формулам (10) и (11) вычисляются коэффициенты An, Bn, Cn и Dn.

В результате получаем аналитическое описание волновых полей вне и внутри упругого тела с помощью выражений (1) и (3)—(5).

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

Используя асимптотическую формулу для цилиндрической функции Ханкеля первого рода при больших значениях аргумента [15] (вг ^ 1)

Hn(er)

2

пвг

exp

ПП П Ч вг - -у - -

из (1) находим

= а/2Г exP

Maz і вг — П

F

где

Vnkrl

^ (-in)An exp[m(p - po)j.

С помощью выражения для амплитуды рассеяния в дальней зоне поля | ^(р) | строятся диаграммы направленности и частотные характеристики рассеянного поля, позволяющие изучить звукоотражающие свойства тела в различных направлениях и резонансное рассеяние звука.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

2

n=—ro

Список литературы

1. Faran J.J. Sound scattering by solid cylinders and spheres // Acoust. Soc. Amer. 1951. V. 23. № 4. P. 405-420.

2. Flax L., Varadan V.K., Varadan V.V. Scattering of an obliquely incident acoustic wave by an infinite cylinder // J. Acoust. Soc. Amer. 1980. V. 68. N 6. P. 1832-1835.

3. Безруков А.В., Приходько В.Ю., Тютекин В.В. Рассеяние звуковых волн упругими радиально-слоистыми цилиндрическими телами // Акуст. журн. 1986. Т. 32. № 6. С. 762-766.

4. Коваленко Г.П. К задаче о дифракции акустической волны на неоднородном твердом теле // Акуст. журн. 1987. Т. 33. № 6. С. 1060-1063.

5. Скобельцын С.А., Толоконников Л.А. Рассеяние звуковых волн трансверсаль-но-изотропным неоднородным цилиндрическим слоем // Акуст. журн. 1995. Т. 41. № 1. С. 134-138.

6. Толоконников Л.А. Дифракция звуковых волн на неоднородном анизотропном полом цилиндре // Оборонная техника. 1998. № 4-5. С. 11-14.

7. Романов А.Г., Толоконников Л.А. Рассеяние плоской звуковой волны неоднородным упругим полым цилиндром в вязкой жидкости // Изв. ТулГУ. Естественные науки. 2009. Вып. 1. С. 63-72.

8. Ларин Н.В., Толоконников Л.А. Дифракция плоской звуковой волны на неоднородном термоупругом цилиндрическом слое, граничащем с невязкими теплопроводными жидкостями // Прикладная математика и механика. 2009. Т. 73. Вып. 3. С. 474-483.

9. Романов А.Г., Толоконников Л.А. Рассеяние звуковых волн цилиндром с неоднородным упругим покрытием // Прикладная математика и механика. 2011. Т. 75. Вып. 5. С. 850-857.

10. Шендеров Е.Л. Волновые задачи гидроакустики. Л.: Судостроение, 1972. 352 с.

11. Иванов Е.А. Дифракция электромагнитных волн на двух телах. Минск: Наука и техника. 1968. 584 с.

12. Седов Л.И. Механика сплошной среды. Т. 2. М.: Наука, 1973. 584 с.

13. Новацкий В. Теория упругости. Т. 2. М.: Мир, 1975. 872 с.

14. Толоконников Л.А., Ларин Н.В. Рассеяние звука неоднородными термоупругими телами. Тула: Изд-во ТулГУ, 2008. 232 с.

15. Справочник по специальным функциям / Под ред. Абрамовица М., Стигана И.

- М.: Наука, 1979. 832 с.

Толоконников Лев Алексеевич (tolokonnikovla@mail.ru), д.ф.-м.н., профессор, кафедра прикладной математики и информатики, Тульский государственный университет.

Scattering of an obliquely incident plane sound wave by an elastic cylinder with a non-uniform covering

L.A. Tolokonnikov

Abstract. The problem of the scattering of an obliquely incident plane sound wave by an elastic cylinder with a radially non-uniform elastic coating is considered. An analytical expression describing the scattered acoustic field is obtained.

Keywords: scattering, sound waves, elastic cylinder, non-uniform elastic coating.

Tolokonnikov Lev (tolokonnikovla@mail.ru), doctor of physical and mathematical sciences, professor, department of applied mathematics and computer science, Tula State University.

Поступила 20.03.2013

i Надоели баннеры? Вы всегда можете отключить рекламу.