Аналитические и численные решения статических и динамических задач несимметричной теории упругости

Аналитические и численные решения статических и динамических задач несимметричной теории упругости

В.В. Корепанов, М.А. Кулеш, В.П. Матвеенко, И.Н. Шардаков

Институт механики сплошных сред УрО РАН, Пермь, 614013, Россия

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

Analytical and numerical solutions for static and dynamic problems of the asymmetric theory of elasticity

V.V. Korepanov, M.A. Kulesh, V.P. Matveenko, and I.N. Shardakov Institute of Continuous Media Mechanics UB RAS, Perm, 614013, Russia

The paper offers a number of new analytical solutions of static and dynamic wave problems for a linear elastic Cosserat continuum. We compare them to corresponding analytical solutions within the classical theory of elasticity and illustrate the manifestations of moment effects using model materials. Concurrently with the derivation of analytical solutions, the finite-element method has been used to develop a numerical solution algorithm for problems of the asymmetric theory of elasticity on static deformation. The algorithms gave rise to numerical results that significantly expand the considered class of analytical solutions. The obtained results provide the basis for the performance of physical experiments that investigate moment behavior of materials under elastic deformation.

1. Введение

В линейной модели Коссера [1] каждая материальная точка среды является малым абсолютно твердым телом и имеет ориентацию (полярная среда), а нагрузка на произвольную область тела передается посредством распределенного усилия р = п -ст и распределенного момента т = п Д, действующих на поверхности с нормалью п (рис. 1, а).

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

произвольная область тела (рис. 1, б) состоит из абсолютно твердых частиц, имеющих, в отличие от классической теории упругости, линейный размер Ь и связанных между собой упруго, но эти связи действуют на внешний край частиц. В этом случае кинематика частицы будет определяться вектором перемещения центра масс и и вектором поворота ю. В случае континуума Коссера эти векторы являются непрерывными функциями, т.е. каждой материальной точке приписываются свойства частицы, имеющей масштаб и ориентацию. Таким образом, в среде появляются дополнительные физические параметры, характеризующие линейный размер и момент инерции частиц, а также свойства связей относительно вращения, а основные соотношения для линейной упругой среды Коссера имеют вид [2]:

© Корепанов В.В., Кулеш М.А., Матвеенко В.П., Шардаков И.Н., 2007

уравнения движения:

