Динамика чувствительных элементов пороговых устройств автоматики. Учебное издание
.pdfr1 0 i 0 и |
r2 0 i 0. |
Таким образом, корни характеристического уравнения (2.10) различные и при этом комплексные, поэтому паре комплексных корней уравнения (2.7) соответствуют два частных решения
q e 0 i 0 t |
и |
q e 0 i 0 t . |
(2.11) |
1 |
|
2 |
|
Вместо этих двух решений рассмотрим некоторые их комбинации, которые также будут решениями уравнения (2.7) вследствие следующих теорем [7] о свойствах частных решений линейных дифференциальных уравнений (для рассматриваемого случая).
Теорема 1. Если функция q1 является решением уравнения (2.7), то и функция Cq1 есть решение этого уравнения, где C – произвольная постоянная.
Теорема 2. Если функции q1 и q2 являются решениями уравнения (2.7), то и функция q1 + q2 есть решение этого уравнения.
На основании данных теорем рассмотрим функции
q |
q1 q2 |
(2.12) |
||
1 |
|
|
2 |
|
и |
|
|
|
|
q1 q2 . |
|
|||
q |
(2.13) |
|||
2 |
|
|
2i |
|
|
|
|
|
|
Применяя известные формулы Эйлера [5] |
||||
cos x |
eix |
e ix |
, |
|
|
2 |
|||
|
|
|
|
|
sin x |
eix |
e ix |
, |
|
|
2i |
|||
|
|
|
|
|
преобразуем выражения (2.12) и (2.13) с учетом (2.11) следующим образом:
31
q |
|
e 0 i 0 t e 0 i 0 t |
e0t |
ei 0t e i 0t |
e0t cos t, |
|
|
||||
1 |
2 |
2 |
0 |
||
|
|
||||
q2 e 0 i 0 t 2ie 0 i 0 t e0t ei 0t 2ie i 0t e0t sin 0t.
Так как e0t 1 , то общее решение однородного уравнения (2.7) запишется в виде
qодн С1 sin 0t C2 cos 0t, |
(2.14) |
где С1 и С2 – произвольные постоянные, определяемые из начальных условий.
Теперь найдем частное решение qчаст. Исходя из вида правой частиуравнения (2.3), частное решение будем искать в виде
qчаст Аt B. |
(2.15) |
|
Так как |
|
|
qчаст A |
и qчаст 0, |
|
то, подставив соответствующие выражения в (2.3), получим уравнение
0 2 |
At B kt j . |
(2.16) |
0 |
0 |
|
Раскрыв скобки и приравняв соответствующие части, найдем
A |
k |
, |
B |
j0 |
. |
2 |
2 |
||||
|
|
|
|
||
|
0 |
|
|
0 |
|
Тогда частное решение запишется в виде
q |
|
kt j0 |
. |
(2.17) |
|
||||
част |
|
2 |
|
|
|
0 |
|
|
|
Общее решение неоднородного дифференциального уравнения (2.3) запишется в виде
q t C |
sin t C |
2 |
cos t kt j0 . |
(2.18) |
|
1 |
0 |
0 |
2 |
|
|
|
|
|
|
0 |
|
|
32 |
|
|
|
|
Для нахождения постоянных коэффициентов С1 и С2 продифференцируем уравнение (2.18), получим
q t C cos t C |
2 |
sin t |
k |
. |
|
(2.19) |
|||
2 |
|||||||||
0 1 |
0 |
0 |
0 |
|
|
|
|||
Введем начальные условия |
|
|
0 |
|
|
|
|||
|
|
|
j0 |
|
|
||||
q t0 q t0 0 |
при |
t t0 , |
где t0 |
. |
|
||||
|
|
||||||||
|
|
|
|
|
|
k |
|
||
Данные начальные условия определяют, что движение инерционного элемента начинается в тот момент, когда выполнено условие j0 = kt0.
Из уравнений (2.18) и (2.19) с учетом начальных условий составим систему уравнений
|
0 C |
sin t |
|
C |
|
cos t |
|
|
kt0 |
j0 |
, |
||||||||||
|
0 |
2 |
0 |
|
|
|
|||||||||||||||
|
1 |
0 |
|
|
|
|
|
|
0 |
|
|
|
2 |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
(2.20) |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
0 |
|
C cos t |
0 |
|
C |
2 |
sin t |
0 |
|
k |
. |
||||||||||
2 |
|||||||||||||||||||||
|
|
0 |
1 |
|
0 |
|
|
0 |
|
|
|
0 |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|||
В результате решения данной системы уравнений находим постоянные коэффициенты С1 и С2:
|
|
|
|
|
|
C |
|
|
k |
cos t |
, |
|
|
|||||
|
|
|
|
|
|
|
3 |
|
|
|||||||||
|
|
|
|
|
|
1 |
|
|
|
|
0 0 |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
||
|
|
|
|
|
|
C |
2 |
|
|
k |
|
sin t |
. |
|
|
|
||
|
|
|
|
|
|
3 |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
0 0 |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
Подставляя С1 и С2 в уравнение (2.18), получим |
|
|||||||||||||||||
q t |
k |
cos t |
0 |
sin t |
|
|
k |
|
sin t |
0 |
cos t kt j0 |
. (2.21) |
||||||
3 |
|
3 |
||||||||||||||||
|
0 |
0 |
|
|
0 |
|
0 |
2 |
|
|||||||||
0 |
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
0 |
|
||
Преобразовав уравнение (2.21) к виду |
|
|
|
|
||||||||||||||
|
|
|
k |
|
|
|
|
|
|
|
|
|
|
|
|
kt j0 |
||
q t |
|
cos 0t0 sin |
0t sin 0t0 cos 0t 2 |
|||||||||||||||
3 |
||||||||||||||||||
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
33 |
|
|
|
|
|
|||
и используя тождество [5] |
sin α sin α cos sin cos α, |
|||||||||||||||||||||||
получим уравнение (2.21) в виде |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
q t |
|
k |
sin |
|
|
|
t |
|
j0 |
|
|
kt j0 |
. |
|
(2.22) |
|||||||||
|
3 |
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
0 |
|
|
|
|
k |
|
|
2 |
|
|
|
||||||
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
||||||
При межконтактном зазоре |
|
|
|
|
время tср срабатывания |
|||||||||||||||||||
датчика можно определить, численно решив уравнение |
|
|||||||||||||||||||||||
|
k |
|
|
|
|
|
|
|
|
|
|
|
|
|
j |
|
|
ktср j0 |
|
|
||||
|
|
|
sin |
|
|
|
t |
ср |
|
|
0 |
|
|
|
|
|
|
. |
(2.23) |
|||||
3 |
|
|
2 |
|
||||||||||||||||||||
|
|
|
|
0 |
|
|
|
|
k |
|
|
|
|
|||||||||||
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|||||
Из уравнения (2.23) можно получить аналитическую зависимость (с некоторым приближением) времени tср срабатывания датчика от параметров j0, k и величины межконтактного зазора при нагружении линейно возрастающим ускорением.
Для этого используем формулу разложения функции
sin(z) в ряд Тейлора sin z z |
|
z3 |
|
... [5] (используя первые |
||||||||||||||||||
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
3! |
|
|
|
|
|
|
|
|
|
|
||
два члена ряда) |
и заменим q t . |
|
Подставив в уравнение |
|||||||||||||||||||
(2.23), получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
3 |
|
|
|
|
j |
3 |
|
|
|
|
||||
|
|
|
|
|
|
|
|
tср |
0 |
|
|
|
|
|
|
|||||||
|
k |
|
|
|
|
j0 |
|
|
0 |
|
|
|
|
|
ktср j0 |
|
|
|||||
|
t |
|
|
|
|
k |
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
(2.24) |
||||||||
|
ср |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
3 |
0 |
|
|
k |
|
|
|
|
|
6 |
|
|
|
|
2 |
|
|
|||||
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
После преобразований имеем |
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
j |
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k |
tср |
0 |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
k |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
, |
|
|
|
|
|
(2.25) |
|||
6
откуда находим искомую зависимость
34
tср |
j0 |
3 |
6 |
. |
(2.26) |
k |
|
||||
|
|
k |
|
||
Формула (2.26) впервые была получена в работах кандидата технических наук В. П. Дмитриева и доктора технических наук Ю. Н. Бухарева и широко применяется в настоящее время для инженерной оценки времени срабатывания пороговых датчиков инерционного типа при их нагружении с различными темпами нарастания ускорения.
Следует подчеркнуть, что формула (2.26) будет давать точный результат при условии, что разница между sin z
|
z3 |
|
|
и z |
|
|
будет достаточно мала. При z ≤ 0,1644 разница |
|
|||
|
3! |
|
|
|
|
|
|
составляет не более 10–6 с. Поэтому применять формулу (2.26) для определения времени срабатывания датчика при действии линейно возрастающего ускорения с точностью порядка единиц микросекунд можно только при выполнении условия
|
6 |
0,1644 |
m |
(2.27) |
3 |
k |
c . |
Определим время tср срабатывания датчика №1, имеюще-
го параметры j0 = 15000 м/с2 (1500 g), m = 1 10–3 кг, с= 6560 H/м,
= 3 10–4 м, при действии линейно возрастающего ускорения j(t) = kt, где k = 108 м/с3 – темп нарастания ускорения, полученный путем решения уравнения (2.23) численным методом. В результате получим tср = 0,000414 c (414 мкс) как время от момента начала действия ускорения на датчик до момента замыкания его электрического контакта или до момента пе-
ремещения чувствительного датчика на величину рабочего хода = 3 10–4 м.
На рис. 2.2 показан график зависимости перемещения чувствительного элемента датчика от времени.
35
Рис. 2.2. Зависимость перемещения чувствительного элемента датчика от времени в соответствии с уравнением (2.22)
2.3. Построение динамической модели порогового датчика инерционного типа с упругим элементом отрицательной жесткости при силовом возмущении
Получим уравнение движения чувствительного элемента датчика, у которого силовая характеристика упругого элемента описывается уравнением F(q) = F0 – cq.
Составим дифференциальное уравнение движения чувствительного элемента датчика в виде
mq mj t F q . |
(2.28) |
Подставив в (2.28) уравнения j(t) = kt и F(q) = F0 – cq, получим уравнение
mq mkt F0 cq. |
(2.29) |
Разделив данное уравнение на m и преобразовав, получим уравнение движения в окончательном виде:
q 02q kt j0. |
(2.30) |
В уравнении (2.30) приняты следующие обозначения: j0 Fm0 ,
36
где j0 – порог по ускорению срабатывания датчика, и
02 mс ,
где 02 – параметр, характеризующий отношение коэффици-
ента жесткости упругого элемента к массе чувствительного элемента датчика, так как рассматриваемая динамическая система уже не является колебательной.
Уравнение (2.30) является неоднородным дифференциальным уравнением. Его решение представляет собой сумму общего решения соответствующего однородного дифференциального уравнения и частного решения:
qобщ qодн qчаст. |
(2.31) |
Однородное дифференциальное уравнение имеет вид |
|
q 02q 0. |
(2.32) |
Найдем решение линейного однородного уравнения (2.32). Будем искать частные решения также в виде q = ert. Так как
q rert |
и q r2ert , |
то, подставив соответствующие выражения в (2.32), получим
r2ert 2ert 0 |
(2.33) |
0 |
|
или, вынеся ert за скобки, |
|
ert r2 02 0. |
(2.34) |
Множитель ert в формуле (2.34) не обращается в нуль ни при каком значении t. Поэтому функция q = ert удовлетворяет линейному однородному дифференциальному уравнению (2.32) тогда и только тогда, когда число r является корнем характеристического уравнения
r2 2 |
0. |
(2.35) |
0 |
|
|
Корни характеристического уравнения (2.35)
37
r1 0 |
и r2 0 |
различные и действительные, поэтому паре действительных корней уравнения (2.32) соответствуют два частных решения
q e 0t |
и q e 0t , |
(2.36) |
1 |
2 |
|
а общее решение однородного уравнения (2.32) запишется в виде
q |
C e 0t C |
e 0t , |
(2.37) |
|
одн |
1 |
2 |
|
|
где С1 и С2 – произвольные постоянные, определяемые из начальных условий.
Теперь найдем частное решение qчаст. Исходя из вида правой части уравнения (2.30), частное решение будем искать также в виде
qчаст At B. |
(2.38) |
|
Так как |
|
|
qчаст A |
и qчаст 0, |
|
то, подставив соответствующие выражения в (2.30), получим уравнение
0 2 |
At B kt j |
(2.39) |
0 |
0. |
|
Раскрыв скобки и приравняв соответствующие части, найдем
A |
k2 , |
B |
j0 |
. |
2 |
||||
|
|
|
|
|
|
0 |
|
0 |
|
Тогда частное решение запишется в виде
q |
|
j0 kt |
. |
|
|
|
(2.40) |
||
|
|
|
|
||||||
част |
|
2 |
|
|
|
|
|||
|
|
|
0 |
|
|
|
|
|
|
Общее решение неоднородного |
дифференциального |
||||||||
уравнения (2.30) запишется в виде |
|
j0 kt |
|
|
|||||
q t C e 0t C |
e 0t |
|
. |
(2.41) |
|||||
|
|||||||||
1 |
2 |
|
|
|
2 |
|
|||
|
|
|
|
|
|
0 |
|
|
|
|
38 |
|
|
|
|
|
|
||
Для нахождения постоянных коэффициентов С1 и С2 продифференцируем уравнение (2.41), в результате чего получим уравнение
q t |
C e 0t |
C e 0t |
k |
. |
|
(2.42) |
|||
2 |
|||||||||
0 |
1 |
0 |
2 |
|
|
|
|
||
|
|
|
|
|
0 |
|
|
|
|
Введем начальные условия |
|
|
|
|
j0 |
|
|
||
q t q t 0 |
при |
t t0 , |
гдеt0 |
. |
|
||||
|
|
||||||||
|
|
|
|
|
|
k |
|
||
Данные начальные условия определяют, что движение инерционного элемента начинается в тот момент, когда выполнено условие j0 = kt0.
Из уравнений (2.41) и (2.42) с учетом начальных условий cоставим систему уравнений
|
0 |
C e 0t0 |
C |
e 0t0 |
j0 kt0 |
, |
|||||
|
|
|
|||||||||
|
1 |
|
2 |
|
|
2 |
|
|
|
||
|
|
|
|
|
|
|
0 |
(2.43) |
|||
|
|
|
|
|
|
|
|
||||
0 |
|
C e 0t0 |
C e 0t0 |
k |
. |
||||||
2 |
|||||||||||
|
|
0 |
1 |
|
0 |
2 |
|
|
|||
|
|
|
|
|
|
|
|
0 |
|||
В результате решения данной системы уравнений находим постоянные коэффициенты С1 и С2:
|
|
|
|
|
1 |
|
|
k |
|
|
|
|
j |
kt |
0 |
|
|
|
||||
C |
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
, |
(2.44) |
|||||
2e 0t0 |
|
|
|
|
|
2 |
|
|
||||||||||||||
|
1 |
|
|
3 |
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
0 |
|
|
|
|
0 |
|
|
|
|
|
||||
|
|
|
|
|
1 |
|
|
kt |
0 |
|
j |
|
k |
|
|
|
||||||
C |
2 |
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
. |
(2.45) |
||||
|
t |
|
|
|
2 |
|
|
3 |
|
|||||||||||||
|
|
2e |
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
0 0 |
|
|
|
|
0 |
|
|
0 |
|
|
|
||||||||
Подставляя С1 и С2 вуравнение (2.41), получим уравнение
39
|
q t |
|
|
1 |
|
|
|
|
k |
|
|
j |
|
kt |
0 |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
e 0t |
|
||||||||
|
|
|
|
|
|
3 |
|
|
|
|
|||||||||||||||
|
|
|
2e 0t0 |
|
|
|
2 |
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
1 |
|
kt |
0 |
j |
|
|
|
k |
|
|
|
|
|
j |
kt |
|
|||||||
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
e 0t |
0 |
|
. |
||||||
|
t |
|
|
|
|
2 |
|
3 |
|
2 |
|||||||||||||||
|
2e |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
0 0 |
|
|
|
|
0 |
|
|
|
|
|
0 |
|
|
|
|
|
0 |
|
||||||
Используя тождества [5] |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
ez e z |
sh z , |
|
|
|
ez e z |
сh z , |
|
||||||||||||||||||
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
||
преобразуем уравнение (2.46) к виду |
|
|
|
|
|
|
|||||||||||||||||||
|
|
q t |
|
1 |
k |
sh t t |
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
0 |
|
|
||||||||||||||||
|
|
|
|
|
|
02 |
|
|
|
|
|
|
|
0 |
|
0 |
|
|
|
||||||
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
kt0 j0 ch 0t 0t0 j0 kt .
(2.46)
(2.47)
Так как kt0 – j0 = 0, то получим окончательный вид урав-
нения (2.47)
q t |
1 |
|
k |
sh |
|
t |
j0 |
|
j |
kt |
. |
(2.48) |
|
2 |
|
|
|
|
|
||||||||
|
0 |
|
0 |
0 |
k |
0 |
|
|
|
||||
|
0 |
|
|
|
|
|
|
|
|
|
|
||
Решение уравнения (2.30) в аналитическом виде было получено также В. П. Гандуриным.
Чтобы определить время tср срабатывания датчика, в уравнении (2.48) необходимо вместо q(t) подставить – величину рабочего хода чувствительного элемента датчика, которая одновременно является величиной межконтактного зазора. Поэтому время срабатывания датчика с учетом времени на перемещение чувствительного элемента датчика на величину межконтактного зазора необходимо определять путем численного решения уравнения
|
1 |
|
k |
sh |
|
t |
|
|
j0 |
|
j |
kt |
. |
(2.49) |
2 |
|
|
|
ср |
|
|
||||||||
|
0 |
|
0 |
0 |
k |
0 |
|
ср |
|
|||||
|
0 |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
40 |
|
|
|
|
|
|
