Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Динамика чувствительных элементов пороговых устройств автоматики. Учебное издание

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
825 Кб
Скачать

r1 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

 

 

 

 

 

 

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]