Колебания неоднородного пороупругого слоя с пустыми порами Текст научной статьи по специальности «Математика»

УДК 539.3

Колебания неоднородного пороупругого слоя с пустыми порами1

А. О. Ватульян, Д. С. Шведов

(Донской государственный технический университет)

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

Ключевые слова: пороупругость, система интегральных уравнений Фредгольма второго рода, функция относительного объёма.

Введение. Теория пористых линейно упругих материалов с пустотами первоначально была разработана Коуином и Нунциато [1]. Она предназначена для моделирования упругих материалов, содержащих распределение малых пор (рис. 1). Главной особенностью данной теории является введение в определяющие соотношения новой переменной, характеризующей относительный объём пор, который берётся в качестве независимой кинематической переменной. Включение новой переменной требует дополнительных сил для обеспечения равновесия объёма пор. Если объём пор обращается в ноль, поведение материала описывается классической теорией упругости.

Рис. 1. Упругая среда с распределёнными порами

Теория Коуина—Нунциато хорошо описывает поведение композитных, керамических, пенистых и пористых материалов. При этом ряд функционально-градиентных материалов отличается тем, что распределение физических характеристик неравномерно, что приводит к их зависимости от координат.

Настоящая работа посвящена исследованию колебаний неоднородного по толщине пороупругого слоя с пустыми порами.

1 Работа выполнена в рамках инициативной НИР.

Постановка задачи. Рассмотрим двумерную задачу об установившихся колебаниях с частотой О. неоднородного изотропного пористого упругого слоя 0 < у< /7, |л1 < оо. Нижняя грань слоя горизонтальна и сцеплена с абсолютно жёстким основанием.

Необходимо найти поля перемещений в слое и (х,у) = {и,у}еи функцию относительного объёма <&(х,у= ,у)е~'м при известных законах неоднородности под действием

нагрузок, приложенных к верхней границе слоя.

Считая режим колебаний установившимся, отделим временной множитель е_,£“. Будем рассматривать амплитудные значения функций.

Уравнения движения имеют вид [2]:

Компоненты тензора напряжений находятся из следующих определяющих соотношений [2]:

где Л = Л(у), ^ = ^(у) — положительные функции координаты у, являющиеся аналогами параметров Ляме, р — массовая плотность материала; а, (3, со, к — положительные физические параметры, связанные с пористостью [1].

Система интегральных уравнений Фредгольма 2 рода. Отметим, что в представленной постановке колебания описываются системой уравнений в частных производных второго порядка с переменными коэффициентами. Наиболее эффективным средством анализа колебаний в слое является использование преобразования Фурье по переменной хи сведение к системам обыкновенных дифференциальных уравнений с переменными коэффициентами [3].

(1)

(2)

оху =м(и; +у’х),

Полевые характеристики задачи должны быть подчинены следующему уравнению [2]:

—г-н-----------------------------------------г-н-Ф —■ —0.

дх ду а

Соответствующие граничные условия для слоя могут быть записаны в виде:

д2Ф д2Ф /СОЙ + рко2 - 5

(3)

у(*,0) = 0, \оху (*,/?) = т(х),

и(х,0) = 0, \оуу(хгЬ) = о(х),

(4)

Функция

СО

—со

называется прямым преобразованием Фурье f(x). В свою очередь, функция

называется обратным преобразованием Фурье и с физической точки зрения это выражение представляет собой интегральную запись принципа суперпозиции гармоник с комплексной амплитудой и волновым числом 5 с непрерывным спектром.

Системы (1), (2) и уравнение (3) в образах Фурье выглядят следующим образом:

+-^°ху +9&й = 0, ~‘5дху+^дуу+р^ = 0-

ахх = -'яй (Л + 2^) + (1 - 2 с2) (Л + 2^) —V + НФ (Л + 2^), ауу = (1 - 2с2) (Л + 2^) (-/5) й + (Л + 2^) ^-\/ + НФ (Л + 2^),

ох=Щ —и-ву \ = \л — и -\JilSV.

с!у

й2 Л Г

-Ф +

С1у2

с!у

/СОФ + рАФ2-?; ^2

^ зГ^ - - -1 п

Ф-— -------V — /5£/ | = 0.

а {&

(5)

(6)

(7)

Если из (5)—(7) выразить функции —о„/,—д1П/,—й,—\/,-^-Ф, исключить й , а также

с1у у с1у " с1у с1у б у

ввести —Ф = Ф, получается каноническая система обыкновенных дифференциальных уравнений ду

1-го порядка с переменными коэффициентами:

с! ~ Л _

—о = /5-а -

Оу ху Л + 2и

рф2 - 45'

р(А + и)

Л + 2^

й + /52 \}НФ,

с1

с1уиуу

и = -0Х¥ + /5Н,

1 ~ Л ~ игк

V =-.---—о,,,, + /5----—и - НФ,

йу

с!у Л + 2\} |_ф = ф, с!у