V-ô + X = pii, àT : É + V (1. + Y = jtü, геометрические соотношения:

у = Vu-É ев, % = Vcb, физические уравнения:

д = 2^у(s) + 2ау(a) + X/1 ( y)ê (1 = 2yx(S) + 2ex(A) + P/j (X)ë.

С учетом данных соотношений уравнения движения для вектора перемещения u и вектора поворота ю примут вид:

(2^ + X) grad div u - (^+а) rot rot u +

+ 2а rot œ + X = pü,

(2y + P) grad div ffl - ( y+e) rot rot œ +

+ 2a rot u - 4aœ + Y = jœ.

В приведенных выше уравнениях X — вектор удельной плотности объемных сил; Y — вектор удельной плотности объемных моментов; u — вектор перемещения; ffl — вектор поворота; у и % — тензоры деформаций и изгиба-кручения; â и Д — тензоры напряжений и моментных напряжений; ц, X — постоянные Ламе; a, Р, у, Е — физические постоянные материала в рамках моментной теории упругости; р — плотность; j — плотность момента инерции (мера инерции среды при вращении); É — тензор Леви-Чивиты третьего ранга;

л (S) —

операция симметрирования;

Л (A) —

операция

альтернирования; У(-) — набла-оператор; /)(•) — первый инвариант тензора. В отличие от классической теории тензоры у и б являются несимметричными. Однако несложно видеть, что в случае равенства нулю новых материальных констант а, Р, у, Е, j соотношения (1) трансформируются в уравнения Ламе классической теории упругости.

Для сравнительного анализа будем также рассматривать модель псевдосреды Коссера [2] или среды со

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

Ю = rot и/2,

а физические соотношения и уравнения движения для этой среды имеют вид:

a = 2^y(A) + Х11(у)ё - -2 V-Д- E,

Д = 2yx(S) + 2ex(A), p,V2u + (^ + X) grad div и +

1 2 + 4 (y + s) rot rot V и + X = pii.

Таким образом, для псевдосреды Коссера имеется одна независимая кинематическая неизвестная — вектор перемещения u и в рассмотрение вводятся несимметричные тензоры напряжений д и моментных напряжений Д. Псевдосреда Коссера является частным случаем среды Коссера при условии а^^. Необходимо отметить, что соотношение (2) не позволяет в трехмерном случае сформулировать необходимое число независимых краевых условий. Например, если на поверхности заданы перемещения, то не удается независимо выразить нормальную составляющую вектора поворота. Аналогично обстоит дело и с условиями в нагрузках.

Известны и другие частные случаи модели среды Коссера, например, редуцированная среда Коссера, для которой Р=у = Е = 0 и, как следствие, отсутствует тензор моментных напряжений. Впервые эта модель была исследована в изотропном случае в [3]. Для описания этой среды используются следующие физические соотношения и уравнения движения:

а = 2^y(S) + 2ay(A) + XI1 (y)e,

(2^, + A) grad divu - (^ + a) rot rot u +

+ 2a rot ffl + X = pii,

2a rot u - 4аю+Y = jo.

1,

Рис. 1. Модель среды Коссера

Однако несложно убедиться, что в статическом случае (<Й = 0) и при отсутствии объемных моментов (Y = = 0) данные уравнения преобразуются к классическим уравнениям Ламе, где вектор поворота будет определяться соотношением (2), поэтому в дальнейшем данная модель нами рассматриваться не будет.

Анализ публикаций по моментной теории свидетельствует о том, что с момента опубликования работы братьев Коссера [1] интерес к этой теории возрастает с периодичностью примерно в 30 лет. В настоящий момент мы наблюдаем очередной пик интереса к этой теме в связи с зарождающейся наукой микромеханикой и нанотехнологиями. Однако следует отметить, что в подавляющем большинстве эти работы имеют сугубо теоретический характер. Одной из главных причин, ограничивающих практические приложения несимметричной теории упругости, является отсутствие надежных методов определения материальных констант. Другим важным обстоятельством является недостаточное количество экспериментов, свидетельствующих о моментных эффектах при упругом деформировании материалов.

Важным условием преодоления этих проблем является получение и анализ новых решений задач несимметричной теории упругости. Цель анализа получаемых решений — определение таких значений параметров (характерные геометрические размеры, их соотношение и т.п.), при которых отклик моментных свойств материала на какую-либо макровеличину (экспериментально измеряемую) максимален. Таким образом, результаты анализа позволят определить наиболее информативные принципиальные схемы возможных экспериментов. Именно проведение подобного анализа и явилось целью данной работы и выполнено на примере решения некоторых статических и динамических задач.

2. Аналитические решения статических задач несимметричной теории упругости

В качестве попытки отыскания моментных эффектов при статическом деформировании авторами был проанализирован класс задач, когда напряженно-деформированное состояние в цилиндрической системе координат (р, ф, z) зависит только от координат р и ф. Данный класс задач характеризуется тем, что решение для него можно представить в виде отрезка ряда Фурье:

T

u(p,ф) = {up(р, ф), иф(р, ф),0| ,

up (р, ф) = U (0)(р) + U (1)(p)cos ф +

N г ^

+ X U(й)(р) cos(^),

n=2

Up (Р, Ф) = V (0)(р) + V (1)(р) sin Ф +

N

X V(й)(р) sin(^),

и(р, ф) = {0,0, rnz (р, ф)}т,

rnz (р, ф) = ю(0)(р) + ю(1)(р) sin ф +

N , ч

+ X  (Р) sin(йф).

п=2

Такая форма записи решения будет соответствовать полному отрезку ряда Фурье, если в соотношениях (5) cos(-) заменить на sin(-), а sin(-) на -cos(-). Считаем, что краевые условия представимы рядами, аналогичными (5). Подставляя векторы перемещений и поворота (5) в уравнения (1), получим систему дифференциальных уравнений относительно коэффициентов Фурье U(n^(р), V(n\р), n\р). В силу того, что рассматривае-

мые области являются телами вращения, данная система распадается на три независимых системы относительно нулевой (n = 0), первой (n = 1) и старших гармоник. Например, для коэффициентов нулевой гармоники система разрешающих уравнений имеет вид:

i!U(0>(р) + — U(0)(р) - \U(0)(р) = 0, dp2 Р dp р2

il v <0>(р) + LL v <%) - -1 v (0)(р) -

dp2 Р dp р2

= 0,

а + ц dp

-dT rn(0)(p) +1-d-rn(0)(p) ——ю(0)(р) + dp2 P dp Y + E

+ _2^_d_v(°)(p) + _^v(0>(p) = 0.

Y + e dp (y + e)p

Первое уравнение данной системы описывает решение классической задачи Ламе [4] и не связано с двумя последующими уравнениями этой системы, которые, в свою очередь, описывают задачу кручения несимметричного упругого материала. Из этого видно, что наипростейшей задачей, для которой уже проявляются мо-ментные свойства среды, является задача кручения.

При построении общих решений получающихся систем используется известный факт о том, что частные решения соответствующих краевых задач симметричной теории являются частными решениями для краевых задач несимметричной [2]. Метод получения решения в вышеприведенном виде подробно описан в [5], поэтому здесь приведем лишь итоговые соотношения.

Для удобства дальнейшего анализа перейдем к безразмерной форме записи решений. Для этого в качестве характерного размера выберем некоторый радиус R и воспользуемся новыми безразмерными величинами:

A = R

B =

ар,

n=2

а + р

а ,

C =

Y-e Y + e

Для и

и (0)(р) = С1(0)р +

V(0) (р) = С3(0)р + -4-^ - с5(0) 11(2Ар) -

2 А

—

2 А2

-Г

ю(0)(р) = С(/’ -

2А

+ -

(0) ВК0(2АР) 2А

Можно также выписать решение для первой гармоники:

и (1)(р) = С1(1) + -2) + С3(1)р2 + С 1п р + р2

+ С

V (1)(р) = --

-6

С

С2

Р2

А + ц.

+ —

,(1) 3^+5ц^2 3 р-А,

3р + А ]

- С5(1) |11(2 Ар) - 2 А10 (2 Ар) I +

+ С61) |К-

ю(-)(Р) = -31) _ -4а) Л + 2^ _

