Научный журнал
Успехи современного естествознания
ISSN 1681-7494
"Перечень" ВАК
ИФ РИНЦ = 0,653

РАСЧЕТ ПЬЕЗОАКТЮАТОРОВ МЕТОДОМ КОНЕЧНЫХ ЭЛЕМЕНТОВ

Ершов В.С. 1 Ивашов Е.Н. 1 Федотов К.Д. 1
1 ФГАОУ ВПО «Московский институт электроники и математики Национального исследовательского университета «Высшая школа экономики»
Даны базовые понятия метода конечных элементов. Выполнен машинный эксперимент с пьезокерамическим актюатором. Проведено сравнение результатов машинного эксперимента с результатами теоретических исследований в данной области.
Метод конечных элементов
механика твердого деформируемого тела
пьезокерамика
пьезоактюатор
пьезоэффект.
Список литературы
1. Розин Л.А. Метод конечных элементов. – СПб.: СПбГТУ. – 2000. – С. 126-127.
2. Стренг. Г., Фикс Дж. Теория метода конечных элементов. – М.: Мир. 1977. – 349 с.
3. Зинкевич О., Морган К. Конечные элементы и аппроксимация. – М.: Мир. 1986. – 318 с.
4. Федотов К.Д. Взаимосвязь пьезоэлектрических и механических свойств материала. – Материалы НТК студентов, аспирантов и молодых специалистов МИЭМ НИУ ВШЭ. – М.: РИО МИЭМ. – 2013.
5. Феодосьев В.И. Сопротивление материалов. – М.: Наука. 1967. – 552 с.

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

В основе МКЭ лежит идея замены задачи отыскания самой функции на задачу отыскания конечного числа ее приближенных значений в произвольно выбранных узлах. В одномерном случае для решения подобной задачи выбирается разбиение отрезка на некоторое количество узлов, между которыми задаются отрезки кусочно-полиномиальных функций, которые позволяют провести дальнейшую аппроксимацию к искомой функции. При достаточно большом количестве данных отрезков можно говорить о том, что выполнена сколь угодно точная аппроксимация к искомой функции

Приведем несколько общеизвестных математических выражений из МКЭ [1].

Для одного элемента задается функция q(x), причем – u'' = q(x), u(0) = u(1) = 0, функция u неизвестна. u(x) можно представить в виде линейной комбинации

02.pdf

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

01.pdf. φj(х) – функции, совпадающие с пробными функциями, участвующими в записи вышеприведенной линейной комбинации. Дальнейшие преобразования приводят к матричным методам счисления и позволяют представить функции-элементы в удобном для вычислительной техники виде.

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

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

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

Конструкция имеет конфигурацию пьезотрубки, толщина подводящих электрическое напряжение электродов условно принимается равной нулю (при создании пьезотрубок толщина подводящих электродов подбирается таким образом, чтобы не оказывала влияние на деформации самой пьезотрубки при подаче напряжения). Высота элемента пьезотрубки 2 мм, всего восемь элементов. Внешний диаметр пьезотрубки 10 мм, внутренний – 4 мм. Материал пьезотрубки – пьезокерамика ЦТС-38.

а) 06.pdfб)06.pdf

Рис. 1. Вид пьезоактюатора сбоку (а) и сверху (б)

Порядок выполнения машинного эксперимента:

1) Создание элемента-кольца с высотой 2 мм, диаметром 10 мм и внутренним диаметром 4 мм.

2) Создание сборки из восьми жестко скрепленных торцами пьезоколец;

3) Задание материала и справочной геометрии для конструкции;

4) Создание исследования на статические нагрузки, создание сетки элементов и закрепление нижнего торца конструкции;

06.pdf

Рис. 2. Сетка конечных элементов

на исследуемом объекте

06.pdf

Рис. 3. Схема приложения деформирующей

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

5) Приложение деформирующей распределенной силы, равной 100 Н к торцу каждого пьезоэлемента;

6) Дублирование исследования и реверс направления деформирующих сил;

а)07.pdf б)07.pdfв)07.pdf

Рис. 4. Результаты работы программы

для случая сжатия пьезоактюатора: деформации (а), напряжения по Мизесу (Па) (б), перемещения (мм) (в)

Как видно из результатов машинного эксперимента, наибольшие деформации и напряжения (5,658×10-5 относительных единиц и 7,634 МПа) возникают в основании (области закрепления образца). Это подтверждается эпюрами аналитических решений задач по прикладной механике.

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

Для того, чтобы проверить точность машинного эксперимента, можно использовать формулу для нахождения абсолютного удлинения одного элемента пьезоактюатора [4].

03.pdf

где ε – диэлектрическая проницаемость пьезокерамики, ε0 – электрическая постоянная, E3 – напряженность электрического поля вдоль направления деформации, l – длина элемента, d33 – пьезомодуль материала вдоль направления деформации, Yz – модуль Юнга вдоль направления деформации.

Напряженность электрического поля можно вычислить по формуле

04.pdf

где Р3 – растягивающая/сжимающая сила (в случае эксперимента 100 Н), F – площадь торцевой грани элемента

F = 66×10-6 м2. Yz = 6.8×1010 Н/м2.

Можно преобразовать формулу абсолютного удлинения.

05.pdf

разница между машинным экспериментом и теорией составляет 356-337 = 19 нм, это абсолютная погрешность. (19/356)×100% = 5,4%. Это свидетельствует об адекватности теоретической модели и машинного эксперимента.

В заключение следует отметить, что метод конечных элементов сложнее известного метода конечных разностей. Однако у МКЭ есть ряд преимуществ, проявляющихся на реальных задачах: произвольная форма обрабатываемой области, сетку можно сделать более редкой в тех местах, где особая точность не нужна.


Библиографическая ссылка

Ершов В.С., Ивашов Е.Н., Федотов К.Д. РАСЧЕТ ПЬЕЗОАКТЮАТОРОВ МЕТОДОМ КОНЕЧНЫХ ЭЛЕМЕНТОВ // Успехи современного естествознания. – 2014. – № 1. – С. 45-47;
URL: http://natural-sciences.ru/ru/article/view?id=33201 (дата обращения: 19.07.2019).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.252