3

Л+ 2^

с1у о{Х + 2^)'уу ' аЛ + 2^

Граничные условия в образах Фурье имеют вид:

■ — 5 +^Н а

Ф.

£7(5,0) = О, V (5,0) = О, с/Ф,

йу

(5,0) = О,

ауу(5,/7) = а(5),

сЮ>

бу

(5,/») = 0.

(8)

Краевая задача для системы функций (8), (9) после введения безразмерных переменных и параметров приведена к следующей безразмерной форме:

После замены для получения вещественных коэффициентов в вырожденной системе (со = 0)

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

Мо = М(/?), Р0 =Р(Л), Т0(5)

М -У\/ Юху ~(2\! Оуу -Qг^

(10)

она выглядит следующим образом:

и, - г/эдДф,

(11)

Граничные условия:

и, (б,0) = 0, Г<?! (^,1) = /т0 (^),

•1/(5,0) = 0, <С2(5,1) = а0(5),

ф(^,0) = 0. |ф(5,1) = 0.

(12)

^2=1 (ед (п) + (п) + С3ф(п))^п - /Л5Т0 (1 - г) + о0,

О

1

<4 = -/(^2 (п) + ВД (П) + 66ф(п))^п + /т0иг (г),

О

1

Ф = | (с7с>2 (п) + ВД (п) + 66Ф (п))^п,

С, =с^,ч\), /'= 1-9-

Ядра системы имеют вид:

где

(г,п) = -52/7

2 л, №

+ 2ц 1 М0 + 2ц 2

б,2(^л) =

453/7:

^1 + ^1) _ р02^/73

Л! + 2ц

и0

/С1 +

Ф2/72

К21

Ф25/73 Л,

М0 Л1 + 2М1

-/С,

(?з (г,п) = -2^2/?2М1А/1АГ1 --—Н.К,,

64 (г, п) = -/75 ' \ /сз + /75 у \ К,,

Л! + 2ц

^5(^п) =

2А.2 Л

Л! + 2^

452/7:

^(л1+^) рф2/? Л! + 2^

/С3 + 5/7

2 А»2 "1

Л! + 2^

/С,

£6 (г, п) = -г^/У^з - /75/7^4,

^7(г'П) =

р/72

а (Л,+2^)

к.*1-**-*'

£8(г,Г|) = -2 5

Р/73 И, а Л! + 2[^

£9(*,п) = -

-/72^2 +М1/У,

К; =К,(г,п), / = 1..6; \ =Л1(п), ^ =^(п),

/С1(г,п) = тах{г,п}-г,

1

К2(г, п)= | р£/<;,

тах{2,Г|} т1п{^,п} 1

*зМ) = |

о М1

/С4(г,п) = г-тт{г,п},

2

Къ{г, п)= | 0(у)*/о,

т1п{2/П}

1

^б(^п)= | 0(у)*/о,

/шФ + рАФ2 - ^

- /7252 |с/(;,

и/

О М1

с^П-

(13)

Численный метод решения системы интегральных уравнений и обращения преобразования Фурье. Осуществим дискретизацию системы интегральных уравнений Фредгольма 2 рода (13), используя квадратурную формулу трапеций [4, 5] и метод коллокаций. Придём к следующей алгебраической системе:

&х = —, х =^(/-1), / =1,...,/7,

л-1 у ’

А1=Ап=*±,А2=А3=... = Ап_1=Ьх1

Оц - X (СЛ + С20и1] + ^Ф; ) = а0 - //?5Т0 (1 - */ ) ,

;=1

<4/ + Е Л (^2; + ^5^1; + ^Ф; ) = /Т0^ ,

;=1

Ф/ “ X Л [С7^2, + ^8 Д; + ^Ф; ) = 0'

(14)

;=1

Значения остальных функций в соответствующих узловых точках находятся из инте-

гральных представлений:

оАг) = ітй-\

К + 2и,

452/72

Мі (Аі + Мі) рФ2/72

К + 2Мі М0

л

ФИ = 1

а(Л1+2м1) а Лх + 2мх ^ а а

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

1 “

/Г(х) = ^1/Г(5)е^5'

2п _

используем интерполяционную формулу третьей степени, являющуюся аналогом «правила трёх восьмых» Ньютона-Котеса [5, 6]. Оно получается при замене функции ^ на каждом из отрезков [0,3/7Х],[3/^,6/^],... алгебраическим многочленом третьей степени, интерполирующим ^ по четырём значениям.

Возьмём четыре точки 5к,5к+1,5к+г,5к+ъ и выполним интерполирование ^ по её значениям в этих точках:

(5-5кл)(5-5к 2)(5-5к з) ц | (5-5,)(5-5,+2)(5-5,+з)^ ^ +

-/7, (-2/7,) (-3/7,)

-/7і(-/7і)(-2/7і)

| (5-5,)(5-5,:1)(5-5,+3)5 ^ | (5-5,)(5-5,+1)(5-5.+2)/г^ | ^

3^ • 2/7, • ^

2/7, - /7, (-/7,)

Умножение этого равенства на е~15Х и интегрирование по отрезку + 3/7Х] с последующим

