Динамика чувствительных элементов пороговых устройств автоматики. Учебное издание
.pdf
Рис. 4.4. Качественный вид скелетной кривой динамической системы пороговых датчиков инерционного типа с упругими элементами отрицательной жесткости
Анализ построенной скелетной кривой для нелинейной консервативной системы с упругими элементами отрицательной жесткости показывает, что данная система не является резонансной, так как у подобных систем с увеличением перемещения инерционного элемента возвращающая сила уменьшается в отличие от систем с упругими элементами положительной жесткости, поэтому отсутствуют частоты, при которых перемещения инерционного элемента стремятся к бесконечно большим значениям.
4.4. Устойчивость пороговых датчиков инерционного типа к вибрационным нагружениям
Для проведения расчетов по оценке устойчивости пороговых датчиков инерционного типа к вибрационным нагружениям необходимо получить уравнение движения их чувствительных элементов при действии ускорения, временная зависимость которого имеет вид
81
j t J sin t, |
(4.18) |
где J [м/с2] – амплитуда колебаний; Ω [рад/с] – угловая частота колебаний.
Частота колебаний ν, Гц определяется по формуле
|
|
. |
(4.19) |
|
|||
|
2 |
|
|
Примем условие, что амплитуда колебаний J всегда превышает j0, т. е. порог по ускорению срабатывания датчиков.
Под потерей устойчивости датчиков будем понимать перемещение их чувствительных элементов при действии вибрации на величину межконтактного зазора . Поэтому при расчетах будем определять амплитуду колебаний J при фиксированных значениях частоты колебаний в диапазоне от 200 до 20000 Гц, при которых выполняется условие q(t) = . Получим уравнение движения чувствительного элемента датчика, у которого силовая характеристика упругого элемента описывается уравнением F(q) = F0 + cq.
По аналогии с уравнением (2.3) запишем уравнение
движения чувствительного элемента датчика в виде |
|
q 02q J sin t j0. |
(4.20) |
Уравнение (4.20) является неоднородным дифференциальным уравнением.
Решение дифференциального уравнения (4.20) представляет собой сумму общего решения соответствующего однородного дифференциального уравнения и частного решения:
qобщ qодн qчаст. |
(4.21) |
Однородное дифференциальное уравнение имеет вид |
|
q 02q 0. |
(4.22) |
Общее решение однородногоуравнения (4.22) запишется ввиде
qодн С1 sin 0t C2 cos 0t, |
(4.23) |
82 |
|
где С1 и С2 – произвольные постоянные, определяемые из начальных условий.
Теперь найдем частное решение qчаст. Исходя из вида правой части уравнения (4.20), частное решение будем искать в виде
qчаст Asin t B cos t C. |
(4.24) |
Так как
qчаст A cos t B sin t
и
qчаст A 2 sin t B 2 cos t,
то, подставив соответствующие выражения в (4.20), получим
A 2 sin t B 2 cos t 2 |
Asin t B cos t C |
|
0 |
|
|
J sin t j0. |
(4.25) |
|
Раскрыв в выражении (4.25) скобки и приравняв соответствующие члены в левой и правой частях, получим выражения для А, В, и С в виде
A 02 J 2 ,
В0,
Сj0 .
02
Выражение для частного решения запишется в виде
q |
|
|
J |
sin t |
j0 |
, |
|
2 |
2 |
2 |
|||||
част |
|
|
|
||||
|
|
0 |
|
|
0 |
|
поэтому решение для уравнения (4.20) запишется в виде
q t C |
sin t C |
|
cos t |
|
J |
sin t |
j0 |
. |
|
2 |
2 |
|
|||||
1 |
0 |
2 |
0 |
|
2 |
|||
|
|
|
|
0 |
|
|
0 |
|
(4.26)
(4.27)
83
Для определения постоянных С1 и С2 из начальных условий продифференцируем уравнение (4.27), в результате чего получим уравнение
q t C |
cos t C |
sin t |
J |
cos t. (4.28) |
||
2 |
2 |
|||||
1 0 |
0 2 0 |
0 |
|
|||
|
|
|
0 |
|
|
|
Рассмотрим начальные условия. Чувствительный элемент датчика начнет перемещаться с момента времени t0, когда ускорение j(t), действующее на датчик, достигнет значения j0, равного порогу по ускорению срабатывания датчика:
q t0 q t0 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
j |
|
|||||||
|
при |
t t0 , |
где |
t0 |
|
|
|
|
|
sin 1 |
0 |
, |
||||||||||||||||||||
|
|
|
J |
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
т. е. находится из уравнения |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
j0 J sin t0. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(4.29) |
|||||||
Из уравнений (4.27) и (4.28) с учетом начальных усло- |
||||||||||||||||||||||||||||||||
вий составим систему уравнений |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
С sin t |
0 |
C |
2 |
cos t |
0 |
|
|
|
J |
|
|
sin t |
0 |
|
|
j0 |
0, |
|
||||||||||||||
|
2 |
|
|
2 |
2 |
|
||||||||||||||||||||||||||
|
1 |
|
0 |
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
0 |
|
|
(4.30) |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
J |
|
|
|
|
|
|
|
|
|
|
|
|
||
С |
cos t |
0 |
C |
|
sin t |
0 |
|
|
|
|
cos |
t |
0 |
0. |
|
|||||||||||||||||
2 |
|
2 |
|
|||||||||||||||||||||||||||||
|
1 |
0 |
|
|
0 |
|
|
2 |
0 |
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
В результате решения данной системы уравнений находим постоянные коэффициенты С1 и С2 и, подставив их в уравнение (4.27), после преобразований получим решение уравнения (4.20) в окончательном виде:
q t |
1 |
|
j0 02 j0 2 j0 2 cos 0 t t0 |
||||||
04 02 2 |
|
||||||||
|
|
|
|
|
|
|
|||
|
J 02 sin t J |
|
|
j2 |
0 sin 0 |
t t0 . (4.31) |
|||
|
1 |
|
0 |
||||||
|
J 2 |
||||||||
|
|
|
|
|
|
|
|
||
|
|
|
|
|
84 |
|
|
||
Для датчика № 1, имеющего параметры j0 = 15000 м/с2 |
|
(1500 g), m = 1 10–3 кг, c = 6560 Н/м, |
= 3 10–4 м, определим |
амплитуду колебаний J при частотах колебаний 500, 1000, 2000, 2500, 5000, 7500, 10000, 12500, 15000, 17500 и 20000 Гц,
при которых происходит срабатывание датчика, т. е. выполняется условие q(t) = . Результаты расчетов приведены в табл. 2.
На рис. 4.5 показан график, соответствующий значениям табл. 2. Видно, что полученная кривая делит всю область рассматриваемых параметров синусоидальных колебаний на две подобласти, расположенные выше и ниже кривой. Это означает, что если параметры синусоидальных колебаний лежат выше кривой, то при таких колебаниях датчик срабатывает, если ниже, то датчик сохраняет свою устойчивость.
Таблица 2
Частота |
Угловая частота |
Амплитуда колебаний |
|
колебаний , Гц |
колебаний Ω, рад/с |
J, м/с2 (g) |
|
500 |
3141,6 |
18900 |
(1890) |
1000 |
6283,2 |
23000 |
(2300) |
2000 |
12566,3 |
33000 |
(3300) |
2500 |
15707,9 |
39000 |
(3900) |
5000 |
31415,9 |
80000 |
(8000) |
7500 |
47123,9 |
135000 |
(13500) |
10000 |
62831,8 |
180000 |
(18000) |
12500 |
78539,8 |
230000 |
(23000) |
15000 |
94247,8 |
275000 |
(27500) |
17500 |
109955,7 |
330000 |
(33000) |
20000 |
125663,7 |
370000 |
(37000) |
85
Амплитуда колебаний J, м/с2
Частота колебаний Ω/2π, Гц
Рис. 4.5. Зависимость амплитуд синусоидальных колебаний корпуса датчика от частот, позволяющая определить параметры колебаний, при которых происходит срабатывание датчика
Получим уравнение движения чувствительного элемента датчика, у которого силовая характеристика упругого элемента описывается уравнением F(q) = F0 – cq, при действии вибрационных нагружений.
По аналогии с уравнением (2.30) запишем уравнение движения чувствительного элемента датчика в виде
q 02q J sin t j0. |
(4.32) |
Уравнение (4.32) является неоднородным дифференциальным уравнением.
Решение дифференциального уравнения (4.32) представляет собой сумму общего решения соответствующего однородного дифференциального уравнения и частного решения:
qобщ qодн qчаст. |
(4.33) |
Однородное дифференциальное уравнение имеет вид
86
q 02q 0. |
(4.34) |
Общее решение однородного уравнения (4.34) запишется как
q |
С е 0t С |
е 0t , |
(4.35) |
|
одн |
1 |
2 |
|
|
где С1 и С2 – произвольные постоянные, определяемые из начальных условий.
Теперь найдем частное решение qчаст. Исходя из вида правой части уравнения (4.32), частное решение будем искать в виде
qчаст Asin t B cos t C. |
(4.36) |
Так как
qчаст A cos t B sin t
и
qчаст A 2 sin t B 2 cos t,
то, подставив соответствующие выражения в (4.32), получим
A 2 sin t B 2 cos t 2 |
Аsin t B cos t C |
|
0 |
|
|
J sin t j0. |
(4.37) |
|
Раскрыв в выражении (4.37) скобки и приравняв соответствующие члены в левой и правой частях, получим выражения для А, В и С в виде
A02 J 2 ,
В0,
Сj0 .
02
Выражение для частного решения запишется в виде
q |
|
|
J |
sin t |
j0 |
, |
(4.38) |
|
2 |
2 |
2 |
||||||
част |
|
|
|
|
||||
|
|
0 |
|
|
0 |
|
|
|
|
|
|
87 |
|
|
|
|
поэтому решение для уравнения (4.32) запишется как
q t C e 0t C |
e 0t |
|
J |
sin t |
j0 |
. (4.39) |
||
2 |
2 |
2 |
||||||
1 |
2 |
|
|
|
||||
|
|
|
0 |
|
|
0 |
|
|
Для определения постоянных С1 и С2 из начальных условий продифференцируем уравнение (4.39), в результате получим уравнение
q t C e 0t C |
e 0t |
J |
cos t. (4.40) |
|||
2 |
2 |
|||||
1 0 |
2 |
0 |
|
|||
|
|
|
0 |
|
|
|
Рассмотрим начальные условия. Чувствительный элемент датчика начнет перемещаться с момента времени t0, когда ускорение j(t), действующее на датчик, достигнет значения j0, равного порогу по ускорению срабатывания датчика, т. е.
q t0 q t0 0 |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
j |
|
||||
|
при |
t t0 , |
где t0 |
|
|
|
sin 1 |
0 |
. |
||||||||||||
|
|
J |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Из уравнений (4.39) и (4.40) с учетом начальных усло- |
|||||||||||||||||||||
вий составим систему уравнений |
|
|
|
|
|
|
|
|
|
|
|
||||||||||
C e 0t |
C |
e 0t |
|
J |
|
|
sin t |
|
j0 |
0, |
|
|
|||||||||
2 |
|
|
2 |
|
2 |
|
|
||||||||||||||
|
1 |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
0 |
|
(4.41) |
||||
|
|
|
|
|
|
|
|
|
|
|
|
J |
|
|
|
|
|
|
|
||
|
C e |
0t |
C |
e |
0t |
|
|
|
cos t 0. |
|
|
||||||||||
|
|
|
|
2 |
|
2 |
|
|
|||||||||||||
1 |
0 |
|
|
2 |
|
0 |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
||
В результате решения данной системы уравнений находим постоянные коэффициенты С1 и С2 и, подставив их в уравнение (4.39), после преобразований получим решение уравнения (4.32) в окончательном виде:
88
|
|
|
1 |
|
|
2 |
|
2 |
|
|
|
q t |
|
|
|
|
j0 0 |
|
|
|
|
||
02 |
02 2 |
|
|
|
|||||||
j0 2 cosh 0 t t0 J 02 sin t |
|
||||||||||
|
|
j2 |
|
|
|
|
|
|
|
||
|
|
|
0 |
|
|
|
|
|
|
|
|
J 1 |
J |
2 0 sinh 0 |
t |
t0 . |
(4.42) |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
Решение уравнения (4.32) в аналитическом виде было |
|||||||||||
получено также В. П. Гандуриным. |
|
|
|
|
|
||||||
Для датчика № 2, имеющего параметры j0 = 15000 м/с2 |
|||||||||||
(1500 g), m = 1,6 10–4 кг, с = 5000 Н/м, |
|
= 3 10–4 м, определим |
|||||||||
амплитуды колебаний J при частотах колебаний 500, 1000, 2000, 2500, 5000, 7500, 10000, 12500, 15000, 17500 и 20000 Гц,
при которых происходит срабатывание датчика, т. е. выполняется условие q(t) = .
Результаты расчетов приведены в табл. 3.
Таблица 3
Частота |
Угловая частота |
Амплитуда колебаний |
|
колебаний , Гц |
колебаний Ω, рад/с |
J, м/с2 (g) |
|
500 |
3141,6 |
16880 |
(1688) |
1000 |
6283,2 |
20400 |
(2040) |
2000 |
12566,3 |
29600 |
(2960) |
2500 |
15707,9 |
35500 |
(3550) |
5000 |
31415,9 |
75000 |
(7500) |
7500 |
47123,9 |
112000 |
(11200) |
10000 |
62831,8 |
151000 |
(15100) |
12500 |
78539,8 |
190000 |
(19000) |
15000 |
94247,8 |
230000 |
(23000) |
17500 |
109955,7 |
270000 |
(27000) |
20000 |
125663,7 |
310000 |
(31000) |
|
89 |
|
|
На рис. 4.6 показан график, соответствующий значениям табл. 3. Видно, что полученная кривая делит всю область рассматриваемых параметров синусоидальных колебаний на две подобласти, расположенные выше и ниже кривой. Это означает, что если параметры синусоидальных колебаний лежат выше кривой, то при таких колебаниях датчик срабатывает, если ниже, то датчик сохраняет свою устойчивость.
Таким образом показано, что с помощью полученных уравнений (4.31) и (4.42) можно определять параметры вибрационных нагружений, при которых пороговые датчики инерционного типа сохраняют свою устойчивость. При этом уравнение (4.31) необходимо применять, если силовая характеристика упругого элемента датчика описывается уравнением F(q) = F0 + cq, а уравнение (4.42) – если силовая характеристика упругого элемента описывается уравнением
F(q) = F0 – cq.
Амплитуда колебаний J, м/с2
Частота колебаний Ω/2π, Гц
Рис. 4.6. Зависимость амплитуд синусоидальных колебаний корпуса датчика от частот, позволяющая определить параметры колебаний, при которых происходит срабатывание датчика
90