3 ц-Л 4 (3ц + А,)р - С5(1) 2 А2 В11 (2 Ар) - с61)2 А2 ВК1 (2 Ар). Решение для старших гармоник будет иметь вид:

г{п) г{п)

и(п )(Р) = -птт + СТл + -3 п)рп- + рп + рп

Лп) Г

+ С4п)рп+1 + 1п (2 Ар) + Кп (2 Ар),

Р Р

V(п) (р) = С-1 + Х(п - 2) +14п - 4) _

рп+1 рп4 Хп + р

(п)пп-1 _ —(п) ^(п + 2) + р(п + 4) пп+1 ^

-С п )р

Хп + р

+С

+ С

'

1 2 А

- 1п (2 Ар)------1п _1(2 Ар) | +

р п

1 2 А

- Кп (2 Ар) +—Кп _1 (2 Ар) р п

п)(р) =

Хп + р

п) Х{2п + 2) +р(4п + 4) ^п —• п

Хп + р

- С5п1п (2Ар) - С6п) —Кп(2 Ар). п п

В выражениях (7)-(9) 1п (р), Кп (р) — модифицированные функции Бесселя первого и второго рода соответственно. Вид частных решений, стоящих при коэффициентах —() - С(), соответствует симметричной

1 4 /-

теории упругости, а при С5 , —6 — несимметричной.

Основная проблема состояла именно в нахождении вида

последних двух частных решений для всех гармоник.

Константы при частных решениях в соотношениях (7)-

(9) определяются из краевых условий соответствующих

задач.

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

Одной из возможных экспериментально-реализуемых схем является задача о деформировании кольца (цилиндра), жестко закрепленного по внешнему контуру р = Rl (рис. 2, а). Деформирование осуществляется за счет жесткого смещения внутреннего контура Р = Ro на величину ^ В этой задаче линейных масштабов два — радиус жесткой шайбы и радиус упругого внешнего кольца. Нами в качестве характерного размера был выбран радиус внешнего контура Rl. Решение этой задачи можно получить, если в рядах (5) удержать одно слагаемое при п = 1. Для этой задачи на рис. 2, б, в представлены распределения по безразмерной радиальной координате р безразмерных компонент тензора напряжений арф и арр. Данные величины позволяют определить осевую силу реакции внутреннего контура Fy, которая может быть измерена экспериментально. На основе данной осевой силы реакции можно ввести величину 5, характеризующую отличия от силы реакции F* в классическом случае:

2л

