Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Пьезокомпозиты и датчики в 3 ч. Ч. 1 Статистическая механика пьезокомпозитов

.pdf
Скачиваний:
6
Добавлен:
12.11.2023
Размер:
4.82 Mб
Скачать

личин uz , ,

,

rz ,

Dr , Br

 

при r = a и условий для осредненных

величин

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E = −v C (1v )C = E* 0 ,

 

 

 

 

 

 

 

 

x

 

 

1

 

(1)

 

 

1

(2)

 

 

x

 

 

 

 

 

 

 

 

 

H

x

= −v E

(1v )E

 

= H * 0 ,

 

 

(2.76)

 

 

 

 

 

 

 

 

1 (1)

 

 

1 (2)

 

x

 

 

 

 

 

 

 

 

 

 

xz

= v1 A(1)

+ (1v1 )A(2) = 0

 

 

 

 

с учетом преобразований компонент

 

 

 

 

 

 

 

 

xz = rz cos z

sin ,

 

Ex

= Er cos E sin ,

Hx

= Hr

cos H sin

В результате получим решение для констант

 

 

 

 

 

 

 

 

 

 

 

A

 

 

= b E*

+b H

* ,

 

 

 

 

 

 

 

 

 

 

 

 

 

(1)

 

 

11 x

 

12

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

 

= b E*

+b H *

,

 

 

 

 

(2.77)

 

 

 

 

 

 

 

 

(1)

 

 

21 x

 

22

x

 

 

 

 

 

 

 

 

 

 

 

 

 

C(1)

 

= b31Ex* +b32 Hx* ,

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b =

 

 

1

,

b =

 

 

2

,

b = b ,

b = b

 

,

11

 

2

 

0

 

12

 

2

 

0

21

1

11

 

22

1 12

 

2

 

 

 

 

 

 

 

b31 = 1b11 2 ,

1

=

e(1)113 (1v1 ) + e(2)113 (1+ v1 )

 

,

 

 

 

 

 

 

(1)11

(1v ) +

(2)11

(1+ v )

 

 

 

 

 

1

1

 

 

1

=

h(1)113 (1v1 ) + h(2)113 (1+ v1 )

,

 

 

 

 

 

(1)11

(1v ) +

(2)11

(1+ v )

 

 

 

 

 

 

1

 

1

 

 

b32

 

= 1b12 ;

 

 

 

 

 

 

2 =

 

 

 

 

 

2 (2)11

 

 

 

,

 

 

(1)11

(1v ) +

(2)11

(1+ v )

 

 

 

 

1

 

1

 

 

2 =

 

 

 

 

 

2 (2)11

 

 

 

,

 

(1)11

(1v ) +

(2)11

(1+ v )

 

 

 

 

1

 

1

 

 

1+ v ,

0 = 3 + 1v1 4 + 1 1 + 1 2

1

101

 

 

 

= e

(1)11

e

 

,

 

 

= h

 

 

 

(1)11

h

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

(1)113

 

(2)11

(2)113

 

 

2

(1)113

 

 

(2)11

(2)113

 

 

 

 

 

 

= C

 

e

e

 

 

h

h

 

 

,

 

= C

 

 

 

e2

 

 

h2

 

+ (1)113

(2)113 +

(1)113 (2)113

 

 

+ (2)113 +

(2)113

 

 

3

(1)1313

 

(2)11

 

 

(2)11

 

 

 

 

 

4

 

 

(2)1313

 

 

(2)11

 

(2)11

 

Осредним оператором ...

 

по объему составного цилиндра

величины

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Dx = Dr

cos D sin ,

Bx = Br

cos B sin ,

xz

 

= rz

cos z sin

и с учетом (2.57), (2.73) – (2.76) получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Dx = (2)11Ex* + v1 (e113 A(1)

 

11C(1) ) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= (2)11Hx* + v1 (

 

113 A(1)

11E(1) ) ,

 

 

