УДК 621.3.06 + 621.01
г
\
БО! 10.21685/2307-4205-2016-4-5
ДИНАМИЧЕСКАЯ МОДЕЛЬ МИКРОПЕРЕКЛЮЧАТЕЛЕЙ С УПРУГИМИ СТЕРЖНЕВЫМИ ЭЛЕМЕНТАМИ
И. Е. Люминарский, Е. С. Люминарская
Введение
В настоящее время микропереключатели широко используются в автомобилях, авиации, космической технике и т.д. Важнейшим требованием, предъявляемым к разрывным электрическим контактам, является безотказность их работы. Отказы в электрических контактах микропереключателей могут происходить из-за тепловых, электрических, химических и механических явлений. Отказы в слаботочных контактах микропереключателей стержневого типа в основном появляются из-за поломки упругого элемента или распорной пружины, залипания подвижного контакта, дребезга контактов при коммутации, внезапного изменения силы прижатия контактов, загрязнения контактных поверхностей и др.
Важнейшим показателем работы микропереключателей является быстродействие, которое определяется временем переключения и дребезгом контактов. Дребезг может вызывать сбои в работе электронной аппаратуры. Поэтому при проектировании микропереключателей необходимо не только определять время переключения, но и исследовать дребезг контактов. В микропереключателях быстродействие обеспечивается рациональным подбором основных параметров.
Для создания надежных миниатюрных переключателей необходимо разрабатывать математические модели, учитывающие основные процессы, протекающие при коммутации, и использовать их для выявления причин и предвестников сбоев и отказов на этапе проектирования.
Существующие методики расчета колебаний контактных систем с упругими стержневыми элементами и упругими ограничениями основываются на использовании метода конечных элементов [1, 2]. В контактных системах упругие ограничители имеют высокую контактную жесткость. При ударе упругого элемента об ограничитель в нем возбуждаются формы колебаний высокой частоты, которые необходимо учитывать при исследовании дребезга контактов. Поэтому при расчете колебаний стержневого упругого элемента с жесткими ограничителями используют малый шаг интегрирования по времени и большое количество конечных элементов, что значительно увеличивает машинное время счета.
Целью данной работы является разработка методики динамического анализа микропереключателей с упругими стержневыми элементами, позволяющей исследовать их работу в момент переключения и дребезга контактов. Методика основана на применении метода Бубнова -Галеркина и использовании собственных форм упругого элемента с подвижным контактом в качестве базисных функций.
В статье предлагается динамическая модель микропереключателя с упругим стержневым элементом, схема которого представлена на рис. 1. В этой модели принимаются следующие допущения: 1) корпус и толкатель являются абсолютно твердыми телами; 2) положение касания контакта бобышки с неподвижными контактами не изменяется; 3) деформации упругого элемента и распорной пружины подчиняются закону Гука; 4) инерция распорной пружины не учитывается; 5) инерция поверхностных слоев контактов не учитывается, а контактные деформации определяются по формулам Герца.
Постановка задачи
Математическая модель
Рис. 1. Схема микропереключателя: 1, 2 - неподвижные контакты; 3 - толкатель; 4 - распорная пружина;
5 - упругий элемент; 6 - бобышка подвижного контакта
Быстродействие микропереключателей обеспечивается подбором таких параметров, при которых все положения равновесия упругого элемента с распорной пружиной без опоры на неподвижные контакты неустойчивы. Неустойчивость положений равновесия обеспечивается малыми углами наклона распорной пружины а (см. рис. 1). В этом случае даже при небольших перемещениях поперечных сечений вдоль оси y продольная составляющая силы со стороны распорной пружины Sx будет оказывать значительное влияние на движение упругого элемента. Поэтому в предлагаемой математической модели микропереключателя учитывается влияние Sx на деформацию упругого элемента.
Схема сил, действующих на упругий элемент, приведена на рис. 2. Продольная составляющая силы со стороны распорной пружины Sx приближенно учитывается поперечной силой P* и изгибающим моментом m*. При таком учете продольной силы Sx нелинейная зависимость внутреннего изгибающего момента от продольной координаты M (x) заменяется линейной M*(x) = -P*-(x4-x)-m* на участке (0, x4). Коэффициенты P* и m* определяются методом наименьших квадратов из условия
x4
Ф = J (M (x) -M*(x))2dx = min, (1)
0
где MM (x) = Sx (y - У4), при x < x4.
Рис. 2. Расчетная схема
После несложных преобразований получается формула для определения P* и m*
О, 5х4 x4
1 0,5x2
р*
т*
(
5
Л
Х4 >"4 - I У^х 0
5
Х4 >4
- х4 I ydx + I у • хёх
0
0
(2)
Для динамического анализа микропереключателя используется дифференциальное уравнение движения балки при изгибе с учетом потерь на внутреннее трение [3]
т.
0 д(2 +Эх2
Ы-т У + Р
I д>
дх2 Г ' дt
= Ч( х, t),
(3)
где у^, х) - функция перемещений поперечных сечений балки; ш0 - погонная масса балки; Е - модуль упругости I рода; J - момент инерции поперечного сечения; Р - линеаризованный параметр диссипации; д( х, t) - интенсивность силы, действующей на упругий элемент. Интенсивность д( х, ^ может быть представлена следующим образом (см. рис. 2):
д(х,t) = -8(х - х4) - Р* ^8(х - х4) - т* •&'(х - х4) + Е^) ^8(х - х1) +
д 2
+ ЕМ) ^8(х - х2) - ЕМ) -8(х - х3) - т —> ^8(х - х2),
дt2
(4)
д 2
где 8(х - х,) - дельта-функция Дирака; т —> ^8(х - х2) - сила инерции бобышки; т - масса бо-
дt2
бышки.
Вертикальная составляющая силы от распорной пружины определяется из выражения 5> = 5 • 8т(а) = к (/° /) • (>4 + а) = г • (>4 + а),
(5)
где к - коэффициент жесткости распорной пружины; 10 - длина недеформированной распорной пружины; I - длина деформированной распорной пружины.
Решение уравнения (4) можно приближенно представить в следующем виде [4, с. 170]:
N
у(/, х) = X С ^) и (х), у=1
(6)
где и у (х) - нормированные по кинетической энергии собственные формы упругого элемента с бобышкой [3]; Су ^) - неизвестные функции; N - число учитываемых собственных форм.
Неизвестные функции Су ^) определяются из обыкновенных дифференциальных уравнений, которые получаются путем подстановки решения (6) в уравнение (4) и последующего скалярного умножения полученного уравнения на нормированные собственные формы и у (х) [3, 5].
Так как скалярное произведение двух различных собственных форм равно нулю, то все дифференциальные уравнения получаются независимыми. Для определения функции Су ^) необходимо
решить одно дифференциальное уравнение второго порядка
СС^) + 2п • С,^) + ю2 • С,(0 = г • (у^,хА) + а) • и,(х4) + + т)и,(х1)-Е2^)иг(х2)-Ез(0и,(хз)-Р* •и,(х4) + т* • иг'(хА).
(7)
где } = 1,...,N ; ю, - собственные частоты упругого элемента с бобышкой; п = 0,5 Рю2 - коэффициент затухания; а - расстояние, показанное на рис. 1.
Собственные формы упругого элемента с бобышкой определялись методом последовательных приближений [4, с. 180; 6, с. 176-189].
Силы реакций F1 и F2 зависят от формы, шероховатости, запыленности и сближения контактных поверхностей. Значения этих сил определяются по методике, изложенной в работе [7].
Сила, действующая на упругий элемент со стороны толкателя, F3 (t), вычисляется по формуле Герца.
Сближения контактных поверхностей зависят от перемещения точек упругого элемента yk
А *
и начальных зазоров Ак
ДДО = Ai* + y(t,xi), A2(t) = Д2 -y(t,x2), Дз(0 = Д3 -y(t,хз).
Для решения дифференциальных уравнений (7) участок интегрирования (0, tk) разбивается на M равных интервалов. На каждом интервале переменные величины r, y(t, x4), F1, F2, F3, P*, m* считаются постоянными и принимаются равными средним значениям на шаге интегрирования.
В начальном положении толкатель расположен в крайнем верхнем положении и упирается в ограничитель, а упругий элемент и бобышка находятся в положении равновесия. Расчет выполняется при постоянной скорости толкателя, который движется из крайнего верхнего в крайнее нижнее положение.
dc
В начальный момент времени производные —'- равны нулю, а функции Ct определяются
dt
из уравнений (7), в которых все производные по времени приравниваются к нулю. Значения величин, входящих в правую часть дифференциальных уравнений (7), определяются из расчета равновесия упругого элемента.
Более подробно последовательность динамического анализа упругой системы, ограниченной упругими односторонними связями, изложена в [5].
Результаты исследования
Для проверки предложенной математической модели были проведены теоретические и экспериментальные исследования работы микропереключателя МП9, показанного на рис. 3.
Рис. 3. Микропереключатель МП9
На рис. 4 показана зависимость силы прижатия контактов от времени в момент коммутации. Из графика видно, что при коммутации происходит дребезг контактов. В момент дребезга наблюдается четыре отрыва подвижного контакта от неподвижного контакта. Первый отрыв происходит в течение 0,12 мс, второй - 0,005 мс, третий - 0,041 мс, четвертый - 0,021 мс. Вычисления при различной жесткости распорной пружины показали, что при дребезге в основном происходит два-три отрыва со временем отрыва более to > 0,02 мс.
о
т. мс ->
а!
5.8 5.9 6 6.1 6.2 б.З 6.4 6.5
Рис. 4. Изменение силы прижатия контактов при коммутации
В исследуемом микропереключателе сила прижатия контактов Е1 в положении равновесия (номинальная) составляет 0,64 Н. При первом взаимодействии максимальное значение этой силы превышает номинальное значение более чем в 10 раз, а время соприкосновения составляет менее 7 мкс (см. рис. 3). За такое короткое время в упругом элементе не успевают появиться большие деформации и напряжения [8, 9].
Расчетные исследования показали зависимость времени переключения от скорости движения толкателя. Поэтому для проверки адекватности математической модели была получена зависимость времени переключения от времени движения толкателя. Время движения толкателя определялось от начала движения толкателя до момента размыкания контактов 2 и 6 (см. рис. 1).
На рис. 5 приведены теоретическая и регрессионная зависимости времени переключения tп от времени движения толкателя ^ . Время переключения tп определяется от начала размыкания контактов 2 и 6 до момента первого касания контактов 1 и 6 (см. рис. 1). При определении экспериментальной зависимости использовалась квадратичная модель регрессии. В исследуемом диапазоне скоростей толкателя отклонение теоретических и экспериментальных значений времени переключения не превышает 5 %.
Рис. 5. Зависимость времени переключения от времени движения толкателя: 1 - теоретическая зависимость; 2 - регрессионная зависимость
Предложенная динамическая модель позволяет исследовать работу микропереключателя в период дребезга контактов. Для проверки адекватности модели в этот период работы исследовалась зависимость времени первого (наибольшего) отрыва при дребезге контактов 1 и 6 to1 от времени переключения tп . На рис. 6 показаны теоретическая и регрессионная зависимости. Для построения регрессионной зависимости также использовалась квадратичная модель регрессии.
Сравнение приведенных зависимостей показывает адекватность модели в фазе дребезга контактов. Расхождение теоретических и экспериментальных зависимостей не превышает 10 %.
Необходимо отметить, что при дребезге появляются два-три больших отрыва (с временем отрыва более 10 мкс). В экспериментальных исследованиях был получен большой разброс значений времени второго отрыва (15...60 мкс). Большой разброс значений может быть связан со случайными вибрациями корпуса, возникающими при коммутации.
Рис. 6. Зависимость времени первого отрыва при дребезге контактов 1 и 6 от времени переключения: 1 - регрессионная зависимость; 2 - теоретическая зависимость
Выводы
1. Предложена математическая модель микропереключателей с упругими стержневыми элементами, позволяющая исследовать их работу в фазах переключения и дребезга контактов.
2. Проведенные теоретические и экспериментальные исследования подтвердили адекватность предложенной модели.
3. Установлено, что время переключения зависит от скорости движения толкателя.
4. Расчетным и экспериментальными путями доказано, что при дребезге контактов происходит три-четыре отрыва со временем отрыва от 10 до 140 мкс.
Список литературы
1. Контактные системы / В. В. Осташавичюс, Б. В. Рудгальвис, В. Л. Рагульскине, Б. П. Бакшис. - Л. : Машиностроение, 1987. - 279 с.
2. Уланов, А. В. Математические модели динамики, алгоритмы и информационно-измерительные средства виброакустической диагностики и неразрушающего контроля контактных устройств авиационных приборов : дис. ... канд. физ.-мат. наук : 05.11.16 и 05.11.13 / Уланов А. В. - Ижевск : РГБ, 2006. - 178 с.
3. Вибрация в технике : справ. : в 6 т. - М. : Машиностроение, 1981. - Т. 6. Защита от вибрации и ударов / под ред. К. В. Фролова. - 456 с.
4. Вибрация в технике : справ. : в 6 т. - М. : Машиностроение, 1978. - Т. 1. Колебания линейных систем / под ред. В. В. Болотина. - 352 с.
5. Люминарский, И. Е. Расчет упругих систем с односторонними связями : моногр. / И. Е. Люминарский. -М. : МГИУ, 2006. - С. 38-99.
6. Коллац, Л. Задачи на собственные значения / Л. Коллац. - М. : Наука : Гл. ред. Физ.-мат. лит., 1968. -504 с.
7. Люминарская, Е. С. Математическая модель слаботочных электрических контактов при сбое / Е. С. Лю-минарская, В. Н. Дианов // Динамика сложных систем - ХХ1 век. - 2014. -Т. 8, № 1.- С. 73-78.
8. Интегродифференциальная кодоимпульсная модуляция в задачах повышения надежности в аппаратуре / В. Н. Дианов, Т. А. Гевондян, И. М. Белоусов, Е. С. Люминарская // Труды Международного симпозиума Надежность и качество. - 2012. - Т. 2. - С. 136-137.
9. Дианов, В. Н. Диагностика интегродифференциальных сбоев в активных элементах интегральной схемотехники / В. Н. Дианов // Труды Международного симпозиума Надежность и качество. - 2014. -Т. 2. - С. 37-40.
Люминарский Игорь Евгеньевич доктор технических наук, профессор, кафедра теории механизмов и машин, Московский государственный технический университет им. Н. Э. Баумана (105005, Россия, г. Москва, 2-я Бауманская ул., 5, стр. 1) Е-шаИ: [email protected]
Люминарская Екатерина Станиславовна
кандидат технических наук, доцент, кафедра электротехники и промышленной электроники, Московский государственный технический университет им. Н. Э. Баумана (105005, Россия, г. Москва, 2-я Бауманская ул., 5, стр. 1) Е-шаП: 1уиштагека| [email protected]
Аннотация. Актуальность и цели. Микропереключатели с упругими стержневыми элементами широко используются в авиационной, космической, автомобильной и другой технике. Отказы в слаботочных контактах микропереключателей стержневого типа в основном появляются из-за поломки упругого элемента или распорной пружины, залипания подвижного контакта, дребезга контактов при коммутации, внезапного изменения силы прижатия контактов, загрязнения контактных поверхностей. В связи с этим появляется необходимость выявления предвестников отказов и сбоев на этапе проектирования и выбора основных параметров микропереключателей с учетом возможных загрязнений. Перспективным направлением в решении поставленной задачи является математическое моделирование, позволяющее исследовать работу микропереключателей при переключении и дребезге контактов в реальных условиях эксплуатации. Материалы и методы. Предложенная методика основана на применении метода Бубнова - Галеркина и использовании собственных форм упругого элемента с подвижным контактом в качестве базисных функций. Результаты. Разработанная динамическая модель микропереключателей позволила установить зависимость времени переключения от скорости движения толкателя и появление при дребезге контактов трех-четырех размыканий со временем отрыва подвижного контакта 10-140 мкс. Адекватность математической модели подтверждается экспериментальными данными. Выводы. Разработка динамических и математических моделей переключателей дает возможность учитывать основные процессы, протекающие в них при коммутации, и использовать их для выявления причин сбоев на этапе проектирования.
Ключевые слова: микропереключатель, дребезг, быстродействие, стержневой упругий элемент, собственные формы.
Lyuminarskji Igor Evgenevich
doctor of technical science, professor,
sub-departmen of theory of mechanisms and machines,
Moscow State Technical University
named after N. E. Bauman
(105005, 2-nd Baumanskaya street, apartment 5,
building 1, Moscow, Russia)
Lyuminarskaja Ekaterina Stanislavovna
candidate of technical science, associate professor,
sub-departmen of electrical engineering
and industrial electronics,
Moscow State Technical University
named after N. E. Bauman
(105005, 2-nd Baumanskaya street, apartment 5,
building 1, Moscow, Russia)
Abstract. Relevance and purpose. Microswitches with elastic core elements are widely used in aerospace, automotive, etc. engineering. The bounce in low-voltage contacts of the microswitch rod type mostly appear because of a breakage of the elastic element or spacer springs, sticking the movable contact to contact bounce when switching, sudden change of the pressing force of the contact, contamination of contact surfaces. In this regard, there is the need to identify precursors of failures and failures at the design stage and selection of basic parameters of switches with consideration of possible contamination. A promising direction in solving this problem is mathematical modeling, which allows to investigate the operation of the switches when switching and contact bounce in real conditions. Materials and methods. The proposed method based on the method of Bubnov-Galerkin and use their own forms of the elastic element with a movable contact as basis functions. Results. Developed the dynamic model of the micro switches has allowed to establish the dependence of the switching time of the speed of movement of a pusher and a contact bounce three or four razmiene with the time separation of the movable contact 10-140 ISS. The adequacy of the mathematical model is confirmed by experimental data. Conclusions. Development of dynamic and mathematical models of switches takes into account basic processes proceeding in them at switching, and use them to identify the causes of failures at the design stage.
Key words: microswitch, bounce, speed, the rod elastic element, their own forms.
УДК 621.3.06 + 621.01
Люминарский, И. Е.
Динамическая модель микропереключателей с упругими стержневыми элементами / И. Е. Люминарский, Е. С. Люминарская // Надежность и качество сложных систем. - 2016. - № 4 (16). - С. 35-42. DOI 10.21685/2307-4205-2016-4-5.