Научная статья на тему 'Передаточная функция момента сопротивления центробежного расходомера сыпучих сельскохозяйственных материалов при осевой загрузке'

Передаточная функция момента сопротивления центробежного расходомера сыпучих сельскохозяйственных материалов при осевой загрузке Текст научной статьи по специальности «Механика и машиностроение»

CC BY
209
30
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ПЕРЕДАТОЧНАЯ ФУНКЦИЯ / ЦЕНТРОБЕЖНЫЙ РАСХОДОМЕР / ДАТЧИК / ТОЧНОСТЬ ИЗМЕРЕНИЯ / РАСХОД ЗЕРНА / ПЕРЕХОДНЫЙ ПРОЦЕСС / ДИФФЕРЕНЦИАЛЬНОЕ УРАВНЕНИЕ / ЯЧЕЙКА / СИЛА / TRANSFER FUNCTION / CENTRIFUGAL FLOWMETER / SENSOR / ACCURACY OF MEASUREMENT / GRAIN FLOW / TRANSITION / DIFFERENTIAL EQUATION / CELL / FORCE

Аннотация научной статьи по механике и машиностроению, автор научной работы — Багаев Андрей Алексеевич, Чернусь Роман Сергеевич

Целью работы является получение системы дифференциальных уравнений для поля центробежных сил и передаточной функции момента сопротивления центробежного расходомера сыпучих сельскохозяйственных материалов при осевой загрузке. Для анализа явлений происходящих в поле центробежных сил используется цилиндрическая система координат. Получены дифференциальные уравнения Эйлера для поля центробежных сил для трех координатных осей цилиндрической системы координат с учетом действия кориолисовых сил. Указанные дифференциальные уравнения позволили получить передаточную функцию момента сопротивления по соответствующим осям. При этом выходной переменной считается масса сыпучего материала, а переменной на выходе момент сопротивления на валу электродвигателя. В среде MatLab была составлена структурная схема для электромеханической системы «асинхронный двигатель крыльчатка» с учетом полученных передаточных функций момента сопротивления, используемая в дальнейшем для имитационного моделирования. Результаты моделирования представлены осциллограммами переходных процессов, анализ которых позволяет сделать следующий вывод: при неизменных параметрах двигателя и момента сопротивления показатели качества переходных процессов при измерении скорости и тока различны. При мгновенном приложении единичной нагрузки, скорость изменяется по апериодическому закону, а процесс изменения тока носит колебательный затухающий характер. Кроме того, время достижения установившегося значения переменной значительно меньше при регистрации тока. Следовательно, меньшей инерционностью и меньшей постоянной времени обладает процесс изменения тока, поэтому в качестве выходной измеряемой переменной целесообразно использовать именно ток статора.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по механике и машиностроению , автор научной работы — Багаев Андрей Алексеевич, Чернусь Роман Сергеевич

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Передаточная функция момента сопротивления центробежного расходомера сыпучих сельскохозяйственных материалов при осевой загрузке»

с 2 X мм

2.5 ~