(2.78)

 

 

 

 

 

Bx

h

 

 

xz = −e(2)113Ex* h(2)113Hx* + v1 (C1313 A(1) + e113C(1) + h113 E(1) ).

Из формул (2.78) с учетом разложений (2.77) и равенств

Dx* = Dx , Bx* = Bx , *xz = xz

следует вид определяющих соотношений на макроуровне композита: а) для индукции электрического поля

 

D* = *

E* + * H

* ,

 

 

(2.79)

 

x

 

11

x

 

11

 

x

 

 

где эффективные модули

 

 

 

 

 

 

 

 

 

 

 

 

 

*

 

 

+ v (e b

 

b ),

 

=

 

 

(2.80)

11

 

(2)11

 

1

113

11

 

11

31

 

 

 

 

 

*

= v (e b

 

b );

 

 

 

 

 

 

 

11

1

113

12

11

32

 

 

б) для индукции магнитного поля

Bx* = 11* Hx* + 11* Ex* ,

(2.81)

102

где

11* = v1 (h113b11 11b21 ),

11* = (2)11 + v1 (h113b12 11b22 );

в) для макронапряжения

*xz = −e113* Ex* h113* Hx* ,

где

e113* = e(2)113 v1 (C1313b11 + h113b21 + e113b31 ), h113* = h(2)113 v1 (C1313b12 + h113b22 + e113b32 ),

(2.82)

(2.83)

(2.84)

тензоры разностей

C = C1 C2 , e = e1 e2 , h = h1 h2 ,

λ = λ1 λ2 , μ = μ1 μ2 ;

определяющие соотношения (2.70), (2.79), (2.81), (2.83) записаны для рассматриваемых частных случаев нагружений композита:

Ex* 0 , Hx* 0 .

На рис. 2.1, а, б, представлены результаты () расчета эффективных коэффициентов электромагнитной связанности 11* , *33 во-

локнистого пьезоэлектромагнетика в зависимости от содержания пьезоэлектрических PVDF-волокон v1 в ферритовой (2.11)–(2.13)

матрице с полидисперсной структурой (см. рис. 2.1, а). Результаты () на рис. 2.1, а, б, соответствуют случаю инверсии свойств 1-й и 2-й фаз (см. рис. 2.1, б), т.е. когда матрица – это пьезоэлектрик PVDF, а волокна – феррит; для наглядности сопоставления графиков здесь по-прежнему через v1 обозначено относительное объем-

ное содержание PVDF в композите. Отметим, что аналитическое решение *33 (2.71) () для композита с полидисперсной структурой

(см. рис. 2.1, а) в точности совпало с решением асимптотического метода осреднения [11]

103

v1

 

v1

 

 

 

 

а

 

 

 

 

 

 

б

 

Рис. 2.1. Коэффициенты электромагнитной связанности

 

*

(а), *

(б) композита в зависимости от объемной доли

 

11

33

 

 

 

 

 

 

 

 

пьезоэлектрической фазы v1

в виде волокон (□) или матрицы (Δ)

 

 

 

v1 (1v1 )e311

 

311

 

 

 

 

*33 =

h

(2.85)

 

 

 

 

 

 

 

 

 

 

k(1)12 v1k12 + G(2)12

 

 

 

 

для идеальной периодической волокнистой структуры. Результаты расчета на рис. 2.1, а, б, позволяют сделать вывод, что инверсия свойств волокон и матрицы композита приводит к значительному увеличению абсолютных значений эффективных коэффициентов

электромагнитной связанности 11* , *33 композита при фиксированных величинах объемного содержания пьезоэлектрической v1 и пьезомагнитной 1v1 фаз. В частности, при объемной доле v1 0,2 для PVDF и 0,8 для феррита абсолютные значения *33 для композита с ферритовыми волокнами в PVDF-матрице более чем в два раза превышают *33 для композита с PVDF-волокнами в фер-