Fy =1 (ОрФ (Ro, Ф^ш Ф+арр (Ro, ф)^ ф) Rodф,

5 =

•100%.

Рис. 2. Параметрический анализ решения задачи о деформировании кольца, Л = 0.02, Я0 = 0.1, R1 = 1. Сплошные линии соответствуют решению для среды Коссера, штрихпунктирные — для псевдосреды Коссера, штриховые — для классической теории упругости

Зависимость данной величины как от внутреннего радиуса Яо, так и от толщины сечения D = Я - Яо приведена на рис. 2, г. Хорошо видно, что отличия от классического случая увеличиваются как по мере уменьшения радиуса Яо (при фиксированном отношении радиусов), так и по мере уменьшения толщины сечения, так как оба этих фактора приводят к возрастанию влияния моментных свойств.

По аналогичному алгоритму были проанализированы две последующие задачи. Задача о вращении кольца показана на рис. 3. В этой задаче, как видно из представленных результатов, степень отклика значительно выше, чем в предыдущей задаче. Другими словами, она наиболее информативна для оценки моментных свойств в предполагаемом эксперименте, где возможным является измерение реактивного момента М — момента, необходимого для поворота шайбы на угол фо. Как видно из рис. 3, г, эта величина возрастает по мере уменьшения как радиуса Яо, так и толщины сечения шайбы и значительно превосходит подобную величину в предыдущей задаче:

2 П

М = | (арф(Яо,Ф)Яо2 + ^(Яо,Ф)Яо^Ф, о

Впервые подробный анализ этих двух задач (в альтернативной форме записи решений без использования специальных функций Бесселя) был проделан авторами

в [6].

Третьей задачей является задача Кирша об одностороннем растяжении пластины, ослабленной круговым отверстием. Эта задача наиболее привлекательна с точки зрения экспериментальной реализации, так как в ней края кругового отверстия свободны от внешних напряжений, а на бесконечности действует растягивающее усилие постоянной интенсивности в направлении горизонтальной оси (рис. 4, а). Решение этой задачи можно получить, если в рядах (5) удержать слагаемые при п = = 0 и п =2. В этой задаче весь сценарий моментного отклика, безусловно, реализуется в окрестности отверстия, где появляются градиенты напряжений и деформаций.

Впервые, в рамках классической теории упругости, эта задача была решена Киршем, позднее, несколько иным путем, Мусхелишвили [7]. Обобщение этой зада-

Рис. 3. Параметрический анализ решения задачи о вращении кольца фо = 0.1, Я = 0.1, Я = 1. Сплошные линии соответствуют решению для среды Коссера, штрихпунктирные — для псевдосреды Коссера, штриховые — для классической теории упругости

\а\

0.2

0.4

0.6

0.8

1.0

Рис. 4. Параметрический анализ решения задачи Кирша. Я0 = 0.1, р = 2. Сплошные линии соответствуют решению для среды Коссера, штрихпунктирные — для псевдосреды Коссера, штриховые — для классической теории упругости

чи на случай псевдосреды Коссера можно найти в работе [8]. Пальмовым [9] была найдена концентрация напряжений вблизи кругового отверстия в рамках несимметричной теории упругости. Следует отметить, что решение, приведенное в [9], не позволяет в полной мере проанализировать напряженно-деформированное состояние в окрестности кругового отверстия, в частности, характеризовать степень искажения отверстия при деформировании. Поэтому данная задача была подробно проанализирована авторами в [10, 11]. Некоторые результаты данного анализа приводятся ниже.

Для этой задачи можно определить два макропараметра: Dl, определяющий степень искажения первоначально круглого отверстия, и , определяющий меру искажения профиля поверхности пластины по контуру отверстия (рис. 4, б, в):

D =

ир

мр

=\е22 (Яо,о) (% V2)^

- О*

8. =

-1 оо %,

Оп

где величина е может быть вычислена в рамках гипотезы о плосконапряженном состоянии:

х-з

Х =

(Р, Ф) =------- (Урр (Р, Ф) + Уфф (Р, Ф)),

Х + 1 Х + 3р

х + р

Параметр привлекателен тем, что в классической

теории он не зависит от механических свойств материала и всегда равен 3. Как и следовало ожидать, на представленном графике (рис. 4, г) меры отклика на момент-ные свойства §1 и §2 возрастают по мере уменьшения линейного масштаба Яо.

Таким образом, рассмотренный класс задач демонстрирует возможность определения макропараметров, которые зависят от моментных свойств и, что очень важно, могут быть экспериментально измерены. Степень отклика введенных макропараметров на момент-

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

3. Численные решения статических задач несимметричной теории упругости

С целью расширения круга решаемых задач авторами был разработан алгоритм метода конечных элементов для решения двумерных статических задач несимметричной теории упругости. Для конечно-элементной реализации используется следующее вариационное уравнение несимметричной теории упругости [2]:

• • 5у + Д -8х)dV- /(X-5и + Y• 5со^V =

V V

= |(р-5и + т,

где р, т — интенсивность внешних поверхностных сил и моментов.

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

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

На основе разработанного конечно-элементного алгоритма получены численные решения ряда задач, в которых моментные эффекты проявляются более ярко, чем в задачах, для которых получены аналитические решения. Например, задача о растяжении пластины, ослабленной несколькими круговыми симметрично расположенными отверстиями (рис. 5, а). Радиусы всех отверстий равны.

Рис. 5. Задача о растяжении пластины с пятью отверстиями, где й = о.25Яо, Ь = 2о

В качестве характеристики различия решений используется величина

5 =

•100%,

где Ыу — решение в рамках несимметричной теории упругости; ы* — решение в рамках симметричной теории упругости.

На рис. 5, б представлена зависимость 5 от изменения радиусов отверстий. Из полученных результатов решения данной задачи следует, что:

1. Во всех случаях мера отклика моментных свойств материала 5 немонотонно возрастает при уменьшении радиуса отверстия Яо. Такая зависимость соответствует характеру этой зависимости и для традиционной задачи Кирша.

2. По мере сближения отверстий (уменьшения параметра Л ) величина меры 5 возрастает, особенно в точках, наиболее близко расположенных к центральному отверстию (рис. 5, б).

3. Увеличение числа отверстий усилило меру отклика моментных свойств материала 5 на выбранный макропараметр. При прочих равных условиях это усиление достигает нескольких раз.

4. Предыдущий вывод позволяет утверждать, что увеличение числа отверстий значительно повышает возможность зафиксировать отклик моментных свойств материала в физическом эксперименте.

Также решена задача об одноосном растяжении пластины с трещиной (рис. 6, а). В качестве макропараметра в данной задаче было принято перемещение Ыу по направлению оси у (направление действия силыр) в середине одного из берегов трещины (рис. 6, б). В физическом эксперименте это перемещение может быть измерено с необходимой степенью точности. Расчеты показывают, что раскрытие берегов трещины, полученное в рамках несимметричной теории упругости, больше, чем в классическом упругом решении. На рис. 6, в приведена зависимость меры отклика 5 от длины трещины Л. Эта мера возрастает при уменьшении длины трещины.

Для всех статических задач, рассмотренных в этом и предыдущем разделах, значения материальных физических констант принимались следующими:

X = 2.096 ГПа, р. = 1.033 ГПа, а = 114.8 МПа, у = 4.1 -106 Н, е = 1.312 105 Н.

4. Аналитические решения волновых задач несимметричной теории упругости

В сейсмических и геофизических исследованиях сегодня применяются (хотя и не очень широко) механические [12] или лазерные [13] сенсоры, позволяющие непосредственно измерять скорости поворотов в трех перпендикулярных направлениях. Это говорит о принципиальной возможности измерений, показывающих характер связи векторов перемещения и поворота. Та-

і і її

' х '612

Рис. 6. Задача о растяжении пластины с трещиной, где L = 20Л

ким образом, волновые эксперименты, особенно в геологических средах, могут являться очень информативными с точки зрения идентификации моделей несимметричных сред. Однако необходимой теоретической базой для подобных экспериментов являются точные аналитические решения для различных типов плоских волн. При этом наиболее подходящим для сравнения с экспериментальными данными, полученными в результате сейсмических измерений, является двумерное немонохроматическое представление волновых решений (рис. 7). Подобное представление характерно тем, что волна, которая распространяется в направлении оси х с фазовой скоростью Ср, не зависит от поперечной координаты у (хотя и имеет поперечную компоненту колебаний Ыу ).

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

При таком представлении общее решение уравнений движения представимо в виде фурье-интегралов относительно всех компонент вектора перемещения и и вектора поворота ю (физический смысл имеют здесь только вещественные части соответствующих компонент):

un(x,z,k(f ),t) = j Un(z) ei(k(f) x+ft) s0(f) df, n ^ ^ k (f),t) = J Wn (z) ei(k(f) x+ft) Sof df,

где п = {х, у, z} — координатный индекс; i — мнимая единица; к (/) — волновое число как функция круговой частоты/; t — время; ип (z) и Жп (z) — амплитудные функции, зависящие от «глубины» г.

Неизвестная функция волнового числа к (/) будет определяться далее из соответствующего алгебраического дисперсионного уравнения, которое получается для каждого типа волны при учете геометрии задачи и граничных условий. Функция 5?о( /) задает комплексный фурье-спектр сигнала-источника, который определяет форму волнового пакета и является начальным условием рассматриваемой динамической задачи.

Оправданным в данном случае является выполнение непрерывного фурье-преобразования уравнений движения (1). Это дает следующую систему относительно фурье-образов искомых компонент векторов перемещения и поворота (учитываем, что массовых сил и моментов нет):

(2р + X) grad div u - (p+a)rot rot u +

+ 2arot to + pf2u = 0,

(2y + P) grad div to - (y + e) rot rot to +

2

+ 2arot ü - (4a- jf )ю = 0.

Фурье-преобразование представления (10) дает: ->T

ü ={UX(z), Uy(z), Uz(z)} eikxS0(f),

Ю = {Wx (z), Wy (z), Wz(z)} eikxS0(f).

Подставляя (12) в (11), в силу линейности исходной системы экспоненциальный член и фурье-образ сигнала-источника сокращаются. Кроме этого, исходная система распадается на две независимые системы линейных дифференциальных уравнений второго порядка, описывающие продольную волну относительно функций Ux (z), Uz(z), Wy (z):

(|А + a)UX'(z) + (pf2 - k2(2p + X))Ux (z) +

+ i k (X + ^-a)U'z (z) -2aWy (z) = 0,

(2ß + X)U'z'( z) + (pf 2 - k 2(ß + a))Uz (z) +

+ ik (X + p-a)UX (z) + 2ikaWy (z) = 0,

(у + e)w;'(z) + (jf 2 - k2(Y + e) - 4a)Wy (y) +

+ 2<xU'x (z) - 2ik aUz (z) = 0 и поперечную волну относительно Uy (z), Wx (z), Wz (z):

(Y + e)W;(z) + (jf2 - k2 (2Y + ß) - 4a) Wx (x) +

+ ik(ß + Y-e)W4 z) - 2uU'y (z) = 0,

+ ik (ß + Y-e)Wx'( z) + 2ikaUy (z) = 0,

(p + a)U'y(z) + (pf2 - k2(p + a))Uy(z) +

+ 2aWx (z) - 2ik aWz (z) = 0.

Рис. 7. Принципиальное представление волновых решений уравнений движения

Системы (13) и (14) допускают три различных типа решений. Конкретизируя виды граничных условий, можно получить решения для затухающей с глубиной волны Рэлея [14, 15] в полупространстве, не затухающей с глубиной волны Лэмба в пластине [16], а также решения для плоских волн с не зависящей от глубины (постоянной по глубине) амплитудой. Первые два типа относятся к классу плоских поверхностных волн, в то время как третий тип соответствует плоским объемным волнам. Далее приведем подробный анализ этих волновых типов, начиная с объемных волн.

Решения для объемных продольных волн получаются из условий:

их (,) = их, и у (z) = о, и, (,) = о,

Жх (,) = Жх, Жу (,) = о, ^ ,) = о.

