Full Text

Введение Впервые модель, состоящая из твердого каркаса с многочисленными связанными между собой порами, была предложена в сороковые годы XX в. Морисом Био (Maurice A. Biot) [1], поэтому такие среды часто называют средами Био. В работе [1], как и в последующих работах цикла [2-8], пористость материала предполагалась постоянной, что не соответствует физической реальности. Несмотря на это допущение, работы М. Био, как пионерские, несомненно, внесли большой вклад в теорию пороупругости, в частности, в них выявлен новый физический эффект - существование двух продольных волн - быстрой и медленной. Формула изменения пористости при малых деформациях впервые предложена Я.И. Френкелем и опубликована в статье [9], являющейся второй классической работой в этой области. В работе Л.Я. Косачевского [10] последовательно изложена линейная теория деформирования двухфазной среды, учитывающая результаты работ М. Био и Я.И. Френкеля, получено выражение для энергии, коэффициенты которого зависят от таких параметров среды, как исходная пористость, коэффициент сжимаемости и др. Развитие основных положений механики и физики пористых материалов содержится в монографиях [11-15]. Вопросам распространения волн в пористых материалах посвящены, в частности, работы [16-25]. С 1999 г. регулярно проводятся крупные многопрофильные международные конференции по механике пористых материалов («Poromechanics») [26-30]. В монографии [31] отмечается, что во многих природных средах (туф, пемза и т.д.), а также в искусственных композитах в сочетании с пористостью встречаются многочисленные, часто заполненные жидкостью, полости. Таким образом, изучение пористо-полостных твердых сред представляет большой интерес с точки зрения теории и практики. 1. Постановка задачи Пусть имеется полубесконечная или в форме слоя среда Био, в которой существуют полости с жидкостью. Считается, что поры цилиндрической формы, а полости - шарообразной. Под действием упругой волны полости будут колебаться, оказывая, в свою очередь, влияние на распространение волны. Уравнения, описывающие деформацию пористой жидконасыщенной среды с полостями, наполненными жидкостью, в лагранжевых координатах приведена в [31]. В этой системе считаем, что вязкость отсутствует, а всеми типами нелинейности пренебрегаем. Рассматриваем плоскую продольную волну, распространяющуюся вдоль оси . Ограничиваясь вторыми производными и линейными членами, получим следующую систему, описывающую распространение плоской продольной волны в среде Био: (1) где , - перемещения твердой и жидкой фаз соответственно вдоль направления оси ; - объем полости, возмущенной волной; , - эффективные начальные плотности массы твердой и жидкой фазы соответственно; , , , , - массы твердой и жидкой фазы соответственно; - присоединенная плотность массы; , - коэффициенты Ламе; , - известные коэффициенты, характеризующие жидконасыщенную среду, их значения даны в [10]; - количество полостей в единице объема; - резонансная частота, ; - радиус полостей; - плотность твердой фазы; - начальное давление в жидкости внутри полости; - показатель адиабаты жидкости. 2. Дисперсионные свойства системы Для получения дисперсионного уравнения подставляем в систему линейных дифференциальных уравнений (1) решение в виде монохроматических волн: группируем слагаемые при амплитудах , , , выписываем соответствующий определитель и приравниваем его нулю. Раскрывая определитель, получаем дисперсионное соотношение (2) где для сокращения записи введены обозначения , , , , , . Выбрав в качестве характерных величин частоты и волнового числа резонансную частоту и отношение резонансной частоты к скорости распространения упругой продольной волны в воде , , где , - безразмерные величины частоты и волнового числа соответственно, получаем дисперсионное соотношение (2) в безразмерном виде (черта опущена): (3) где ( ) - безразмерные коэффициенты: , , , , . Уравнение (3) имеет шестой порядок по и четвертый порядок по . Дисперсионные кривые имеют три ветви, две из которых исходят из начала координат и одна из точек . Для безразмерных частоты и волнового числа получены следующие характеристики при различных значениях параметра уравнения (3) (рис. 1). Из рисунков видно, что внутренние ветви, теряя гладкость, стремятся к нулю. В начале координат происходит четкое разделение дисперсионных кривых на три пары, из которых одна начинает вытягиваться вдоль оси , а другая меняет форму, принимая вид ломаной, касается третьей. В точке касания происходит разрыв дисперсионной кривой. При дальнейшем увеличении параметра одна ветвь кривой исчезает с действительной плоскости, а две другие, разделившись на две общие части, продолжают отдаляться друг от друга. Рис. 1. Динамика дисперсионной кривой при изменении параметра , Fig. 1. Dynamics of disperse curve depending on alterations of parameter , Дисперсионная кривая, заданная уравнением (3), имеет две горизонтальные и четыре наклонные проходящие через начало координат асимптоты где . Вертикальные асимптоты у кривой отсутствуют. Подставляя в (3), находим зависимость фазовой скорости и волнового числа: . (4) При одна ветвь кривой уходит в бесконечность, а две другие принимают значения . У кривых фазовых скоростей, определяемых зависимостью (4), наклонные асимптоты отсутствуют. Но есть шесть горизонтальных асимптот, две из которых нулевые и две - вертикальные нулевые асимптоты: , , . Дифференцируя (3) по и подставляя , получаем связь групповой скорости и волнового числа: . (5) В (5) определяется из дисперсионного уравнения (3). При одна ветвь кривой групповой скорости устремляется в начало координат, а две другие принимают соответствующие значения фазовых скоростей. Кривые групповых скоростей, определяемые зависимостью (5), имеют шесть горизонтальных асимптот, из которых две нулевые: , . Наклонные и вертикальные асимптоты отсутствуют. Рассмотрим варианты расположения трех веток дисперсионной кривой, которые показывают характерное поведение системы. В силу симметрии графиков функций (3), (4), (5) относительно осей координат достаточно их изображения в первой четверти. Для удобства анализа представим зависимости , , в одних осях (рис. 2). На рисунке сплошной линией показана зависимость , длинным пунктиром - , штрихпунктиром - , коротким пунктиром отмечены асимптоты кривых. Из рисунка видно, что первая ветвь дисперсионной кривой зарождается на частоте . При небольших волновых числах ( ), т.е. для длинных волн ( ), частота стремится к значению , на бесконечности частота стремится к значению . Фазовая и групповая скорости соответствующих ветвей стремятся к значению с разных сторон. Вторая ветвь дисперсионной кривой зарождается на частоте . На бесконечности частота стремится к значению . Соответствующие ветви кривых фазовой и групповой скорости стремятся к общей асимптоте с разных сторон. Третья ветвь дисперсионной кривой зарождается на частоте . На бесконечности частота стремится к значению , значения фазовой и групповой скоростей стремятся к нулю. На промежутке групповая скорость меньше фазовой скорости, что говорит о нормальной дисперсии при выбранных значениях параметров системы. Рис. 2. Зависимости , , Fig. 2. Dependences of , , Рис. 3. Плотность спектра частот Fig. 3. Density of spectral frequencies Плотность спектрального распределения [32] позволяет определить промежутки частот, на которых проявляется явление дисперсии волн. Вычисляется плотность спектрального распределения по следующей формуле: где константу можно выбрать равной единице, так как на качественный вид кривой она не повлияет. С точностью до множителя плотность является величиной обратной к групповой скорости. Из графика плотности спектра частот несложно определить, при приближении к каким частотам спектр становится непрерывным. Плотность спектрального распределения для каждой ветви дисперсионной кривой изображена на рис. 3. Из рисунка видно, что спектр частот сгущается и становится непрерывным при . Спектр собственных частот . Значит, для медленной волны (на рис. 1 - третья ветка) дисперсия проявляется на высоких частотах. Также из рис. 3 видно, что спектр частот становится непрерывным при . Спектр собственных частот лежит в диапазоне . Значит, у быстрой волны, которая возникает из-за полостности среды, дисперсия проявляется на низких частотах. Небольшое увеличение плотности спектрального распределения у первой ветки дисперсионной кривой, соответствующей быстрой волне, возникшей в результате пористости среды, соответствует незначительной дисперсии. Спектр собственных частот . При изменении одного из параметров дисперсионная кривая претерпевает разрыв. Значение параметра , , при котором происходит разрыв дисперсионной кривой. Значение координат точки разрыва можно найти из условия пересечения кривой (3) со своей асимптотой или из условия пересечения двух веток дисперсионной кривой. Дисперсионная кривая после разрыва представлена на рис. 4. Дисперсионная кривая (рис. 4) по-прежнему имеет три ветви, но их вид качественно изменился. Появилась ветка с двумя частотами зарождения и «свободная» ветка. На бесконечности кривая сохранила линейную асимптотику. Участки первой и второй веток дисперсионной кривой «перепрыгнули» асимптоту . Кривая групповой скорости приобретает асимптоты, аналитический вид которых не удается отыскать. Дисперсия приобрела аномальный характер (рис. 4) на тех участках дисперсионной кривой, которые сменили расположение относительно асимптоты . Рис. 4. Зависимости (сплошная); (длинный пунктир); (штрихпунктир) Fig. 4. Dependences of (solid line); (longdashedline); (dot-dash line) Рис. 5. Плотность спектра частот Fig. 5. Density of spectral frequencies Спектр частот (рис. 5) становится непрерывным при стремлении частоты к четырем критическим частотам: , и частотам, вблизи которых дисперсионные кривые приобретают пики. Заключение Таким образом, в статье показано следующее: 1. В среде распространяются три продольные волны, две - быстрая и медленная за счет жидконасыщенной пористой среды и одна (на рис. 2. вторая ветка) - за счет полостности среды. При отсутствии полостей в среде дисперсионная кривая имела бы две ветви (как в среде Био): первая, которая всегда присутствует на плоскости , и третья, которая появляется вследствие пористости среды. 2. Волны, распространяющиеся в пористой среде с полостями, обладают значительной дисперсией. Причем быстрые волны обладают дисперсией в области низких частот, а медленные волны - в области высоких частот. 3. При фиксированных параметрах в системе может наблюдаться аномальная дисперсия на двух участках дисперсионной кривой. Существенная дисперсия проявляется вблизи четырех частот. У медленной волны на высоких частотах. У быстрой волны, которой соответствует «свободная» ветка дисперсионной кривой, дисперсия проявляется на низких частотах, а у быстрой волны, которой соответствует ветка с двумя частотами зарождения, - на частотах в области пика и в области ненулевой частоты зарождения.

About the authors

S M Aizikovich

Don State Technical University

V I Erofeev

Institute of Problems of Mechanical Engineering of the Russian Academy of Sciences; Lobachevsky State University of Nizhny Novgorod

A V Leonteva

Don State Technical University

References

Statistics

Views

Abstract - 226

PDF (Russian) - 69