ритовой матрице, т.е. PVDF предпочтительнее использовать в качестве матрицы композита.

Отметим, что приведенные в источнике [22] решения для эффективных модулей гибридных однонаправленных волокнистых полидисперсных композитов нельзя считать точными решениями в

104

рамках традиционных полидисперсных моделей (см. [25]), так как в гибридной полидисперсной структуре [22] наличие различных типов составных ячеек с различными по свойствам волокнами и с едиными напряжениями на внешних поверхностях всех ячеек приводит к невыполнению совместности деформаций с выполнением уравнений равновесия для напряжений в представительной области V нагруженного композита; нарушение сплошности в окрестностях внешних границ ячеек полидисперсной среды происходит изза отсутствия однородности в области V осредненных по разнотипным ячейкам деформаций.

Таким образом, получены новые аналитические решения (см. (2.71), (2.80), (2.82), (2.84)) и точные соотношения (см. (2.72)) для

эффективных констант: k*

,

С*

 

,

С*

, e*

, e*

, h*

, h*

, *

,

 

 

 

 

 

 

 

 

 

 

12

 

3333

 

1133

333

311

311

333

11

 

*

,

*

,

*

,

*

,

*

,

*

,

*

однонаправленного волокнистого

33

 

11

 

33

 

11

 

33

 

11

 

33

 

 

 

 

 

 

 

 

 

композита с полидисперсной структурой. Проведен численный расчет и анализ влияния на эффективные константы композита с пьезоэлектрическими волокнами PVDF в ферритовой пьезомагнитной матрице величины объемного содержания волокон v1 и инвер-

сии свойств фаз в сравнении с известным решением [11] асимптотического метода осреднения для *33 идеальной периодической

структуры. Полученное аналитическое решение *33 (2.71) для

композита с полидисперсной структурой в точности совпало с решением асимптотического метода осреднения [11] для идеальной периодической волокнистой структуры. В волокнистом пьезоэлектромагнетике пьезоэлектрик PVDF предпочтительнее использовать в качестве матрицы, а феррит – волокон композита, что позволяет значительно увеличить абсолютные значения эффективных коэф-

фициентов электромагнитной связанности 11* , *33 по сравнению со случаем инверсии свойств у фаз композита (см. рис. 2.1, а, б).

2.4.2. Самосогласованные решения*

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

* Принимая во внимание данные исследования [47].

105

взаимодействия между включениями композита через замену неоднородной среды, окружающей произвольное включение, например, с прилегающей к нему прослойкой матрицы [25] или с фрагментом близлежащих соседних включений, как в самосогласованном варианте метода локального приближения [85], однородной анизотропной средой с искомыми эффективными свойствами композита. В исследованиях [33; 140] предложен обобщенный метод самосогласования, согласно которому расчет эффективных свойств композита локальным осреднением сводится к рассмотрению расчетной схемы самосогласования – это одиночное включение, окруженное некоторым неоднородным переходным слоем, в эффективной среде. На большом удалении от включения в эффективной среде заданы однородные поля, в частности макродеформаций. Отметим, что переходный слой появляется лишь в локальноосредненной задаче, и свойства этого слоя рассчитываются через корреляционные функции случайной структуры композита, и размер слоя соизмерим с радиусом корреляции структуры. В частных случаях для полидисперсных структур (см. рис. 1.1, б–е) композитов расчетные схемы самосогласования: одиночное включение в эффективной среде (см. рис. 1.1, г, е) и одиночное включение с прослойкой матрицы в эффективной среде (см. рис. 1.1, б, д).

Волокно с прослойкой матрицы в эффективной среде.

Рассмотрим обобщенную модель полидисперсной структуры (см. рис. 2.2, а) из полидисперсных частиц двух типов: с волокном (1-й фазой) и без волокна; a , b – внутренний и наружный радиусы произвольной частицы 1-го типа, варьируемый параметр