После подстановки данных условий в уравнения (13) и (14) получаем два независимых дисперсионных уравнения, первое из которых соответствует продольной волне перемещений, а второе — продольной волне поворотов:

(р/2 - к2(Х + 2^))их (,) = о,

(]/2 - к2(р + 2у) - 4а)Жх (,) = о.

Оправданным в данном случае является переход к безразмерной форме решений с использованием безразмерных параметров С) и С5, что дает две дисперсионные кривые:

Ш) = /, к2(/) = ,

2 _ Х + 2^ г2 _ Р + 2У .

С1 ’ С5 __о „о ’ к0 ~ 2.

рХ о/

jXо/

а

Ф + 2у

где X о — некоторый характерный размер; /0 — некоторая характерная частота. Таким образом, дополнительно к скорости продольной волны Сі необходимо ввести параметр С5, связанный со скоростью продольной волны поворотов. Кроме этого, продольная волна поворотов имеет запрещенную частотную зону, характе-

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

По аналогии можно получить решение для плоской поперечной волны:

кз(/) = 4Ар (/), к4(/) = 4Лт (/),

-*р ,т

С32 + < /2 _ 2л2 ±

2с32с4:

"(Сз2 - С42)2 /4 _

4С34С4

2Л2(С22Сз2 -2Сз2С42 + СС) /2 +2Л4

