М ОДЕЛИРОВАНИЕ КОЛЕБАНИЙ РЕЗОНАТОРОВ НА ОСНОВЕ УГЛЕРОДНЫХ ВИСКЕРОВ Ванюшкина Валентина Научный руководитель: к. ф.-м. н., доц. каф. «Теоретическая механика» И.Е. Беринский 1 Санкт-Петербургский Политехнический Университет Петра Великого Кафедра Теоретической Механики 2015

У ГЛЕРОДНЫЕ НАНОВИСКЕРЫ Камертон Одиночный вискер 200 nm 2 Скальпель Вилка СПб НИУ ИТМО, Исследование параметров роста и механических свойств металл-углеродных многомерных вискерных наноструктур, Лукашенко С.Ю.Мухин И.С..Голубок А.О

П РИМЕНЕНИЕ 1. Скальпель 3 2. Камертон 3. Нановесы Z.L.Wang, P.Poncharal, W.A. de Heer. Measuring physical and mechanical properties of individual carbon nanotubes by in situ TEM. J. Phys. Chem. Solids, 2000, 61(7), pp

Ц ЕЛИ И ЗАДАЧИ Цель работы: исследование динамики углеродных вискеров с целью их оптимального применения в НЭМС(наноэлектромеханические системы). Задачи: 1. Моделирование вискеров, как систем с двумя степенями свободы 2. Исследование модели нано весов. 3. Конечно-элементное моделирование вискеров сложной формы. 4

И ССЛЕДОВАНИЕ ПРОДОЛЬНЫХ КОЛЕБАНИЙ СИСТЕМЫ, СОСТОЯЩЕЙ ИЗ ДВУХ МАСС И ДВУХ ПРУЖИН 5 Уравнения движения системы: Начальные условия:

М ЕТОД ИНТЕГРИРОВАНИЯ Эта и следующая задачи решены с помощью интегрирования методом конечных разностей (методом Эйлера). Аналитически получены выражения для вторых производных координат по времени(ускорений) и, с помощью интегрирования на языке программирования JavaScript, найдены значения текущих координат. 6

Р ЕЗУЛЬТАТЫ МОДЕЛИРОВАНИЯ 7 7 Одинаковые массы шаров и жесткости пружин: x t x v Зависимость координаты от времени Фазовая плоскость

8 Масса синего шара в 10 раз меньше массы красного: Масса синего шара в 10 раз больше массы красного: x tx v Зависимость координаты от времени Фазовая плоскость Зависимость координаты от времени Фазовая плоскость x x t v

9 Жесткость второй пружины в 10 раз больше жесткости первой: x tx v Зависимость координаты от времени Фазовая плоскость Жесткость первой пружины в 10 раз больше жесткости второй: Зависимость координаты от времени x t Фазовая плоскость v x

И ССЛЕДОВАНИЕ КОЛЕБАНИЙ СИСТЕМЫ, СОСТОЯЩЕЙ ИЗ ДВУХ СТЕРЖНЕЙ И ДВУХ УГЛОВЫХ ПРУЖИН 10 Начальные условия:

Для данной системы движение имеет более сложный характер, чем для предыдущей. 11 x x x x v v t t Р ЕЗУЛЬТАТЫ МОДЕЛИРОВАНИЯ Длина синего стержня в 5 раз меньше длины красного, а жесткость первой пружины в 10 раз меньше жесткости второй Длина синего стержня в 5 раз меньше длины красного, а масса первого стержня в 100 раз меньше массы второго

З АДАЧА О КОЛЕБАНИИ ВИСКЕРА С МАССОЙ НА КОНЦЕ 12 Параметры задачи: Цель данной задачи: исследовать зависимость собственных частот от массы на конце балки и найти соотношения, описывающие эту зависимость. Слева на вискер ставятся условия заделки, справа конец свободен и к нему прикреплена точечная масса

13 Таблица значений, полученных в результате моделирования Точечная масса, кг% от массы балки Первая собственная частота, МГц Вторая собственная частота, МГц 1001,38198, ,595211,35518, ,97651,26097, ,952101,16647, ,904201,02687, ,76500,792386, ,664700,7036, ,521000,611986, , ,567836, ,042000,455136,2261