суммированием по значениям к = 0,+3,±6,... приводит к представлению:

(У3+/б3)

2п

2п

У Р е ,3

Зк-1

К(х)

+ — У Рзке~і3~ +-п 2п

(15)

где

в - Т, )3 Е & - т,) -1^ & - 2) & - 3) (1 - т, )3 Е (1 - т,) +

+1^ & -1) в - 3) (2 - т, )3 Е (2 - т,) - и & -:1) & -,2) (3 ■- т, )3 ± ^4> (зк + /»1Т1) е-'^‘ +л^ >*,

к=-со

^(*)=

1,5 > О, 1 2‘

0,5 < 0.

-,5=0,

Коэффициенты а3,у3,53 имеют вид:

0 = хЬх,

/7, ха3 = ^0 2 - 0 4 + [е 4 “ ^9 2)СОБ30 + 0 3 ап 30,

= 30 4 - 30-2 +3^0 2 -0^ак30-40-3 вт30,

Ь~% =50 3 +3^0 2 -0 4^51п30 + 40 3соэ30.

Стоит отметить, что численный эксперимент проводился для п = 30 узловых значений при точности е = 0,001.

Расчёт полей смещений и функции относительного объёма при некоторых законах неоднородности. В качестве примеров реализации предложенной схемы рассмотрим возрастающие и убывающие функции, характеризующие законы неоднородности в случае сосредоточенной и распределённой нагрузки.

Для сосредоточенной нагрузки находим:

СО

а0 (х) = б(х); а0 (5)= | 5(х)е/5Хс1х = 1,

—со

СО