0/2

/~*2 /~*2 /~*2 С2 С3 С4

С2 -С2 - „2

С32 =

р + а

рхо^

С2 -С4 -

у + е

рХоо /о

С использованием (15), (16) приведем численный пример зависимости фазовых скоростей объемных волн от частоты, используя следующие гипотетические значения материальных параметров: X = 28 ГПа, р = 4 ГПа, р = 105 кг/м3, а = 2 ГПа, Р = 108 Н, у = 1.936 • 108 Н, е = = 3.0464 • 109 Н, j = 104 кг/м.

На рис. 8, а приведены фазовые скорости Срд(/) = = //к1(/) и Ср,2(/) = //к2(/). Помимо известной дисперсионной кривой Ср1(/) продольной волны перемещений в среде возникает независимая волна поворотов с фазовой скоростью Срд(/), обладающая дисперсией и имеющая нижнюю частоту = 2^ а/}. Вид-

но, что скорость С5 является предельной для волны поворотов.

Решение (16) имеет следующую интерпретацию. Во-первых, в случае среды Коссера у поперечной волны имеется две волновые моды с волновыми числами кз (/) и к4 (/) (в силу изотропии среды горизонтально и вертикально поляризованные поперечные волны неразличимы, поэтому по две волновые моды имеют каждая из них). Это отличает решения (16) от классического случая, когда имеется лишь одна волновая мода. Во-вторых,

Рис. 8. Численный пример, иллюстрирующий поведение объемных продольных (а) и поперечных (б) волн в среде Коссера

обе волновые моды обладают дисперсией, что видно из рис. 8, б. Одна из мод обладает нижней критической частотой, но в данном случае эта частота определяется из уравнений (16) и не равна критической частоте и продольной волны. Также можно сделать вывод о том, что безразмерные скоростные параметры С3 и С4, используемые при записи решений (16), являются асимптотическими скоростями при / объемных поперечных волновых мод.

Основным отличительным признаком поверхностных волн является то, что амплитуда волны изменяется, например затухает, с глубиной. Поэтому в отличие от предыдущего случая в решении для поверхностных волн присутствуют все шесть компонент перемещений, связанных с функциями затухания Vх (,), и у (,), V, (,), Жх (,), Жу (,), (,). Однако, как и в предыдущем слу-

чае, решение получается из независимых систем (13), (14). Поэтому оно состоит из двух частей, одна из которых описывает продольно-поперечную эллиптическую волну, а вторая соответствует горизонтально-поперечной волне, затухающей с глубиной.

Используя метод, подробно описанный в [14, 15] применительно к (13), (14), можем получить общее безразмерное решение уравнений движения (1). Необходимо только отметить, что в отличие от поверхностной волны, решение для которой рассмотрено в [14, 15], в общем случае необходимо оставить все частные решения, а не только затухающие с глубиной. Таким образом, общее решение в перемещениях имеет вид [16]:

- для продольно-поперечной эллиптической волны:

ых (х, ,,г) = | |о!ке + £2У тв +

+ Озу3е~^ + О^ке1^ -О5у2е^ --О6} ег(к+/1 ] 5о(/) /,

ы, (х, ,,г) = | ^ + О^ке +

+ 0з1ке+ D4v1eVlZ + 051ке ^ +

+Б^ке^ } еК кх+/г} ¿0( /) d/,

(

ю.

/

Сз2

2 А

+ °3,1

/ /2 Л

АР г2

Сз

з

,+vзZ I а‘

Яо(/) 4/,

для горизонтально-поперечной волны:

ыу (x, z, г) =

/2 4

А — ±__________|_____

т с42 р

+£2, , в Ъ2 +

+ Ез,6 х е

/2 4Л

Ар----------2 ^—

Р С42 р

V 4 у

ю

ю

^^о(/) У.

х(х,,,г) = | ^Е^ке+ Е2^2в +

+ Е^е^ + Е41ке^ - Е&2е^ --Е6^} е*к+/г) §о

