Напряжения в трехмерной электроупругой среде от сосредоточенного источника Текст научной статьи по специальности «Физика»

МЕХАНИКА

УДК 539.3

НАПРЯЖЕНИЯ В ТРЕХМЕРНОЙ ЭЛЕКТРОУПРУГОЙ СРЕДЕ ОТ СОСРЕДОТОЧЕННОГО ИСТОЧНИКА

© 2010 г. Л.А. Игумнов, В.П. Пазин, А.Н. Петров

НИИ механики Нижегородского госуниверситета им. Н.И. Лобачевского

[email protected]

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

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

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

Введение

Построение перемещений и напряжений в трехмерной электроупругой среде от сосредоточенного источника позволяет классическим методом граничных интегральных уравнений (ГИУ) получить решение краевой задачи трехмерной электроупругой статики, а методом ГИУ с двойным применением теоремы взаимности получить решение начально-краевой задачи трехмерной динамической теории электроупругости [1]. Ссылки на первые работы по детальному обсуждению соответствующих упругостатических решений для анизотропного случая можно найти в обзоре Bacon и др. (1978), а также в работах Mura (1987) и Ting (1996). Исследования таких упругостатических матриц электроупругости для бесконечных и полубес-конечных трансверсальных изотропных тел в замкнутой форме были проведены Wang и Zheng (1995), Dunn и Wienecke (1996, 1998), Ding и др. (1996 - 2000) и Pan (2002), а также подробно представлены в [2].

1. Постановка задач

Базовые уравнения трехмерной линейной теории электроупругости имеют следующий вид [3, 4]:

аЛ,] + Г1 = 0 А,» - Q = 0 *> У = 13

ау = Су1тг1т - еку Ек, А = еукУ]к + %]Е],

г, у, к, I, т = 1,3,

ег/ = ] + и],г )> Ег = -ф,г > 7 = !Д

где а у, Di - соответственно компоненты тензора напряжения и вектора электрического смещения; ^, Q - компоненты плотности объемной силы и плотности электрического заряда; - компоненты тензора деформации,

Еу - компоненты вектора электрического поля;

Су1т , еук и - упругие м°дули, электр°Уп-ругие коэффициенты, диэлектрические константы; ц и ф - упругое смещение и электрический потенциал соответственно. Здесь и в дальнейшем повторение индексов означает суммирование. Запятая обозначает частное дифференцирование.

Для того чтобы задачи анизотропной теории упругости и электроупругости рассматривать с единых позиций, сгруппируем базовые компоненты [5]:

л к., I = 1, 2, 3, л Г£у, г = 1, 2, 3,

