Рассмотрим собственные колебания в упругой среде, содержащей двухслойный цилиндр с жидкостью в цилиндрической системе координат (r,z,θ). Обозначим через Vpi, Vsi, i, μi, i (i =0,1,2,3) соответственно скорость продольной и поперечной волны, плотность и модуль упругости. Рассмотрим задачу о распространении свободных волн, возникающих в такой системе. Уравнения движения среды для продольных φi и поперечных ψi потенциалов представляется в виде:
(1)
В уравнениях (1) следует положить 
=0, если в затрубном пространстве находится жидкость. Соответствующие напряжения rr, ir и смещения u r, u z, определяются через потенциалы , равенствами
(2)
На границах раздела упругой среды с жидкостью выполняются граничные условия непрерывности нормальных составляющих смещений и напряжений, а также равенство нулю касательных напряжений в твердом теле:
Если жидкость заменена упругой средой, то на контакте двух сред ставятся следующие условия:
Решения уравнений (1), удовлетворяющие условию конечности среды на оси r=0 и условиям убывания на бесконечности, выражаются через модифицированные функции Бесселя.
На границе контакта слоев (r = r2) ставится условие скользящего контакта (непрерывны нормальные составляющие напряжений и смещений, отсутствуют касательные составляющие напряжений) и контакт между внешним слоем и окружающей средой жесткий (непрерывны нормальные и тангенциальные составляющие напряжений и смещений). Дисперсионное уравнение записывается в виде:
(к, ) = 0. (4)
Оно представимо в виде определителя, у которого элементы ij (
), отличные от нуля, имеют следующий вид:
Остальные элементы также записываются в аналогичном виде. Здесь в двухслойном цилиндре V определяет скорость обобщенной волны по двухслойному цилиндру.
Фазовая скорость волны определятся величиной реальной части корня, то есть V = ReV; величины мнимой части корня связаны с затуханием на единице расстояния зависимостью
.
Дисперсионное уравнение (4) решается методом Мюллера. Значение левой части на каждой итерации метода Мюллера определяется методом Гаусса с выделением главного элемента. Нами были составлены программы и проведены расчеты дисперсии и затухания волны Лэмба для моделей скважин, описываемых граничными условиями (3а) и (3b). Исходя из физической постановки задачи, будем считать, что поглощением обладают буровая жидкость, цемент, тампонажная смесь, поглощением же в материале колонны и в окружающей среде будем пренебрегать. Переход к системе с поглощением был сделан посредством введения комплексных параметров сред. Численные результаты получены при следующих значениях параметров:
Vp0 =1500 м/c; Vp0 =1500 м/c; Vp2 =1500 м/c;
Vp3 =5300 м/c; Vs1 =2900 м/c; Vs2 =2000 м/c;
Vs3 =2000 м/c; ρ0 =1 г/cм3; ρ1 =8 г/cм3;
Ρ2 =3 г/cм3; ρ4 =4 г/cм3; r1 =0,05 м; r2 =0,06 м; r3 =0,067 м.
Результаты расчетов представлены в табл. 1. Видно, что фазовая скорость слабо зависит от волнового числа.
Таблица 1
Изменение фазовой скорости V (м/c) в зависимости от волнового числа α/a.
|
α/a |
1 |
2 |
3 |
4 |
5 |
|
1 |
1450,121 |
1450,024 |
1451,01 |
1456,28 |
1450,82 |
|
2 |
1456,3 |
1456,09 |
1456,013 |
1456,72 |
1455,3 |
|
3 |
1462,24 |
1462,50 |
1462,03 |
1462,4 |
1461,9 |
Из анализа значений фазовой скорости выявлено, что разница между скоростями осе симметричных и не осе симметричных волн первой моды мала для всех значений волнового числа, кроме близких к нулю (область очень длинных волн), а минимумы части первой моды для всех значений n совпадают, так что и в данном случае первая резонансная скорость может быть определена из решения соответствующей осе симметричной задачи. Так, при γ>250 получено С=0,31 (3ρ/G)
. Как показали расчеты,
Рис. 1. Зависимости фазовой скорости от волнового числа
Из результатов выясняется, что сжимающие контактные напряжения имеют место в некоторой окрестности приложения каждой силы. С удалением от точки приложения силы по окружности напряжения для всех рассмотренных случаев меняют знак. Это является следствием предположения двухсторонних характеров связи между оболочкой и заполнителем.
Литература:
1. Сафаров И. И. Колебания и волны в диссипативно неоднородных средах и конструкциях.Ташкент; Фан.1992 г.250 с.