г (х,,, г) = | ^-Е^е + Е21ке +

+ Езгке + Е^е^ + Е51ке ^ +

+Е61ке& } еК кх+/г} §о(/) /,

где величины А, В, С определяются соотношением (6), константы Еу, Б у должны быть определены из граничных условий для каждого типа волны отдельно, а экспоненциальные показатели амплитудных функций определяются выражениями:

V, =

.11

С2

^1 =.

С52 РС2

= \1к - Ат ,р , р =

В-1 А2 В'

Приведенные выше соотношения (17), (18) описывают в частном случае хорошо изученные решения для поверхностных волн в упругом полупространстве. Данные решения являются затухающими с глубиной, то есть константы при экспонентах с положительными показателями равны нулю: Е4 5 б = 045б = о. Граничные условия определяют отсутствие сил и моментов на поверхности полупространства (г = 0) и в безразмерном случае имеют вид: а2Х = ° огу = 0, ^ = 0 ^ = 0, Мгу = 0, ^= 0

Используя данные краевые условия последовательно для (17) и (18), получаем два дисперсионных уравнения. Первое дисперсионное уравнение

det(M г (У1, У2, Уз)) = о,

Мг (Pl, P2, Р3) =

2к2 -/2/С22 ~21кр2 -21крз

21кр1 2к2 - / 2 а 2к2 - /21с22

о Р2

определяет поверхностную волну с компонентами ых, ыг, Юу. Этой волне соответствует классическая волна Рэлея.

Второе дисперсионное уравнение det(Mг (^l, ^2, ^з)) = 0,

мг(Рь Р2, Рз) =

2гк 1 - В

1к (1 + С) р\

\

2 +

/ 9

г С2 - С -1

с4

V 4 J

,2 2 С52

к - Р1~^2 С

АтС| - /2 Л

2 А2С42

к С + р 22

2гкр2

Р2

2 +

АрС| - /2

к 2С + р 2

21крз

Рз

определяет поверхностную волну с компонентами Ыу , юх, ю,. Эту волну можно было бы назвать поперечной, так как компонента Ыу перпендикулярна направлению распространения волны. Геометрически эта волна подобна волне Лява, однако в классической теории упругости существование волны Лява как поверхностной волны определяется наличием слоя на полупространстве, при стремлении толщины слоя к нулю волна Лява переходит в объемную. Таким образом, в среде Коссера была обнаружена качественно новая волновая мода, аналогов которой в классической теории упругости нет.

Численный анализ полученных дисперсионных уравнений позволил установить графические зависимости волновых чисел и нормированных скоростей от частоты (рис. 9, 10). Сплошные линии обозначают решения для среды Коссера, точечные относятся к решениям, полученным в рамках классического упругого конти-

нуума, а пунктирные соответствуют скорости С2 упругой среды. Как и для объемных волн, данные решения говорят о наличии дисперсии, что является значимым отличием от классической теории упругости.

На рис. 11 приведены зависимости решений (17),

(18) от глубины, отнесенной к длине волны. Сплошная линия соответствует новой горизонтально-поперечной волновой моде, пунктирная — волне Рэлея в среде Коссера, точечная — волне Рэлея в классической упругой среде. Данная зависимость иллюстрирует, что исследуемые волны являются типично поверхностными, вся энергия которых сосредоточена в тонком приповерхностном слое, причем толщина слоя локализации зависит от частоты.

Полученное решение позволяет также проследить эволюцию исходного сигнала в виде волнового пакета. В качестве сигнала-источника 50 ( /) был использован вейвлет, физический образ которого представлен на

^ рад/с

^ рад/с

Рис. 9. Численный пример, иллюстрирующий дисперсионные свойства волны Рэлея в среде Коссера

0 1000 2000 3000 4000

Ь рад/с рад/с

Рис. 10. Численный пример, иллюстрирующий дисперсионные свойства горизонтально-поперечной волны в среде Коссера

кг/

Рис. 11. Численный пример, иллюстрирующий затухание с глубиной поверхностных волн в среде Коссера

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

Рис. 13. Трехкомпонентная сейсмограмма перемещений

рис. 12. На серии графиков, приведенных на рис. 13, представлен вид изменения исходного сигнала по мере его продвижения по поверхности упругого полупространства для трех компонент перемещения (17), (18).

5. Заключение

В данной работе получены следующие новые результаты.

Предложены подходы к построению точных аналитических решений некоторых одномерных и двумерных статических и волновых краевых задач в рамках несимметричной теории упругости для линейной среды Кос-сера.

На основе метода конечных элементов реализован алгоритм решения задач несимметричной теории упругости о статическом деформировании.

Получены аналитические и численные решения ряда статических и динамических задач.

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

Качественный и численный анализ полученных решений позволяет сделать ряд выводов о поведении упругой среды в рамках модели Коссера.

Безразмерная форма записи аналитических решений позволяет наглядно установить принципиальное разли-

