УДК 671.1.052.2
DOI: 10.17238/issn0536-1036.2016.6.166
МОЩНОСТЬ ПРИВОДА НАСОСА ЗАГРУЗОЧНОЙ ЦИРКУЛЯЦИИ ВАРОЧНОГО КОТЛА, РАСХОДУЕМАЯ НА ДЕФОРМАЦИЮ СУСПЕНЗИИ В ПИТАТЕЛЕ
В.П. Сиваков, д-р техн. наук, проф.
А.П. Панычев, канд. техн. наук, доц.
И.А. Партин, ст. преп.
И.К. Кашапов, асп.
Уральский государственный лесотехнический университет, Сибирский тракт, д. 37, г. Екатеринбург, Свердловская область, Россия, 620100; e-mail: [email protected]
Исследовано влияние факторов режима загрузки на деформацию суспензии «щепа-щелок» в питателе загрузочной циркуляции варочного котла. Показано, что на деформацию суспензии при загрузке расходуется мощность, потребляемая насосом высокого давления. В ходе эксперимента решались следующие задачи: установить по экспериментальным данным зависимость деформации суспензии «щепа-щелок», выполняемой насосом высокого давления в каналах ротора питателя загрузочной циркуляции варочного котла; получить уравнения и численные значения масштабов подобия параметров для деформации суспензии в каналах ротора питателя и экспериментальной модели; определить мощность, расходуемую насосом высокого давления при деформации суспензии «щепа-щелок» в канале ротора питателя загрузочной циркуляции. Использованы следующие методы: полного факторного эксперимента; подобия и анализа размерностей. На основе методов подобия и анализа размерностей составлены уравнения для расчета и определены численные значения следующих масштабов подобия типового питателя и модели: геометрического, коэффициента Пуассона, сил деформации и трения, гидравлических радиусов, работы сил деформации, скорости, деформации и времени нагружения. Зависимость деформации суспензии при подаче из зоны низкого в зону высокого давления загрузочной циркуляции исследована на модели. Установлено, что относительная деформация возрастает пропорционально времени нагружения и давлению и обратно пропорционально жидкостному модулю. Максимальная мощность, расходуемая насосом на деформацию сырья в канале ротора питателя, определена по расчетным параметрам процесса деформации. Показано, что при расчетах привода насоса высокого давления необходимо учитывать расход мощности не только на транспортирование щепы в трубопроводе, но и на деформацию суспензии в питателе загрузочной циркуляции варочного котла.
Ключевые слова: мощность, ротор, питатель высокого давления, суспензия.
Для цитирования: Сиваков В.П., Панычев А.П., Партин И.А., Кашапов И.К. Мощность привода насоса загрузочной циркуляции варочного котла, расходуемая на деформацию суспензии в питателе // Лесн. журн. 2016. № 6. С. 166-174. (Изв. высш. учеб. заведений). DOI: 10.17238/issn0536-1036.2016.6.166
В известных работах [2, 3, 5-8, 10] при расчетах мощности приводов насосов высокого давления учитывалась только работа, затрачиваемая на транспортирование суспензии «щепа-щелок» загрузочной циркуляции установок непрерывной варки целлюлозы с вертикальными варочными котлами. Расход мощности на деформацию суспензии в питателе высокого давления (ПВД) не учитывался. Суспензия при перегрузке из зоны низкого (0,2 МПа) давления в зону высокого (1,2...1,4 МПа) давления нагнетательной линии загрузочной циркуляции рассматривалась как недеформируемая. Ограниченное количество исследований технологических и динамических процессов подачи суспензии из зоны низкого в зону высокого давления и пренебрежение расходом мощности на деформацию и изменение плотности суспензии проявились в перегрузке приводов насосов высокого давления. Для компенсации неточностей метода расчета мощности вводились поправочные коэффициенты. В частности Г.А. Тордуа [6] отмечает, что для трактов загрузочной циркуляции котлов установок непрерывной варки с учетом опыта их эксплуатации следует увеличивать давление от 1,1 до 1,6 МПа.
Установлено, что длина суспензии в канале ротора под действием избыточного давления при выгрузке уменьшается. Деформация (сжатие) суспензии происходит в режиме кратковременных (до 10 с) затухающих колебаний. Сжатие приводит к уменьшению длины суспензии на величину относительной деформации е.
Цель данной работы - определение мощности насоса, затрачиваемой на деформацию (сжатие) суспензии в ПВД.
Исследование деформации суспензии при сжатии проводили на модели, представленной на рис. 1, с применением масштабов подобия. Основные положения применения масштабов подобия сил и мощности при переходе от моделей к расчетам натуральных машин изложены в [1].
F
^6=
<
<
7ZZZZZ2Z2 //////.
П 1
г ill
Рис. 1. Физическая модель
Длина пробки суспензии в канале модели под действием внешней силы F уменьшается. Суспензию в канале модели рассматриваем как жесткую систему, деформация в которой мала по сравнению с размерами канала.
Исследование деформации суспензии в канале ротора выполняли в режиме трехфакторного эксперимента [2]. Относительную деформацию суспензии St в канале модели определяли при варьировании времени нагружения z\ = 10...20 с; статического давления z2 = 1,0...1,2 МПа; жидкостного модуля z3 = 4,5...7,7 дм3/кг. Пределы изменения этих факторов соответствовали режиму работы насоса в цикле выгрузки суспензии из ПВД в трубопровод загрузочной циркуляции [1].
Запишем уравнение линейной регрессии для исследуемого процесса:
st = 7,87 + 0,38xj + 0,63х2 - 0,87х3 . (1)
Относительная деформация суспензии возрастает пропорционально времени нагружения и давлению и обратно пропорционально жидкостному модулю. Для расчета численных значений масштабов подобия типового ротора ПВД с объемом каналов 0,45 м3 относительно модели сводим характеристики типового ротора и модели в таблицу.
Рабочие характеристики типового ротора ПВД и модели
Характеристика Значение характеристики
Ротор Модель
Длина канала l, м 0,480 -
Ширина канала X, м 0,280 -
Диаметр канала dm, м - 0,08
Высота канала 1к , lm , м 0,844 0,156
Время нагружения суспензии давлением t, tm, с 10...20 10.20
Постоянное давление на суспензию P, Pm, МПа 1,2 1,2
Технологическая щепа из древесины Ель Ель
Жидкостной модуль суспензии Д, Дт, дм3/кг 4,5...7,7 4,5...7,7
Схема типового ротора ПВД с объемом каналов ротора 0,45 м3 приведена на рис. 2.
Для выполнения условий подобия линейно-упругих свойств суспензии в модели и реальном канале ротора ПВД выполнен ряд условий [4, 6]. Суспензия в модели и в канале ротора ПВД имеет одинаковые безразмерные зависимости между напряжениями и деформациями, а также одинаковую схему нагружения.
Выразим все компоненты деформации s суспензии в точке xj [1]:
s = f (l, К, К X, F, E, и), (2)
где l, 1к, X - соответственно длина, высота и ширина поперечного сечения канала;
xj - произвольная координата суспензии;
F - внешняя сила;
E - модуль упругости суспензии;
и - коэффициент Пуассона суспензии.
Рис. 2. Схема типового ротора ПВД (объем каналов ротора 0,45 м3): I, III, V, VII - окна ротора с углом сектора 38о; II, IV, VI, VIII - участки корпуса ротора
между окнами
Уточним, что модели, спроектированные на основе подобия для линейно-упругих конструкций, допустимо использовать для исследования конструкций из нелинейного материала, включая пластическую область деформации [4]. Применив метод анализа размерностей, уравнение (2) запишем в виде
s =
F ,х. X. I. и..
-----ф(_1 —L _L —L)
E l ХФ( l ’ X ’ l ’ и ,
(3)
xt X i l и где ф(~т ,-Ь,—)
l X l и
функция п-переменных, включающая только безразмер-
ные члены.
Уравнение (3) выразим через масштабы переменных суспензии в канале ротора и модели:
F
подобия функций
s = -
E l X
I'
i=1
п-
(4)
где n - масштабы подобия;
k - число масштабов подобия.
Авторами работы [2] в качестве передаточных функций применены следующие масштабы подобия рассматриваемой задачи: геометрический (формы и координат), коэффициент Пуассона, сил и времени периода нагружения конструкции. Численные значения масштабов подобия получаем по расчетным формулам, определенным в [1], и экспериментальным данным исследования деформации щепы на модели [5].
Масштаб подобия площади поперечного сечения определяем по фор-
муле
n
a
4 l X
% d2
4•480•280 % 802
= 26,7,
(5)
где l, X, dm - линейные размеры поперечного сечения каналов модели и ротора приведены в таблице.
Масштаб коэффициента Пуассона при одинаковой технологической щепе принимаем одинаковым:
v
= 1.
(6)
Тогда графики зависимости относительных напряжений от деформации (а/Е от е) так же будут одинаковыми. Деформации в пластической области не являются однозначной функцией напряжения. Для выполнения условий подобия (Е = Ет) режим нагружения суспензии в модели должен соответствовать режиму нагружения суспензии в канале ротора. Отсюда следует, что
4 F F
% d2
__i_
X l
где Fm, F\ - соответственно силы при деформации суспензии в каналах моде-
ли и ротора,
77 % d2m
F = p------—
m
6 •lO3 = 1,2 •lO6 %-0^08
= 6 •lO3 Н:
F = F
4 X l % d2
= 6 •lO-
4 • 0,28 • 0,48 3,14 • 0,082
= 161,5 •103H;
l, X, dm - приведены в таблице.
Масштаб подобия сил определяем из соотношения
4 X l
% d
= n = 26,7 .
(7)
(8)
(9)
При деформации суспензии в каналах модели и ротора затрачивается сила F2 на трение по внутренним стенкам каналов. Влияние формы поперечного сечения на силу нагружения определяем из соотношения гидравлических радиусов.
Гидравлический радиус поперечного сечения канала ротора
R = ^ l 7ч = ^ 0,48'0,28 , = 0,09 м. (10)
2(X +1) 2(0,48 + 0,28)
Гидравлический радиус поперечного сечения модели
Кт =
% d2
d_
4
0,08
= 0,02 м.
(11)
% d 4 4
Масштаб подобия гидравлических радиусов пг и силу F2 находим по следующим формулам:
nv =
U
m
= ^ = 0,088 = 4 .
г R 0,02
гт 5
(12)
F = Fm n = 6-103 • 4,4 = 26,4-103Н .
Рассчитываем силу, расходуемую на деформацию и трение суспензии в канале ротора ПВД при p = 1,2• 106Па;dm = 0,08 м :
F = F + F = 161,5 •Ю3 + 26,4 •Ю3 = 187,9 •Ю3 Н, Определяем работу сил деформации сырья в канале ротора:
Q = F l — = 187,9 403 0,156— = 5,5 •103Дж . Q ’ ’ 5,4 , ^
Масштаб подобия высоты lm находим по формуле
4 0,844
п = — =-------
в 1т 0,156
= 5,4,
где 1к - средняя высота канала ротора, 1к = 0,844 ;
1т - средняя высота канала модели, lm = 0,156.
Максимальная скорость деформации щепы в канале модели:
= 8т = 0,016 = 0,0016 м/с,
m t 10 , ,
где гт =Мт lK = 0,1 • 0,156 = 0,016 м;
1mm - максимальное время нагружения (см. таблицу);
Alm - максимальная деформация суспензии по длине 1т модели,
Ы = 0,11
Максимальная скорость деформации щепы в канале ротора:
а = -
0,0844
= 0,0084 м/с,
t 10
min
где ек= 0,1 4 = 0,10,844 = 0,0844 м .
Масштаб подобия времени нагружения определяем по формуле
а 0,0084
а 0,0016
■ = 5,25 .
8
(13)
(14)
(15)
(16)
(17)
(18)
(19)
(20)
Находим максимальную мощность, расходуемую на деформацию сырья в канале ротора:
Q
N =^ п
t
mui
5,5-103 10
5,25 = 2,9 кВт .
(21)
Применение масштабов подобия при переходе от модели к типовому ротору ПВД позволяет определить мощность насоса, расходуемую на деформацию суспензии. Расход мощности на деформацию суспензии увеличивается из-за плотности щепы, загружаемой в варочный котел, что положительно сказывается на технологии варки [9].
СПИСОК ЛИТЕРАТУРЫ
1. Кабаяси А. Экспериментальная механика: пер. с англ. Т. 2. М.: Мир, 1990.
552 с.
2. Кашапов И.К., Партин И.А., Сиваков В.П. К вопросу о расчете мощности привода насоса высокого давления загрузочной циркуляции варочного котла // Перспективы развития техники и технологий в целлюлозно-бумажной промышленности: материалы II Всерос. отраслевой науч.-практ. конф. Пермь: ЦНИТИ, 2014. 297 с.
3. Малков С.Ю. Модернизация котлов непрерывной варки по технологии Lo-Solids™ // Целлюлоза, бумага и картон. 2004. № 7. С. 42 - 49.
4. Партин И.А., Сиваков В.П. Динамические процессы в питателе высокого давления при выгрузке // Лесн. журн. 2011. № 1. С. 117-125. (Изв. высш учеб. заведений).
5. Сиваков В.П., Партин И.А. Критерии подобия модели деформации суспензии в питателе высокого давления // Лесн. журн. 2013. № 4. С. 106-111. (Изв. высш. учеб. заведений).
6. Тордуа Г.А. Машины и аппараты целлюлозного производства // Машины и аппараты целлюлозного производства: учеб. пособие для вузов. М.: Лесн. пром-сть, 1986. С. 440.
7. Gullichsen J., Sundqvist H. On the Importance of Impregnation and Chip Dimensions on the Homogeneity of Kraft Pulping. Proc. 1995 Pulping Conf. Chicago, 1995, pp. 227-234.
8. Malkov S., Tikka P., Kuzmin V., Baltakhinov V. Efficiency of Chip Presteaming -Result of Heating and Air Escape Processes. Nordic Pulp and Paper Research Journal, 2002, vol. 17, no. 4, pp. 420-426.
9. Malkov S., Leavitt A., Stromberg B. Improved Understanding of Chip Steaming and Impregnation. Proc. Engineering, Pulping andPCE&I Conference. US, Atlanta, 2004.
10. Marcoccia B., Johanson J.R., Williams G., Bruce P. Proc. Pacific Paper’95. Canada, Vancouver B.C., 1995, pp. 8-14.
Поступила 17.04.16
UDC 671.1.052.2
DOI: 10.17238/issn0536-1036.2016.6.166
Driving Power of the Digester Load Circulation Pump for the Slurry Deformation in a Feeder
V.P. Sivakov, Doctor of Engineering Sciences, Professor
A.P. Panychev, Candidate of Engineering Sciences, Associate Professor
I.A. Partin, Senior Lecturer
I.K. Kashapov, Postgraduate Student
Ural State Forest Engineering University, Sibirskiy trakt, 37, Yekaterinburg, Sverdlovsk region, 620100, Russian Federation; e-mail: [email protected]
For citation: Sivakov V.P., Panychev A.P., Partin I.A., Kashapov I.K. Driving Power of the Digester Load Circulation Pump for the Slurry Deformation in a Feeder. Lesnoy zhurnal, 2016, no. 6, pp. 166-174. DOI: 10.17238/issn0536-1036.2016.6.166
This paper investigates the influence of loading mode factors on the “chips-liquor” slurry deformation in the load circulation feeder of a digester. A power consumed by the high pressure pump is also consumed for the slurry deformation at loading. The experiment objectives are: to establish a “chips-liquor” slurry deformation dependence, performed by the high-pressure pump in the feeder rotor channels of the digester load circulation; to obtain the equations and numerical values of the scales of similarity parameters for the slurry deformation in the feeder rotor channels and an experimental model; to determine the power consumed by a high-pressure pump at the “chips-liquor” slurry deformation in the rotor channel of the load circulation feeder. We use the methods of the full factorial experiment and of similarity and dimensional analysis. On the basis of the methods of similarity and dimensional analysis we compile the design equations and determine numerical values of the following similarity scales of a standard feeder and a model: geometric, Poisson's ratio, strain energy and friction force, hydraulic radius, work of strain energy, velocity, loading time and deformation. The dependence of the slurry deformation during supply from the low-pressure zone to the high-pressure zone of load circulation is examined on a model. The relative deformation increases in proportion to the load time and pressure and is inversely proportional to a liquid module. The maximum power, consumed by the pump to deform raw material in the feeder rotor channel, is determined by the design parameters of the deformation process. We should take into account the power consumption not only for the transportation of chips in the pipeline, but also for the slurry deformation in the load circulation feeder of a digester to calculate a high-pressure pump drive.
Keywords: power, rotor, high-pressure feeder, slurry.
REFERENCES
1. Kobayashi A.S. Handbook on Experimental Mechanics. US, New York, 1987.
1002 p.
2. Kashapov I.K., Partin I.A., Sivakov V.P. K voprosu o raschete moshchnosti privo-da nasosa vysokogo davleniya zagruzochnoy tsirkulyatsii varochnogo kotla [On the Driving Power Calculating of the High Pressure Pump of the Digester Load Circulation]. Perspek-tivy razvitiya tekhniki i tekhnologiy v tsellyulozno-bumazhnoy promyshlennosti: materialy II Vseros. otraslevoy nauch.-prakt. konf [Prospects of Development of Techniques and Technologies in the Pulp and Paper Industry: Proc. 2nd All-Russ. Sci. Prac. Ind. Conf.]. Perm, 2014. 297 p.
3. Malkov S.Yu. Modernizatsiya kotlov nepreryvnoy varki po tekhnologii Lo-SolidsTM [Modernization of Continuous Cookers by the Lo-SolidsTM Technology]. Tsel-lyuloza, bumaga i karton [Pulp. Paper. Board], 2004, no. 7, pp. 42-49.
4. Partin I.A., Sivakov V.P. Dinamicheskie protsessy v pitatele vysokogo davleniya pri vygruzke [Dynamic Processes in the High Pressure Feeder when Unloading]. Lesnoy zhurnal, 2011, no. 1, pp. 117-125.
5. Sivakov V.P., Partin I.A. Kriterii podobiya modeli deformatsii suspenzii v pitatele vysokogo davleniya [Criteria of Similarity of a Deformation Model of Suspension in a High Pressure Feeder]. Lesnoy zhurnal, 2013, no. 4, pp. 106-111.
6. Tordua G.A. Mashiny i apparaty tsellyuloznogo proizvodstva [Machines and Devices of the Pulp Production]. Moscow, 1986. 440 p.
7. Gullichsen J., Sundqvist H. On the Importance of Impregnation and Chip Dimensions on the Homogeneity of Kraft Pulping. Proc. 1995 Pulping Conf. Chicago, 1995, pp. 227-234.
8. Malkov S., Tikka P., Kuzmin V., Baltakhinov V. Efficiency of Chip Presteaming -Result of Heating and Air Escape Processes. Nordic Pulp and Paper Research Journal, 2002, vol. 17, no. 4, pp. 420-426.
9. Malkov S., Leavitt A., Stromberg B. Improved Understanding of Chip Steaming and Impregnation. Proc. Engineering, Pulping andPCE&I Conference. US, Atlanta, 2004.
10. Marcoccia B., Johanson J.R., Williams G., Bruce P. Proc. Pacific Paper’95. Canada, Vancouver B.C., 1995, pp. 8-14.
Received on April 17, 2016