Scientific journal
Modern high technologies
ISSN 1812-7320
"Перечень" ВАК
ИФ РИНЦ = 0,899

FLAP VALVE MODEL DEVELOPMENT FOR DAMPING OPTIMIZATION AT CLOSING

Chubarov F.L. 1 Akimenko D.A. 1 Sizov A.N. 2 Nikitin A.V. 3
1 Federal State Budgetary Institution of Higher Education Kaluga branch Moscow State Technical university named after Bauman (National Research University)
2 OJSC «Kaluga turbine plant»
3 Public Autonomous Vocational Educational Institution Kaluga region Kaluga Technical College
A mathematical model of a quick-closing flap valve was built and researched for the heat-extraction selection of a steam turbine operating in the forced-shutdown mode. The rotary-type flap valve is a steam check valve equipped with a hydraulic damper and actuator that is unloaded in working condition. The design, main details and operating principle of the Dy400 slam is shown. When triggered, it performs a protective function for the turbine, prevents the reverse flow of steam from the selections into the turbine when its emergency protection is triggered. The model is applied to get some recommendations for the design and adjustment of the hydraulic damper in the composition of the flap, based on the condition of the fastest and at the same time possible unstressed fit on the saddle, in order to avoid its early destruction. Setting the closing speed is possible due to the fact that the design of the water and damper has the means to adjust. A mathematical model of a valve-collapse with a damper has been developed, its operation has been simulated for various initial vapor pressure drops. The model parameters are determined both by real construction and by the results of parametric identification of its field tests. Also, when creating the model, a number of assumptions were made. A comparison of the results of modeling on the angle of opening of the valve-collapse with the results obtained experimentally on the stand, the optimal coefficient of viscous friction was determined. Scheduled further let research.
check valve
flap valve
steam extraction
mathematical model
hydraulic damper
closing speed
unstressed fit
dynamic study
load
matching criterion
restrictions
parametric identification

Стационарные паровые турбины типа «ПТ», как правило, имеют два регулируемых отбора пара: производственный и теплофикационный. Для штатного подключения и отключения линий отборов на паропроводы устанавливаются электроприводные задвижки, скорость перемещения которых ограничена и составляет обычно несколько десятков секунд. Однако во время работы турбины возможно возникновение ряда аварийных ситуаций, при которой вместе с остановом самой турбины требуется быстро закрыть и её отборы для предотвращения обратного потока пара в турбину из коллектора питаемой тепловой сети. Для решений этой задачи в составе паровых турбин ОАО «КТЗ» [1] зачастую используются обратные клапаны поворотного типа c условными диаметрами Dy от 150 до 500 мм. Такой клапан представляет собой поворотную заслонку, подвешенную в трубопроводе отбора на горизонтальной оси. В нерабочем состоянии захлопка закрыта и плотно прижата к своему седлу пружиной в составе привода. При необходимости подачи пара в отбор на привод заслонки подаётся команда на взвод, в результате которой действие пружины преодолевается гидравликой и заслонка становится разгруженной, благодаря чему и открывается потоком пара (для начала открытия, как правило, достаточно перепада давления всего ~15 кПа). При возникновении обратного потока пара этот же перепад давления и закрывает заслонку, действуя вместе с собственной силой тяжести заслонки. Таким образом, реализована быстродействующая защита. Проблемой при этом является то, что под воздействием больших перепадов давления в паропроводе с разных сторон захлопки происходит слишком быстрое ее закрытие, и в совокупности с большими габаритами и массой этой захлопки возникает сильный удар о седло, что ведет к его разрушению и заклиниванию самой захлопки в седле. Во избежание подобных дефектов применяются ряд конструктивных приспособлений [1, 2], таких как:

– ограничительные упоры;

– демпфер пружинного типа;

– демпфер гидравлического типа.