Р ЕЗУЛЬТАТЫ МОДЕЛИРОВАНИЯ Слева балка заделана, а справа приложена точечная масса. 14 Тогда первая собственная форма будет выглядеть так: А вторая так:

15 f – частота колебаний. d – диаметр балки M – масса на конце балки l – длина балки m – масса балки 0.227*m + M – «эффективная» масса системы. «Dynamik der Baukonstruktionen», Christian Petersen st edition. 722 Получены формулы для связи первой собственной частоты и массы, приложенной к концу балки. При сравнение аналитических результатов с результатами моделирования наибольшая погрешность не превысила 3%

16 Расстояние от заделки до места приложения точечной массы, мкм Отношение расстояния от заделки до места приложения точеной массы к длине балки Первая собственная частота, MГц Вторая собственная частота, MГц 1001,38198, ,720,11,38188, ,440,21,38088, ,160,31,37688, ,880,41,36757, ,60,51,35097, ,320,61,32648, ,040,71,29438, ,760,81,25598, ,480,91,21288, ,211,16647,5956 З АДАЧА О КОЛЕБАНИИ ВИСКЕРА С МАССОЙ, ПРИКРЕПЛЕННОЙ НА РАЗНОМ РАССТОЯНИИ Постановка задачи: вискер заделан с левого края, правый конец свободен. Масса прикладывается к различным точкам.

Получены графики для частот, с помощью моделирования в ANSYS. Видим, что в данной задаче сложно вывести теоретическую зависимость. Но для второй собственной частоты видно, что характер графика повторяет по форме вторую собственную форму колебаний, 17

И ССЛЕДОВАНИЕ СОБСТВЕННЫХ И ВЫНУЖДЕННЫХ КОЛЕБАНИЙ КОНСТРУКЦИЙ ИЗ НАНОВИСКЕРОВ 18 Для исследования собственных колебаний модель вискера была заделана на торце(синяя стрелка) по всем направлениям движения. Для исследования вынужденных колебаний модели было придано гармоническое смещение по оси Х на торце(желтая стрелка)

19 Первая собственная форма Вторая собственная форма Напряжения при гармонических смещениях правого края, равных 0,1*sin(9*time)

nm Для исследования собственных колебаний модель камертона была заделана на торце(синяя стрелка) по всем направлениям движения, а для исследования вынужденных колебаний модели было придано гармоническое смещение по оси Z на торце(желтая стрелка)

21 Первая собственная форма Напряжения при гармонических смещениях нижнего края, равных 0,1*sin(9*time) Вторая собственная форма

22 Для исследования собственных колебаний модель вилки была заделана на торце(синяя стрелка) по всем направлениям движения, а для исследования вынужденных колебаний модели было придано гармоническое смещение по оси Z на торце(желтая стрелка)

23 Первая собственная форма Напряжения при гармонических смещениях нижнего края, равных 0,1*sin(9*time) Вторая собственная форма

24 Для исследования собственных колебаний модель скальпеля была заделана на торце(синяя стрелка) по всем направлениям движения, а для исследования вынужденных колебаний модели было придано гармоническое смещение по оси Х на торце(желтая стрелка)

Первая собственная форма Вторая собственная форма Напряжения при гармонических смещениях нижнего края, равных 0,1*sin(9*time) 25

В ЫВОДЫ В ходе работы были исследованы и смоделированы системы с двумя степенями свободы, что позволяет оценить колебания вискера и иглы. Была исследована модель нано весов и найдена зависимость массы, взвешиваемого объекта, от получившейся частоты колебаний. Получены деформации и напряжения при вынужденных и собственных колебаниях сложных углеродных структур. В дальнейшем, планируется найти теоретические соотношения для определения местоположения массы на весах и получить результаты для большего количества сложных нано конструкций. 26

С ПАСИБО ЗА ВНИМАНИЕ ! 27

Э КСПЕРИМЕНТАЛЬНАЯ УСТАНОВКА С ПЬЕЗОТРУБКОЙ

29

30

31