v

(a / b)2 = v / p

,

v v

1,

(2.86)

0

1

0

 

1

0

 

 

где p0 (v1;1) – относительное число частиц 1-го типа в полидисперсной структуре. В предельном случае при p0 1 (см. рис. 1.1, б):

v0 v1 , a b

v1

(2.87)

и в случае p0 v1 (см. рис. 1.1, г):

 

v0 1, a b

(2.88)

для всех значений v1 (0;1) .

 

106

y

x

а

б

Рис. 2.2. Фрагмент обобщенной полидисперсной структуры (а)

и расчетная схема (б): 1 – волокно, 2 – прослойка матрицы, 3 – эффективная среда

Ось z цилиндрической системы координат r , , z совместим с осью волокна (см. рис. 2.2, б); оси электрической и магнитной поляризации обеих трансверсально-изотропных пьезоактивных фаз: волокон (1-я фаза) и матрицы (2-я фаза) совпадают с координатной осью z и с направлением ориентации осей волокон. Определяющие соотношения (2.57) для обеих фаз связывают компонен-

ты тензора напряжений σ , векторов индукций электрического D и

магнитного B полей с компонентами тензора деформаций ɛ, век-

торов напряженностей электрического E и магнитного H полей

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

f =

 

компо-

1,2

ненты тензоров упругих C f , пьезоэлектрических

e f и пьезомаг-

нитных h f свойств, диэлектрических λ f

и магнитных μ f прони-

цаемостей; углы сдвига z = 2 z , rz = 2 rz

, r = 2 r .

Рассмотрим случай осесимметричного нагружения компо-

зита, когда на макроуровне заданы: Ez* 0 , Hz* 0 ,

* 0 , *zz 0 ,

0 , остальные компоненты векторов напряженностей электрического E* и магнитного H* полей и тензора деформаций ε* равны нулю, * = *rr + * – относительное изменение объема при

плоской деформации. Для такого случая нагружения определяющие соотношения (2.70) на макроуровне композита и вычисление

эффективных констант: k*

,

C*

, C*

, e*

, e*

, h*

, h*

,

*

,

12

 

3333

1133

333

311

333

311

 

33

 

 

 

 

 

 

 

 

 

 

107

*33 , *33 , *33 , 11* , *33 , *3 , *3 в (2.70) на основе принципа самосогласования сводятся к решению осесимметричной задачи (см. рис. 2.2, б)

с учетом в формуле (2.57) дополнительных

равенств Ez = Ez* ,

Hz = Hz* , zz = *zz .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для такой задачи уравнения равновесия

 

 

 

 

d rr

 

+

rr

= 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dr

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

соотношения Коши

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rr =

dur

, =

ur

 

 

 

 

(2.89)

 

dr

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и отличными от нуля перемещениями в плоскости xy

будут лишь

радиальные перемещения ur

; электрический и магнитный потен-

циалы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −Ez*z ,

 

= −Hz*z .

 

 

 

Общие решения для радиальных перемещений в волокне при

r (0;a)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(1)r = A(1)r ,

 

 

 

 

 

(2.90)

в матрице при r (a;b)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

= A

 

r +

B(2)

 

 

 

 

 

(2.91)

 

 

 

 

 

 

 

 

 

 

(2)r

(2)

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и в эффективной среде при r (b; )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= A*r +

B*

.

 

 

 

 

 

(2.92)

 

ur

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

Константы интегрирования

A

, A

,

 

B

,

A* , B*

в формулах

 

 

 

 

 

 

 

 

 

(1)

(2)

 

 

(2)

 

 

 

(2.90)–(2.92) определим из условий непрерывности радиальных перемещений ur и напряжений rr на границах при r = a

108

u(1)r = u(2)r ,

(1)rr

= (2)rr ,

(2.93)

при r = b

 

 

 

u(2)r =

ur ,

(2)rr

= rr

(2.94)