Целью настоящей работы является исследование и выработка рекомендаций по конструкции и настройке последнего типа, как наиболее функционального и эффективного. Такой демпфер способен с переменной скоростью замедлять закрытие захлопки при подходе к седлу и тем самым уменьшать силу удара о него. Одним из возможных путей решения данной проблемы является проведение динамических исследований демпфера при различных условиях работы захлопки, некоторые результаты по которым имеются [3, 4]. Однако целесообразным представляется и динамическое исследование с использованием адекватной математической модели объекта, что значительно упрощает сам процесс исследования, сокращает временные и материальные затраты на проведение самих исследований и снимает риски аварийных ситуаций. Для проведения такого исследования необходима разработка математической модели реальной захлопки Dy400 Py80 (рис. 1) с гидравлическим демпфером, представленная ниже в первой части работы.

chubar1a.tif

Рис. 1. Конструкция обратного клапана-захлопки Dy400 Pу80 с гидравлическим демпфером, Основные детали: 1 – фланец нажимной, 2 – сальник, 3 – втулка, 4 – крышка, 5 – корпус, 8 – клапан, 10 – рычаг клапана, 12 – рычаг, 13 – рычаг привода клапана, 15 – шток привода, 17 – цилиндр, 18 – пружина, 19 – поршень, 20 – седло ускорителя, 21 – поршень ускорителя, 23 – пружина ускорителя, 25 – корпус ускорителя, 26 – установка концевых выключателей, 28 – седло

Построение адекватной математической модели является сложной и трудоёмкой задачей, поскольку энергетические объекты, в частности турбины, относятся к классу существенно нелинейных объектов, и при построении модели достаточно трудно учесть все внешние и внутренние факторы, влияющие на работу рассматриваемого объекта [5]. Даже если и удаётся построить модель, учитывающую наиболее важные факторы, остается проблема с определением числовых значений параметров модели, поскольку большинство параметров имеет справочные значения, определенные при стандартных условиях, что существенно отличается от условий эксплуатации [6]. Для решения данной проблемы можно воспользоваться методами параметрической идентификации, которые позволяют уточнить числовые значения по результатам натурных испытаний [7].

Ниже рассматривается решение задачи построения математической модели клапана-захлопки в режиме форсированного закрытия на основе методов исследования динамики и моделей гидро- и пневмоприводов [8].

На рис. 1 представлена конструкция обратного клапана-захлопки Ду 400, а на рис. 2 представлены экспериментальные данные по скорости ее закрытия в зависимости от перепада давления пара на ней.

chubar2.tif

Рис. 2. Экспериментальный график зависимости степени закрытия захлопки от времени при различных перепадах давления пара

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

Построение математической модели

При построении математической модели будем учитывать наличие гидравлического демпфера, который вступает в работу на конечном этапе закрытия клапана-захлопки. Кроме того, необходимо учитывать процесс истечения рабочей жидкости из нижней полости гидравлического цилиндра в дренажную линию [8], что также может усложнить модель рассматриваемого объекта. Однако, исходя из полученных экспериментальных данных, процесс закрытия протекает достаточно медленно и, следовательно, динамикой изменения давления в демпфере и в нижней полости гидроцилиндра можно пренебречь. Также будем пренебрегать аэродинамическими процессами обтекания паром захлопки ввиду их квазистационарного характера при малых углах открытия. Следовательно, для построения динамической модели клапана в первом приближении достаточно рассмотреть баланс сил, действующих на поршень гидравлического цилиндра привода.

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

chubar3.tif

Рис. 3. Схема определения нагрузки со стороны захлопки

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

С учётом этого и расчётной схемы, представленной на рис. 3, можно записать следующее уравнение:

chubar01.wmf

chubar02.wmf

Здесь α(t) – текущий угол открытия захлопки (причём chubar03.wmf); Δpпара – перепад давления пара на захлопке; Dз – диаметр захлопки; lз – расстояние от центра клапана до оси вращения; mз – его масса; g – ускорение свободного падения; Cпр – жесткость возвратной пружины; хпр – начальное поджатие возвратной пружины; rс – радиус сектора реечной передачи; mп – масса подвижных частей гидроцилиндра; Kвт и Kст – соответственно коэффициенты вязкого и сухого трения подвижных частей гидроцилиндра; Jз – момент инерции заслонки относительно оси вращения; Kзвт – момент коэффициента вязкого трения в оси заслонки.