т0 (х) = б(х); т0(5)= | 6(х)еЬх(Ьс = 1.

— со

Для равномерно распределённой на отрезке нагрузки получим: к(,)=я ( ^

[ст0 (л-) = 0, \х\ > а у 5

(т0(х) = т0, |*|<а апф)

Ьм=о. М>* ' ° 5 '

где а — полудлина области приложения нагрузки.

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

= |м1(г)сз!? =1, р(г) = 1. (16)

о о

Согласно физическому смыслу параметры Ляме должны удовлетворять условию:

0<

К (2) 1 г и

Л Vz е [0; 1].

2(Л1(г) + м1(г)) 2

(17)

Рассмотрим некоторый закон, удовлетворяющий условиям (16), (17):

A1(z) = z + 0,5, |j1(z) = l,5-z.

Пусть h =0,5; а = 0,01; 3 = 10; со = 0,01; ^ = 15; к =0,01; Q = 100 , а также

а0 (s) = 1; т0 (s) = 0.

Ниже приведены графики полей перемещения и функции относительного объёма. Сплошной линией обозначена действительная часть, пунктирной — мнимая часть.

Re(U) —■— Im(U) | |------- R.e(V)------------Im(V) |

Re (Phi)--------Im(Phi)

Рис. 2. Поля перемещений в слое и функция относительного объёма

Заключение. Разработан метод решения краевой задачи для неоднородного пороупругого слоя, позволяющий находить поля внутри и на поверхности для широкого класса неоднородностей. Библиографический список

1. Cowin, S. С. Linear elastic materials with voids / S. C. Cowin, J. W. Nunziato // J. Elasticity. — 1983. — V. 13, № 2. — Pp. 125—147.

2. Скалия, А. Гармонические колебания жёсткого штампа на пористо упругом слое /

А. Скалия // Прикладная математика и механика. — 2009. — Т. 73, вып. 3 — С. 484—491.

3. Бабешко, В. А. Динамика неоднородных линейно-упругих сред / В. А. Бабешко,

Е. В. Глушков, Ж. Ф. Зинченко. — Москва : Наука, 1989. — 343 с.

4. Манжиров, А. В. Методы решения интегральных уравнений: справочник / А. В. Манжи-ров, А. Д. Полянин. — Москва : Факториал, 1999. — 272 с.

5. Крылов, В. И. Приближённое вычисление интегралов / В. И. Крылов. — Москва : Наука, 1967. - 500 с.

6. Крылов, В. И. Методы приближённого преобразования Фурье и обращения преобразования Лапласа / В. И. Крылов, Н. С. Скобля. — Москва : Наука, 1974. — 600 с.

Материал поступил в редакцию 28.11.2012.

References

1. Cowin, S. С., Nunziato, J. W. Linear elastic materials with voids. J. Elasticity, 1983, vol. 13, no. 2, pp. 125-147.

2. Skaliya, A. Garmonicheskie kolebaniya zhestkogo shtampa na poristo uprugom sloe. [Harmonic vibrations of a rigid stamp on porous elastic layer.] Prikladnaya matematika i mexanika, 2009, vol. 73, iss. 3, pp. 484-491 (in Russian).

3. Babeshko, V. A., Glushkov, E. V., Zinchenko, Z. F. Dinamika neodnorodny'x linejno-uprugix sred. [Nonhomogenous linear-elastic media dynamics] Moscow : Nauka, 1989, 343 p. (in Russian).

4. Manzhirov, A. V., Polyanin, A. D. Metody' resheniya integral'ny'x uravnenij : spravochnik. [Integral equation solution methods: reference guide.] Moscow : Faktorial, 1999, 272 p. (in Russian).

5. Krylov, V. I. Priblizhennoe vy'chislenie integralov. [Integral approximate computation.] Moscow : Nauka, 1967, 500 p. (in Russian).

6. Krylov, V. I., Skoblya, N. S. Metody' priblizhennogo preobrazovaniya Fur'e i obrashheniya preobrazovaniya Laplasa. [Fourier approximate transformation methods and Laplace transformation inversions.] Moscow : Nauka, 1974, 600 p. (in Russian).

OSCILLATIONS OF INHOMOGENEOUS POROELASTIC LAYER WITH VOIDS1

A. O. Vatulyan, D. S. Shvedov

(Don State Technical University)

As part of the plane strain, the problem on steady oscillations of the inhomogeneous through-thickness isotropic poroelastic layer with voids is considered. The layer bottom is connected with the perfectly rigid foundation, the oscillations are caused by the surface stress on the top face. Under the known inhomogeneity laws for Lame parameter analogs - positive functions of the vertical coordinate - the layer oscillations are described by the system of three partial differential second-order equations with variable coefficients. Using the Fourier transform and some access statements, the problem is reduced to the system of three Fredholm integral equations of the second kind with continuous kernels. The numerical technique is proposed for identifying the transforms by the collocation method. The displacement vector and the relative volume are calculated using numerical inverse Fourier transform. Keywords: poroelasticity, system of Fredholm integral equations of the second kind, relative volume function.

1 The research is done within the frame of the independent R&D.