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

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

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

ячейками: коррелированными полидисперсными (см. рис. 1.18, б) (случайное взаимное расположение центров сечений волокон которых совпадает с расположением центров сечений волокон в начальной полидисперсной структуре III, (см. рис. 1.18, а) или монодисперсными (с квазипериодическим расположением в узлах гексагональной решетки, см. рис. 1.18, г) составными ячейками заполнено некоррелированными по размерам и взаимному расположению такими же, но, например, бесконечно малыми составными ячейками. Величины, относящиеся к коррелированным полидисперсным структурам (см. рис. 1.18, б, г), обозначим верхним индексом « I' ».

Радиусы «привязанных» составных ячеек в коррелированных полидисперсных структурах (см. рис. 1.18, б, г) попарно совпада-

ют, а радиусы расположенных в них волокон rI'

rIII . Коэффици-

ент корреляции коррелированных

полидисперсных структур III

(см. рис. 1.18, а) и I(см. рис. 1.18, б)

 

 

 

 

p =

p0 v1

,

 

 

(1.105)

1v

 

 

 

 

 

 

 

1

 

 

 

 

так как для рассматриваемых структур v

= p v с учетом rI'

rIII .

 

 

11

0 1

 

 

а

б

в г

Рис. 1.18. Фрагменты реализаций коррелированных полидисперсных структур III (а) и I(б), квазипериодических монодисперсной (в)

и коррелированной полидисперсной I(г) структур

61

Квазипериодическая полидисперсная структура (рис. 1.18, г) образована внесением корреляций или «привязок» к заданной монодисперсной квазипериодической структуре. Считаем, что в поперечной плоскости изотропии над каждым круговым сечением волокна радиуса r монодисперсной квазипериодической структу-

ры мысленно расположено концентрическое круговое сечение волокна радиуса r I' некоторой составной ячейки с внешним радиу-

сом RI' полидисперсной структуры; центр каждой гексагональной

ячейки с волокном квазипериодической монодисперсной структуры совпадает с центром соответствующей составной ячейки полидисперсной структуры. Так как в заданной квазипериодической монодисперсной структуре радиусы всех волокон одинаковые и равны r , поэтому в «привязанных» к ним волокнах полидисперс-

ной структуры радиусы тоже считаем одинаковыми и равными r I' и в общем – r rI' . В полидисперсной структуре пространство

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

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

ты корреляции различных пар: квазипериодических монодисперсной (см. рис. 1.18, а) и полидисперсной (см. рис. 1.18, б), монодисперсных квазипериодической и периодической, периодических монодисперсной и полидисперсной случайных структур.

Рассмотрим расчет коэффициента корреляции p1 квазипери-

одических монодисперсной и полидисперсной структур. Коэффициент корреляции

iI'/i/

p1 = 1I' 1 (1.106)

1 1

62

случайных статистически однородных индикаторных функций 1-й фазы (волокон) двухфазных квазипериодических монодисперсной

i (r) и полидисперсной

iI' (r) структур

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

i (r) =

1,

r V

iI' (r) =

1,

r VI'

(1.107)

0,

 

1

0,

1

1

 

r V

1

 

r VI'

 

 

 

 

 

 

1

 

 

 

1

 

рассчитывается

 

 

 

 

 

 

 

 

 

 

 

p1

=

 

v11 v12

,

v11

= i1I'i1

 

(1.108)

v (1v )

 

 

 

 

 

1

1

 

 

 

 

 

 

через одноточечный начальный смешанный момент v11 , где пуль-

сации i1I'/ (r) = i1I' (r) i1I' , i1/ (r) = i1 (r) i1

, осредненные значения

iI'

 

= v ,

среднеквадратические

отклонения

I' = =

 

,

= i

D

1

1

1

 

 

 

1

1

11

 

дисперсия D

(1.15), области 1-й V ,

VI'

и 2-й V

, VI' фаз обеих

 

 

11

1

1

2

 

2

 

 

структур в представительной области

V = V1 V2 .

Для тривиаль-

ного случая, когда полидисперсная некоррелированная структура произвольным образом «накладывается» на монодисперсную ква-

зипериодическую структуру имеем v11 = v12 , p1 = 0 .

Из формулы для объемной доли волокон в монодисперсной квазипериодической структуре

v = p (r

/ r

)2 v

(1.109)

1

 

max

1max

 

следует, что радиус вписанной в гексагональную ячейку окружности

rmax = r

p v1max / v1

.

(1.110)

Объемная доля волокон в монодисперсной квазипериодической структуре (1.109) равна объемной доле волокон в полидис-

персной структуре и в каждой ее составной ячейке v1 = (rI' / RI' )2 , где внутренний r I' и внешний RI' радиусы однотипной составной ячейки полидисперсной структуры. Минимальная гарантированная

63

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

= 2r (

p v1max / v1

1).

(1.111)

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

RI' rmax :

0 x

p v1max

 

1

(1.112)

с учетом формулы (1.110) и неравенств:

p 1, v1max

1, где мак-

симально возможная объемная доля волокон при гексагональной

укладке

 

 

v1max =

3

/ 6 0,906,

радиус составной

ячейки

RI' = rI'

/

 

 

, отношение радиусов связанных волокон

x = rI' / r .

v

 

 

1

 

 

 

 

 

 

Из формулы (1.112) следует, что всегда rI' r при требовании

RI' rmax .

В частных случаях (1.112):

1) если RI' r

