книги / Пьезокомпозиты и датчики в 3 ч. Ч. 1 Статистическая механика пьезокомпозитов
.pdfличин uz , , |
, |
rz , |
Dr , Br |
|
при r = a и условий для осредненных |
|||||||||||||||||
величин |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
E = −v C − (1− v )C = E* 0 , |
|
|
|
||||||||||||||
|
|
|
|
|
x |
|
|
1 |
|
(1) |
|
|
1 |
(2) |
|
|
x |
|
|
|
|
|
|
|
|
|
|
H |
x |
= −v E |
− (1− v )E |
|
= H * 0 , |
|
|
(2.76) |
|||||||||
|
|
|
|
|
|
|
|
1 (1) |
|
|
1 (2) |
|
x |
|
|
|
|
|||||
|
|
|
|
|
|
xz |
= v1 A(1) |
+ (1− v1 )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 (1− v1 ) + e(2)113 (1+ v1 ) |
|
, |
|||||||
|
|
||||||||||
|
|
|
|
(1)11 |
(1− v ) + |
(2)11 |
(1+ v ) |
|
|||
|
|
|
|
1 |
1 |
|
|
||||
1 |
= |
h(1)113 (1− v1 ) + h(2)113 (1+ v1 ) |
, |
||||||||
|
|||||||||||
|
|
|
|
(1)11 |
(1− v ) + |
(2)11 |
(1+ v ) |
|
|||
|
|
|
|
|
1 |
|
1 |
|
|
b32 |
|
= 1b12 ; |
|
|
|
|
|
|
||||
2 = |
|
|
|
|
|
2 (2)11 |
|
|
|
, |
||
|
|
(1)11 |
(1− v ) + |
(2)11 |
(1+ v ) |
|||||||
|
|
|
|
1 |
|
1 |
|
|
||||
2 = |
|
|
|
|
|
2 (2)11 |
|
|
|
, |
||
|
(1)11 |
(1− v ) + |
(2)11 |
(1+ v ) |
||||||||
|
|
|
|
1 |
|
1 |
|
|
1+ v ,
0 = 3 + 1− v1 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 (1− v1 )e311 |
|
311 |
|
|
||
|
|
*33 = |
h |
(2.85) |
|||||
|
|
|
|
|
|
|
|
||
|
|
k(1)12 − v1k12 + G(2)12 |
|||||||
|
|
|
|
для идеальной периодической волокнистой структуры. Результаты расчета на рис. 2.1, а, б, позволяют сделать вывод, что инверсия свойств волокон и матрицы композита приводит к значительному увеличению абсолютных значений эффективных коэффициентов
электромагнитной связанности 11* , *33 композита при фиксированных величинах объемного содержания пьезоэлектрической v1 и пьезомагнитной 1− v1 фаз. В частности, при объемной доле 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