Поскольку динамикой процессов, протекающих в демпфере и при истечении жидкости из нижней полости гидроцилиндра, пренебрегаем, то демпфер можно рассматривать как дополнительное сопротивление, и, следовательно, коэффициент Kвт можно рассматривать как «составной», учитывающий вязкое трение, сопротивление истечение жидкости из отверстия и демпфер, который «включается» на конечном участке закрытия клапана; т.е. Kвт = Kзвт и можно определить как: Kзвт = Kвт, если chubar04.wmf, и chubar05.wmf, если chubar06.wmf, где lдемп – длина хода демпфера, и chubar07.wmf значительно больше Kвт.

Пренебрегая сухим трением, окончательно получим

chubar08.wmf

chubar09.wmf,

или

chubar10.wmf

chubar11.wmf (1)

Последнему уравнению будет соответствовать структурная схема, представленная на рис. 4.

chubar4.tif

Рис. 4. Структурная схема математической модели захлопки

Кроме того, на значение угла накладываются ограничения chubar12.wmf и начальными условиями для решения последнего уравнения являются chubar13.wmf и chubar14.wmf.

Параметрическая идентификация

В полученной математической модели клапана-захлопки большинство параметров являются конструктивными, значения которых определяются геометрией и характеристиками материала, из которых они изготовлены. Уточнению подлежит параметр Kвт = Kвт(α), поскольку в соответствии с допущениями данный параметр включает коэффициент вязкого трения, зависящего от условий эксплуатации и характеристик используемой рабочей жидкости, и учитывает сопротивление, возникающее при истечении рабочей жидкости из рабочей полости гидроцилиндра в дренажную линию и сопротивление со стороны демпфера. Поэтому числовое значение данного параметра необходимо уточнить по результатам эксперимента.

В качестве критерия совпадения угла открытия захлопки, полученного при решении дифференциального уравнения модели, с экспериментальными данными, воспользуемся критерием [3, 7]

chubar15.wmf (2)

Здесь αэкс(t) – экспериментально полученная зависимость угла открытия захлопки от времени; α(t, Kвт) – решение уравнения (1) модели захлопки; T – интервал исследования.

Задача уточнения числового значения коэффициента Kвт решается в два этапа. На первом этапе подбирается значение, соответствующее участку траектории, на котором не задействован демпфер. Минимизируя функционал (2) при начальных условиях chubar16.wmf и chubar17.wmf решения уравнения (1), находим значение коэффициента Kвт, которое будет учитывать коэффициент вязкого трения и сопротивление истечению жидкости из полости гидроцилиндра в дренажную линию. В результате решения для данного участка получено следующее оптимальное значение: K*вт = 13605. На этом же этапе определяется скорость закрытия захлопки и значение угла на момент времени, соответствующий началу демпфирования, которые будут являться начальными значениями для нахождения решения уравнения (1) при минимизации функционала (2) на втором участке. Для данного участка K*вт = 2102062. Таким образом,

chubar18.wmf (3)

На рис. 5 представлены результаты моделирования для (3) при различных перепадах давления пара на захлопке, из которых видно, что построенная математическая модель для сделанных допущениях с приемлемой точностью описывает протекающие процессы при закрытии клапана-захлопки.

chubar5.tif

Рис. 5. Графики изменения угла открытия захлопки от времени

Заключение

Построенная математическая модель с учётом сделанных допущений является упрощенной и поэтому может быть использована для предварительной оценки времени закрытия клапана-захлопки, и тем самым можно оценить соответствие рассматриваемого объекта предъявляемым в техническом задании требованиям. В дальнейшем планируется повысить адекватность построенной модели путем учёта дополнительных факторов, в частности моделей демпфера и истечения жидкости из полости гидроцилиндра в дренажную линию, с последующим уточнением числовых значений параметров модели по экспериментальным данным. Конечной целью является исследование влияния и оптимизация размеров отверстия и поршня демпфера с точки зрения обеспечения наискорейшего и безударного закрытия (посадки на седло). Оно будет проведено на последующих этапах работы.