или 0

x

 

, имеем

 

 

v

 

 

 

 

 

 

 

 

1

 

 

 

 

 

v

= p

 

rI'2

+

(r2 RI'2 ) v

 

= v2

,

p = 0 ;

11

 

 

S

 

S

1

 

1

 

1

 

 

 

 

 

 

 

 

 

 

2) если RI'

 

или

 

 

 

 

 

 

 

 

,

имеем

 

r

 

v

 

x v

 

/ p

 

 

 

 

 

1

 

 

 

1max

 

 

 

 

 

 

 

 

rI'2

+ 0

 

= v x2 ,

 

 

x2

v

,

(1.113)

v

= p

 

 

 

p =

 

 

1

 

 

1v

11

 

S

 

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

где площадь гексагональной ячейки S. Решение (1.113) для одното-

чечного момента v11

и коэффициента корреляции

p1

можно пред-

ставить в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

64

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

= p v v

, p =

v1

,

(1.114)

1v

11

1

1max

1

 

 

 

 

 

 

1

 

 

где параметр = ' (1+ ' / 4)

рассчитывается через заданное значе-

ние фиксированной относительной толщины ' = / r минималь-

ной гарантированной прослойки между волокнами в монодисперсной квазипериодической структуре; в частном случае, когда = 0,

имеем p1 = 0 (1.114). При постоянстве прослойки при варьировании объемной доли волокон v1 (0;v1max ) вероятность наличия монодисперсных волокон квазипериодической структуры в узлах

гексагональной решетки изменяется p0 = v1 / v1max ,

максимальное

значение объемной доли волокон

 

 

 

 

 

 

 

v1max =

3

.

(1.115)