чие безразмерных моментных и классических решений, а именно: безразмерное моментное решение зависит от характерного геометрического размера, а классическое — нет. Одним из интересных направлений дальнейших исследований может быть определение этих зависимостей для волновых задач.

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

Наибольшая степень возрастания имеет место в задаче о кручении. Однако с точки зрения простоты экспериментальной схемы и надежности измерений более привлекательными являются задачи о растяжении пластины с отверстиями (рис. 4, 5), так как в них не требуется задавать граничные условия для вектора поворота.

Задачи о распространении объемных волн (рис. 8) являются достаточно интересными с точки зрения интерпретации новых параметров, возникающих при записи волновых решений в безразмерном виде. Однако с точки зрения экспериментальной реализации данные задачи не представляются нам перспективными.

В отличие от объемных волн поверхностные волны для упругого полупространства несут в себе значительный экспериментальный потенциал. Во-первых, дисперсионный анализ экспериментальных трехкомпонентных сейсмограмм (полученных по аналогии с сейсмограммой, приведенной на рис. 13) позволяет получить экспериментальные дисперсионные кривые и сравнить их с теоретическими кривыми (рис. 9, 10). Во-вторых, с использованием ротационных сенсоров [12, 13] можно было бы проследить связь компонент перемещений и поворотов в соответствии с соотношениями (17), (18). В-третьих, наличие в среде Коссера поверхностной поперечной волны также может быть объектом экспериментальных исследований с применением внедренных на различные глубины сенсоров.

Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (проекты 07-01-96029-р-Урал-а, 07-01-97613-р-офи), а

также при поддержке Американского фонда гражданских исследований и развития (CRDF) в рамках грантов молодым ученым Y2-P-09-04, Y4-P-09-04.

Литература

1. Cesserai E., Cesserai F. Théorie des corps déformables. - Paris: A. Hermann et fils, 1909. - 226 p.

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

3. Grekova E.F., Herman G.C. Wave Propagation in Rocks Modeled as Reduced Cosserat Continuum // Proc. of 66th EAGE (European Association of Geoscientists and Engineers) Conf. - Paris, 2004. - P. 098.

4. Лурье А.И. Теория упругости. - М.: Наука, 1970. - 939 с.

5. Кулеш М.А., Шардаков И.Н. Построение и анализ некоторых точных аналитических решений двумерных упругих задач в рамках континуума Коссера // Математическое моделирование систем и процессов. - Пермь: ПГТУ, 2001. - № 9.- С. 187-201.

6. Кулеш М.А., Матвеенко В.П., Шардаков И.Н. Построение аналити-

ческих решений некоторых двумерных задач моментной теории упругости // Изв. РАН. Механика твердого тела. - М.: Наука, 2002.- № 5. - С. 69-82.

7. Мусхелишвили Н.И. Некоторые основные задачи математической теории упругости. - М.: Наука, 1966. - 707 с.

8. Mindlin R.D. Influence of couple-stress on stress concentrations // Experimental Mechanics. - 1963. - V. 3. - No. 1. - P. 1-7.

9. Пальмов В.А. Плоская задача теории несимметричной упругости // Прикладная математика и механика. - 1964. - Т. 28. - Вып. 6. -C. 1117-1120.

10. Kulesh M.A., Maiveenko VP., Shardakov I.N. Parametric analysis of analytical solutions to one- and two-dimensional problems in couple-stress theory of elasticity // Z. Angew. Math. Mech. - 2003. - V. 83. -No. 4. - P. 238-248.

11. Кулеш М.А., Матвеенко В.П., Шардаков И.Н. Построение и анализ точного аналитического решения задачи Кирша в рамках континуума и псевдоконтинуума Коссера // Прикладная механика и техническая физика. - 2001. - Т. 42. - № 4. - С. 145154.

12. Nigbor R.L. Six-degree-of-freedom ground-motion measurement // Bulletin of the Seismological Society of America. - 1994. - V. 84. -No. 5. - P. 1665-1669.

13. Igel H., Schreiber U., Flaws A., Schuberth B., Velikoselisev A., Cochard A. Rotational motions induced by the M8.1 Tokachi-oki earthquake, September 25, 2003 // Geophysical Research Letters. -2005. - V. 32. - P. L08309.

14. Кулеш М.А., Матвеенко В.П., Шардаков И.Н. Построение и анализ аналитического решения для поверхностной волны Рэлея в рамках континуума Коссера // Прикладная механика и техническая физика. - 2005. - Т. 46. - № 4. - С. 116-124.

15. Кулеш М.А., Матвеенко В.П., Шардаков И.Н. О распространении упругих поверхностных волн в среде Коссера // Акустический журнал. - 2006. - Т. 52. - № 2. - C. 227-235.

16. Кулеш М. А., Матвеенко В.П., Шардаков И.Н. Построение аналитического решения волны Лэмба в рамках континуума Коссера // Прикладная механика и техническая физика. - 2007. - Т. 48. -№ 1. - С. 143-150.

Поступила в редакцию 21.06.2007 г.