Научная статья на тему 'Математическое моделирование высокоскоростных МНЛЗ с электрогидравлическим приводом'

Математическое моделирование высокоскоростных МНЛЗ с электрогидравлическим приводом Текст научной статьи по специальности «Механика и машиностроение»

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

In the article is given the mathematical model of the table of swinging with hydroelectric drive of the crystallizer of the slugs continuous casting machine. The mathematical model is received by methods of the theory of automatic operating with outputting in dimensionless form of transfer functions of the system several elements.

Текст научной работы на тему «Математическое моделирование высокоскоростных МНЛЗ с электрогидравлическим приводом»

I /лггттгп гс ктьглрггл

II 1 (29). 2004-

In the article is given the mathematical model of the table of swinging with hydroelectric drive of the crystallizer of the slugs continuous casting machine. The mathematical model is received by methods of the theory of automatic operating with outputting in dimensionless form of transfer functions of the system several elements.

J

В. В. БОЛОТСКИЙ, ГГТУим. П. О. Сухого, Д. Н. АНДРИАНОВ, БИТУ

УДК 621.74.047

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ВЫСОКОСКОРОСТНЫХ МНЛЗ С ЭЛЕКТРОГИДРАВЛИЧЕСКИМ ПРИВОДОМ

В зависимости от скорости разливки сортовых заготовок существующие электромеханические привода качалки машин непрерывного литья заготовок (МНЛЗ) обеспечивают частоту колебаний от 0,8 до 4,2 Гц и амплитуду перемещений от 3 до 6 мм [1]. Одним из возможных путей повышения скорости разливки является повышение частоты колебаний и изменение амплитуды перемещений кристаллизатора МНЛЗ. Решение этой задачи может быть достигнуто применением электрогидравлического следящего привода (ЭГСП), использование которого позволит осуществить возможность реализации выходным звеном привода различных по профилю периодических сигналов и регулирования частоты и амплитуды колебаний в режиме реального времени, а также обеспечит повышение качества сортового металла и экономической эффективности использования МНЛЗ за счет значительного увеличения скорости разливки.

Режимы функционирования МНЛЗ предъявляют комплекс требований к динамическим характеристикам привода и в связи с этим, несомненно, актуальными являются исследования динамики привода, выполняемые на основе алгоритма, который сводится к построению амплитудной частотной характеристики по перемещению выходного звена привода Z Согласно теории автоматического управления, ЭГСП МНЛЗ можно представить в виде структурной схемы (см. рисунок). Выражения отдельных передаточных функций в виде изображений по Лапласу по этой схеме будут иметь вид [2-6]:

звено усилителя сигнала рассогласования:

щг{(3) = НМ = к19 (1)

и/ и^) 1

где 1/{ и иг - соответственно напряжение на входе и выходе усилителя сигнала рассогласования,

В; кх — коэффициент усиления усилителя сигнала рассогласования [1]; звено усиления по силе тяги ЭМП:

иъ{з) Т23 + \

где /г — тяговое усилие ЭМП без нагрузки, Н; иъ — результирующее напряжение с учетом его

потерь, В; к2 — коэффициент усиления по силе тяги ЭМП, Н/В; Т2 — электромагнитная постоянная времени обмотки управления, с;

звено преобразования перемещения золотника управляющего каскада в перемещение золотника распределительного каскада:

(3)

где к3 — коэффициент усиления по перемещению золотника распределительного каскада [1]; звено усиления гидравлического исполнительного механизма (ГИМ) по расходу:

ЯГ4Ь)=Ш = к4, (4)

Я*)

/;гггггГ; IÏ гстштшо!

-1 (29). 2004

и»

W„

и- и, ^-U2 л -.F F„ ^-^ h- Fa Pw_х ,--

—-ffl—w, |—н£)-н w2 h-ffl-H w, [H wz, l—ФН w,7 Ы w3

"f F- П^Мщ-

W18

w19

ufl

w9 w8

U,„

w„ W10

Структурная схема ЭГСГТ МНЛЗ

W5

W12

W1:

wt

/31

Wfi

где 20 ~~ расход жидкости через щель без нагрузки, м3/с; у - перемещение распределительного золотника, м; кА — коэффициент усиления по расходу, м2/с;

звено преобразования полезного расхода ()Р в перемещение выходного звена привода 2 •

QP(s) Sns

(5)

где Sn — площадь плунжера, м2; звено нагружения ГИМ:

W6(s) =

pm(s) ms +hTps + C1

(6)

где ртц — перепад давления нагрузки, Н/м2; тп - приведенная масса подвижных элементов на

выходном звене привода, кг; ^ - коэффициент демпфирования выходного звена привода,

Н с/м; — суммарная жесткость упругих элементов, Н/м;

звено преобразования перемещения заслонки в давление в диагонали гидравлического мостика ГУ без учета его потерь:

wM-

Pa(s) к7

h(s) 7]6s + l

где ра — давление в диагонали гидравлического мостика ГУ без учета его потерь, Н/м2; к7

коэффициент усиления по давлению, Н/м3.

звено датчика перемещения золотника распределительного каскада:

t/sW

(7)

fVg(s)=-

y(s)

' - 9

(8)

где С/8 — напряжение на выходе датчика перемещения, В; к8 — коэффициент усиления датчика, В/м;

звено усилителя сигнала датчика перемещения золотника распределительного каскада:

и9 (¿) к9

W9(s)=•

(9)

£/8(5) Г95 +1

где и9 - напряжение на выходе усилителя, В; к9 - коэффициент усиления усилителя сигнала

датчика перемещения [1]; Т9 — постоянная времени усилителя сигнала датчика перемещения, с; звено датчика перемещения выходного звена привода:

104 г{8)

■ - ^ю >

(Ю)

у

22 //^ГГгП ГГ ГСЩМ/^РГТО

I 1 (29). 2004 -

где 1/{0 — напряжение на выходе датчика перемещения, В; к10 — коэффициент усиления датчика, В/м;

звено усилителя сигнала датчика перемещения выходного звена привода:

где £/и — напряжение на выходе усилителя, В; кп — коэффициент усиления усилителя сигнала

датчика перемещения выходного звена привода [1]; Тп — постоянная времени усилителя, с; звено потерь расхода на сжимаемость жидкости:

= (12)

АцСО Сж

где Й2 ~ потери расхода, м3/с; Сж - жесткость жидкости, Н/м; звено утечек через неуплотненные элементы:

(13)

Ргц(з)

где - потери расхода, м3/с; кп - коэффициент утечек, м5/(Н с); звено потерь расхода в ГУ:

= (14)

Р ГЦ СО

где ()[4 - потери расхода, м3/с; кы - коэффициент внутренней проводимости, м3/(Н-с); звено потерь в ЭМП от влияния скоростной э.д.с.:

= (15)

где и15 - потери напряжения в ЭМП, В; - коэффициент скоростной э.д.с., В/м; звено потерь давления в ГУ:

= (16)

где рос - потери давления в ГУ, Н/м2; х - перемещение золотника управляющего каскада, м; к16

— коэффициент усиления по давлению звена обратной связи, Н/м3; Г16 — постоянная времени, с;

звено преобразования давления в диагонали гидравлического мостика ГУ в перемещение золотника управляющего каскада:

^Ь^-Л1^ (17)

где /?ду - давление в диагонали гидравлического мостика ГУ, Н/м; кХ1 - коэффициент усиления

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

по перемещению золотника управляющего каскада, м3/Н; Т{7, Т[7 — постоянные времени, с; ¡¡17 —

коэффициент демпфирования [1];

звено потерь тягового усилия ЭМП от действия гидродинамических сил на заслонку:

= (18)

где ^ - потери тягового усилия ЭМП от действия гидродинамических сил на заслонку, Н; -

коэффициент усиления звена обратной связи, м2;

звено потерь тягового усилия ЭМП от действия упругой механической обратной связи на заслонку:

К19(*)=^ = к19, (19)

■Ф')

/;птггг г: г^шглтгп? / оо

--1 (29). 2004/ иб

где ~ потери тягового усилия ЭМП от действия упругой механической обратной связи на заслонку, Н; к[Ч - коэффициент усиления звена обратной связи, Н/м;

звено преобразования результирующей силы тяги в перемещение заслонки:

h(s) _ к20

Fh(s) (r20V+2i20^ + l)' (20)

где Fh - результирующая сила тяги, действующая на заслонку, Н; h ~ перемещение заслонки, м;

к20 - коэффициент усиления по перемещению заслонки, м/Н; Г20 — постоянная времени, с; ¡¡20 - коэффициент демпфирования [1].

На основании выражений (1)—(20) и теоремы эквивалентного преобразования структурных схем запишем выражение общей передаточной функции с выходным сигналом по перемещению выходного звена привода:

W L =_Wmr(s)W4(s)W5(s)_

°6UlW £/„(*) 1 + ^ + + ^14 + ^пэг (^>^4 (^)^IO ! '

где

W (J) = W20 № № fr frfr1 (^3 № + W16(S)W[7(S)] Lx (s) + W20 (s)WX72(s)W2 (s)WX5 (s) + L2 (s) ;

Lx (s) = [1 + Wx 6{s)VX7 {s)+ W20 (s)W7 (s)(fVX8 (s) + WX9(s)WX7(s))][^WX6{s)VX7{s)\;

L20) = W2o(s)W7(s)W2(s)WX7(s)Wx(s)JV3(s)W8(s)W90)[1 + WXb(syVX7(s)] . Используя приведенную математическую модель, можно провести аналитические исследования динамических режимов нагружения механизма качания MHJI3 с целью обоснования оптимальных рабочих режимов для получения необходимого качества в зависимости от марки разливаемого сортового металла.

Литература

1. Кель Н., Л он ар д и Э., Ш ер ген П. Гильзовый кристаллизатор высокоскоростных MHJI3 со встроенным механизмом гидравлического качания Вибромолд // Сталь. 1999, №4. С. 23-26

2. Гамынин Н.С., Каменир Я.А., Коробочкин Б.Л. и др. Гидравлический следящий привод. М.: Машиностроение, 1968.

3. Гамынин Н.С. Гидравлический привод систем управления. М.: Машиностроение, 1972.

4. Лещенко В.А. Гидравлические следящие приводы станков с программным управлением. М.: Машиностроение, 1975.

5. Попов Д.Н. Динамика и регулирования гидро- и пневмосистем. М.: Машиностроение, 1976.

6. Шагинян A.C. Динамика сейсмических вибраторов с электрогидравлическим сервоприводом // Исследование Земли невзрывными сейсмическими источниками. М.: Наука, 1981.

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