4 ={ф, I = 4, ^ = {- Бу, г = 4,

а У

К-, У = 1, 2, 3,

ІА, у = 4,

Сук/

С.

ук/ ’

І, к = 1, 2, 3,

І = 1, 2, 3, к = 4,

е/к/> І = 4 к = 1 2 3

-%, і',к = 4-

е/у ’

где

\РГ і = 1,2,3,

і

=^т / (Мк (г,. (Ф)))-ч

°* = ,А

м7к = С/*/Г *1 >

где Г - расстояние между произвольной и исходной точкой.

Для получения матрицы Неймана найдем первую производную от представления из (2) [1]:

і 2п ТГ I(-Г0(МТк )-1 - ^)#

^к = (МП )-1 (М|к )-1,

по определению, исследуются свойства особых точек подынтегрального выражения и на основе теории вычетов матрица Грина представляется в виде:

С]к(х)=

А}к (р+Сщд)

-, (4)

т=1

В новых переменных физические соотношения и уравнения равновесия принимают вид:

а9 (Ст-Ст ) П(Ст-Ск )(Ст-Ск )

к=1,к&п

■е, 7 = 4.

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

2. Получение формул представлений

Для построения матрицы Неймана необходимо знать матрицу Грина. Математически матрица Грина может быть определена следующим образом:

Аук/^кт,// ( х) = -8 ут8(х), г, I, т = 1, 3; у, к = 1, 4.

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

Матрица Оук имеет следующий вид [1]:

2п

где - корни В(р + Сд) = 0 - многочлена восьмой степени от ^ ; Сш - сопряженный к ^т

о

корень; а9 - коэффициент при ^ ; 1т ^т > 0, т = 1, 4.

Для произвольной точки X = (Хі, Х2, Х3 ) производные тензора Грина О ж по координатам имеют вид:

дС 1

~РК ~^-[СРК(х1 + х2,х3) -дх-[ 2А

- СРК(Х1 - Ах25х3)]

дС

РК

дх0

2А

[СРК (х15 х2 + А х3 ) -

(5)

- СРК (х15 х2 - А х3 )]

(1)

дСРК \

дхт

(2)

(3)

Щ'

0 0 В = Мг 2 + м2Г

^п# -/Г1-пд 1 -/Г1-п£ 5

00

0 2 Г'' 0

М,к = Сук12 у г! > М,к = Сук1 Г г! ■

Второй способ опирается на интерполяцию Лагранжа [4]: обратная матрица записывается

— [СРК (х1, х2, х3 + А) -

2А

- СРК(х15х25х3 - А)]-Величина шага И определяется из соотно-

_3

шения: к = г -10 .

3. Численные эксперименты

Для удобства тензор упругих модулей С,-^ запишем через константы сар (а = 1, 6; в = 1, 6) [6]. Электроупругие константы могут быть

записаны с двумя индексами екр(к = 1,3;р = 1, 6) следующим образом:

ек1 = ек\\; ек2 = ек22; ек3 = ек33; ек 4 = ек 23 = ек 32; ек 5 = ек13 = ек 31; ек 6 = ек12 = ек 21В качестве примера рассмотрим трансвер-сальный изотропный электроупругий материал, для которого известны точные решения в замкнутой форме [7]. Предположим, что ось симметрии материала параллельна оси х3, тогда отличные от нуля элементы констант материала следующие:

4

0

0

с11, с22 (= С11), с33’ с13’ с23 (= с13 )’ с44, с55(= с44)> с66> с12(= С11 - 2с66) > е31,е33,е15,е32(= е31),е24(= е15) ,

£11, е22(= Е11), Е33 ■

Для трансверсально изотропной пьезокерамики (Р2Т-4) [8] соответствующие константы материала приведены в табл. 1, упругие константы сар даны с точностью до множителя в 109 Н/м2, пьезоэлектрические коэффициенты

е^р - Кл/м2 и диэлектрические постоянные - 10-9 Кл/Вм.

Для точки х = (1,1,1) компоненты элек-

троупругих матриц Неймана ОрК 1, ОрК 2 ,

ОрК з представлены соответственно в табл.

2-4 в сравнении с точным решением [7]. В табл. 3 также представлено сравнение с численным решением Б. С. Вигош и А. 8аег

[9].

Таблица 1

С11 с13 с12 3 3 О О с66 е31 3 3 е15 811 833

139.0 74.3 77.8 115.0 25.6 30.6 -5.2 15.1 12.7 6.4605 5.61975

Таблица 2

(Р, к) Точное решение Численное решение

1,1 - 9.26183802345 х 10-14 - 9.26211971502 х 10-14

1,2; 2,1 1.65932816546 х 10-14 1.65938011817 х 10-14

1,3; 3,1 1.82113322728 х 10-14 1.82111361861 х 10-14

1,4; 4,1 6.70402857710 х 10-6 6.703 86239193 х 10-6

2,2 - 4.81188998237 х 10-13 - 4.81190985218 х 10-13

2,3; 3,2 -1.54203110547 х 10-13 -1.54203820596 х 10-13

2,4; 4,2 - 2.01554138734 х 10-4 - 2.01554974395 х 10-4

3,3 - 3.30447677887 х 10-13 - 3.3 0448786793 х 10-13

3,4; 4,3 6.20021679052 х 10-4 - 6.20023799701 х 10-4

4,4 1.18596300519 х 106 1.18597183363 х 106

Таблица 3

(Р, к) Точное решение Численное решение Вигош Численное решение

1,1 - 4.81188998237 х 10-13 - 4.81189369756 х 10-13 - 4.81164650629 х 10-13

1,2; 2,1 1.65932816546 х 10-14 1.65929102010 х 10-14 1.65912699071 х 10-14

1,3; 3,1 -1.54203110547 х 10-13 -1.54205935560 х 10-13 -1.54199292809 х 10-13

1,4; 4,1 - 2.01554138734 х 10-4 - 2.01546078340 х 10-4 - 2.01549777817 х 10-4

2,2 - 9.26183802345 х 10-14 - 9.26187902596 х 10-14 - 9.25985855923 х 10-14

2,3; 3,2 1.82113322728 х 10-14 1.82108618017 х 10-14 1.82182962166 х 10-14

2,4; 4,2 6.70402857710 х 10-6 6.70401912700 х 10-6 6.71225201843 х 10-6

3,3 - 3.30447677887 х 10-13 - 3.30458722395 х 10-13 - 3.30436364465 х 10-13

3,4; 4,3 - 6.20021679062 х 10-4 - 6.19996630452 х 10-4 - 6.20000274462 х 10-4

4,4 1.18596300519 х 106 1.18591339858 х 106 1.18588086747 х 106

Таблица 4

(р, К) Точное решение Численное решение

1,1 - 5.76929903985 х 10-13 - 5.76904645572 х 10-13

1,2; 2,1 - 2.27471872311 х 10-13 - 2.27474197093 х 10-13

1,3; 3,1 - 3.64226645455 х 10-14 - 3.64183817407 х 10-14

1,4; 4,1 -1.347080571542 х 10-5 -1.34032007654 х 10-5

2,2 - 5.76929903985 х 10-13 - 5.76909046396 х 10-13

2,3; 3,2 - 3.64226645455 х 10-14 - 3.64152995171 х 10-14

2,4; 4,2 -1.347080571542 х 10-5 -1.33994666571 х 10-5

3,3 -1.36396808876 х 10-13 -1.36385029260 х 10-13

3,4; 4,3 - 2.83889208291 х 10-4 - 2.83867035216 х 10-4

4,4 2.01966683615 х 106 2.01958092858 х 106

4. Визуализация

Опишем построение матриц Неймана в виде двумерных поверхностей на основе интерполяционной схемы вычисления. Функции Оук і (х) зависят от трех пространственных переменных Хі,Х2 и Х3. Перейдем к сферическим координатам (г, 0І5 02). Тогда на единичной сфере ; (®1, ) - функции только двух перемен-

ных (полярного угла 0 <0! < п и азимутного угла 0 <02 < 2п), поэтому функция Од і от

конкретных 0^ и 02 дает соответствующие значения во всем электроупругом трехмерном пространстве с точностью до значения множителя г. Из дискретных значений с помощью интерполяции Лагранжа строится вариант непрерывного поведения искомой матрицы [1].

На рис. 1-12 представлены компоненты электроупругой матрицы Неймана для транс-версального изотропного пьезокерамика (Р2Т-4). Для графито-эпоксидного моноклинического материала параметры материала приведены в табл. 5, а на рис. 13-18 представлены компоненты анизотропной матрицы Неймана.

Таблица 5

С11 СЧ СЧ О СП СП О с12 с13 3 2 О О О с45 с66 с16 с26 с36

95.5 25.9 16.3 28.9 4.03 4.65 4.40 6.45 -1.78 32.7 44.7 15.6 0.54

Рис. 3. Компонента Оц 3

Рис. 5. Компонента 2

•12 ................ :

и ю ...--г"'' : | '

о о

Рис. 7. Компонента 6^2 1

■12

о о

Рис. 9. Компонента 6^2 3

Рис. 4. Компонента 1

Рис. 6. Компонента 6^2 3

•12

х 10

э

о о

Рис. 8. Компонента 0^22 2

7

к 10

О о

Рис. 10. Компонента 644 1

Рис. 11. Компонента G44 2

Рис. 13. Компонента G^ 1

Рис. 15. Компонента Оц 3

Рис. 17. Компонента G22 2

Рис. 12. Компонента G44 3

Рис. 14. Компонента Gn 2

о о

Рис. 16. Компонента G22 і

Рис. 18. Компонента G22 3

Заключение

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

Работа выполнена в рамках реализации ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009-2013 годы и гранта РФФИ № 10-08-01017-а.

Список литературы

1. Баженов В.Г., Игумнов Л.А. Методы граничных интегральных уравнений и граничных элемен-

тов в решении задач трехмерной динамической теории упругости с сопряженными полями. М.: Физмат-лит, 2008. 352 с.

2. Qing-Hua Qin. Green’s function and boundary elements of multifield materials. Elsevier, 2007. 254 p.

3. Партон В.З., Кудрявцев Б.А. Электромагнито-упругость пьезоэлектрических и электропроводных тел. М.: Наука, 1988. 470 с.

4. Pan E., Tonon F. Three-dimensional Green’s functions in anisotropic piezoelectric solids // International Journal of Solids Structures. 2000. 37. P. 943-958.

5. Li X., Wang M. Three-dimensional Green’s functions for infinite anisotropic piezoelectric media // International Journal of Solids Structures, 2007. 44. P. 16801684.

6. Ting T.C.T. Anisotropic elasticity - theory and applications // Oxford University Press. New York, 1996.

7. Dunn M.L., Wienecke H.A. Green's functions for transversely isotropic piezoelectric solids // International Journal of Solids Structures. 1996. 33. P. 4571-4581.

8. Dunn M.L., Taya M. An analysis of piezoelectric composite materials containing ellipsoidal in homogeneities // Proc. R. Soc. Lond. 1993. A443. P. 265-287.

9. Buroni F.C., Saez A. Three-dimensional Green's function and its derivative for materials with general anisotropic magneto-electro-elastic coupling // Proc. R. Soc. 2010. A466. P. 515-537.

STRESSES IN A THREE-DIMENSIONAL ELECTRIC ELASTIC MEDIUM CAUSED

BY A LUMPED SOURCE

L.A. Igumnov, V.P. Pazin, A.N. Petrov

Infinite 3D anisotropic and electric elastic media are considered. Two approaches to obtain matrix solutions on the basis of a unified formalization have been presented. The first one permits to obtain matrix solutions in the form of an integral representation, whereas the second one, in the form of an analytical formula using root sets of functional solution matrices. The components of the obtained representations of matrix solutions have been numerically built for specific material properties. The results have been compared with those available elsewhere. Visualization of matrix solutions in the form of surfaces has been presented.

Keywords: Neumann matrix, 3D formula, anisotropic elasticity, electric elasticity, numerical simulation (modeling).