Влияние начальных напряжений на квазилэмбовские моды в гидроупругих волноводах

В данной работе для исследования распространения квазилэмбовских мод в гид роупругих волноводах, состоящих из полупространства жидкости и упругого слоя, а также из упругого полупространства и жидкого слоя, привлекаются модели предварительно напряженного тела и идеальной сжимаемой жидкости. При этом...

Full description

Saved in:
Bibliographic Details
Published in:Прикладная механика
Date:2020
Main Authors: Гузь, А.Н., Багно, А.М.
Format: Article
Language:Russian
Published: Інститут механіки ім. С.П. Тимошенка НАН України 2020
Online Access:https://nasplib.isofts.kiev.ua/handle/123456789/188211
Tags: Add Tag
No Tags, Be the first to tag this record!
Journal Title:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Cite this:Влияние начальных напряжений на квазилэмбовские моды в гидроупругих волноводах / А.Н. Гузь, А.М. Багно // Прикладная механика. — 2020. — Т. 56, № 1. — С. 4-22. — Бібліогр.: 35 назв. — рос.

Institution

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id nasplib_isofts_kiev_ua-123456789-188211
record_format dspace
spelling Гузь, А.Н.
Багно, А.М.
2023-02-16T17:04:45Z
2023-02-16T17:04:45Z
2020
Влияние начальных напряжений на квазилэмбовские моды в гидроупругих волноводах / А.Н. Гузь, А.М. Багно // Прикладная механика. — 2020. — Т. 56, № 1. — С. 4-22. — Бібліогр.: 35 назв. — рос.
0032-8243
https://nasplib.isofts.kiev.ua/handle/123456789/188211
В данной работе для исследования распространения квазилэмбовских мод в гид роупругих волноводах, состоящих из полупространства жидкости и упругого слоя, а также из упругого полупространства и жидкого слоя, привлекаются модели предварительно напряженного тела и идеальной сжимаемой жидкости. При этом используются трехмерные линеаризованные уравнения теории упругости конечных деформаций для упругого тела и трехмерные линеаризованные уравнения Эйлера для жидкости. Предполагается, что жидкость находится в состоянии покоя. В качестве подхода выбраны постановки задач и метод, основанные на применении представлений общих решений уравнений движения идеальной сжимаемой жидкости и предварительно напряженного сжимаемого упругого тела.
Досліджено поширення акустичних хвиль у попередньо деформованому стисливому пружному шарі, що взаємодіє з шаром в'язкої стисливої рідини. Результати одержано на основі тривимірних рівнянь лінеаризованої теорії пружності скінченних деформацій для пружного шару та тривимірних лінеаризованих рівнянь Нав'є - Стокса для в'язкої стисливої рідини. Застосовано постановку задачі та підхід, які базуються на використанні представлень загальних розв'язків лінеаризованих рівнянь для пружного та рідкого шарів. Одержано дисперсійне рівняння, яке описує поширення гармонічних хвиль у гідропружній системі у широкому діапазоні частот для товстого рідкого шару. Проаналізовано вплив початкових напружень, а також товщини шарів пружного тіла і рідини на фазові швидкості та коефіцієнти загасання мод Лемба. Показано, що в'язкість рідини призводить до зменшення глибини проникнення нижчої квазілембової моди в рідкий шар.
The problem of quasi-Lamb waves propagation in a pre-deformed compressible elastic layer that interacts with a half-space of an ideal compressible fluid and in a predeformed compressible elastic half-space that interacts with a layer of an ideal compressible fluids is studied. The study is carried out basing on the three-dimensional linearized equations of the theory of elasticity of finite deformations for the compressible elastic solid and the three-dimensional linearized Euler equations for the ideal compressible fluid. A problem statement and an approach, based on the representations of general solutions of the linearized equations for elastic solid and fluid are used. The dispersion equations, which describe a propagation of quasi-Lamb waves in the hydroelastic systems are obtained in the wide frequency range. An effect of the initial stresses as well as the thickness of elastic layer and layer of ideal compressible fluid on the phase velocities of quasi-Lamb modes are analyzed. A criterion of existence of the quasi-Lamb waves in the hydroelastic waveguides is proposed. The developed approach and the findings make it possible to establish the limits of applicability of the models of wave processes, based on different versions of the theory of small initial deformations, as well on the classical theory of elasticity. The numerical results are presented in the form of graphs, and their analysis is given.
ru
Інститут механіки ім. С.П. Тимошенка НАН України
Прикладная механика
Влияние начальных напряжений на квазилэмбовские моды в гидроупругих волноводах
Effect of Initial Stresses on Quasi-Lamb Modes in Hydroelastic Waveguides
Article
published earlier
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
title Влияние начальных напряжений на квазилэмбовские моды в гидроупругих волноводах
spellingShingle Влияние начальных напряжений на квазилэмбовские моды в гидроупругих волноводах
Гузь, А.Н.
Багно, А.М.
title_short Влияние начальных напряжений на квазилэмбовские моды в гидроупругих волноводах
title_full Влияние начальных напряжений на квазилэмбовские моды в гидроупругих волноводах
title_fullStr Влияние начальных напряжений на квазилэмбовские моды в гидроупругих волноводах
title_full_unstemmed Влияние начальных напряжений на квазилэмбовские моды в гидроупругих волноводах
title_sort влияние начальных напряжений на квазилэмбовские моды в гидроупругих волноводах
author Гузь, А.Н.
Багно, А.М.
author_facet Гузь, А.Н.
Багно, А.М.
publishDate 2020
language Russian
container_title Прикладная механика
publisher Інститут механіки ім. С.П. Тимошенка НАН України
format Article
title_alt Effect of Initial Stresses on Quasi-Lamb Modes in Hydroelastic Waveguides
description В данной работе для исследования распространения квазилэмбовских мод в гид роупругих волноводах, состоящих из полупространства жидкости и упругого слоя, а также из упругого полупространства и жидкого слоя, привлекаются модели предварительно напряженного тела и идеальной сжимаемой жидкости. При этом используются трехмерные линеаризованные уравнения теории упругости конечных деформаций для упругого тела и трехмерные линеаризованные уравнения Эйлера для жидкости. Предполагается, что жидкость находится в состоянии покоя. В качестве подхода выбраны постановки задач и метод, основанные на применении представлений общих решений уравнений движения идеальной сжимаемой жидкости и предварительно напряженного сжимаемого упругого тела. Досліджено поширення акустичних хвиль у попередньо деформованому стисливому пружному шарі, що взаємодіє з шаром в'язкої стисливої рідини. Результати одержано на основі тривимірних рівнянь лінеаризованої теорії пружності скінченних деформацій для пружного шару та тривимірних лінеаризованих рівнянь Нав'є - Стокса для в'язкої стисливої рідини. Застосовано постановку задачі та підхід, які базуються на використанні представлень загальних розв'язків лінеаризованих рівнянь для пружного та рідкого шарів. Одержано дисперсійне рівняння, яке описує поширення гармонічних хвиль у гідропружній системі у широкому діапазоні частот для товстого рідкого шару. Проаналізовано вплив початкових напружень, а також товщини шарів пружного тіла і рідини на фазові швидкості та коефіцієнти загасання мод Лемба. Показано, що в'язкість рідини призводить до зменшення глибини проникнення нижчої квазілембової моди в рідкий шар. The problem of quasi-Lamb waves propagation in a pre-deformed compressible elastic layer that interacts with a half-space of an ideal compressible fluid and in a predeformed compressible elastic half-space that interacts with a layer of an ideal compressible fluids is studied. The study is carried out basing on the three-dimensional linearized equations of the theory of elasticity of finite deformations for the compressible elastic solid and the three-dimensional linearized Euler equations for the ideal compressible fluid. A problem statement and an approach, based on the representations of general solutions of the linearized equations for elastic solid and fluid are used. The dispersion equations, which describe a propagation of quasi-Lamb waves in the hydroelastic systems are obtained in the wide frequency range. An effect of the initial stresses as well as the thickness of elastic layer and layer of ideal compressible fluid on the phase velocities of quasi-Lamb modes are analyzed. A criterion of existence of the quasi-Lamb waves in the hydroelastic waveguides is proposed. The developed approach and the findings make it possible to establish the limits of applicability of the models of wave processes, based on different versions of the theory of small initial deformations, as well on the classical theory of elasticity. The numerical results are presented in the form of graphs, and their analysis is given.
issn 0032-8243
url https://nasplib.isofts.kiev.ua/handle/123456789/188211
citation_txt Влияние начальных напряжений на квазилэмбовские моды в гидроупругих волноводах / А.Н. Гузь, А.М. Багно // Прикладная механика. — 2020. — Т. 56, № 1. — С. 4-22. — Бібліогр.: 35 назв. — рос.
work_keys_str_mv AT guzʹan vliânienačalʹnyhnaprâženiinakvazilémbovskiemodyvgidrouprugihvolnovodah
AT bagnoam vliânienačalʹnyhnaprâženiinakvazilémbovskiemodyvgidrouprugihvolnovodah
AT guzʹan effectofinitialstressesonquasilambmodesinhydroelasticwaveguides
AT bagnoam effectofinitialstressesonquasilambmodesinhydroelasticwaveguides
first_indexed 2025-11-27T00:52:30Z
last_indexed 2025-11-27T00:52:30Z
_version_ 1850789578751868928
fulltext 2020 П Р И К Л А Д Н А Я М Е Х А Н И К А Том 56, № 1 4 ISSN0032–8243. Прикл. механика, 2020, 56, № 1 А . Н . Г у з ь , А . М . Б а г н о ВЛИЯНИЕ НАЧАЛЬНЫХ НАПРЯЖЕНИЙ НА КВАЗИЛЭМБОВСКИЕ МОДЫ В ГИДРОУПРУГИХ ВОЛНОВОДАХ Институт механики им. С.П.Тимошенко НАНУ, ул Нестерова, 3, 03057, Киев, Украина; e-mail: alexbag2016@gmail.com Abstract. The problem of quasi-Lamb waves propagation in a pre-deformed compress- ible elastic layer that interacts with a half-space of an ideal compressible fluid and in a pre- deformed compressible elastic half-space that interacts with a layer of an ideal compressible fluids is studied. The study is carried out basing on the three-dimensional linearized equa- tions of the theory of elasticity of finite deformations for the compressible elastic solid and the three-dimensional linearized Euler equations for the ideal compressible fluid. A problem statement and an approach, based on the representations of general solutions of the linear- ized equations for elastic solid and fluid are used. The dispersion equations, which describe a propagation of quasi-Lamb waves in the hydroelastic systems are obtained in the wide frequency range. An effect of the initial stresses as well as the thickness of elastic layer and layer of ideal compressible fluid on the phase velocities of quasi-Lamb modes are analyzed. A criterion of existence of the quasi-Lamb waves in the hydroelastic waveguides is pro- posed. The developed approach and the findings make it possible to establish the limits of applicability of the models of wave processes, based on different versions of the theory of small initial deformations, as well on the classical theory of elasticity. The numerical results are presented in the form of graphs, and their analysis is given. Key words: compressible elastic layer, compressible elastic half-space, half-space of ideal compressible fluid, layer of ideal compressible fluid, initial stresses, quasi-Lamb waves. Введение. Закономерности распространения поверхностных акустических волн широко ис- пользуются в сейсмологии, сейсморазведке, гидроакустике, а также при разработке ультразвуковых неразрушающих методов выявления дефектов и определения напря- жений в материалах и элементах конструкций [1, 16, 17, 21, 22, 27 – 31, 33 – 35]. Зна- чительное практическое использование поверхностных волн ставит задачу более пол- ного учета свойств реальных упругих сред. К числу таких факторов принадлежат начальные напряжения в твердом теле. Созданные целенаправленно или возникшие в результате технологических операций при изготовлении, они оказывают существен- ное влияние на волновые процессы. Рассмотренные задачи и результаты, полученные с учетом в упругих телах начальных напряжений, приведены в [3, 4, 6 – 8, 10 – 12, 14, 15, 18, 19, 23 – 26, 32]. В данной работе для исследования распространения квазилэмбовских мод в гид- роупругих волноводах, состоящих из полупространства жидкости и упругого слоя, а также из упругого полупространства и жидкого слоя, привлекаются модели предвари- тельно напряженного тела и идеальной сжимаемой жидкости. При этом используются трехмерные линеаризованные уравнения теории упругости конечных деформаций для упругого тела и трехмерные линеаризованные уравнения Эйлера для жидкости. Пред- полагается, что жидкость находится в состоянии покоя. В качестве подхода выбраны постановки задач и метод, основанные на применении представлений общих решений 5 уравнений движения идеальной сжимаемой жидкости и предварительно напряженно- го сжимаемого упругого тела, предложенные в работах [3 – 9, 18 – 20]. §1. Постановка задачи. Основные уравнения. Далее рассматриваем такие динамические процессы в гидроупругой системе, при которых возникающие дополнительные деформации, т.е. возмущения деформаций, значительно меньше начальных. Исследуем гармонические волновые процессы малой амплитуды. При этом примем, что упругое тело находится в начальном состоянии. Заметим, что в отличие от упругих тел, соотношения которых записаны в лагранже- вых координатах, равенства для жидкой среды записываются в эйлеровых координа- тах, введенных в естественном состоянии жидкости. Следует подчеркнуть, что начальное состояние упругого тела при рассмотрении гидроупругой задачи является естественным состоянием по отношению к жидкости и системе в целом. Поскольку в дальнейшем исследуем распространение малых возмущений, то, как известно, в этом случае подходы Эйлера и Лагранжа в описании поведения сред совпадают. Поэтому в дальнейшем не делаем различий между лагранжевыми и эйлеровыми координатами и характерные для нелинейных задач трудности при записи граничных условий при указанных двух подходах не возникают. Далее предположим, что изотропное нелинейно-упругое твердое тело, упругий потенциал которого является произвольной дважды непрерывно-дифференцируемой функцией компонент тензора деформаций Грина, занимает объем ( 1 ,z    2 0,h z   3z    ) и контактирует с полупространством идеальной сжимаемой жидкости, заполняющей объем: ( 1 ,z    20 ,z   3z    ). В случае вто- рой гидроупругой системы жидкий слой заполняет объем ( 1 ,z    20 ,z h  3z    ) и взаимодействует с упругим полупространством, занимающим объем ( 1 ,z    2 0,z   3z    ). Примем, что внешние силы, действующие на указанные среды, распределены равномерно вдоль оси 3.Oz В этом случае задача яв- ляется плоской и можно ограничиться изучением процесса распространения возму- щений в плоскости 1 2.Oz z Для упругого слоя, взаимодействующего с полупространством жидкости, задача сводится к решению системы уравнений движения упругого тела и жидкости при сле- дующих динамических и кинематических граничных условиях: 21 0 0;zQ   2 22 0 2 0 ;z zQ P   21 0;z hQ   22 0;z hQ   2 2 2 2 0 0 .z z u v t     (1.1) В случае упругого полупространства, контактирующего со слоем жидкости, зада- ча сводится к решению системы уравнений гидроупругости при таких граничных условиях: 21 0 0;zQ   2 22 0 2 0 ;z zQ P   22 0;z hP   2 2 2 2 0 0 .z z u v t     (1.2) Здесь iQ и iP – составляющие напряжений, соответственно, в упругом теле и жидкости. §2. Методика решения. Воспользуемся постановками задач гидроупругости для тел с начальными напря- жениями и идеальной жидкости, а также представлениями общих решений, предло- женными в работах [3 – 9, 18 – 20]. В дальнейшем исследуем волновые процессы в предварительно деформированных сжимаемых упругих телах, взаимодействующих с идеальной жидкостью, начальное состояние которых является однородным. В случае однородного напряженно-деформированного состояния для плоского случая общие решения имеют вид [3 – 9, 18 – 20]: 1) для упругого тела из сжимаемого материала – 6 2 1 1 1 2 ;u z z                2 0 2 2 02 2 2 1 11 11 2 1 12 22 2 12 2 2 22 2 0 2 2 0 2 12 12 1 21 1 11 11 1 1 11 11 ; a s s u a z z ta s a s                           2) для идеальной сжимаемой жидкости – 2 2 1 1 ; z t      2 2 2 2 ,v z t     где введенные функции i являются решениями следующих уравнений: 1) для упругого тела из сжимаемого материала –           2 2 0 2 2 02 2 2 2 2 2 1 12 22 2 2 22 22 2 2 2 2 22 2 0 2 2 0 2 2 0 1 2 1 21 1 11 11 1 1 11 11 1 2 12 11 s a s z z t z za s a s s                                          242 4 2 12 12 12 2 22 2 0 2 0 2 0 1 21 1 12 11 1 11 11 2 12 11 0; a t z zs a s s                     2) для идеальной сжимаемой жидкости – 2 2 2 22 2 2 2 1 2 0 1 0. z z a t                Здесь и выше введены такие обозначения: iu – компоненты вектора смещений упру- гого тела u ;  – плотность материала упругого тела; ija и ij – величины, определя- емые из уравнений состояния и зависящие от вида упругого потенциала [3, 4, 6 – 8, 10]; і – удлинения упругого тела в направлениях координатных осей; 0 ii 0 0 2 1 2 3( ( ) )іі ii is      – начальные напряжения; iv – составляющие вектора возмуще- ний скорости жидкости v относительно состояния покоя; 0 и 0a – плотность и ско- рость звука в жидкости в состоянии покоя. Заметим, что для упругих тел приведенные соотношения относительно возмуще- ний являются линейными, однако величины начального состояния, входящие в них, определяются из общих нелинейных уравнений. В связи с этим, несмотря на то, что основные соотношения приведены в координатах начального деформированного со- стояния iz и все величины отнесены к размерам тела в этом состоянии, общая поста- новка задач гидроупругости для тел с начальными напряжениями в координатах iz по форме аналогична формулировке линейных задач классической теории гидроупруго- сти; есть и существенные различия, которые относятся к структуре уравнений и гра- ничных условий. Для анализа распространения возмущений, гармонически изменяющихся во вре- мени, решения системы уравнений определяем в классе бегущих волн 2( )j jX z    1exp ( )i k z t  ( 1, 2),j  где k – волновое число;  – круговая частота; i – мни- мая единица ( 1i ). Заметим, что выбранный в данной работе класс гармонических волн, являясь наиболее простым и удобным в теоретических исследованиях, не ограничивает общ- ности полученных результатов, поскольку линейная волна произвольной формы, как известно, может быть представлена набором гармонических составляющих. Далее применяя метод Фурье, приходим к двум задачам о собственных значениях для урав- нений движения упругого тела и жидкости. Решая их, определяем соответствующие 7 собственные функции. После подстановки полученных общих решений в граничные условия (1.1) и (1.2) получаем однородные системы линейных алгебраических урав- нений относительно произвольных постоянных. Исходя из условия существования нетривиальных решений этих систем, получаем дисперсионные уравнения. Для упру- го-жидкостной системы, состоящей из упругого слоя и жидкого полупространства, дисперсионное соотношение имеет вид  0 0 0det , , , , , , , , , , 0lm i ij ij ii sc a s а h c         ( , 1, 5).l m  (2.1) Для упруго-жидкостной системы, жидкий слой которой взаимодействует с упру- гим полупространством, аналогичное уравнение будет таким:  0 0 0det , , , , , , , , , , 0lm i ij ij ii sc a s а h c         ( , 1, 4),l m  (2.2) где с – фазовая скорость нормальных волн в гидроупругой системе с предварительно напряженным упругим телом; h – толщина упругого слоя для первой системы или толщина слоя жидкости в случае второго гидроупругого волновода; sc ( 2 sc   ) – скорость волны сдвига в ненапряженном упругом теле;  – модуль сдвига материала упругого тела. Как известно, в неограниченном сжимаемом упругом теле существуют продоль- ная и сдвиговая волны. В идеальной сжимаемой жидкой среде распространяется толь- ко продольная волна. Именно эти волны, взаимодействуя между собой на свободных граничных поверхностях, а также на поверхностях контакта сред, порождают слож- ные волновые поля в гидроупругих системах. Отметим, что дисперсионные уравнения (2.1) и (2.2) не зависят от формы упруго- го потенциала и получены для сжимаемых упругих тел, подверженных большим (ко- нечным) начальным деформациям. Они являются наиболее общими и из них можно получить соотношения для ряда частных случаев [1, 4, 5, 10, 11, 13 – 15, 26]. Если принять 0 0,ii  то получим равенства для основательно исследованных в рамках классической теории упругости волн Рэлея, Стоунли – Шольте и Лэмба [1]. §3. Числовые результаты и их анализ. В дальнейшем дисперсионные уравнения (2.1) и (2.2) решаем численно. При этом расчеты проводим для двух гидроупругих систем, состоящих из реальных упругих тел и жидкости, а также для четырех систем с гипотетическими жидкостями. Первая сис- тема состоит из органического стекла и воды. Она характеризуется следующими па- раметрами: упругое тело – 31160 кг/м ,  93,96 10 Па,   91,86 10 Па;   жидкость – 3 0 1000 кг/м ,  0 1459,5 м/с,а  0 0 1,152595.sa a c  У этого волновода материал упругого тела (оргстекло) является жестким. Вторая система представляет собой волновод из стали марки 09Г2С и воды. При этом параметры выбираем такими: упругое тело – 37800 кг/м ,  109,26 10 Па,   107,75 10 Па;   жидкость – 3 0 1000 кг/м ,  0 1459,5 м/с,а  0 0,463021a  . Этот волновод отличается тем, что материал упругого тела (сталь) относится к разряду более жестких, чем оргстекло. Гипотетические жид- кости выбирали со скоростью распространения звуковой волны 0 0,9a  и 0 0,97a  при взаимодействии с оргстеклом, а также со скоростью волны звука 0 0,97a  и 0 1,05a  для пары сталь – вода. Такой выбор позволяет проиллюстрировать справед- ливость предложенных на основании аналитического анализа дисперсионных уравне- ний и результатов численного их решения критериев расщепления единственной по- верхностной волны на множество мод, распространяющихся в гидроупругих системах. Заметим, что уравнения (2.1) и (2.2) выведены без введения каких-либо дополни- тельных ограничений к виду функции упругого потенциала, поэтому они справедли- вы для упругих потенциалов произвольной формы. В данной работе при численном 8 решении этих уравнений для описания упругих свойств оргстекла и стали использо- вался потенциал Мурнагана, зависящий от трех алгебраических инвариантов тензора деформаций Грина [3 – 8, 10, 19, 22]. При рассмотрении конкретных примеров и чис- ленного решения уравнений (2.1) и (2.2) учитывалось то обстоятельство, что оргстек- ло и сталь, не разрушаясь, не допускают больших деформаций и поэтому коэффици- енты уравнений состояния ija и ij определены в рамках линейного акустического приближения [3 – 8, 10, 19, 22]. Результаты вычислений в виде графиков представлены на рис. 1 – 19. На рис. 1 – 10 приведены результаты численных расчетов для упруго-жидкостной системы, состоящей из упругого слоя и жидкого полупространства. Графики, полу- ченные для гидроупругого волновода, компонентами которого являются жидкий слой и упругое полупространство, представлены на рис. 11 – 19. В частности, на рис. 1 – 3 представлены графики для упругого слоя из органического стекла (жесткий материал) и воды ( 0 0 1,152595sa a c  ), на рис. 4 и 5 – для оргстекла и гипотетической жидко- сти, скорость распространения звуковой волны в которой равна 0 0,9a  . На рис. 6 – 8 приведены графики для упругого слоя из стали марки 09Г2С (более жесткий матери- ал) и воды, а на рис. 9 – 10 – для стали и гипотетической жидкости со скоростью звука в ней 0 0,97a  . На рис. 1, 6 для упругого слоя, невзаимодействующего с жидкостью, приведены зависимости безразмерных величин фазовых скоростей нормальных волн Лэмба c ( sc c c ) от безразмерной величины толщины упругого слоя (частоты) h ( schh  ) [3 – 8, 10] при отсутствии начальных напряжений ( 0 11 0  ). На этих ри- сунках для наглядности штриховыми линиями отмечены асимптотики, к которым стремятся фазовые скорости первой (нулевой антисимметричной) и второй (нулевой симметричной) мод при возрастании толщины упругого слоя (частоты) h . Из графиков, полученных для чисто упругого слоя из оргстекла [3 – 8, 10] и пред- ставленных на рис. 1, следует, что скорость первой (нулевой антисимметричной) мо- ды Лэмба при росте частоты или толщины упругого слоя h стремится к скорости волны Рэлея Rc ( 0,9335596R R sc c c  ) снизу, а скорость второй (нулевой симмет- ричной) моды стремится к скорости волны Рэлея Rc ( 0,9335596Rc  ) сверху. Скорости всех мод Лэмба высокого порядка при увеличении толщины упругого слоя (частоты) h стремятся к скорости волны сдвига в материале упругого тела sc [1, 3 – 8, 10]. Рис. 1 9 На рис. 2, 4, 7 и 9 изображены дисперсионные кривые для гидроупругих волново- дов, отражающие зависимости безразмерных величин фазовых скоростей квазилэм- бовских мод c от безразмерной величины толщины упругого слоя (частоты) h при отсутствии начальных напряжений. На этих рисунках для наглядности штриховыми линиями отмечены асимптотики, к которым стремятся фазовые скорости первой и второй мод при возрастании частоты (толщины упругого слоя) h . Рис. 2 Графики для гидроупругой системы, приведенные на рис. 2, показывают, что при росте толщины упругого слоя (частоты) h скорость первой моды стремится к скоро- сти волны Стоунли stc ( 0,7717101st st sc c c  ) снизу, а скорость второй моды – к скорости волны Рэлея Rc ( 0,933558Rc  ) сверху. Моды высокого порядка распро- страняются в упругом слое в его толще с фазовыми скоростями, стремящимися с ро- стом частоты (толщины упругого слоя) к скорости волны сдвига в материале упругого тела sc . При стремлении толщин упругого слоя или частот h к частотам зарождения мод ( kph ) скорости мод стремятся к величине скорости волны звука в жидком полу- пространстве 0a ( 0 1,152595a  ). В работе [13] показано, что движения моды 1 при 20h в случае упругого слоя из оргстекла происходят как в жидкости, так и в упру- гом теле. При этом глубина проникновения движений этой моды в упругий слой пре- вышает глубину проникновения в жидкость. Рис. 3 10 Характер влияния предварительного растяжения ( 0 11 0,004  ) на скорости нор- мальных волн в упруго-жидкостных системах иллюстрируют графики на рис. 3, 5, 8, 10, отображающие зависимости относительных изменений величин фазовых скоро- стей с [ ( )с с с c   , с – фазовая скорость мод в системе с предварительно напря- женным слоем, с – фазовая скорость нормальных волн в системе с упругим слоем, не- подверженным начальным деформациям] от толщины упругого слоя (частоты) h . Приведенные на рис. 3 графики позволяют заключить, что для рассматриваемого диапазона частотного спектра начальное растяжение ( 0 11 0,004  ) упругого слоя из менее жесткого материала приводит к повышению фазовых скоростей всех мод. При этом наибольшее влияние начальные напряжения оказывают на скорость моды 1 в окрестности частоты запирания. Графический материал для гидроупругой системы с гипотетической жидкостью ( 0 0,9a  ), приведенный на рис. 4, показывает, что при росте толщины упругого слоя (частоты) h скорость, возникающей единственной моды, стремится снизу к скорости волны Стоунли stc ( 0,732467stc  ), которая меньше скорости волны звука в жидкой среде 0a ( 0 0,9a  ). Рис. 4 Рис. 5 11 Анализ графика, полученного для гидроупругой системы и представленного на рис. 5, позволяет заключить, что начальное растяжение ( 0 11 0,004  ) упругого слоя оказывает существенное влияние на величину фазовой скорости квазилэмбовской моды, в основном, в окрестности частоты ее запирания. В дальнейшем с ростом ча- стоты влияние предварительных деформаций на скорость этой квазиповерхностной волны (волна типа Стоунли) ослабевает. Рис. 6 Из графиков, полученных для упругого слоя из стали [1, 3 – 8, 10], представлен- ных на рис. 6, следует, что скорость первой (нулевой антисимметричной) моды Лэмба при росте частоты или толщины упругого слоя h стремится к скорости волны Рэлея Rc ( 0,923008Rc  ) снизу, а скорость второй (нулевой симметричной) моды стремит- ся к скорости волны Рэлея Rc ( 0,923008Rc  ) сверху. Скорости всех мод Лэмба вы- сокого порядка при увеличении частоты или толщины упругого слоя h стремятся к скорости волны сдвига в материале упругого тела sc . Рис. 7 График для гидроупругой системы, приведенный на рис. 7, показывает, что при росте толщины упругого слоя (частоты) h скорость существующей единственной 12 моды 1 стремится снизу к скорости волны Стоунли stc ( 0, 462886stc  ), которая не- сколько меньше скорости волны звука в жидкой среде 0a ( 0 0, 463021a  ). В работе [13] показано, что глубина проникновения этой моды в жидкость значительно больше глубины проникновения в упругий слой, т.е. движения моды 1 при 20h  в случае более жесткого материала (сталь) происходят в отличие от менее жесткого упругого слоя из оргстекла, преимущественно, в жидкости. Рис. 8 Анализ графика, полученного для гидроупругой системы и представленного на рис. 8, позволяет заключить, что начальное растяжение ( 0 11 0,004  ) упругого слоя из стали (более жесткий материал) оказывает существенное влияние на величину фазо- вой скорости моды, в основном, в окрестности частоты ее запирания. В дальнейшем с ростом частоты влияние предварительных деформаций на скорость квазиповерхност- ной волны (волна типа Стоунли) ослабевает. Рис. 9 Графический материал для гидроупругой системы с гипотетической жидкостью ( 0 0,97a  ), приведенный на рис. 9, показывает, что при росте толщины упругого слоя (частоты) h скорость первой моды стремится снизу к скорости волны Стоунли 13 stc ( 0,878214stc  ), которая меньше скорости волны звука в жидкой среде 0a ( 0 0,97a  ), а скорость второй моды – к скорости волны Рэлея Rc ( 0,923007Rc  ) сверху. Скорость второй моды при стремлении частоты к частоте ее зарождения ( kph ) стремится к скорости волны звука в жидкости 0a ( 0 0,97a  ). Анализ графиков, полученных для этой гидроупругой системы и представленных на рис. 10, показывает, что начальное растяжение ( 0 11 0,004  ) упругого слоя оказы- вает существенное влияние на величину фазовой скорости моды 1, в основном, в окрестности частоты ее запирания. В дальнейшем с ростом частоты предварительные деформации оказывают влияние на фазовые скорости обеих мод на всем частотном интервале. Рис. 10 Рис. 11 На рис. 11 – 19 приведены результаты численных расчетов для упруго- жидкостной системы, состоящей из упругого полупространства и жидкого слоя. На рис. 11, 13, 15, 16 и 18 изображены дисперсионные кривые для гидроупругих волноводов, отражающие зависимости безразмерных величин фазовых скоростей ква- зилэмбовских мод c от безразмерной величины толщины слоя жидкости (частоты) h 14 при отсутствии начальных напряжений. В частности, на рис. 11 представлены графи- ки для упругого полупространства из органического стекла (менее жесткий материал) и воды ( 0 0 1,152595sa a c  ), на рис. 13 – для оргстекла и гипотетической жидкости, скорость распространения звуковой волны в которой равна 0 0,97a  . На рис. 15, 16 приведены графики для упругого полупространства из стали марки 09Г2С (более жесткий материал) и воды, а на рис. 18 – для стали и гипотетической жидкости со скоростью звука в ней 0 1,05a  . Штриховыми линиями на рисунках отмечены асимптотики, к которым стремятся фазовые скорости первой моды при возрастании толщины (частоты) жидкого слоя h . Нетрудно видеть (рис. 11), что скорость единственной существующей поверх- ностной волны в упругом полупространстве из оргстекла (менее жесткий материал), взаимодействующем со слоем воды, изменяется от скорости волны Рэлея Rc ( 0,933557Rc  ) при 0h  до скорости волны Стоунли stc ( 0,7717101stc  ) при h . Отметим, что, как известно [2], фазовая скорость и структура волны Стоунли при взаимодействии твердого и жидкого полупространств зависят от механических параметров гидроупругой системы и определяются соотношением между скоростью волны звука в жидкости и скоростью волны Рэлея в твердом полупространстве. В рассматриваемом случае механические параметры гидроупругой системы «оргстекло – вода» таковы, что скорость распространения звуковой волны в жидкости 0a ( 0 1,152595a  ) больше скорости рэлеевской волны Rc ( 0,933557Rc  ). Как следует из результатов, полученных в работе [24], это приводит к тому, что в коротковолно- вой части спектра мода 1, распространяясь вдоль границы раздела сред, локализуется в приконтактных областях как жидкости, так и упругого полупространства. При этом глубина проникновения движений этой поверхностной волны (волна типа Стоунли) в упругое тело немного превышает глубину проникновения их в жидкость. Начальное растяжение ( 0 11 0,004  ), как следует из графика на рис. 12, оказывает влияние на фазовую скорость поверхностной волны во всем частотном интервале. При этом происходит повышение величины фазовой скорости этой волны во всем диапазоне частот. Наибольшее влияние предварительные деформации оказывают на скорость поверхностной волны в окрестности частоты ее зарождения. Рис. 12 15 Рис. 13 В случае гипотетической жидкости ( 0 0,97a  ) (см. рис. 13) в гидроупругом вол- новоде распространяется множество мод. При этом скорость первой моды изменяется от скорости волны Рэлея Rc ( 0,933557Rc  ) при 0h до скорости волны Стоунли stc ( 0,747764stc  ) при h . С ростом толщины слоя жидкости мода 1 трансфор- мируется в поверхностную волну Стоунли. Нетрудно видеть, что при зарождении мод высокого порядка их скорости равны скорости волны сдвига в материале упругого полупространства sс . Далее с ростом частоты (толщины жидкого слоя) h их скоро- сти стремятся к скорости волны звука в жидкости 0a ( 0 0,97a  ). Учитывая результа- ты, полученные в работе [24], в данной гидроупругой системе при 0 0,97 0,933557Ra c   в коротковолновой части спектра мода 1, распространяясь вдоль границы раздела сред, локализуется в приконтактных областях не только жид- кости, но и упругого полупространства. Рис. 14 Из графиков, приведенных на рис. 14, следует, что начальное растяжение 0 11(  0,004) оказывает влияние, в основном, на фазовую скорость квазиповерхностной 16 моды 1 (волна типа Стоунли). Предварительные деформации оказывают небольшое влияние на скорости мод высокого порядка, преимущественно, в окрестности их кри- тических частот ( kph ). Далее с ростом толщины слоя жидкости воздействие началь- ных напряжений заметно ослабевает. На рис. 15, 16 приведены дисперсионные кривые для системы «упругое полупро- странство из стали (более жесткий материал) – слой воды». Как видно из графиков, фазовая скорость первой моды изменяется от скорости волны Рэлея Rc ( 0,923007Rc  ) при 0h до скорости волны Стоунли stc ( 0, 462886stc  ) при h . Скорости мод высокого порядка при частотах их зарож- дения ( kph ) равны скорости волны сдвига в материале упругого полупространства sc и в дальнейшем с увеличением толщины жидкого слоя h стремятся к скорости волны звука в жидкости 0 0( 0,463021).a a  Отметим, что при 0 0, 463021 0,923007Ra c   , как показано в работе [24], при увеличении частоты (толщины слоя жидкости) низшая мода, распространяясь вдоль границы раздела сред, локализуется в приконтактной области жидкого слоя. Моды высокого порядка с ростом частоты также локализуются в жидком слое. При этом движения в модах высокого порядка при h  смещаются от поверхностей слоя жидкости в его толщу. Рис. 15 Рис. 16 17 Рис. 17 Из графиков, приведенных на рис. 17, следует, что начальное растяжение 0 11(  0,004) оказывает влияние на фазовые скорости всех мод, в основном, в окрестности критических частот ( kph ). Далее с ростом толщины слоя жидкости h воздействие предварительных деформаций значительно ослабевает. Рис. 18 В случае слоя гипотетической жидкости ( 0 1,05a  ), как следует из графика на рис. 18, в гидроупругой системе распространяется лишь одна волна, скорость которой изменяется от скорости волны Рэлея ( 0,923007)R Rc c  при 0h  до скорости волны Стоунли ( 0,887755)st stc c  при .h  В рассматриваемом случае механические параметры гидроупругой системы таковы, что скорость распространения звуковой волны в жидкости 0 0( 1,05)a a  больше скорости рэлеевской волны (R Rc c  0,923007). Учитывая результаты, полученные в работе [24], в данной гидроупругой системе это приводит к тому, что в коротковолновой части спектра поверхностная мода 1 (волна типа Стоунли), распространяясь вдоль границы раздела сред, локализуется в приконтактных областях как слоя жидкости, так и упругого полупространства. 18 Рис. 19 Начальное растяжение 0 11( 0,004),  как следует из рис. 19, приводит к повыше- нию фазовой скорости этой волны во всем интервале частот. Наибольшее влияние предварительные деформации оказывают на скорость поверхностной волны в окрест- ности частоты ее зарождения. §4. Критерии существования квазилэмбовских мод в гидроупругих волноводах. Анализ дисперсионных уравнений, проведенные дополнительные расчеты и по- лученные в данной работе числовые результаты, позволяют установить критерии су- ществования квазилэмбовских мод в упруго-жидкостных волноводах. В системе, состоящей из упругого слоя и полупространства жидкости, такими критериями могут служить соотношения между скоростями волны Рэлея в упругом слое Rc , волны звука в жидкости 0a и волны сдвига в материале упругого слоя sc . В случае выполнения условия 0 Ra c идеальная жидкость не способствует генерирова- нию продольной и сдвиговой волн и взаимодействию их на поверхностях упругого слоя. Поэтому в упругом слое не формируются незатухающие квазилэмбовские моды высокого порядка. В гидроупругом волноводе возникает лишь одна низшая первая мода, которая распространяется без демпфирования вдоль границы контакта сред. Ее скорость изменяется от нуля (при 0h ) до скорости волны Стоунли stc (при ).h  Это подтверждают графики, представленные на рис. 4 (оргстекло – гипоте- тическая жидкость 0 0,9 0,933558)Ra c   и рис. 7 (сталь – вода 0 0,463021a   0,923007).Rc  Если 0R sc a c  , то в волноводе формируются две моды [рис. 9 (сталь – гипотетическая жидкость 00,923007 0,97 1)].R sc a c     При этом ско- рость первой моды изменяется от нуля (при 0h ) до скорости волны Стоунли stc (при h .), а скорость второй моды – от скорости волны звука 0a (при крh h ) к скорости волны Рэлея Rc (при h  ). В случае выполнения условия 0 sa c , как видно из графиков, представленных на рис. 2 (оргстекло – вода 0 1,152595 1sa c   ), идеальная сжимаемая жидкость не препятствует взаимодействию продольной и сдви- говой волн на поверхностях упругого слоя. Вследствие этого в упругом слое возника- ет полный спектр незатухающих квазилэмбовских мод высокого порядка. При этом скорость первой моды изменяется от нуля при 0h  до скорости волны Стоунли stc 19 при h . Скорость второй моды при 0h стремится к скорости волны звука в жидкости 0a , а при h  – к скорости волны Рэлея Rc . Скорости мод высокого по- рядка изменяются от скорости звуковой волны в жидкости 0a до скорости волны сдвига в материале упругого тела .sc Случай 0 sa c на примере упругого слоя из резины рассмотрен и подробно изложен в работе [24]. В ней показано, что дисперси- онная картина волнового поля гидроупругой системы, несмотря на ряд различий, по- хожа на картину волн Лэмба в упругом слое, не взаимодействующем с жидкостью. Здесь отметим лишь одно из отличий, состоящее в том, что в гидроупругом волноводе фазовые скорости мод высокого порядка, начиная с третьей моды, при стремлении частот h к их частотам зарождения крh стремятся к скорости волны звука в жидкости 0 0( ).ia c a В чисто упругом же слое, не взаимодействующем с жидкостью, при крh h фазовые скорости мод Лэмба высокого порядка неограниченно возрастают ( ).ic  В системе, компонентами которой является упругое полупространство и слой жидкости, таким критерием может служить соотношение между скоростью волны звука в жидком слое 0a и скоростью волны сдвига в материале упругого тела sc . При 0 sa c , как видно из графиков, представленных на рис. 11 (оргстекло – вода 0 1,152595 1sa c   ) и рис. 17 (сталь – гипотетическая жидкость 0 1,05 1sa c   ), в гидроупругой системе распространяется лишь одна поверхностная волна, скорость которой изменяется от скорости волны Рэлея Rc (при 0h ) до скорости волны Сто- унли stc (при h ). В случае выполнения условия 0 sa c в гидроупругом волно- воде распространяется множество квазилэмбовских мод. При этом, как следует из графиков на рис. 13, в системе «оргстекло – гипотетическая жидкость» с параметрами 0 0,97 1sa c   на интервале изменения частоты h от 0 до 50 генерируется четыре моды, а в системе «сталь – вода» при 0 0, 463021 1sa c   (рис. 15, 16) на частотном интервале 0 – 10 распространяется без демпфирования 7 мод. При этом скорость пер- вой моды изменяется от скорости волны Рэлея Rc при 0h до скорости волны Сто- унли stc при h . Скорости мод высокого порядка изменяются от скорости волны сдвига в материале упругого полупространства sc при kph h до скорости волны звука в жидкости 0a при .h  §5. Особенности влияния начальных напряжений на квазилэмбовские волны в гидроупругих системах. В гидроупругой системе, состоящей из упругого слоя и полупространства жидко- сти, начальные напряжения вызывают изменение частот зарождения нормальных волн и смещение их дисперсионных кривых. Как видно из графиков, представленных на рис. 3, 5, 8 и 10, начальное растяжение 0 11( 0,004)  приводит к сдвигу критиче- ских частот и дисперсионных кривых мод в коротковолновую часть частотного спек- тра. Следствием этого является повышение величин фазовых скоростей квазилэмбов- ских волн. При этом наибольшее влияние начальные напряжения оказывают на ско- рости низших мод в окрестности их частот запирания. Характер и степень воздей- ствия начальных напряжений на фазовые скорости мод в данной гидроупругой систе- ме тесно связаны со свойствами их локализации. При 0 ,Ra c как отмечалось, низшие моды, распространяясь вдоль границы раздела сред, проникают в упругий слой [13]. Этим и обусловлено влияние упругого тела и предварительных деформаций на фазо- вые скорости этих мод во всем интервале частот. Если между механическими пара- 20 метрами гидроупругой системы выполняется условие 0 Ra c , то низшая мода локали- зуется в жидкости [13]. Поэтому упругое тело и начальные напряжения оказывают су- щественное влияние на величину фазовой скорости этой моды, в основном, в окрестнос- ти частоты ее зарождения. В дальнейшем с ростом частоты их влияние ослабевает. В случае жидкого слоя, взаимодействующего с упругим полупространством, начальные напряжения также вызывают изменение частот зарождения мод и смеще- ние их дисперсионных кривых. Как видно из графиков, приведенных на рис. 12, 14, 17 и 19, начальное растяжение 0 11( 0,004)  смещает критические частоты и дисперси- онные кривые мод в высокочастотную часть спектра. Это приводит к повышению величин фазовых скоростей всех нормальных волн. При этом начальное растяжение оказывает значительное влияние на величины фазовых скоростей низших мод, в ос- новном, в окрестности частот их зарождения. Различия в закономерностях воздей- ствия предварительных деформаций на скорости первых мод для гидроупругих си- стем при 0 Ra c и 0 Ra c связаны с локализацией этих мод в разных средах. При 0 Ra c первая мода распространяется как в жидкости, так и в упругом теле. Этим обусловлено влияние начальных напряжений на скорость этой моды во всем частот- ном диапазоне. При 0 Ra c , как показано в работе [24], с ростом частоты мода 1 ло- кализуется в жидком слое. Поэтому упругое тело и начальные напряжения не оказы- вают значительного влияния на величину фазовой скорости этой моды. В случае многомодовых волноводов (рис. 14 и 17) наибольшее влияние начальные напряжения оказывают на фазовые скорости мод высокого порядка в окрестности критических частот. Характер влияния начальных напряжений на скорости этих мод при увеличении толщины слоя жидкости (частоты) вызван тем, что при ,h  дви- жения в модах высокого порядка, начиная со второй (рис. 14 и 17), смещаются от по- верхностей слоя жидкости в его толщу. Этим и объясняется незначительное влияние упругого полупространства и начальных напряжений на фазовые скорости, а также дисперсию этих мод в коротковолновой части частотного спектра. Заключение. В рамках трехмерных уравнений линеаризованной теории упругости конечных деформаций для упругого тела и трехмерных линеаризованных уравнений Эйлера для идеальной сжимаемой жидкости даны постановки и решение задач о распространении квазилэмбовских нормальных волн в системах, состоящих из идеального сжимаемого жидкого полупространства и предварительно напряженного сжимаемого упругого слоя, а также из слоя жидкости и предварительно деформированного сжимаемого упругого полупространства. С использованием представлений общих решений полу- чены характеристические уравнения и построены дисперсионные кривые мод в ши- роком диапазоне частот. Анализ полученных результатов показал, что основными критериями существо- вания незатухающих квазилэмбовских нормальных волн в гидроупругой системе «идеальное сжимаемое жидкое полупространство – упругий слой» являются соотно- шения между величинами скоростей волны звука в жидкости, сдвига в материале упру- гого тела и квазирэлеевской волны, распространяющейся вдоль свободной поверхности упругого слоя. В случае системы, компонентами которой являются слой жидкости и упругое полупространство, таким критерием может служить соотношение между ско- ростью волны звука в жидкости и скоростью волны сдвига в материале упругого тела. Показано, что в гидроупругих волноводах начальное растяжение вызывает измене- ние частот зарождения квазилэмбовских мод, конфигурации их дисперсионных кривых, а также смещение их в коротковолновую часть спектра. Характер влияния начальных напряжений на скорости мод обусловлен их свойствами локализации в средах. Развитый подход и полученные результаты позволяют установить пределы при- менения моделей волновых процессов, основанных на различных вариантах теории малых начальных деформаций и классической теории упругости. 21 РЕЗЮМЕ. Досліджено поширення квазілембових хвиль у попередньо деформованому стисли- вому пружному шарі, що взаємодіє з півпростором ідеальної стисливої рідини, а також у шарі ріди- ни, що контактує з попередньо напруженим стисливим пружним півпростором. Результати отримано на основі тривимірних рівнянь лінеаризованої теорії пружності скінченних деформацій для стисливо- го пружного тіла та тривимірних лінеаризованих рівнянь Ейлера для ідеальної стисливої рідини. Застосовано постановку задачі та підхід, основані на використанні представлень загальних розв'язків лінеаризованих рівнянь для пружного тіла та рідини. Отримано дисперсійні рівняння, які описують поширення квазілембових хвиль у гідропружних системах у широкому діапазоні частот. Проаналізо- вано вплив початкових напружень, а також товщини пружного шару та шару ідеальної стисливої рідини на фазові швидкості квазілембових мод. Запропоновано критерії існування квазілембових мод в гідропружних хвилеводах. Розвинутий підхід та отримані результати дозволяють встановити гра- ниці застосування моделей хвильових процесів, oснованих на різних варіантах теорії малих початко- вих деформацій, а також на класичній теорії пружності. Числові результати наведено у вигляді гра- фіків та дано їх аналіз. 1. Викторов И.А. Звуковые поверхностные волны в твердых телах. – М.: Наука, 1981. – 288 с. 2. Волькенштейн М.М., Левин В.М. Структура волны Стоунли на границе вязкой жидкости и твердо- го тела // Акуст. журнал. – 1988. – 34, № 4. – С. 608 – 615. 3. Гузь А.Н. Упругие волны в телах с начальными напряжениями: в 2-х томах. Т.1. Общие вопросы – К.: Наук. думка, 1986. – 374 с. 4. Гузь А.Н. Упругие волны в телах с начальными напряжениями: в 2-х томах. Т.2. Закономерности распространения. – К.: Наук. думка, 1986. – 536 с. 5. Гузь А.Н. Динамика сжимаемой вязкой жидкости. – К.: А.С.К., 1998. – 350 с. 6. Гузь А.Н. Упругие волны в телах с начальными (остаточными) напряжениями. – К.: А.С.К., 2004. – 672 с. 7. Гузь А. Упругие волны в телах с начальными (остаточными) напряжениями: в 2-х частях. Ч. 1. Общие вопросы. Волны в бесконечных телах и поверхностные волны. – Saarbrucken: LAP LAM- BERT Academic Publishing, 2016. – 501 с. 8. Гузь А. Упругие волны в телах с начальными (остаточными) напряжениями: в 2-х частях. Ч. 2. Волны в частичноограниченных телах. – Saarbrucken: LAP LAMBERT Academic Publishing, 2016. – 505 с. 9. Гузь А.Н. Введение в динамику сжимаемой вязкой жидкости. – Saarbrucken: LAP LAMBERT Academic Publishing RU, 2017. – 244 с. 10. Гузь А.Н., Жук А.П., Махорт Ф.Г. Волны в слое с начальными напряжениями. – К.: Наук. думка, 1976. – 104 с. 11. Жук А.П. Волны Стонли в среде с начальными напряжениями // Прикл. механика. – 1980. – 16, № 1. – С. 113 – 116. 12. Bagno A.M. Effect of Prestresses on the Dispersion of Quasi-Lamb Waves in the System Consisting of an Ideal Liquid Layer and a Compressible Elastic Layer // Int. Appl. Mech. – 2017. – 53, N 2. – P. 139 – 148. 13. Bagno A.M. Dispersion Properties of Lamb Waves in an Elastic Layer–Ideal Liquid Half-Space System // Int. Appl. Mech. – 2017. – 53, N 6. – P. 609 – 616. 14. Bagno A.M., Guz A.N. Elastic Waves in Pre-stressed Bodies Interacting with a Fluid (survey) // Int. Appl. Mech. – 1997. – 33, N 6. – P. 435 – 463. 15. Bagno A.M., Guz A.N. Effect of Prestresses on the Dispersion of Waves in a System Consisting of a Viscous Liquid Layer and Compressible Elastic Layer // Int. Appl. Mech. – 2016. – 52, N 4. – P. 333 – 341. 16. Drinkwater B.W., Wilcox P.D. Ultrasonic arrays for non-destructive evaluation: A review // NDT and E Int. – 2006. – 39, N 7. – P. 525 – 541. 17. Gibson A., Popovics J. Lamb wave basis for impact-echo method analysis // J. Engng. Mech. – 2005. – 131, N 4. – P. 438 – 443. 18. Guz A.N. Aerohydroelasticity Problems for Bodies with Initial Stresses // Int. Appl. Mech. – 1980. – 16, N 3. – P. 175 – 190. 19. Guz A.N. Elastic Waves in Bodies with Initial (Residual) Stresses // Int. Appl. Mech. – 2002. – 38, N 1. – P. 23 – 59. 20. Guz A.N. Dynamics of Compressible Viscous Fluid. – Cambridge: Cambridge Scientific Publishers, 2009. – 428 p. 21. Guz A.N. On the foundations of the ultrasonic non-destructive determination of stresses in near-the- surface layers of materials. Review // J. Phys. Science and Application. – 2011. – 1, N 1. – P. 1 – 15. 22 22. Guz A.N. Ultrasonic Nondestructive Method for Stress Analysis of Structural Members and Near-Surface Layers of Materials: Focus on Ukrainian Research (Review) // Int. Appl. Mech. – 2014. – 50, N 3. – P. 231 – 252. 23. Guz A.N., Bagno A.M. Effect of Prestresses on the Dispersion of Lamb Waves in a System Consisting of a Viscous Liquid Layer and a Compressible Elastic Layer // Int. Appl. Mech. – 2018. – 54, N 3. – P. 249 – 258. 24. Guz A.N., Bagno A.M. Effect of Prestresses on Lamb Waves in a System Consisting of an Ideal Liquid Half-Space and an Elastic Layer // Int. Appl. Mech. – 2018. – 54, N 5. – P. 495 – 505. 25. Guz A.N., Bagno A.M. Effect of Prestresses on Quasi-Lamb Waves in a System Consisting of a Com- pressible Viscous Liquid Half-Space and an Elastic Layer // Int. Appl. Mech. – 2018. – 54, N 6. – P. 617 – 627. 26. Guz A.N., Zhuk A.P., Bagno A.M. Dynamics of Elastic Bodies, Solid Particles, and Fluid Parcels in a Compressible Viscous Fluid (Review) // Int. Appl. Mech. – 2016. – 52, N 5. – P. 449 – 507. 27. Jhang K.Y. Nonlinear ultrasonic techniques for nondestructive assessment of micro damage in material: a review // Int. J. Precision Engng. and Manufacturing. – 2009. – 10, N 1. – P. 123 – 135. 28. Kessler S.S., Spearing S.M., Soutis C. Damage detection in composite materials using Lamb wave meth- ods // Smart Materials and Structures. – 2002. – 11, N 2. – P. 269 – 279. 29. Kobayashi M., Tang S., Miura S., Iwabuchi K., Oomori S., Fujiki H. Ultrasonic nondestructive material evaluation method and study on texture and cross slip effects under simple and pure shear states // Int. J. Plasticity. – 2003. – 19, N 6. – P. 771 – 804. 30. Leonard K.R., Malyarenko E.V., Hinders M.K. Ultrasonic Lamb wave tomography // Inverse Problems. – 2002. – 18, N 6. – P. 1795 – 1808. 31. Liu L., Ju Y. A high-efficiency nondestructive method for remote detection and quantitative evaluation of pipe wall thinning using microwaves // NDT and E Int. – 2011. – 44, N 1. – P. 106 – 110. 32. Ottenio M., Destrade M., Ogden R.W. Acoustic waves at the interface of a pre-stressed incompressible elastic solid and a viscous fluid // Int. J. Non-Lin. Mech. – 2007. – 42, N 2. – P. 310 – 320. 33. Ramadas C., Balasubramaniam K., Joshi M., Krishnamurthy C.V. Interaction of the primary anti- symmetric Lamb mode (Ao) with symmetric delaminations: numerical and experimental studies // Smart Materials and Structures. – 2009. – 18, N 8. – P. 1 – 7. 34. Rossini N.S., Dassisti M., Benyounis K.Y., Olabi A.G. Methods of measuring residual stresses in compo- nents // Materials and Design. – 2012. – 35, March. – P. 572 – 588. 35. Spies M. Analytical methods for modeling of ultrasonic nondestructive testing of anisotropic media // Ultrasonics. – 2004. – 42, N 1 – 9. – P. 213 – 219. Поступила 17.01.2018 Утверждена в печать 05.11.2019