CC BY
12
2019 Математика и механика № 58
УДК 539.3; 539.5
DOI 10.17223/19988621/58/11
А.С. Скалиух
ФУНКЦИОНАЛЬНАЯ ЗАВИСИМОСТЬ ФИЗИЧЕСКИХ ХАРАКТЕРИСТИК ОТ НЕОБРАТИМЫХ ПАРАМЕТРОВ ПРИ ЭЛЕКТРОМЕХАНИЧЕСКОМ ВОЗДЕЙСТВИИ НА СЕГНЕТОЭЛЕКТРИЧЕСКИЕ КЕРАМИКИ1
Для необратимых процессов деформирования и поляризации устанавливается линейная функциональная зависимость упругих, пьезоэлектрических и диэлектрических характеристик деформированного и частично поляризованного поликристаллического сегнетоэлектрического тела от необратимых параметров, к которым относятся остаточная деформация и поляризация, когда направление вектора электрического поля не совпадает ни с одной из главных осей тензора напряжения.
Ключевые слова: определяющие соотношения, деформирование, поляризация, сегнетоэлектрики, сегнетоэластики, индуцированные и остаточные параметры, необратимый процесс.
Обозначения и операции тензорной алгебры и анализа
В данной статье приняты следующие обозначения и операции тензорной алгебры и анализа: Тп - линейное пространство тензоров п-го ранга;
{ет >т=1,(ет >т=1 - векторы основного и взаимного базисов; gmn = ет • еп - метрические коэффициенты; о, £ еТ2 - тензоры напряжений и деформации; Е, Б еТ1 - векторы электрического поля и электрической индукции; £0 еТ2 -тензор остаточной деформации; Р0 еТ1 - вектор остаточной поляризации; 8 еТ4, d еТ3, э еТ2 - тензор упругих податливостей, тензор пьезоэлектрических модулей и тензор диэлектрических проницаемостей соответственно. Простое умножение d • Е = dm'гkEsgьетеп; полное умножение тензоров одного ранга о: £ = стmn6ksgmkgns = <зтггетп , полное умножение тензоров разных рангов 8: о = 8тпЬетеп = Бтпкаьетеп ; если ранг тензора первого множителя меньше ранга второго множителя, то о :8 = а']8тпк^]^текец =зтпк<зтпекеж.
Операция транспонирования 8Т((12),(34)) = Бтпкцекепеп . Набла - оператор Гад
мильтона V = ет-. Производная тензорной функции по тензорному аргументу
дхт
д8 дБтпк* г
8,£ = етепеке*ее1 и Т.д.
д£ дбу
1 Работа поддержана РФФИ, грант 17-08-00860-а.
Введение
Стремительное внедрение сегнетоэлектрических материалов в технику за последние два десятка лет обусловлено их замечательными пьезоэлектрическими и диэлектрическими свойствами. Большую роль в этом играют поликристаллические материалы, или керамики, технология изготовления которых позволяет получать рабочие элементы сенсоров и актуаторов различных форм и размеров. Для придания рабочих свойств керамику подвергают предварительной поляризации в сильном электрическом поле. В результате этого изменяется структура материала, изменяется класс анизотропии и появляются пьезоэлектрические свойства. Кроме электрического поля на процесс поляризации влияют также интенсивные механические напряжения. В отличие от пластических материалов здесь при интенсивных нагрузках появляется не только остаточная деформация, но и остаточная поляризация. В дальнейшем поляризованная керамика используется в малых электрических и механических полях, когда ее структура уже не изменяется. Однако сами процессы поляризации электрическим полем и деполяризации механическими напряжениями являются необратимыми процессами и играют важную роль в первоначальном формировании структуры материала. Математическое моделирование необратимых процессов деформирования и поляризации (квазистатический процесс) осуществляется на основе общих законов механики сплошных сред [1]. Полная система уравнений включает в себя полевые уравнения (уравнения равновесия и статики диэлектриков), геометрические соотношения (связь тензора деформаций с вектором перемещения и электрического поля с электрическим потенциалом) и определяющие соотношения (дифференциальные или интегральные операторы гистерезисного типа). Вывод определяющих соотношений связан с учетом многих факторов и должен опираться на физику явления, учитывать внутреннюю структуру материала и достаточно точно согласовываться с экспериментальными данными. К настоящему времени накоплен определенный опыт моделирования и можно выделить несколько основных направлений исследования, принципиально отличающиеся друг от друга. Это феноменологические методы [2-8], типа методов пластичности, методы двухуровневой сплошной среды [9-14], опирающиеся на методы микромеханических переключений и некоторые другие [15]. Эксперименты показывают, что в процессе поляризации меняется класс анизотропии материала, появляются пьезоэлектрические свойства, соответственно упругие, диэлектрические и пьезоэлектрические свойства также меняются. В отмеченных работах для упрощения делаются предположения о том, что упругие и диэлектрические свойства изменяются столь незначительно, что их можно взять такими, какие они в неполяризованном состоянии. Основной акцент делается на пьезоэлектрические модули, причем учитывается только влияние остаточной поляризации посредством введения множителя, в виде отношения модуля вектора остаточной поляризации на величину поляризации насыщения. Делается это безо всякого доказательства, исключительно с помощью интуитивных соображений. Однако при таком подходе остается открытым вопрос о механическом деформировании поликристаллического сегнетоэлектрика интенсивными механическими напряжениями, когда появляется остаточная деформация, меняется класс анизотропии материала, а описать эту анизотропию невозможно. Исследований, посвященным изучению влияния остаточной деформации на упругие модули в этом случае, автором не обнаружено. В [16] рассмотрена модель Джила - Атерто-на, являющаяся разновидностью микромеханических моделей. В ней авторы для
учета обратимых параметров вводят один скалярный коэффициент, да и тот приходится подбирать при сравнении с экспериментальными результатами. В [17] авторы предложили конечно-элементный подход, в рамках которого проводят усреднение и получают эффективные модули поляризованной керамики. Но и здесь функциональная зависимость от остаточных параметров не раскрыта.
В настоящей работе с помощью термодинамического подхода установлена функциональная зависимость для тензоров упругих, пьезоэлектрических и диэлектрических модулей от вектора остаточной поляризации и тензора остаточной деформации, когда протекает необратимый процесс поляризации и деформирования в поликристаллическом сегнетоэлектрике при нагрузках, не приводящих к разрушению. Предметом исследования являются поликристаллические сегнето-электрические материалы. Из общей модели, описывающей меняющиеся упругие и необратимые параметры, целенаправленно выбраны определяющие соотношения для обратимых частей поляризации и деформирования. В отличие от многих задач пластических материалов, где упругие деформации подчиняются закону Гу-ка с модулем упругости и коэффициентом Пуассона изотропного тела, здесь тело претерпевает фазовый переход «твердое тело - твердое тело», переходя из изотропного состояния в анизотропное, при котором физические модули изменяются в процессе развития пластической деформаций и остаточной поляризации.
1. Определяющие соотношения для обратимых параметров
Для определения параметров рассматриваемых процессов обратимся к основным положениям внутреннего строения поликристаллических сегнетоэлектриков [1]. Сегнетоэлектрическая керамика состоит из огромного числа кристаллитов, размерами порядка 10-6 - 10-5 м. На линейном размере кристаллита оказывается в среднем от одного до сотни доменов, каждый из которых способен содержать сотни и более атомных ячеек. Поэтому как в мелкозернистых, так и в крупнозернистых керамиках количество атомных ячеек по линейному размеру кристаллита составляет число порядка 104. В каждой атомной ячейке в низкотемпературной сегнетоэлектрической фазе центр положительных и отрицательных зарядов не совпадают, что характеризуется объемной плотностью электрических диполей в виде вектора спонтанной поляризации р, и тензора спонтанной деформации е;, одна из главных осей которого совпадает с направлением вектора р;. Пусть для определенности рассматриваются сегнетоэлектрики типа перовскита, типичным представителем которых является титанат бария. Под представительным объемом понимаем объем, линейные размеры которого значительно меньше линейных размеров конечного тела, но содержащий огромное количество кристаллитов, для сохранения физических свойств тела. Для наглядности на рис. 1 показан представительный объем, в котором каждый кристаллит просто совпадает с одним доменом. В деполяризованном состоянии векторы спонтанной поляризации, показанные стрелками, располагаются произвольным образом. Взяв некоторую точку за точку приведения, можно поставить в соответствие каждому вектору спонтанной поляризации коллинеарный ему единичный вектор. Деполяризованное состояние характеризуется равномерным распределением таких векторов на единичной сфере, как показано на том же рисунке.
Если приложить интенсивное электрическое поле, то векторы спонтанной поляризации поворачиваются наиболее близко к направлению поля, подчиняясь кристаллографическому строению сегнетоэлектрика, что и показано на рис. 2.
Рис. 1. Представительный объем в деполяризованном состоянии и распределение осей доменов
Fig. 1. Representative volume in a depolarized state anddistribution of the domain axes
Рис. 2. Представительный объем в поляризованном электрическим полем состоянии и распределение осей доменов Fig. 2. Representative volume in a polarized state provided by electric field and distribution of the domain axes
Изображающие единичные векторы расположатся в конусе, показанном на рис. 3, а. После снятия электрического поля векторы спонтанной поляризации остаются в таком положении и в начальное состояние не возвращаются, что характеризует изменившуюся структуру. Однако и механические напряжения могут поворачивать векторы спонтанной поляризации, при достижении ими пороговых значений. В частности, для сжимающих напряжений оси доменов располагаются в торе треугольного сечения, а для растягивающих - в двух одноосных конусах с общей вершиной, как показано на рис. 3, Ь и с.
Рис. 3. Распределение осей векторов спонтанной поляризации после воздействия: а - электрического поля; b - сжимающих напряжений; c - растягивающих напряжений Fig. 3. Distribution of the axes of spontaneous polarization vectors after the action of: (a) electric field, (b) compressive stresses, and (c) tensile stresses
В этом проявляются отличительные особенности сегнетоэлектриков от ферромагнетиков, у которых векторы спонтанной намагниченности могут поворачиваться в направлении магнитного поля на любой угол.
В представительном объеме, содержащем N элементарных ячеек, определяется вектор остаточной поляризации и тензор остаточной деформации путем операции осреднения:
1 N 1 N
ро = N ^(р4 )к, £° = N ^)к.
Л к=1 Л к=1
Пусть к представительному объему поликристаллической сегнетоэлектриче-ской среды приложено внешнее электрическое поле Е и механическое напряже-
ние ст. Считается, что эти поля неизменны во всех точках представительного объема, что позволяет считать их определяющими параметрами. Под их воздействием доменные стенки деформируются. Эти явления описываются результирующими эффектами в виде обратимых составляющих вектора поляризации и тензора деформации Ре, ее, которые исчезают при исчезновении нагрузок. Как только механические или электрические воздействия достигают пороговых значений, начинаются повороты векторов спонтанной поляризации, и связанные с ними изменениями направлений главных осей тензоров спонтанной деформации. Доменные стенки передвигаются, что приводит к появлению вышеопределенных необратимых параметров Р0, е0. В связи с вышесказанным вектор поляризации Р и тензор деформации е в необратимых процессах подразделяют на обратимые Ре, ее и необратимые Р0, е0 части: Р = Ре + Р0, £ = £е + £0, и относят к искомым параметрам. Обратимые части являются параметрами состояния и связаны с внутренней энергией, а необратимые части являются параметрами процесса, и для их определения приходится использовать дополнительные соотношения, наступающие в процессе достижения нагрузками пороговых значений.
Целью нашего исследования является построение определяющих соотношений обратимых составляющих и установление функциональной зависимости физических характеристик материала от остаточных параметров. Поскольку речь идет о характеристиках представительного объема, воспользуемся элементами термодинамики необратимых процессов. Известно, что физические характеристики термически деполяризованной среды в отличие от среды, поляризованной до состояния насыщения, отличаются между собой. Поэтому в случае частичной поляризации они будут зависеть от остаточных параметров Р0, е0, и будут меняться с изменением этих параметров. Найдем уравнения, связывающие обратимые искомые параметры Ре, ее с определяющими параметрами Е и ст. С этой целью воспользуемся уравнениями первого и второго начала термодинамики необратимых процессов в дифференциальной форме [1]:
где р, и, о, £, Е, Б, q, ю,Т, V - плотность массы, массовая плотность внутренней энергии, тензор напряжений, тензор деформаций, вектор электрического поля, вектор электрической индукции, вектор потока тепла, мощность внутренних источников тепла, массовая плотность энтропии, абсолютная температура, набла оператор Гамильтона соответственно. Точка сверху обозначает субстанциональную производную по времени. Но так как все деформации малы, то справедливы гипотезы геометрически линейной механики, откуда следует, что субстанциональная производная совпадает с частной производной.
Внутренняя энергия является функцией внешних параметров и энтропии и (£е, Бе, я), Бе = Ре +е* Е, где е* - диэлектрическая проницаемость вакуума. Чтобы определяющие параметры стали независимыми, введем термодинамическую функцию Гиббса
и проведем стандартные действия, связанные с исключением из соотношений (1)
ри = о: £ + Е• I)- V-q + рю,
(1)
О = и - Тя -1 о: £ -1Е • I Р Р
мощности внутренних источников тепла. Получим неравенство диссипации, в котором для случая процесса холодной поляризации можно пренебречь термическими слагаемыми (q, я, Т), и которое можно записать в виде
-(рО,о +£): о - (рО,Е + Б) • Е - рО,£о: £0 - рО,Ро ^ = $ > 0. (2)
Динамика переключения диполей составляет менее 10-11 с. Процесс переключения домена имеет временной масштаб порядка 10-8 - 10-5 с. Поэтому, когда напряжение и электрическое поле изменяются в пределах 10-1 - 10-2 с, имеем квазистатический процесс и малые отклонения от термодинамического равновесия. Тогда функцию скорости диссипации $ можно рассматривать как линейную относительно скоростей своих параметров, т.е.
$ = Х$ : о + X Б • Е + х$0 : £ о + хР0 • Ро. (3)
Обобщенные диссипативные силы х$, х), х$0, хР0 в общем случае представляют собой функции, которые зависят от набора параметров о, Е, £0, Р0, равно как и от скоростей этих параметров о, Е, £0, Р0, причем они должны обеспечивать выполнение неравенства скорости диссипации (2).
Для керамических материалов обратимые части деформации и поляризации £е, Ре достаточно малы (приблизительно на порядок меньше остаточных параметров), в силу чего определяющие соотношения естественно строить в виде линейных соотношений. С этой целью конкретизируем функцию Гиббса, выбрав ее в квадратичном виде:
О = --(о : £0 + Е• Р0)-±[о : 8(^,Р0): о + 2Е• Р0): о + Е• э(^,Р0)• Е]. (4) Р 2Р
Входящие сюда упругие податливости, пьезоэлектрические модули и диэлектрические проницаемости являются функциями остаточных параметров 8 = 8(£0, Р0), d = d(£0, Р0), э = э(£0,Р0) и обладают следующими свойствами
симметрии 8 = 8т(12) = 8т(34) = 8т((12),(34)), а = аТ(23), э = эт(12). Воспользовавшись независимостью обобщенных скоростей о, Е, £0, Р0 и проведя стандартные математические операции в (2) на основании (4), получаем две пары соотношений. Первая пара связывает тензор деформаций и вектор электрической индукции с тензором механических напряжений, вектором электрического поля и первой парой диссипативных сил. Вторая пара уравнений связывает тензор механических напряжений, вектор электрического поля и производные модулей физических характеристик по остаточным параметрам со второй парой диссипативных сил. Избегая чрезмерного обобщения, будем пренебрегать вязкоупругими свойствами, для чего надо положить х$ = 0, XI = 0 . Тогда первую пару уравнений можно представить в виде
£ - £0 = 8(£0, Р0): о + аТ(£0, Р0) • Е, (5)
Б -Р0 = Р0): о + э(£0,Р0) • Е,
где для тензора пьезомодулей аТ = ¿Т(1,(23)). Вторую пару уравнений можно записать так:
а+1 в: г*. 1т((34),(5б)): в+Е г*.1Т((23),(45)): в+1Е г*.Г(34)) Е *
° + 2 ° : ) : ° + \3£°) : ° + 2 \д£° ) ' Хе°
1 ( ЯС \Т((34),5) , \Т((23),4) , \ Т(2,3)
Е +1 а : М : а + Е-И : а +1Е-М • Е = хР .
2 Г5Р°1 ГдР°1 2 Г5Р° J АР°
(б)
Учитывая, что Б - Р° = Бе, £ - £° = ге, можно сказать, что (5) представляют
собой определяющие соотношения для обратимых параметров. Остается только выяснить характер зависимости физических характеристик от остаточных параметров. С этой целью рассмотрим обобщенные диссипативные силы х^0, хРо,
входящие в систему (б). Для необратимых процессов мы не вправе считать, что они являются аналитическими функциями своих аргументов. Сделаем предположение, что они состоят как из аналитической, так и неаналитической частей, входящих аддитивным образом, с обязательным выполнением неравенства (2).
Обозначим аналитические части через х5(ап), XР^""), а неаналитические -
хЕ°( % Xр°( % тогда
х 5 = х 5 ("" ) + х 5 ("" ) х 5 = х 5 ("") + х 5 ("") >4 ' ЛР° Лр°
Очевидно, что аналитические части не могут зависеть от остаточных параметров, поэтому они будут функциями только внешних параметров:
X5°("") = X5°("") (а, Е), хр°("") = Xр°("") (а, Е).
Надо также отметить, что при отсутствии внешних нагрузок эти функции равны нулю. Все это позволяет разложить их в ряд Тейлора в окрестности ненагру-женного состояния. Выполнив такое разложение с точностью до квадратичных составляющих и введя для удобства записи дальнейших формул соответствующие (б) операции транспонирования, для тензоров такого разложения получим
X5("") = А: а + В - Е + а : Кт((34),(5б)) : а + Е - Мт(2,(34)) - Е + Е - Nт((23),(45)) : а, Xр(а") = А1 : а + В1 - Е + а : к/«34«56» : а + Е - М^^34» - Е + Е - ^т((23),(45)) : а.
Тензоры В) еТ2; В,А1,М1 еТ3; А,М,N1 еТ4; N,К) еТ5, КеТб не зависят от остаточных £°, Р° и от определяющих параметров а, Е . Введенные аналитические части диссипативных сил никак не изменяют неравенство диссипации (2), они лишь исключают из функции скорости диссипации слагаемые, которые не влияют на диссипацию энергии. Подставляя эти разложения в (б) и учитывая независимость определяющих параметров, получаем систему равенств, которые связывают эти параметры с неаналитическими частями диссипативных сил
X5(""), Xр^"") и с производными физических характеристик, т.е. имеем
а = А: а + В - Е + у^"а\ Е = А1 : а + В1 - Е + хр(""), (7)
58 М _т дэ 58 _т 5э -= К,-= N-= М, -= К1,-= -= М1. (8)
д£° д £° д£° дР° дР° дР°
Соотношения (7) описывают изменение необратимых параметров. Однако построение неаналитических частей обобщенных диссипативных сил представляет собой очень сложную задачу, которая, по-видимому, даже в простейших случаях не решена. Вместо этого для нахождения необратимых параметров применяются другие методы и подходы [18], с условием удовлетворения неравенству (2).
Однако для наших целей важны соотношения (8), которые позволяют сформулировать линейную зависимость тензоров 8(£0, Р0), Д(£0,Р0), э(£0, Р0) от остаточных параметров с точностью до постоянных тензоров 80, а0, э0:
8(£0,Р0) = 80 + к: £0 + К • Р0,
Р0) = ^ + N : £0 + К! • Р0, (9)
э(£0, Р0) = Э0 + М : £0 + М1 • Р0.
Очевидно, что в полученных соотношениях 80, э0 - тензоры упругих податли-востей и диэлектрических проницаемостей для термически деполяризованной керамики, когда £0 = 0, Р0 = 0. В деполяризованном состоянии керамика не обладает пьезоэлектрическими свойствами, поэтому а0 = 0. Простые рассуждения позволяют еще больше упростить эти соотношения, рассмотрев ситуации с изменением направления вектора остаточной поляризации. Действительно, если направление вектора остаточной поляризации изменить на противоположное, то тензор остаточной деформации не изменится. Соответственно тензор упругих податли-востей и тензор диэлектрических проницаемостей также не изменятся, в отличие от тензора пьезоэлектрических модулей, который поменяет знак на противоположный. Отсюда сразу же вытекают следующие равенства:
к1 = 0, М1 = 0, N = 0 . Следовательно, вместо (9) имеем:
8(£0) = 80 + К: £0,
од = N • Р0, (10)
э(£0) = Э0 + М : £0.
Тензоры упругих податливостей и диэлектрических проницаемостей для термически деполяризованной керамики 80,э0 должны быть определены до начала
процесса поляризации. Соотношения (10) составляют одну из главных частей исследования. Остается лишь конкретизировать компоненты входящих в них тензоров.
2. Обсуждение результатов
Тензоры К, N2, М методами термодинамики определить уже невозможно, здесь, как и при определении физических характеристик материала, необходимы экспериментальные методы. Заметим, что эти тензоры, будучи связанными с физическими модулями, должны обладать такими же свойствами симметрии, как и тензоры физических модулей 8, а, э . Кроме того, К, М должны быть симметричными по последней паре индексов, что вытекает из (8), в силу дифференцирования физических характеристик по симметричному тензору остаточной деформации. Это обстоятельство позволяет вместо тензорного представления воспользоваться матричным представлением Фойхта. Используя декартову прямоуголь-
ную систему координат, введем векторы
£0 = (е011, е022, 8033, 2е023, 2е013,
Тензорам К е Т6, ],М,8,80 е Т4, а е Т3. матрицы К, ], М, 8,8 0 ,(1,э,э0 по следующему правилу:
2е012 Р0 = {Р01,, Р02,, Р03 } •
э, э0 е Т2 поставим в соответствие
К
(тп)(кз)( р1)
•К
ару, М1 т(пк )з
•м
1 таз >
$ та
м
тп(кз)
>м
"(тп)(кз) ^ "ар,
(тп)(кз) 0аР' т(пк) та' тп тп' тп 0 тп'
где в скобках выделены пары латинских индексов, которые заменяются на один греческий, а а, р, у определяются по одному и тому же правилу, например:
(т, т = п;
9 - (т + п), т Ф п.
Вначале определяются компоненты матриц 80,э0 • С этой целью отметим [1], что, термически деполяризованная керамика является изотропным телом. Поэтому, проведя два упругих и один электрический эксперимент, можно полностью определить компоненты этих матриц. Действительно, среди упругих экспериментов наиболее простыми являются статические эксперименты по определению модуля Юнга Е и коэффициента Пуассона V, а в качестве электрического эксперимента - эксперимент по определению диэлектрической проницаемости пластинки 80. Тогда компоненты матриц 80,э0 определяются как
5011 = 5011 = 5011 =
"044 = "0>55 = "066 =
1
Е
2(1 + V)
Е
5012 = "013 = "023 =
Е
(11)
эт 1 = эт 1 = эт 1 = £п
Для определения остальных матриц воспользуемся также известными значениями упругих податливостей и диэлектрических проницаемостей поляризованной до насыщения керамики 8ка1, Й8а1, э^ . Известно [1], что поляризованная
электрическим полем до насыщения керамика относится к классу трансверсально-изотропного тела и ее матрицы констант в системе координат, когда ось 07 направлена по вектору остаточной поляризации, имеют следующий вид:
«г =
0 А 0 0 0 0
0 0 0 0 0 2(511 - 512)
("11 "12 "13 0 0
"11 "13 0 0
"13 "13 "33 0 0
0 0 0 "44 0
0 0 0 0 "44
( 0 0
0
0 ё15 0А
0
ч ^31 ^31 ^33 0
0 0 0 0
Э11 0
0 Э 0
11 0 э
0 А 0
(12)
Одна из главных осей тензора деформации совпадает с осью Ох, а две другие лежат в плоскости Oxy• Поэтому будем определять компоненты матриц тензоров
К, N1, М в этих осях. Обратимся вначале к матрице К. С этой целью рассмотрим матрицу К • £0. В состоянии, когда керамика поляризована до состояния насыщения, компоненты е01, е02, е03 являются главными значениями тензора остаточной деформации, удовлетворяющими условию е01 < е02 < е03. Как следует из экспе-
риментальных данных, тензор остаточных деформаций является несжимаемым, тогда из матричного равенства = £50 + К • £0 следует, что:
^ар = ^0ар + + Кар2802 + Кар3803, 801 = 802 = 803 = ^ ,
где 8м: - осевая деформация вдоль направления поляризации в состоянии насыщения. Естественно считать, что всем З^р = 0 будут отвечать КГару = 0, Уу. Так
как поворот осей вокруг оси 02 не влияет на все главные значения, а главные значения деформаций в плоскости Оху равны между собой, можно считать, что
^ар1 = Кар2 = 1^ ККарз = тар^ар1 = тар /СТ0 , где ст0 - некоторый множИтель, размерность которого совпадает с размерностью механического напряжения. Тогда ^р = ^0ар +8м: (-1 + тар)/ ст0. Отсюда определяем тар, после чего находим
л л ~ ~ ~ (^ар — *^0ар К
^ар (80) = *^0ар + (801 +802 + тар803)/ст0, тар = 1 + .
В этих соотношениях 801 + 802 + 803 = 0 , так как это первый инвариант и выполняются условия несжимаемости. В силу этого неопределенный множитель ст0 уходит, а выражение принимает вид
*^ар (80) = *^0ар + 803. (12)
88а1
Отметим, что сюда входит лишь наибольшее главное значение тензора остаточной деформации, а все компоненты матриц в правой части по (11), (12) полностью известны. Аналогичные рассуждения позволяют найти матрицу М и определить изменяющиеся диэлектрические проницаемости:
л , ч л (Этп Э0тп) ^
Этп (80) = Э0 тп +-803 . (13)
88а1
Остается найти матрицу ], после чего определить пезоэлектрические модули. С этой целью рассмотрим ту же систему координат и состояние полной поляризации, тогда компоненты матричного равенства (1 ш1 = ] • Р0:
¿та = ^~\та\Р0\ + ^1та2Р02 + ^1та3Р03, ^^01 = ^^01 = 0, Р03 = .Psat ,
где - поляризация в состоянии насыщения. Опуская схожие рассуждения, можем окончательно записать:
¿та (Р0) = — ¿та . (14)
Аа1
Соотношения (12) - (14) полностью определяют физические характеристики материала в необратимом процессе деформирования и поляризации.
Соотношения (12), (13) надо рассматривать в локальной системе координат, оси которой совпадают с главными осями тензора остаточной деформации, когда ось 02 направлена по третьему главному направлению и направлена по вектору остаточной поляризации.
3. Заключение
Необратимые процессы в поликристаллических сегнетоэлектрических средах связаны с внутренней перестройкой структуры материала, что приводит к изменению и механических и электрических свойств материала. Изменяются не только пьезоэлектрические модули, но и упругие податливости и диэлектрические проницаемости. Если моделируется необратимый процесс поляризации и деформирования, то необходимо знать функциональную зависимость этих модулей от изменяющихся остаточных деформаций, равно как и от остаточной поляризации. Очевидно, что физические материальные характеристики должны определяться для представительного объема, включающего в себя множество микроструктурных особенностей в виде кристаллитов и доменов. Для таких объемов можно использовать элементы термодинамики необратимых процессов, что и было сделано в настоящей работе. Основным моментом в таком подходе явилось аддитивное представление диссипативных сил при скоростях остаточных параметров в виде аналитической и неаналитической частей. Именно это позволило конкретизировать аналитические части и построить линейные тензорные соотношения относительно остаточных параметров для физических характеристик материала. С помощью простого анализа проведена конкретизация этих зависимостей: в частности, выяснилось, что упругие податливости и диэлектрические проницаемости зависят от деформаций, в то время как пьезомодули зависят от поляризации, что наглядно продемонстрировано соотношениями (10). Входящие в них тензоры с постоянными компонентами имеют высокие ранги, что оказывает существенное влияние на их определение. Очевидно, что методами термодинамики определить эти коэффициенты не представляется возможным. Поэтому были задействованы данные двух состояний керамики, когда она деполяризована и поляризована до состояния насыщения. И в том и в другом случае физические модули материала известны и их можно найти в справочной литературе для конкретного типа керамики. Реже встречаются данные для деполяризованного состояния. В этом случае предложено находить модули экспериментальным путем. При определении неизвестных величин были сделаны очевидные предположения о нулевых значениях некоторых компонентов тензоров из простых соображений анизотропии. Еще одна порция предположений сделана на основе имеющихся экспериментальных данных об осевой симметрии, несжимаемости и взаимовлияний полей. В итоге были получены соотношения (12) - (14), представляющие собой искомые функциональные зависимости. Входящие в них характеристики материала в деполяризованном состоянии и поляризованном до состояния насыщения материала определяются по (11), (12).
Подводя итоги, можно сказать, что, во-первых, построены определяющие соотношения (5), которые необходимы при моделировании необратимых процессов. Они представляют собой линейные тензорные соотношения, связывающие Бе, £е с определяющими параметрами Е и о. В них входят тензоры физических характеристик Б(е0), й(Р0), э(е0), меняющиеся в процессе эволюции остаточной деформации и поляризации. Интересно отметить, что в определяющих соотношениях линейная зависимость сохраняется при любых значениях и ориентации остаточных параметров е0, Р0. Во-вторых, исследованы функциональные зависимости физических модулей от остаточных параметров, которые представлены соотношениями (10). Этими соотношениями доказано, что физические характеристики частично поляризованного или деполяризованного тела, равно как и в случае не-
однородной поляризации, зависят линейным образом от остаточной деформации и остаточной поляризации. И наконец, выбраны компоненты постоянных тензоров по известным данным двух состояний материала. Интересно также отметить, что соотношения (14) полностью совпадают с данными работ [5, 9], а соотношения (12), (13) получены впервые.
ЛИТЕРАТУРА
1. Белоконь А.В., Скалиух А.С. Математическое моделирование необратимых процессов поляризации. М.: Физматлит, 2010.
2. Huber J.E., Fleck N.A. Multi-axial electrical switching of a ferroelectric: theory versus experiment // J. Mechanics and Physics of Solids. 2001. No. 49. P. 785-811. DOI: 10.1016/S0022-5096(00)00052-1.
3. Kamlah M., Bohle U. Finite element analysis of piezoceramic components taking into account ferroelectric hysteresis behavior // Int. J. Solids and Structures. 2001. No. 38. P. 605-633. DOI: 10.1016/S0020-7683(00)00055-X.
4. McMeeking R.M., Landis C.M. A phenomenological multi-axial constitutive law for switching in polycrystalline ferroelectric ceramics // Int. J. Engineering Science. 2002. No. 40. P. 1553-1577. DOI: 10.1016/S0020-7225(02)00033-2.
5. Landis C.M. Fully coupled, multi-axial, symmetric constitutive laws for polycrystalline ferroelectric ceramics // J. Mechanics and Physics of Solids. 2002. No. 50. P. 127-152. DOI: 10.1016/S0022-5096(01)00021-7.
6. Haug A., Knoblauch V., McMeeking R.M. Combined isotropic and kinematic hardening in phenomenological switching models for ferroelectric ceramics // Int. J. Engineering Science. 2003. No. 41. P. 867-901. DOI: 10.1016/S0020-7225(02)00320-8.
7. Elhadrouz M., Zineb T.B., Patoor E. Constitutive Law for Ferroelastic and Ferroelectric Pie-zoceramics // J. Intelligent Material Systems and Structures. 2005. V. 16. P. 221-236. DOI: 10.1177/1045389X05049650.
8. Muller R., Schroder J., Lupascu D.C. Thermodynamic consistent modelling of defects and microstructures in ferroelectrics // GAMM-Mitt. 31. 2008. No. 2. P. 133-150. DOI: 10.1002/gamm.200890007.
9. Huber J.E., Fleck N.A. Ferroelectric switching: a micromechanics model versus measured behaviour // European J. Mechanics A/Solids. 2004. No. 23. P. 203-217. DOI: 10.1016/ j.euromechsol.2003.11.006.
10. Kamlah M., Liskowsky A.C., McMeeking R.M., Balke H. Finite element simulation of a poly-crystalline ferroelectric based on a multidomain single crystal switching model // Int. J. Solids and Structures. 2005. No. 42. P. 2949-2964. DOI: 10.1016/j.ijsolstr. 2004.09.045.
11. Haug A., Onck P.R., Giessen E.V. Development of inter- and intragranular stresses during switching of ferroelectric polycrystals // Int. J. Solids and Structures. 2007. No. 44. P. 20662078. DOI: 10.1016/j.ijsolstr.2006.07.024.
12. Pane I., Fleck N.A., Chu D.P., Huber J.E. The influence of mechanical constraint upon the switching of a ferroelectric memory capacitor // European J. Mechanics A/Solids. 2009. No. 28. P. 195-201. DOI: 10.1016/j.euromechsol.2008.09.002.
13. Jayabal K., Menzel A., Arockiarajan A., Srinivasan S.M. Micromechanical modelling of switching phenomena in polycrystalline piezoceramics. Application of a polygonal finite element approach // Computational Mechanics. 2011. V. 48. No. 4. P. 421-435. DOI: 10.1007/s00466-011-0595-4.
14. Daniel L., HallD.A., Withers P.J. A multiscale model for reversible ferroelectric behavior of polycrystalline ceramics // Mechanics of Materials. 2014. No. 71. P. 85-100. DOI: 10.1016/j.mechmat.2014.01.006.
15. Aleshin V.I., Luchaninov A.G. Modeling of domain processes in piezoceramic materials // Ferroelectrics. 2002. V. 266. P. 111-124. https://doi.org/10.1080/00150190211447. DOI: 10.1080/00150190211447.
16. Smith R.C., Ounaies Z. A Domain Wall Model for Hysteresis in Piezoelectric Materials // J. Intelligent Material Systems and Structures. 2000. V. 11. No. 1. P. 62-79. DOI: 10.1106/HPHJ-UJ4D-E9D0-2MDY.
17. Осипова Н.Г., Семёнов А.С. Моделирование нелинейного поведения пьезокерамики тетрагональной структуры методами конечно-элементной гомогенизации // Научно-технические ведомости СПбГПУ, физико-математические науки. Санкт-Петербург. 2011. № 4(134). C. 56-64.
18. Skaliukh A. About Mathematical Models of Irreversible Polarization Processes of a Ferroelectric and Ferroelastic Polycrystals // Ferroelectrics and Their Applications / Husein Irzaman (ed). IntechOpen. DOI: 10.5772/intechopen.78262. URL: https://www.intechopen.com/books/ ferroelectrics-and-their-applications/about-mathematical-models-of-irreversible-polarization-processes-of-a-ferroelectric-and-ferroelastic
Статья поступила 10.10.18
Skaliukh A.S. (2019) FUNCTIONAL DEPENDENCE OF THE PHYSICAL CHARACTERISTICS ON IRREVERSIBLE PARAMETERS UNDER ELECTROMECHANICAL ACTION ON THE FERROELECTRIC CERAMICS. Vestnik Tomskogo gosudarstvennogo universiteta. Matematika i mekhanika [Tomsk State University Journal of Mathematics and Mechanics]. 58. pp. 128-141
DOI 10.17223/19988621/58/11
Keywords: constitutive equations, strain, polarization, ferroelectrics, ferroelastics, induced and residual parameters, irreversible process.
The main purpose of the work is to study the functional dependences of the material physical characteristics on the residual parameters of strain and polarization for irreversible process of polarization and strain. The dissipative forces in the inequality for dissipation are presented as a sum of analytical and non-analytical parts. Expansion of the analytical parts in powers of determining parameters provides the differential equations, which, being integrated, yield the unknown functional dependences. Integration constants represent the physical modules of a depolarized state. The tensor components of expansion in powers are determined by using elastic compliance modules, piezoelectric modules, and dielectric constants of two states: thermally depolarized ceramics and ceramics polarized up to saturation. As a result, the constitutive equations for reversible parameters were obtained as linear tensor equations with material elastic, piezoelectric, and dielectric properties representing tensor functions studied above.
Financial support.This study was supported by Russian Foundation for Basic Research (RFBR) under the research project № 17-08-00860-а.
SKALIUKH Alexander Sergeevich (Candidate of Physics and Mathematics, Institute of Mathematics, Mechanics and Computer Sciences named after I.I. Vorovich, Southern Federal University, Rostov-on-Don, Russian Federation)
REFERENCES
1. Belokon' A.V., Skaliukh A.S. (2010) Matematicheskoe modelirovanie neobratimykh protses-sovpolyarizatsii [Mathematical modeling of irreversible processes of polarization]. Moscow: Fizmatlit.
2. Huber J.E., Fleck N.A. (2001) Multi-axial electrical switching of a ferroelectric: theory versus experiment. Journal of the Mechanics and Physics of Solids. 49. pp. 785-811. DOI: 10.1016/S0022-5096(00)00052-1.
3. Kamlah M., Bohle U. (2001) Finite element analysis of piezoceramic components taking into account ferroelectric hysteresis behavior. International Journal of Solids and Structures. 38. pp. 605-633. DOI: 10.1016/S0020-7683(00)00055-X.
4. McMeeking R.M., Landis C.M. (2002) A phenomenological multi-axial constitutive law for switching in polycrystalline ferroelectric ceramics. International Journal of Engineering Science. 40. pp. 1553-1577. DOI: 10.1016/S0020-7225(02)00033-2.
5. Landis C.M. (2002) Fully coupled, multi-axial, symmetric constitutive laws for polycrystalline ferroelectric ceramics. Journal of the Mechanics and Physics of Solids. 50. pp. 127-152. DOI: 10.1016/S0022-5096(01)00021-7.
6. Haug A., Knoblauch V., McMeeking R.M. (2003) Combined isotropic and kinematic hardening in phenomenological switching models for ferroelectric ceramics. International Journal of Engineering Science. 41. pp. 887-901. DOI: 10.1016/S0020-7225(02)00320-8.
7. Elhadrouz M., Zineb T.B., Patoor E. (2005) Constitutive law for ferroelastic and ferroelectric piezoceramics. Journal of Intelligent Material Systems and Structures. 16. pp. 221-236. DOI: 10.1177/1045389X05049650.
8. Muller R., Schroder J., Lupascu D.C. (2008) Thermodynamic consistent modelling of defects and microstructures in ferroelectrics. GAMM-Mitteilungen. 31(2). pp. 133-150. DOI: 10.1002/gamm.200890007.
9. Huber J.E., Fleck N.A. (2004) Ferroelectric switching: a micromechanics model versus measured behavior. European Journal of Mechanics A/Solids. 23. pp. 203-217. DOI: 10.1016/j.euromechsol.2003.11.006.
10. Kamlah M., Liskowsky A.C., McMeeking R.M., Balke H. (2005) Finite element simulation of a polycrystalline ferroelectric based on a multidomain single crystal switching model. International Journal of Solids and Structures. 42. pp. 2949-2964. DOI: 10.1016/ j.ijsolstr.2004.09.045.
11. Haug A., Onck P.R., Giessen E. (2007) Development of inter- and intragranular stresses during switching of ferroelectric polycrystals. International Journal of Solids and Structures. 44. pp. 2066-2078. DOI: 10.1016/j.ijsolstr.2006.07.024.
12. Pane I., Fleck N.A., Chu D.P., Huber J.E. (2009) The influence of mechanical constraint upon the switching of a ferroelectric memory capacitor. European Journal of Mechanics A/Solids. 28. pp. 195-201. DOI: 10.1016/j.euromechsol.2008.09.002.
13. Jayabal K., Menzel A., Arockiarajan A., Srinivasan S.M. (2011) Micromechanical modelling of switching phenomena in polycrystalline piezoceramics. Application of a polygonal finite element approach. Computational Mechanics. 48(4). pp. 499-514. DOI: 10.1007/s00466-011-0595-4.
14. Daniel L., Hall D.A., Withers P.J. (2014) A multiscale model for reversible ferroelectric behaviour of polycrystalline ceramics. Mechanics of Materials. 71. pp. 85-100. DOI: 10.1016/ j.mechmat.2014.01.006.
15. Aleshin V.I., Luchaninov A.G. (2002) Modeling of domain processes in piezoceramic materials. Ferroelectrics. 266. pp. 111-124. DOI: 10.1080/00150190211447.
16. Smith R.C., Ounaies Z. (2000) A domain wall model for hysteresis in piezoelectric materials. Journal of Intelligent Material Systems and Structures. 11(1). pp. 62-79. DOI: 10.1106/ HPHJ-UJ4D-E9D0-2MDY.
17. Osipova N.G., Semionov A.S. (2011) Modelirovanie nelineynogo povedeniya p'ezokeramiki tetragonal'noy struktury metodami konechno-elementnoy gomogenizatsii [Simulation of the nonlinear behavior of ferroelectric ceramics with tetragonal structure by finite-element homogenizing]. Nauchno-tekhnicheskie vedomosti SPbGPU. Fiziko-matematicheskie nauki - St. Petersburg Polytechnic University Journal - Physics and Mathematics. 4(134). pp. 56-64.
18. Skaliukh A. (2018) About Mathematical Models of Irreversible Polarization Processes of a Ferroelectric and Ferroelastic Polycrystals. In «Ferroelectrics and Their Applications». IntechOpen. DOI: 10.5772/intechopen.78262. URL: https://www.intechopen.com/books/ ferroelectrics-and-their-applications/about-mathematical-models-of-irreversible-polarization-processes-of-a-ferroelectric-and-ferroelastic.
Received: October 10, 2018
