CC BY
22
УДК 621.01
DOI: 10.17586/0021-3454-2016-59-10-843-846
РАЗРАБОТКА МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ ДИНАМИКИ УПРУГИХ МЕМБРАННЫХ ЧУВСТВИТЕЛЬНЫХ ЭЛЕМЕНТОВ
СИСТЕМ УПРАВЛЕНИЯ
А. С. Козлов, Р. Я. Лабковская, О. И. Пирожникова, В. Л. Ткалич
Университет ИТМО, 197101, Санкт-Петербург, Россия E-mail: [email protected]
Представлены математические модели динамики упругих мембранных чувствительных элементов, входящих в состав тензорезистивных микродатчиков давления. Моделирование по методу конечных элементов осуществлено в условиях комплексных механических и температурных нагрузок на тензорезистивный микродатчик давления. При построении моделей учитывались силы вязкого трения, связанные с наличием контакта чувствительного элемента с вязкоупру-гим клеевым слоем. Адекватность предложенных моделей подтверждена сравнением расчетных значений с данными экспериментальных исследований. Разработанные математические модели динамики упругих мембранных чувствительных элементов позволяют обеспечить высокие показатели надежности и хорошую конкурентоспособность этой элементной базы систем управления. Созданный пакет прикладных программ, предназначенный для работы в современных компьютерных средах, позволяет визуализировать результаты, полученные при математическом моделировании динамики упругих мембранных чувствительных элементов.
Ключевые слова: упругие мембранные чувствительные элементы, метод конечных элементов, комплексные воздействия, математические модели динамики, тензорезистивные микродатчики давления
В состав тензорезистивных микродатчиков давления входят упругие мембранные чувствительные элементы (УМЧЭ), реагирующие на изменение механических и температурных нагрузок. Эти элементы находятся в непосредственном контакте с вязкоупругим клеевым слоем, что приводит к необходимости учета их нелинейных свойств. Вязкий клеевой слой обеспечивает передачу УМЧЭ не только полезного сигнала (измеряемое давление), но и вызванных механическими и температурными нагрузками помех, которые приводят к деформационным изменениям тензорезисторов и приращению их сопротивлений [1].
При моделировании динамики УМЧЭ целесообразно применять метод конечных элементов, базирующийся на ансамбле конечных элементов треугольной или четырехугольной формы с последующим разложением матрицы деформаций по главным формам свободных колебаний [2—4].
Анализ математических моделей динамики УМЧЭ с использованием метода конечных элементов осуществлялся с учетом сил вязкого трения. Приведем уравнение динамики в общем виде [1]:
MU + CU + K1U = F(t, U), (1)
где F(t, U) = Q(t) + K^U, KH — матрица дополнительной жесткости, K1 — матрица изгиб-ных жесткостей УМЧЭ, М — матрица масс, С — диссипативная матрица, U — вектор обобщенных перемещений, U — вектор скоростей, U — вектор ускорений, Q(t) — вектор обобщенных силовых нагрузок.
Диссипативная матрица С может быть описана одним из следующих уравнений:
C1 =p1M + P2K;
С2 = Мф(М_1К) ; С = /уК,
где Р1 и Р2 — коэффициенты демпфирования Рэлея, ф(М _1К) — скалярная функция матрицы М" К, параметры которой подбираются экспериментально, К — матрица жесткости; у = у/ 2л, / = , у — коэффициент удельного рассеивания энергии в материале.
Собственные значения матрицы М_1К обозначим через у/ (/ = 1, 2,..., п), а координаты смещения каждого конечного элемента — через и1 (^). Тогда вектор перемещения представляется выражением [3]:
п
и(*) = ^уи, (I).
/=1
Координаты смещения и1 (^) каждого конечного элемента задаются уравнением:
С2 ё 2
-у + фе-- + Шг.
Л2 сИ
и/ (') = Ч/ (') =
где Ш/ — собственная частота свободных колебаний УМЧЭ, / = 1,2,..., п ; е е {1;2;3} — индекс скалярной функции матрицы М _1К в уравнении (1).
Передаточная функция может быть представлена выражением:
Ще (з) =
( о о2 У1
1 $ Фе $
1Л—— Л--
2 2
Ш Ш/ у
Ч/ (з )= Р (*) Ч/.
где ч1 — постоянные коэффициенты, $ — площадь срединной поверхности конечного элемента.
Тогда смещение конечного элемента УМЧЭ задается выражением:
и, (з) = Ж/ (з)ч/ш72р (з). При этом вектор перемещений УМЧЭ принимает вид:
и (з ) = ]Г^ (з )Ч/ Шг"2У гР (з ). /=1
Результирующий вектор деформации УМЧЭ представляется выражением:
Е( х, з) = £ Ер (х)Щ (з) р(з ),
/=1
где Е;р (х) — парциальные деформации УМЧЭ в точке X, вызванные единичным статическим давлением.
Вычислить векторы Е;р(х) можно с помощью алгоритма [3], который базируется на
спектральном расчете динамики УМЧЭ. Алгоритм позволяет проанализировать динамику УМЧЭ в случае воздействия импульса при монотонно возрастающем действии давления, которое воспринимается микродатчиком в качестве статической нагрузки.
Экстраполируя приведенные выражения для варианта комплексного нагружения УМЧЭ, запишем:
Е(х, з) = ¿Ж (з )[Егр (х)р(з) + Е/а (х)а(з) + ЕгТ (х)Т (з)], /=1
где p(s) — давление; a(s) — ускорение; T(s)— температурное поле; Eia (x), EiT (x) — парциальные деформации УМЧЭ, вызванные ускорением и единичным приращением температуры.
Рассчитаем отклик УМЧЭ на воздействие длительной импульсной объемной нагрузки р( х)Ф( x), где р( x) — плотность материала УМЧЭ, Ф( x) — вектор-функция, связанная с определением ускорений a (x, t) = Ф( x)a(t) по объему УМЧЭ.
Математическая модель динамики, которая предусматривает учет таких физических характеристик материалов УМЧЭ, как модули упругости и коэффициенты теплового расширения, а также удельные коэффициенты рассеивания энергии вязкоупругого клеевого слоя, определяемые для циклического одноосного нагружения образцов клея, тензорезисторных микродатчиков давления, описывается выражением:
E(x, s) = Wp (x, s) p(s) + Wa (x, s)a(s) + Wt (x, s )T (s),
n
где Wa (x, s) = (s)Eia (x); a = P, a, T, Wp (x, s), Wa (x, s), WT (x, s) — комплексные переда-
i=1
точные функции УМЧЭ для давления, ускорения и температуры. Далее рассчитываются механические напряжения a(x, t) [4]:
a=Es,
где ст — вектор механических напряжений, е — вектор перемещений, E — матрица упруго-стей УМЧЭ.
Таким образом, разработанная математическая модель динамики УМЧЭ тензорезистив-ных микродатчиков давления представляется общей передаточной функцией в комплексной форме, что позволяет перейти к оценке динамических погрешностей измеряемого давления в случае воздействия переменных механических и температурных нагрузок.
Полученные в работе модели и алгоритмы УМЧЭ реализованы в пакетах прикладных программ, работающих в современных компьютерных средах, что позволяет визуализировать полученные результаты, сделав предлагаемый математический аппарат удобным и наглядным для пользователя.
список литературы
1. Лабковская Р. Я., Пирожникова О. И., Ткалич В. Л. Метод повышения надежности упругих чувствительных элементов систем управления и автоматики // Научно-технический вестник СПбГУ ИТМО. 2011. № 1(71). С. 136—138.
2. Лабковская Р. Я., Ткалич В. Л., Пирожникова О. И. Разработка библиотеки конечных элементов для САПР упругих конструкций герконов // Изв. вузов. Приборостроение. 2013. Т. 56, № 3. С. 21—24.
3. Лабковская Р. Я., Ткалич В. Л. Библиотека конечных элементов в приложении к упругим чувствительным элементам пластин и мембран датчиков систем управления // Научная перспектива. Уфа, 2010. № 3—4. С. 86—89.
4. Лабковская Р. Я. Исследование статики и динамики мембранных и пластинчатых электромеханических элементов систем управления // Сб. докл. конгр. молодых ученых. СПб: НИУ ИТМО, 2014. Вып. 1. С. 202—203.
Сведения об авторах
Алексей Сергеевич Козлов — аспирант; Университет ИТМО, кафедра проектирования и безопас-
ности компьютерных систем; E-mail: [email protected] Римма Яновна Лабковская — канд. техн. наук; Университет ИТМО, кафедра проектирования и
безопасности компьютерных систем; E-mail: [email protected] Ольга Игоревна Пирожникова — канд. техн. наук; Университет ИТМО, кафедра проектирования и
безопасности компьютерных систем; E-mail: [email protected]
Вера Леонидовна Ткалич — д-р техн. наук, профессор; Университет ИТМО, кафедра проектиро-
вания и безопасности компьютерных систем; E-mail: [email protected]
Рекомендована кафедрой Поступила в редакцию
проектирования и безопасности 26.02.16 г.
компьютерных систем
Ссылка для цитирования: Козлов А. С., Лабковская Р. Я., Пирожникова О. И., Ткалич В. Л. Разработка математических моделей динамики упругих мембранных чувствительных элементов систем управления // Изв. вузов. Приборостроение. 2016. Т. 59, № 10. С. 843—846.
DEVELOPMENT OF MATHEMATICAL MODELS OF DYNAMICS OF ELASTIC MEMBRANE SENSITIVE ELEMENTS OF CONTROL SYSTEMS
A. S. Kozlov, R. Ya. Labkovskaya, O. I. Pirozhnikova, V. L. Tkalich
ITMO University, 197101, St. Petersburg, Russia E-mail: [email protected]
Mathematical models of dynamics of elastic membrane sensitive elements in piezoresistive pressure microsensors are proposed. Modelling is carried out with the use of the method of finite elements by representing the membrane elastic sensitive element by an ensemble of triangular or quadrangular forms. The models account for viscous friction forces associated with the contact between the sensitive element and a viscoelastic adhesive layer. To confirm the proposed models adequacy, calculated model data are compared to results of experimental study. The developed mathematical models for the dynamics of elastic membrane sensitive elements are reported to provide high reliability and good competitiveness of this element base management systems. Created software package designed to work in modern computer environments allows to visualize the results obtained by mathematical modeling of the dynamics of elastic membrane of sensitive elements.
Keywords: membrane elastic sensitive elements, finite element method, complex influences, mathematical models of dynamics, piezoresistive pressure microsensors
Data on authors
Aleksei S. Kozlov — Post-Graduate Student; ITMO University, Department of Computer
Systems Design; E-mail: [email protected] Rimma Ya. Labkovskaya — PhD; ITMO University, Department of Computer Systems Design;
E-mail: [email protected] Olga I. Pirozhnikova — PhD; ITMO University, Department of Computer Systems Design;
E-mail: [email protected] Vera L. Tkalich — Dr. Sci., Professor; ITMO University, Department of Computer Systems
Design; E-mail: [email protected]
For citation: Kozlov A. S., Labkovskaya R. Ya., Pirozhnikova O. I., Tkalich V. L. Development of mathematical models of dynamics of elastic membrane sensitive elements of control systems // Izv. vuzov. Pri-borostroenie. 2016. Vol. 59, N 10. P. 843—846 (in Russian).
DOI: 10.17586/0021-3454-2016-59-10-843-846