1—1—[—1—I—1—I—1—I—1—I—т~мпн

Ю 20 за ¿О 50 ■ Рис. 3. Результаты испытаний на износ различных масляных композиций при нагрузке 5 Н: ° — чистое масло М-8В; Ал — масло М-8В + Деста; □ — масло М-8В + нанокобальт 0,2%

+ +

Библиографический список

1. Мышкин Н.К., Петроковец М.И. Трение, смазка, износ. Физические основы и технические приложения трибологии. - М.: ФИЗМАТЛИТ, 2007. - 368 с.

2. Выгодский М.Я. Справочник по элементарной математике. — М.: АСТ Аст-рель, 2006. — 509 с.

3. Сакович Г.В., Комаров В.Ф., Петров Е.А. Синтез, свойства, применение и производство наноразмерных синтетических алмазов. Часть 2. Применение и производство / / Сверхтвердые материалы. — 2002. — № 4. — С. 8-23.

+

УДК 621.9Т А.А. Багаев,

Р.С. Чернусь

ПЕРЕДАТОЧНАЯ ФУНКЦИЯ МОМЕНТА СОПРОТИВЛЕНИЯ ЦЕНТРОБЕЖНОГО РАСХОДОМЕРА СЫПУЧИХ СЕЛЬСКОХОЗЯЙСТВЕННЫХ МАТЕРИАЛОВ ПРИ ОСЕВОЙ ЗАГРУЗКЕ

Ключевые слова: передаточная функция, центробежный расходомер, датчик, точность измерения, расход зерна, переходный процесс, дифференциальное уравнение, ячейка, сила.

Перспективным средством измерения расхода сыпучих сельскохозяйственных материалов является центробежный расходомер при условии повышения его точности. Повышение точности невозможно без анализа динамических свойств системы «двигатель-крыльчатка». Динамические свойства двигателей достаточно хорошо изучены [1]. Вместе с тем основной трудностью является отсутствие сведений о динамическом поведении воспринимающего органа — крыльчатки.

Известно, что момент сопротивления центробежного расходомера зависит от конструктивных параметров крыльчатки и физических свойств сыпучего материала [2].

В результате реализации численного эксперимента было получено уравнение регрессии момента сопротивления в виде

полинома второго порядка [3]. Переменными факторами являются: угловая скорость крыльчатки о2, расход материала Q, угол изгиба лопасти турбинки относительно радиального направления ал, коэффициенты трения скольжения f1 и восстановления частиц е. Полученные выражения не позволяют исследовать динамические процессы в электромеханической системе «электрический двигатель — редуктор — крыльчатка расходомера», поскольку эти уравнения не являются дифференциальными, не позволяют получить передаточную функцию момента сопротивления и осуществить имитационное моделирование переходных процессов в центробежном расходомере. В связи с этим целью работы является получение системы дифференциальных уравнений для поля центробежных сил и передаточной функции момента сопротивления.

Для получения дифференциальных уравнений воспользуемся подходом, изложенным в работе [4]. При рассмотрении явлений, происходящих в поле цен-

тробежных сил, удобно использовать цилиндрическую систему координат.

От г

Рис. 1. Цилиндрическая система координат

На рисунке 1 изображен элемент объёма сыпучего материала, находящегося между лопастями крыльчатки и ограниченный тремя парами смежных координатных поверхностей. Выделенный элементарный объём ячейки сыпучего материала с ребрами dr, dz, ^ф имеет объём

V = Мф ■ dr ■ dz . При этом будем считать, что загрузка осуществляется по оси z.

На материал, находящийся в состоянии покоя относительно лопастей крыльчатки, действует ряд сил: F1 — сила, обусловленная действием давления на любую из граней выделенного объёма; F2 — центробежная сила; F3 — сила тяжести. При движении сыпучего материала в окружном или радиальном направлении появляются Кориолисовы силы инерции F4. Давление и, соответственно, сила есть функция всех трех координат

^ = /(г, ф^) = р ■ гйф ■ Ф ■ dz ■ а ,

где р — плотность; а — ускорение.

Тогда на основе общих принципов механики сумма проекций сил на соответствующие оси, действующих на движущийся элемент объёма сыпучего материала, равна произведению массы материала М на её ускорение.

На основании этого запишем:

F = F + F + F . r 1r 2 4p'

F = F + F . p lp 4r '

F = F + F z 1z 3 '

(1) (2) (3)

F F

где Г^г ' 12 , ^фф — проекции силы F1 на оси цилиндрической системы координат;

^ф , — проекции силы F4 на оси цилиндрической системы координат.

Или в дифференциальной форме:

р • rdp • dr • dz • a = - dp dr • rdp • dr • dz + r

pa'2r • rdp • dr • dz + 2av р • rdp • dr • dz ^^

P

p • rdp • dr • dz • a = -1 r • dp dp • P

; (5)

• rdp • dr • dz - 2av p • rdp • dr • dz

p • rdp • dr • dz • a = - dp/dz • rdp • z

, (6)

• dr • dz - pg • rdp • dr • dz

где аг , аф ,az — ускорения по соответствующим осям координат; ш — угловая скорость; vr ; vf — проекции линейной скорости на соответствующие оси цилиндрической системы координат.

С учётом M = pV = p • rdp • dr • dz уравнения (4)-(6) принимают вид:

M • a =-cpjdr • V + a2r • M + 2av ■ M . (7)

p

M • a =- l/r •dp/dp^ V - 2av • M ;

p

(8)

М ■ а = - ф/dz ■ V - рgV . (9)

Для получения уравнений баланса моментов домножим обе части уравнения (7) на г, уравнения (8) — на ^ф, уравнения (9) — на L=z. Тогда получим:

М ■ а ■ г = - с!р/дг ■ V ■

2 2 ;

• r + a r • M + 2av M • r p

M • a rdp = -1 r • dp dp • V • rdp p

- 2av • M • rdp

(10)

; (11)

М ■ а L = - др/dz ■ VL - рgLV . (12)

Рассмотрим уравнения (10)-( 12) отдельно. Момент сопротивления Мс, проекции линейной скорости Vг и vф, угловая скорость ш и ускорение аг могут быть определены следующим образом:

I dr dф

М = др дг ■ V ■ г V = — с =-

' г dt ' dt

cr

dp

v = r-= ra

p dt

r

r

dv

a = r

dt

r

d2r

dt

2

- ra

Тогда уравнение (10) после соответствующих преобразований примет вид:

M

cr .22 =4-a r -r

d2r

M

dt

2

(13)

Передаточная функция момента сопротивления по оси r равна:

^ M ^

W = 1

M V У

. 2 2 2 2 = 4-a r -r p .

(14)

В уравнении (11) введем следующие

обозначения a = r — M = dp/dp ■ V

p dt cp и запишем его в виде:

M

cy 2 da dr

-=-r -- 2ra —.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

M dt dt

(15)

Передаточная функция момента сопротивления по оси ф:

^ M ^

W = 2

M V У

2

= -3pa - r .

Ф

уравнении

М =ф/dz - VL , cL

(12)

d2L

a = z

dtz

образований получим:

ML

M

= -L-

d2L

dt

2

- gL

(16) обозначим и после пре-

(17)

Передаточная функция момента сопротивления по оси ф:

К 1 2 2

Ж = - =-р L - gL. (18)

3 М (18)

V / L

Для проведения имитационного моделирования в среде MatLab составлена структурная схема электромеханической системы «двигатель — крыльчатка» (рис. 2) с учетом полученных выше передаточных функций момента сопротивления (14), (16) и (18).

Параметры двигателя следующие: мощность асинхронного двигателя — 1,1 кВт, момент инерции суммарный — 0,026 кг*м, номинальная частота напряжения питания — 50 Гц, активное сопро-

тивление обмотки статора

R, = 0,118 Ом

, напряжение на обмотке

статора

U1 = 220В , угловая скорость

вращения поля статора — т0 = 157рад / с .

Результаты моделирования представлены на рисунке 3 а, б, где отчетливо определяются время пуска и момент приложения мгновенной единичной нагрузки.

Анализ осциллограмм переходных процессов позволяет сделать следующий вывод: при неизменных параметрах двигателя и момента сопротивления показатели качества переходных процессов при измерении скорости и тока различны. При мгновенном приложении единичной нагрузки скорость изменяется по апериодическому закону, а процесс изменения тока носит колебательный затухающий характер. Кроме того, время достижения установившегося значения переменной значительно меньше при регистрации тока. Следовательно, меньшей инерционностью и меньшей постоянной времени обладает процесс изменения тока, поэтому в качестве выходной измеряемой переменной целесообразно использовать именно ток статора.

Выводы

Получены дифференциальные уравнения момента сопротивления центробежного расходомера в цилиндрической системе координат и передаточные функции момента сопротивления по координатным осям. Синтезирована структурная схема электромеханической системы «двигатель — крыльчатка» и осуществлено имитационное моделирование. Установлено, что в качестве выходной измеряемой переменной целесообразно использовать ток статора.

Библиографический список

1. Соколовский Г.Г. Электроприводы переменного тока с частотным регулированием. — 2-е изд., испр. — М.: Издательский центр «Академия», 2007. — 272 с.

2. Багаев А.А., Лукьянов В.Г., Чер-нусь Р.С. Передаточная функция центробежного расходомера сыпучих сельскохозяйственных продуктов // Вестник Алтайского государственного аграрного университета. — 2010. — № 1 (63). — С. 71-75.

3. Багаев А.А., Чернусь Р.С. Уравнение регрессии момента сопротивления центробежного расходомера сыпучих

2

V

2

r

c

r

c

В

сельскохозяйственных материалов // 4. Соколов В.И. Центрифугирование. —

Вестник Алтайского государственного аг- М.: Химки, 1976. — 408 с. рарного университета. — 2010. — № 6

Slider

Рис. 2. Схема расходомера в среде Ма^аЬ

ш, 1/с 100

50

1 2 3 4 5 6 7 (, с

I, А

3 2 1

0 1 2 3 4 5 6 7

а

О), 1/с 100

50

° 1 2 3 4 5 6 7 1, С

I, А

3

2

1

О 1 2 3 4 5 6 7 I, С

б

Рис. 3. Осциллограммы скорости двигателя и тока статора при пуске двигателя и приложении единичной нагрузки: а — длина крыльчатки L=0,1 м; радиус крыльчатки R=0,1 м; б — L=0,2 м; R=0,2 м

i Надоели баннеры? Вы всегда можете отключить рекламу.