6(1+ ' / 2)2

Из формулы (1.113) видно, что максимальное значение коэффициент корреляции p1 принимает при наибольшем допусти-

мом значении параметра

x =

p v1max

.

(1.116)

Например, при v1 0,65

и p = v1max 0,906 из

формулы

(1.116) следует, что относительный радиус привязанного волокна полидисперсной структуры x 0,906 и коэффициент корреляции

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

p1 0,49. (1.117)

Отметим, что минимальная гарантированная прослойка между волокнами для квазипериодической монодисперсной структуры (1.111) в рассматриваемом случае 0,248r .

Для сравнения рассчитаем коэффициент корреляции p2 связанных монодисперсных квазипериодической ( p (0;1)) и привя-

65

занной периодической ( p =1) структур; центры и границы ячеек

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

r = rmax

v1 / ( p v1max )

(1.118)

и для идеальной периодической структуры

 

 

(1.119)

r p = r p

с учетом формул (1.109), (1.110) и равенства объемных долей волокон v1p = v1 в обеих структурах

 

 

v = (r p )2 /

S = p r2 / S.

 

 

(1.120)

 

 

1

 

 

 

 

 

 

 

Далее, так как в выражении (1.119) всегда

r p

r ,

поэтому

v

= p (r p )2

/ S, и в результате

 

 

 

 

 

 

11

 

 

 

 

 

 

 

 

 

 

 

v

= p v ,

p =

p v1

,

 

 

(1.121)

 

 

 

 

 

 

 

11

1

2

1v1

 

 

 

 

 

 

 

 

 

 

 

и для рассматриваемого случая, когда v1 0,65 и

p = 1max

0,906,

имеем

 

 

 

 

 

 

 

 

 

 

 

p2 0,73.

 

 

 

 

(1.122)

Аналогично формулам (1.118) – (1.122) рассчитаем коэффициент корреляции p3 связанных периодических монодисперсной и

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

 

rI''

= r

 

 

 

(1.123)

 

 

v

 

 

max

1

 

 

и монодисперсной

 

 

 

 

 

 

 

 

 

(1.124)

r p = r

 

v / v

 

max

1

1max

 

структур с учетом выражений (1.109), (1.110) и равенства объемных долей волокон v1p = v1 в обеих структурах

66

 

 

v

 

= (r p / r

)2 v

 

= (rI''

/ r

)2.

 

 

 

(1.125)

 

 

1

 

max

1max

 

 

 

max

 

 

 

 

 

Далее, так как из (1.123), (1.124) следует rI'' = r p

 

 

 

и не-

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1max

 

 

равенство

rI''

r p

,

поэтому

v

= (rI'' )2

/ S = x2 (r p )2

/ S, и в ре-

 

 

 

 

 

 

 

11

 

 

 

 

 

 

 

 

 

 

 

 

зультате

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= x2v

,

 

 

x2

v

,

 

 

 

(1.126)

 

 

 

 

v

p =

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11

1

 

3

 

1v1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где x rI'' / r p =

 

,

и для

рассматриваемого

случая,

 

когда

v

 

 

 

 

 

1max

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v1 0,65 и v1max 0,906, имеем коэффициент корреляции

 

 

 

 

 

 

 

 

p3 0,73.

 

 

 

 

 

 

 

(1.127)

Таким

образом,

для

рассматриваемого

случая,

 

когда

v1 0,65,

p = v1max

0,906, коэффициенты корреляции связанных

структур имеют значения: 0,49 для квазипериодических монодисперсной и полидисперсной (1.117), 0,73 для монодисперсных квазипериодической и периодической (1.122) и также для периодических монодисперсной и полидисперсной (1.127). Отметим, что численное значение коэффициента корреляции p 0,5 (1.102) хо-

рошо согласуется со значением 0,49, найденным из условия наилучшей полидисперсной аппроксимации корреляционной функции реальной структуры (1.117), и со значением для «приведенного коэффициента корреляции» p= p2 p3 0,53 [55].

Двухточечная корреляционная функция квазипериодической полидисперсной структуры. Рассмотрим вычисление

нормированной двухточечной корреляционной функции k11I' (ρ) (1.67) индикаторной функции 1-й фазы i1I' (r) (1.107) для двухфаз-

ной квазипериодической полидисперсной структуры (рис. 1.18, б) в представительной области V, где вектор ρ (1.52), дисперсия D11

(1.15), начальное значение k11I' (0) =1. Считаем, что в формуле (1.67) точки: r V , r1 V , где V – область рассматриваемой составной цилиндрической ячейки с радиусом кругового поперечного сечения

67

волокна r в представительной области V полидисперсной структуры. Нормированная корреляционная функция k11I' (ρ) (1.67) рас-

считывается (1.73) через условную вероятность PI' (r,r ) – вероят-

11 1

ность появления в точке r1 V 1-й фазы полидисперсной структуры при условии, что точка r V 1 находится 1-й фазе V 1 V рассматриваемой составной цилиндрической ячейки V . Так как вероятность попадания точки r V в 1-ю фазу составной ячейки V равна объемной доле v1 1-й фазы в этой ячейке, поэтому двух-

точечный момент iI' (r)iI' (r )

= v PI'

(r,r ) выражается через услов-

 

1

1

1

1

11

1

ную вероятность PI'

(r,r )

(1.73).

 

 

11

 

1

 

 

 

 

k11

 

 

 

 

 

k I'

 

 

 

 

 

 

11

/Rp

/Rp

а

б

k11

 

/Rp

в

Рис. 1.19. Нормированные корреляционные функции традиционных полидисперсных (а), квазипериодических полидисперсных (б), квазипериодических и периодической (в) структур

68

Для монодисперсных периодической и квазипериодических структур точки: r V p , r1 V , где V p – область ячейки перио-

дичности или квазипериодичности этих структур. На рис. 1.19 представлены результаты расчета нормированных корреляционных

функций k11(ρ) различных структур при объемной доле волокон v1 = 0,65, расстояние ρ = ρ (1.52), (1.68) между точками r и r1 . На рис. 1.19, а, изображены графики функций k11(ρ) для традицион-

ных полидисперсных структур: (□) для структуры на рис. 1.18, а,

(Δ) для структуры на рис. 1.18, б, с отсутствием корреляций размеров и типов ячеек в каждой из этих структур; в сравнении с кривыми (□), (Δ) на рис. 1.19, а, приведен график нормированной корреляционной функции (●) для заданной квазипериодической монодисперсной структуры (рис. 1.18, в) при значении величины

минимальной относительной толщины прослойки ' 0,25 между волокнами (1.111), ' Rp . На рис. 1.19, б – графики k11I' (ρ) для

квазипериодических полидисперсных структур при различных значениях вероятности расположения монодисперсных волокон в уз-

лах гексагональной решетки: p0 = 0 (□); 0,2 (◊); 0,4 (Δ); 0,6 (○); 0,8 (); в предельном случае, когда p0 = 0, такая квазипериодическая

структура вырождается в традиционную некоррелированную полидисперсную структуру (см. рис. 1.18, а; (□) рис. 1.19, а). На рис. 1.19,

в – графики k11(ρ) для коррелированных квазипериодических мо-

нодисперсной

(●)

(рис. 1.18,

в),

полидисперсной

(○)

( k11I' (ρ) ) (рис. 1.18, г) структур и, дополнительно, монодисперсной периодической структуры с гексагональной укладкой волокон (◊). Отметим наличие периодичности у «хвостов» функций k11(ρ) для

квазипериодических (см. рис. 1.18, в, г) и периодической структур, обусловленной периодичностью в расположении узлов гексагональной решетки.

69

Глава 2. ПОСТАНОВКА И МЕТОДЫ РЕШЕНИЯ СТОХАСТИЧЕСКОЙ СВЯЗАННОЙ КРАЕВОЙ ЗАДАЧИ МАГНИТОЭЛЕКТРОУПРУГОСТИ

Прогнозирование эффективных физико-механических свойств структурно-неоднородных, в том числе пьезоактивных [9; 19; 34; 76; 102], сред является одной из центральных задач механики композитов. В частности, на основе решения задачи протекания для микронеоднородной области могут быть определены эффективные, например диэлектрические, проницаемости композита [2; 9; 25; 102; 104]. Без учета пьезоактивности фаз известны различные решения [9; 104] для эффективных диэлектрических проницаемостей композита, наиболее известные из которых – это решения Фойгта – Рейсса, Хашина – Штрикмана и Максвелла – Гарнета.

Решения Фойгта

λ*F = λ

(2.1)

и Рейсса

λ*R1 = λ1

образуют вариационную «вилку», внутри которой лежат действительные λ* значения диэлектрических проницаемостей композита

λ*R λ* λ*F .

(2.2)

Решения, или вилка Хашина – Штрикмана λ*(1) , λ*(2) [9, 104] лежат внутри границ решения Фойгта – Рейсса (2.2)

λ*R λ*(1) λ* λ*(2) λ*F .

Для однонаправленного двухфазного волокнистого композита тензоры λ*(1) , λ*(2) могут быть получены [104] из решения задачи

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

70

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