и условия заданной однородной макродеформации эффективной среды при r

 

 

 

 

 

= = * / 2 ,

 

 

 

 

 

(2.95)

 

 

 

rr

 

θθ

 

 

 

 

 

 

 

 

 

 

из которого с учетом соотношений (2.89), (2.92) для

 

 

 

 

 

 

 

dur

 

 

*

B*

 

 

 

 

 

 

 

 

 

 

 

= A

 

 

 

 

 

 

 

 

 

 

dr

 

r2

 

 

 

 

 

 

rr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

из условия (2.95) следует решение для константы

 

 

 

 

 

 

 

 

A* = * / 2.

 

 

 

 

 

(2.96)

Напряжения в волокне

 

 

 

 

 

 

 

 

 

 

 

 

 

(1)rr = (1) = 2k(1)12 A(1) + p(1) ,

 

(2.97)

 

 

 

(1) zz = 2C(1)1133 A(1)

+ q(1) ,

 

 

 

 

 

 

в матрице

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 2k

 

 

A

2G(2)12 B(2)

 

+ p

,

 

 

 

 

 

 

r2

 

 

(2)rr

(2)12

(2)

 

 

 

(2)

 

 

 

 

= 2k

 

 

A

+

 

2G(2)12 B(2)

+ p

,

(2.98)

 

 

 

 

r2

 

 

(2)θ

 

(2)12

(2)

 

 

 

(2)

 

 

(2)zz = 2C(2)1133 A(2) + q(2)

и в эффективной среде

 

 

= 2k*

A*

2G*

B*

,

 

12

+ p

 

rr

12

 

r2

(3)

 

109

= 2k*

A* +

2G*

B*

+ p

,

(2.99)

12

 

θθ

 

12

 

r2

 

 

(3)

 

 

zz

= 2C*

A* + q

 

,

 

 

 

 

1133

 

(3)

 

 

 

где

p(1) = C(1)1133 *zz e(1)311Ez* h(1)311Hz* (1)11 ,

p(2) = C(2)1133 *zz e(2)311Ez* h(2)311Hz* (2)11 ,

p

= C*

*

e*

E* h*

H * * ;

(2.100)

(3)

1133

zz

311

z 311

z 11

 

q(1) = C(1)3333 *zz e(1)333Ez* h(1)333Hz* (1)33 ,

q(2) = C(2)3333 *zz e(2)333Ez* h(2)333Hz* (2)33 ,

q(3) = C3333* *zz e333* Ez* h333* Hz* *33 ,

объемный модуль плоской деформации k12 = (C1111 + C1122 ) / 2 и мо-

дуль сдвига G12 = (C1111 C1122 ) / 2 в плоскости изотропии xy волокна, матрицы и эффективной среды.

В результате, записав условия (2.93), (2.94) в виде системы четырех алгебраических уравнений относительно искомых кон-

стант A(1) , A(2) , B(2) , B* с учетом (2.96), получим решение для константы

A

=

1 (a *

+ a *

+ a E* + a H *

+ a )

,

(2.101)

(1)

 

2

11

 

 

12

zz

 

 

13

 

z

 

14

 

 

z

 

15

 

 

 

 

где коэффициенты

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(k*

+ G*

),

 

 

 

 

 

*

 

 

 

 

 

 

 

 

a =

 

a =

C

 

 

+ C ,

 

 

 

11

3

 

 

12

 

12

 

 

 

12

 

3

 

1133

 

 

4

1133

 

 

 

= −

e *

 

 

,

 

 

 

= −

 

*

 

 

 

(2.102)

a

e

 

 

a

 

h

h

 

13

 

3

 

311

 

4

311

 

14

 

 

3

 

311

 

4

311

 

 

 

 

 

 

 

 

= −

 

*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

15

 

 

 

3

11

 

4

11

 

 

 

 

 

 

 

 

 

 

110

Соседние файлы в папке книги