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

кандидата технических наук
Соколов, Евгений Вячеславович
город
Екатеринбург
год
2008
специальность ВАК РФ
05.21.03
цена
450 рублей
Диссертация по технологии, машинам и оборудованию лесозаготовок, лесного хозяйства, деревопереработки и химической переработки биомассы дерева на тему «Моделирование и исследование динамических и гидродинамических процессов в центробежных насосах массоподводящих систем бумагоделательных машин»

Автореферат диссертации по теме "Моделирование и исследование динамических и гидродинамических процессов в центробежных насосах массоподводящих систем бумагоделательных машин"

На правах рукописи

СОКОЛОВ ЕВГЕНИЙ ВЯЧЕСЛАВОВИЧ

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

05 21 03 - технология и оборудование химической переработки биомассы дерева, химия древесины

Автореферат

диссертации на соискание ученой степени кандидата технических наук

ООЗ167952

Санкг-Петербург-2008

003167952

Работа выполнена на кафедре сопротивления материалов и теоретической механики Уральского государственного лесотехнического университета

Научный руководитель:

Официальные оппоненты:

доктор технических наук, профессор Анкудинов Дмитрий Тимофеевич

доктор технических наук, профессор Кокушин Николай Николаевич

кандидат технических наук, доцент Куцубина Нелли Валерьевна

Ведущая организация: Пермский государственный технический

университет

Защита состоится « 2(Эу> /и а Я 2008 года ъ^ часов на заседании диссертационного совета Д 212.231.01 при ГОУ ВПО Санкт-Петербургский государственный технологический университет растительных полимеров по адресу 198095, Санкт-Петербург, ул Ивана Черных, д 4

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

Ваши отзывы на автореферат в двух экземплярах, заверенные печатью организации, просим направлять по адресу 198095, Санкт-Петербург, ул Ивана Черных, д 4 СПбГТУРП, ученому секретарю диссертационного совета

Автореферат разослан «

Ученый секретарь диссертационного совета

Швецов ЮН

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

ЦН специального назначения являются специализированным химическим оборудованием ЦБП, при этом большая их часть входит в состав поточных линий производства бумаги, картона, целлюлозы, древесной массы Непрерывная и безаварийная работа ЦН влияет на производительность и эффективность работы БМ. Работа ЦН в нестационарном режиме в МПС влияет на качество процесса формирования капиллярно-пористой структуры бумажного полотна. Продолжительная работа в нестационарном режиме снижает долговечность и приводит к аварийным отказам с разрушением основных деталей ЦН При эксплуатации ЦН с помощью систем мониторинга и диагностики важно обеспечить контроль рабочих процессов ЦН, поэтому нужны соответствующие научно-обоснованные рекомендации

Наиболее актуальными в таком случае представляются исследования динамических и гидродинамических процессов в ЦН Данные исследования позволят разработать методы анализа динамических характеристик ротора, получить картину течения жидкости в проточной части насосов, определить параметры напряженно-деформированного состояния (НДС) основных деталей насосов и установить причины их разрушения в эксплуатации. Применение разработанных и экспериментально проверенных методик расчета на практике позволяет обоснованно вносить изменения в конструкцию ЦН для повышения эффективности эксплуатации в массоподводящих системах БМ и других подобных сложных гидравлических системах.

Целью работы является повышение эффективности работы массоподводящих систем бумагоделательных машин путем решения следующих задач.

- моделирования динамических процессов нагружения деталей в насосах,

- выявления собственных форм и частот колебаний ротора и определения амплитудно-частотных характеристик,

- моделирования гидродинамических процессов в насосах и определение параметров потока жидкости,

- определения параметров напряженно-деформированного состояния ротора и спирального отвода насоса,

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

Методы исследования. Теоретические исследования проведены с привлечением основных положений динамики машин, теории колебаний, методов конечных элементов и контрольных объемов, гидродинамики и теории лопастных машин Экспериментальные исследования выполнены на испытательном стенде завода ОАО «Уралгидромаш», на месте эксплуатации ЦН в массоподводящей системе БМ №14 ОАО «МБП Сыктывкарский ЛПК» При экспериментальных исследованиях использован анализатор вибрации «Топаз» фирмы «ДИАМЕХ», а также лицензионное программное обеспечение Геометрические построения выполнены в автоматизированной системе трехмерного моделирования Unigraphics Численные расчеты проведены в пакетах CAE и CFD анализа ANSYS, CFX

Предмет исследования. Динамические и гидродинамические процессы в центробежных насосах

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

Достоверность научных результатов исследования обуславливается использованием при моделировании динамических и гидродинамических процессов в ЦН развитых теорий колебаний, методов конечных и контрольных объемов, известных программных средств При этом значения резонансных частот ротора по формам крутильных и изгибных колебаний экспериментально подтверждены путем идентификации гармоник в частотном спектре вибрации Интегральные параметры потока на выходе и входе конечно-объемной модели подтверждены экспериментально путем сравнения с параметрами центробежного насоса полученных на испытательном стенде ОАО «Уралгидромаш» При этом погрешность полного давления с экспериментом составила для режима оптимальной подачи 1,2 %, нулевой подачи - 2,2 % и для режима максимальной подачи насоса -1,6% Апробация расчета проведена на серийно выпускаемом для ЦБП центробежном насосе консольного типа

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

Результаты работы следует использовать при проектировании ЦН на машиностроительных заводах, выпускающих насосы специального назначения для ЦБП, а также при модернизации насосов на предприятиях ЦБП. Использование разработанных методик расчета динамики и гидродинамики ЦН подтверждены актом внедрения работы на ОАО «Уралщдромаш» и ОАО «МБП Сыктывкарский ЛПК» Результаты исследования используются в учебном процессе для подготовки инженеров-механиков по специальности 150504 в курсе «Теория и конструкция машин и оборудования отрасли», а также в разделе диагностика потокообра-зующих и потокопроводящих систем курса «Диагностика машин и оборудования» Апробация работы. Основные положения диссертационной работы доложены на научно-технической конференции студентов и аспирантов Уральского государственного лесотехнического университета в апреле 2004г (УГЛТУ, Екатеринбург), всероссийской научно-технической конференции, посвященной 90-летию Уральской государственной горно-геологической академии «Математическое моделирование механических явлений» 28-29 мая 2004г (УГГТА, Екатеринбург), П всероссийской научно-технической конференции студентов и аспирантов 17-26 апреля 2006г (УГЛТУ, Екатеринбург), евразийском симпозиуме «Деревообработка технологии, оборудование, менеджмент XXI века» 20-21 сентября 2006г (УГЛТУ, Екатеринбург)

Публикации. По теме диссертации опубликовано 9 печатных работ Работа выполнялась в рамках госбюджетной темы по единому наряд-заказу Минобразования РФ

Научные положения выносимые на защиту:

- результаты моделирования собственных и вынужденных колебаний ротора с учетом упругих свойств опор, под воздействием центробежных сил инерции неуравновешенных масс ротора, присоединенной массы волокнистой суспензии и гироскопического момента;

- результаты моделирования трехмерного течения жидкости в проточной часш ЦН,

- методика определения параметров напряженно-деформированного состояния ротора и спирального отвода ЦН,

- результаты расчета колебаний спирального отвода на лопастной частоте Объём и структура диссертационной работы Диссертационная

работа состоит из введения, четырех разделов, заключения, списка литературы и двух приложений Она включает 176 страниц машинописного текста, 95 ил , 35 табл , 134 наименования использованных источников, в том числе 18 иностранных

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ Во введении обоснована актуальность темы диссертационной работы, определены ее цель и задачи, положения выносимые на защиту, показаны научная новизна и практическая значимость работы

1. Анализ состояния вопроса и постановка задачи В разделе приведены общие сведения о массоподводящих системах назначение и принцип их действия Исследованы массоподводящие системы БМ №2, №3 ОАО «Соликамскбумпром», БМ №14, №15 ОАО «МБП Сыктывкарский ЛПК». Описан процесс массоподвода на примере БМ №14 ОАО «МБП Сыктывкарский ЛПК» Рассмотрена номенклатура и назначение эксплуатируемых в МПС насосов консольного типа БМ, АХ, X, ХРО, АРР, ЕРА Сделан анализ конструкций, выявлены отличия и обоснован тип исследуемого ЦН в качестве примера для построения динамической модели ЦН консольного типа. Рассмотрен теоретический процесс работы и основные характеристики ЦН Показано, каким образом осуществляют выбор насоса для гидравлической системы Описана методика гидравлического расчета массопроводов, а также проанализированы методы проектирования гидравлических систем на предприятиях ЦБП Сделан вывод, что основой расчета не всегда являются инженерные методы, поэтому выбор ЦН в таких случаях не оптимальный, что приводит к неоправданным затратам Это подтверждается опытом эксплуатации насосов на многих предприятиях ЦБП, когда по истечению некоторого времени, требуется вносить изменения в режим работы ЦН, то есть корректировать его подачу и напор

Показано, что насосы в МПС работают не в оптимальном режиме, так как при их выборе необходимо учитывать реологические свойства водоволокни-стой суспензии и возможность изменения гидравлической характеристики МПС При этом отмечено, что вероятность появления нестационарного режима работы ЦН в частности кавитации и помпажа увеличивается, в том случае, когда ЦН работает вне области рабочей части напорной характеристики. Таким образом, в нестационарном режиме работы снижается долговечность основных деталей насосов и происходит их разрушение. Анализ эксплуатационной информации на предприятиях ЦБП позволил установить, что свыше 50 % отказов ЦН происходит из-за узлов уплотнения и подшипников При этом наиболее опасные отказы сопровождаются разрушением вала, элементов РК и спирального отвода, что приводит к аварийному останову ЦН и значительному экономическому ущербу

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

Проведен аналитический обзор работ российских и зарубежных ученых по исследованию динамических и гидродинамических процессов в обо-

рудовании МПС, что позволило выявить актуальные проблемы, решение которых направлено на снижение числа аварийных отказов насосов и повышение эффективности работы МПС.

2. Моделирование и исследование динамических процессов в центробежных насосах В разделе рассмотрены основные методы исследования колебания ротора ЦН как механической системы. Построена принципиальная схема ротора ЦН типа ХРОЮОО/34 (рис. 1). Ротор представлен в виде балки переменной жесткости на упругих опорах с массивным диском на консоли.

На рис. 1 обозначено: Е, р тЛ„ - модуль упругости и

плотность материала вала;

1), - площадь и момент инерции сечения вала; ш, .1га - масса и приведенный момент инерции РК; С/«, Сду, СВх, Сву, Св2 И Ьдх, ЬАу, Ьв« ЬВу, ЬВ2 - коэффициен-Рис. 1. Принципиальная схема ротора ЦН ты жесткости и неупругих

сопротивлений опоры А и В вдоль осей X, У, Ъ\ 1к - расстояние между опорой А и центром масс РК, 1п -расстояние между опорами А и В; 1„ - расстояние между опорой В и концом вала. Начало системы координат О - точка пересечения средней плоскости РК с осью вала ОХ.

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

[мщ}+[вт+[ст={т), о)

где [М], [В], [С] - матрицы масс (инерции), демпфирования и жесткости

системы; {(?}, {д}, {д} -

2 обобщенные перемещения

МВВштшШ узлов и их производные;

// /\ {Р(*)} - обобщенные силы.

МЙЙЁйь. Л//, \ Построение матема-

Иг'^'^У ж *** \\ тической модели ротора

^ЩдВЙ^ 7ч ^^ ^ сводиться к представлению

у 1 принципиальной схемы ро-

^^ V тора (рис. 1) в виде конечно-

элементной модели (рис. 2). При этом выбраны типы конечных элементов (КЭ), для моделирования упругих

Рис. 2. Конечно-элементная модель ротора ЦН: 1 - рабочее колесо; 2 - вал; 3 - подшипники

свойств ротора и упругодемпфирующих свойств опор. Определены места установки КЭ опор качения, для соответствия упругим свойствам реального ротора вращающемуся в подшипниках качения.

Рассмотрены основные динамические воздействия на ротор. Получены графики изменения радиальной и осевой силы, а также крутящего момента от подачи ЦН. Текущие значения подачи и напора насоса определены путем проведения параметрических испытаний на стенде ОАО «Уралги-дромаш». Динамические нагрузки на ротор ЦН приведены в табл. 1.

Таблица 1

Динамические нагрузки для трех режимов работы ЦН_

Режим Подача Напор Радиаль- Осевая Крутящий

подачи ЦН <2, м3/ч Я, м ная сила сила, момент

/V, кН Л, кН Мкр, кНм

Нулевой 0 34,0 7,06 18,6 0,63

Оптимальный 999 32,4 0,01 17,8 1,14

Максимальный 1294 27,4 3,73 15,0 1,25

Для каждого режима (табл. 1) составлена расчетная схема нагружения ротора как балки на шарнирно-жестких опорах. При этом получены картины напряженного состояния и определены места максимальных эквивалентных напряжений в роторе (табл. 2).

Таблица 2

Режим Составляющие реакций Результи- Экв. на-

подачи ЦН опор вдоль осей X и У, кН рующая реак- пряжения

ций опор, кН кв В

опасном

К* ** К, А В сечении, МПа

Нулевой -2,55 -12,53 1,32 7,10 12,8 7,22 28,6

Оптимальный — 1,87 — -0,35 1,87 0,35 30,2

Максимальный 3,87 8,57 -2,01 -3,82 9,4 4,32 36,2

Выявлено, что наиболее опасным режимом для ротора оцениваемым по эквивалентным напряжениям является режим максимальной подачи ЦН

^^ теристики ротора"

о 4,о 8,1 12,1 16,1 20,1 24,2 28,2 32,2 36,2 МП а Определены свобод-Рис. 3. Картина напряженного состояния ротора ЦН ные колебания ротора

с учетом упругоподатливых свойств опор при действии центробежных сил инерции неуравновешенных масс ротора, а также присоединенной массы водоволок-нисгой суспензии и гироскопического момента (табл 3)

Таблица 3

Свободные колебания ротора__

Форма колебаний ротора Значения частот свободных колебаний, Гц Действительный запас устойчивости

жесткие опоры упругие опоры упругие, с присоединенной массой жидкости

без гиро- скопическ момента с гироско-пическ моментом

изгибная, вдоль оси X 81,9 72,9 63,6 64,3 4,0

изгибная, вдоль оси У 81,9 75,6 66,1 68,9 4,3

крутильная, вокруг оси Ъ 111,1 111,1 100,1 100,1 6,3

осевая, вдоль оси Ъ 806,3 263,4 250,4 250,4 15,6

Для подтверждения достаточной жесткости ротора выполнен расчет собственных частот колебаний при условии изменения скорости вращения ротора с 16 с"1 до 24,1 с"1 Собственные частоты ротора изменяются не существенно в пределах 3 % Вынужденные колебания ротора определены с учетом диссипативных сил введением коэффициента неупругих сопротивлений Резонансные состояния ротора определены путем определения источников его возбуждения (табл 4) при измерении и анализе вибрации подшипниковых опор в трех взаимно перпендикулярных направления

Таблица 4

Частоты возбуждения колебаний ротора_

Оборотная £,б, Г'ц Лопастная Гц

Номер гармоники, к Номер гармоники, к

4 5 6 7 8 9 10 1 2

66 82,5 99 115,5 132 148,5 165 99 198

Для низкочастотной области ротора, которой соответствует диапазон до 200 Гц, построена резонансная диаграмма Кэмпбэлла (рис 4) Найдены возможные резонансные и критические режимы ротора Определены наиболее вероятные резонансы на рабочей частоте вращения ротора (103,7 рад/с) по формам, соответствующим частотам 64,2 Гц и 100,1 Гц

Наличие резонансной частоты отмечено точкой 1 и критической точкой 4, что соответствует изгибной форме колебания ротора вдоль оси X и крутильной форме вокруг оси Ъ Путем нагружения КЭ модели ротора последовательностью единичных импульсов получены отклики на опорах в виде

амплитудно-частотных характеристик соответственно на частотах 64,2 Гц и

100,1 Гц Анализ вибрации подтип-

¿Г

Диаграмма Кэмпбелла

£

а юод

а о

--4

ников подтвердил наличие данных частот в спектре при возбуждении колебаний ротора высшими гар-

§ д § мониками ло-$ п^-ТГ,-^ л $ пастной и оборотной частоты (рис 5) Установлено,

что наибольшую опасность представляет резонанс ротора по форме крутильных колебаний на частоте 104,6 Гц Так как крутильные колебания с малым демпфированием в сочетании с поперечными колебаниями ротора „„,„ значительно уве-

м^с СтячаНасосы Агр*»Т45Ю 70ЧМН21П) ЗмирДвп10ЛЭО5 время №1237

ОДцуйуромнь. 2^422 Рюброс 0.1418

198

165 ю 148,5 § 132 * 115,5 | 99 8 82,5 1

66 "8

0 «

1

Скорость вращения ротора, рад/с Рис 4 Частотная диаграмма для ротора ЦН

личивают механические напряжения в вале, что приводит к появлению трещин При эксплуатации выявлено несколько случаев аварийных остановов насосов по причине разрушения (среза) вала под РК, в

том числе и у смесительных насосов Причинами внезапного разрушения валов являются статические напряжения, которые увеличиваются до уровня опасных при резонансных колебаниях ротора В результате исследования определены свободные колебания и резонансные режимы работы ротора

3. Моделирование гидродинамических процессов и определение параметров напряженно-деформированного состояния основных

деталей центробежных насосов В разделе приведено моделирование гидродинамических процессов в ЦН дня определения параметров трехмерного потока жидкости в проточной части и

» (О N

О)

см

(О 5 с с ; £

О N <о СО ® СО (о я ч о —* <0

{ ~<о ц 9 7 А 1 —

№ 1 1 * * »Дм а | с е- .=3 ....... (О &

Рис 5 Спектр вибрации корпуса приводного подшипника

пазухах ЦН, а также методика расчета параметров объемного напряженно-деформированного состояния ротора и спирального отвода насоса и расчет уровней колебания СО насоса на лопастной частоте.

Волокнистая суспензия представляет модель сплошной жидкой среды в области проточной части ЦН. Для определения картины ее внутреннего течения построена конечно-объемная модель сплошной жидкой среды, которая представлена в виде расчетной области с назначенными граничными условиями. Построение конечно-объемной модели выполнено на основе натурной геометрии проточной части рабочего колеса и спирального отвода ЦН. Для создания расчетной области течения использована гибридная сетка с общим количеством ячеек 916,5 тысяч.

Для получения картины течения жидкости использован метод контрольных объемов. При анализе работ по гидродинамике показано, что при турбулентном режиме движения используют криволинейный участок полной реологической характеристики волокнистой суспензии с низкой концентрацией волокна до 2 %. При этом значение реологического параметра волокнистой суспензии и чистой воды отличаются не существенно. Так как для ЦН характерен турбулентный режим движения потока, при расчете использована диссипативная двухпараметри-ческая «к-е» модель турбулентности. Для режимов нулевой, оптимальной и максимальной подач граничные условия и параметры «к-а> модели турбулентности определены по результатам параметрических испытаний ЦН.

Сходимость результатов расчета достигнута по параметрам средней скорости потока й, кинетической энергии турбулентных пульсаций к, скорости диссипации энергии е среднеквадратической невязки (остатка) уровня 0,001. По разнице полного давления на выходе и входе потока (рис. 6) определено полное давление в конечно-объемной модели, которое сравнивалось с давлением ЦН при испытаниях для каждого режима работы. При этом значения статического

давления в проточ-¡р,к=2.43е+04.Па ной части использо-

ваны в качестве основной нагрузки на детали ЦН. Одновременно получена картина распределения статического давления на периферии РК (рис. 7). По графику

„ **** п видно, что при про-

Полиое давление, пи г г

Рис. 6. Картина течения в режиме максимальной подачи хождении лопасти

РК мимо языка спирального отвода наблюдается скачок (импульс) статического давления. При вращении РК периодическая последовательность таких скачков генерирует пульсацию потока жидкости и вибрацию ЦН на лопастной частоте. Для режима макси-

3.12е+05,/Л

мальной подачи амплитуда равна 0,05 МПа, что составляет 17,6 % относительно давления ЦН. Для режима оптимальной подачи амплитуда равна 0,04 МПа, что составляет относительно давления насоса 12 %. Для режима нулевой подачи амплитуда импульсов статического давления равна 0,055 МПа, что составляет 15,4 % относительно давления насоса.

i.eot.ss

I.3U-05 i.eie.05 7.)0е.04

М4«'И - I.E3.-04

-i.ei.-fu -7.79«*04

-i.iet*os

ш -l.J7o.05

Ш -l.S7.-05

Н -J.S7.-05

■ -J.J7.-05

И -2.57. *05

Статическое давление, Па

3.2Se>s5 3.00е*05

и" 2.80e'0s

Я

Я 2.69«'05

ч

2.»8е-05 3 2.2»c«S g г.Мс-вб

, 1.60е-«5 |3 Ы0е.05 О !.J5«.»5

Длина наружного диаметра РК, м Рис. 7. Распределение статического давления жидкости на периферии РК

Дня нагружения наружных поверхностей основного и покрывающего дисков РК получена картина распределения статического давления в пазухах ЦН. Адекватность расчета подтверждена значением расхода жидкости через пазухи, который составляет на оптимальном режиме 2,4 % от расхода ЦН. Для определения картины напряженного состояния ротора и спирального отвода построены их КЭ модели с помощью регулярных сеток.

Напряженное состояние ротора определено при действии статического давления. Установлены максимальные напряжения, которые возникают на валу, где расположено РК и в местах соединения лопасти с покрывающим диском, а также на кромках лопасти. Аналогичный расчет проведен при одновременном действии статического давления и крутящего момента со стороны привода При этом напряжения в опасном сечении на валу увеличиваются в 4,5 раза (рис. 8). Параметры НДС ротора определены для режимов нулевой, оптимальной и максимальной подач насоса. Установлено, что ротор наиболее нагружен в режиме максимальной подачи насоса. При нагружении спирального отвода статическим давлением в соответствии с картиной течения определены параметры НДС. Установлено, что в режиме нулевой подачи напряжения в языке составляют 6,1 МПа, что больше на 20 % относительно напряжений в режиме оптимальной подачи ЦН. В

МПа

Рис. 8. Напряженное состояние ротора ЦН

режиме максимальной подачи ЦН напряжения составляют 3,7 МПа, что соответственно меньше на 28 % относительно напряжений в режиме оптимальной подачи

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

Т = 2я1саг, (2)

где ъ - число лопастей рабочего колеса, а - угловая скорость вращения РК

Импульсы силы статического давления определены на основе амплитуды давления на лопастях РК дня соответствующих режимов работы насоса. Колебания определены на поверхности СО (рис 9, а) в месте, где располагают точку измерения параметров вибрации Полученный уровень колебаний (рис 9,6) сравнивался с виброперемещением СО определенного с помощью спектра виброскорости

Рис 9 Колебания СО во временной (а) и частотной (б) области

4. Экспериментальное определение параметров центробежных насосов и исследование вибрации В разделе приведена методика измерения параметров ЦН в соответствии с ГОСТ 6134-87 при проведении натурных испытаний на стенде ОАО «Уралгидро-маш» Испытания проведены для подтверждения результатов численного расчет интегральных параметров конечно-объемной модели сплошной жидкой среды (габл 5)

Таблица 5

Сравнение параметров центробежного насоса _

Режим подачи ЦН Давление конечно-объемной модели жидкой среды, МПа Давление насоса при испытаниях рн, МПа Погреш ность, %

Оптимальный 0,34 0,336 1,2

Нулевой 0,365 0,357 2,2

Максимальный 0,284 0,283 1,6

Экспериментальное определение параметров вибрации ЦН проведено в системе массоподвода БМ №14 ОАО «МБП Сыктывкарский ЛПК» При этом использована аппаратура и программное обеспечение фирмы «ДИАМЕХ 2000» Измерение вибрации проведены в соответствии с требованиями ГОСТ 26493-85, пьезоак-селерометры крепились в соответствии с ГОСТ Р ИСО 5348-99, точки для измерения вибрации выбирались в соответствии с рекомендациями ГОСТ ИСО 10816-1-97

С помощью разработанного маршрута выполнен необходимый сбор данных для анализа вибрации, который проведен с помощью спектра и спектра огибающей сигнала При этом определены подшипниковые частоты первых пяти гармоник для анализа возможных дефектов Замеры параметров вибрации выполнены в трех взаимно перпендикулярных направлениях на корпусе подшипников ЦН, а также на входном патрубке и спиральном отводе ЦН Результаты экспериментальных исследований вибрации подтверждают наличие в частотном спектре гармоник указывающих на резонанс ротора по форме изгибных и крутильных колебаний (табл 6)

Таблица 6

Сравнение частот вынужденных колебаний_

Форма колебаний ротора Значения частот, Гц

Численный расчет Спектр вибрации

Изгибная 64,2 62,8

Крутильная 100,1 104,6

В результате моделирования взаимодействия поверхности ротора ЦН с потоком жидкости получена его объемная картина напряженного состояния, которая подтверждается анализом случаев разрушения вала в месте крепления РК Подтверждение численного расчета уровня колебания СО на лопастной частоте получено путем сопоставления с виброперемещением спирального отвода определенного с помощью спектра виброскорости

Выводы

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

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

3 Расчет динамических характеристик ротора на основе конечно-элементной модели показал, что с учетом упругих свойств опор, центробежных сил инерции неуравновешенных масс ротора, а также присоединенной массы водоволокнистой суспензии происходит снижение собственных частот изгибных колебаний вдоль оси X на 22,3 %, вдоль оси У на 19,2 % При этом частота крутильных колебаний уменьшаются на 9,9 %, осевых колебаний ротора на 68,9 % относительно собственных частот неподвижного ротора на шарнирно-жестких опорах Гироскопический

момент на рабочем колесе приводит к увеличению собственной частоты изгибной формы колебаний вдоль оси X на 1,1 % и вдоль оси У на 4,2 %.

4 Вынужденные колебания определены с учетом коэффициента неупругих сопротивлений в виде частотной диаграммы ротора. Резонансные и критические режимы ротора в низкочастотной области получены путем построения частотной диаграммы, в которой основные частоты воздействия равны и кратны оборотной и лопастной Установлено, что наиболее опасный резонанс ротора происходит по форме 1футильных колебаний, частота которого экспериментально подтверждена путем идентификации гармоник в частотном спеюре вибрации

5 Гидродинамические процессы в насосах обуславливают нагруженность основных его деталей. На основе конечно-объемной модели сплошной жидкой среды получены параметры потока и установлены картины трехмерного течения чистой воды, реологическая характеристика которой мало отличается от полной реологической характеристики волокнистой суспензии с низкой концентрацией волокна до 2 % в диспергированном турбулентном режиме движения Интегральные параметры потока на выходе и входе конечно-объемной модели подтверждены экспериментально путем сравнения с параметрами центробежного насоса полученных на испытательном стенде ОАО «Уралгидромаш» При этом погрешность полного давления с экспериментом составила для режима оптимальной подачи 1,2 %, для режима нулевой подачи 2,2 % и для режима максимальной подачи 1,6 %.

6. По результатам моделирования трехмерного течения в проточной части насоса определена амплитуда импульсов статического давления на лопастях рабочего колеса Для режима оптимальной подачи насоса амплитуда равна 0,04 МПа, что составляет 12 % относительно его давления; для режима нулевой подачи 0,055 МПа, что составляет 15,4 % относительно давления насоса; для режима максимальной подачи 0,05 МПа, что составляет 17,6 % относительно давления насоса.

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

8 Разработанные методы расчета собственных и вынужденных колебаний пригодны при анализе вибрации и диагностике рабочих режимов центробежных насосов массопододящих систем, а также при проектировании с целью снижения их общей виброактивности Методика расчета параметров напряженно-деформировашого состояния применима в целях снижения механических напряжений в роторе и спиральном отводе Методика моделирования трехмерного течения жидкости в проточной части и определения амплитуды импульсов статического давления пригодна для массных и смесительных насосов массоподводящих сис-

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

Основное содержание диссертации изложено в публикациях:

1 Соколов Е В Диагностическая имитационная модель смесительного насоса в массоподводящей системе [Текст] / Е.В Соколов II Сборник материалов международной научно-технической конференции «Социально-экономические и экологические проблемы лесного комплекса»/Урал гос лесотехн ун-т-Екатеринбург, 2003 -С 166-168

2 Соколов Е В Особенности эксплуатации центробежных насосов на предприятиях целлюлозно-бумажной промышленности [Текст] /ЕВ Соколов, Д Т Ан-кудинов // Материалы научно-технической конференции студентов и аспирантов Уральского государственного лесотехнического университета апрель 2004г / Урал гос лесотехн ун-т -Екатеринбург, 2004 -С 119-122

3 Соколов Е В Анализ вибрационного состояния насосного агрегата при контрольных испытаниях как способ оценки его надежности [Текст] /ЕВ Соколов, Д Т Анкудинов, А В Феофанов // Материалы всероссийской научно-технической конференции «Математическое моделирование механических явлений» 28-29 мая 2004г / Урал гос горно-геолог акад - Екатеринбург, 2004 - С 68-72

4 Соколов Е В Применение численного метода для определения собственных частот колебаний ротора центробежного химического насоса [Текст] /ЕВ Соколов, ДТ Анкудинов // Материалы II всероссийской научно-технической конференции студентов и аспирантов 17-26 апреля 200бг / Матер науч -техн конф Ч 1 - Урал гос лесотехн ун-т-Екатеринбург, 2006 -С 188-191

5 Соколов Е В Динамические процессы нагружения деталей центробежных химических насосов [Текст] / ЕВ Соколов, ДТ Анкудинов, А В Феофанов // «Насосы и оборудование» -2006 -№2 -С 22-24

6 Соколов Е В Моделирование течения жидкости в центробежных насосах [Текст] / Е В Соколов, Д Т Анкудинов // Труды Евразийского симпозиума «Деревообработка технологии, оборудование, менеджмент XXI века» 20-21 сентября 2006г Урал гос лесотехн ун-т - Екатеринбург, 2006 - С 217-221

7 Соколов ЕВ Моделирование гидродинамических процессов в гидротранспорте [Текст] /ЕВ Соколов, Д Т Анкудинов // Известия вузов Горный журнал -2007 -№1 -С 84-88

8 Соколов Е В Численное моделирование рабочих процессов центробежных насосов как способ повышения эффективности его работы [Текст] / Е В Соколов, Д Т Анкудинов // Вестник Московск государ универ леса Лесной вестник -2007 - №8 - С 165-169

9 Соколов Е В Повышение эффективности работы центробежных насосов массоподводящих систем [Текст] /ЕВ Соколов // Целлюлоза Бумага Картон -2008 -№2 -С 70-73

Подписано в печать 15 04 2008 Формат 60x84/16 Уел печ л 1,0 Тираж 100 экз Заказ № 131

Размножено с готового оригинал-макета в типографии «Уральский центр академического обслуживания» 620219, г Екатеринбург, ул Первомайская, 91

Оглавление автор диссертации — кандидата технических наук Соколов, Евгений Вячеславович

ПЕРЕЧЕНЬ СОКРАЩЕНИЙ.

ВВЕДЕНИЕ.

1. Анализ состояния вопроса и постановка задачи.

1.1 Насосное оборудование массоподводящих систем.

1.2 Назначение, устройство и классификация центробежных насосов.

1.3 Теоретические представления о работе центробежных насосов.

1.4 Отказы центробежных насосов массоподводящих систем

1.5 Анализ работ по динамике, гидродинамике и прочности центробежных насосов.

1.6 Выводы и постановка задачи исследования.

2. Моделирование и исследование динамических процессов в центробежных насосах.

2.1 Постановка задачи.

2.2 Динамическая модель центробежных насосов консольного типа.

2.3 Динамические воздействия.

2.3.1 Поперечные силы.

2.3.2 Осевые силы.

2.3.3 Крутящий момент.

2.4 Прочность деталей центробежных насосов. Статическая прочность ротора.

2.5 Колебания ротора центробежных насосов консольного типа.

2.5.1 Свободные колебания ротора.

2.5.2 Вынужденные колебания ротора.

2.6 Выводы по разделу.

3. Моделирование гидродинамических процессов и определение параметров напряженно-деформированного состояния основных деталей центробежных насосов.

3.1 Постановка задачи.

3.2 Описание модели сплошной жидкой среды.

3.2.1 Построение геометрической модели рабочего колеса.

3.2.2 Построение геометрической модели спирального отвода

3.2.3 Построение конечно-объемной модели сплошной жидкой среды.

3.3 Моделирование гидродинамических процессов в центробежных насосах консольного типа.

3.3.1 Определение трехмерного течения жидкости.

3.3.2 Определение плоского течения жидкости.

3.4 Определение параметров напряженно-деформированного состояния ротора и спирального отвода.

3.5 Выводы по разделу.:.

4. Экспериментальное определение параметров центробежных насосов и исследование вибрации.

4.1 Постановка задачи.

4.2 Методика определения параметров центробежных насосов.

4.3 Методика измерения вибрации центробежных насосов на месте эксплуатации.

4.4 Результаты экспериментальных исследований.

4.5 Выводы по разделу.

ВЫВОДЫ.

Введение 2008 год, диссертация по технологии, машинам и оборудованию лесозаготовок, лесного хозяйства, деревопереработки и химической переработки биомассы дерева, Соколов, Евгений Вячеславович

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

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

Непрерывная и безаварийная работа предприятий ЦБП зависит от многих факторов. Одними из главных является надёжная работа узлов и механизмов БМ. Кроме этого необходимо чтобы системы, обеспечивающие непрерывную работу БМ, также работали надежно. К таким системам, в частности, относят массоподводящую, которая предназначена для непрерывного подвода бумажной массы к напорному ящику с последующим напуском ее на сетку формующей части БМ [1,2].

Подвод бумажной массы к напорному ящику включает в себя несколько технологических стадий. На каждой стадии количественные и качественные изменения свойства бумажной массы изменяют с помощью технологических агрегатов и установок входящих в систему. Бумажная масса проходит стадии разбавления, очистки, деаэрации, сортирования, после чего поступает в напорный ящик БМ. Для подачи бумажной массы используют центробежные насосы (ЦН) специального назначения [3].

ЦН эксплуатируют в цехе БМ, подготовки древесной массы, варочном, регенерации щелока. Таким образом, ЦН являются одними из основных технологических агрегатов входящих в состав поточных линий производства бумаги, картона, целлюлозы и древесной массы [4]. Непрерывная и безаварийная работа ЦН влияет на производительность и эффективность работы БМ. Работа ЦН в нестационарном режиме в системах массоподвода влияет на процесс формирования капиллярно-пористой структуры бумажного полотна и на качество продукции в целом [5]. Продолжительная работа в нестационарном режиме снижает долговечность и приводит к аварийным отказам с разрушением основных деталей ЦН [6,7], поэтому для обеспечения нормальной эксплуатации ЦН при использовании систем мониторинга и диагностики необходимы научно-обоснованные рекомендации.

Наиболее актуальными в таком случае представляются исследования динамических и гидродинамических процессов в ЦН. Данные исследования позволят разработать методы анализа динамических характеристик ротора, получить картину течения жидкости в проточной части насосов, определить параметры напряженно-деформированного состояния (НДС) основных деталей насосов и установить причины их разрушения в эксплуатации.

Применение разработанных и экспериментально проверенных методик расчета на практике позволят обоснованно вносить изменения в конструкцию ЦН для повышения эффективности эксплуатации в массоподводящих системах БМ и других подобных сложных гидравлических системах.

Целью работы является повышение эффективности работы массоподводящих систем бумагоделательных машин путем решения следующих задач:

- моделирования динамических процессов нагружения деталей в насосах;

- выявления собственных форм и частот колебаний ротора и определения амплитудно-частотных характеристик;

- моделирования гидродинамических процессов в насосах и определение параметров потока жидкости;

- определения параметров напряженно-деформированного состояния ротора и спирального отвода насоса;

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

Методы исследования. Теоретические исследования проведены с привлечением основных положений динамики машин, теории колебаний, методов конечных элементов и контрольных объемов, гидродинамики и теории лопастных машин. Экспериментальные исследования выполнены на испытательном стенде завода ОАО «Уралгидромаш», на месте эксплуатации ЦН в массопод-водящей системе БМ №14 ОАО «МБП Сыктывкарский ЛПК». При экспериментальных исследованиях использован анализатор вибрации «Топаз» фирмы «ДИАМЕХ», а также лицензионное программное обеспечение. Геометрические построения выполнены в автоматизированной системе трехмерного моделирования Unigraphics. Численные расчеты проведены в пакетах CAE и CFD анализа: ANSYS, CFX.

Предмет исследования. Динамические и гидродинамические процессы в центробежных насосах.

Объект исследования. Центробежные насосы массоподводящих систем.

Научная новизна работы. Реализовано численное моделирование динамики ротора ЦН с распределенными параметрами. Определены формы пространственных колебаний ротора с учетом упругих свойств опор, под воздействием центробежных сил инерции неуравновешенных масс ротора, присоединенной массы водоволокнистой суспензии и гироскопического момента. Получены значения резонансов и критических режимов ротора в низкочастотной области. Реализовано численное моделирование гидродинамических процессов в ЦН. Определены параметры трехмерного потока жидкости в проточной части для режимов оптимальной, нулевой и максимальной подач насоса. Получены значения амплитуд импульсов статического давления на лопастях рабочего колеса. Установлена картина распределения статического давления жидкости в пазухах насоса. Определены параметры объемного напряженно-деформированного состояния ротора и спирального отвода, а также уровни колебания спирального отвода на лопастной частоте.

Достоверность научных результатов исследования обуславливается использованием при моделировании динамических и гидродинамических процессов в ЦН развитых теорий колебаний, методов конечных и контрольных объемов, известных программных средств. При этом значения резонансных частот ротора по формам крутильных и изгибных колебаний экспериментально подтверждены путем идентификации гармоник в частотном спектре вибрации. Интегральные параметры потока на выходе и входе конечно-объемной модели подтверждены экспериментально путем сравнения с параметрами центробежного насоса полученных на испытательном стенде ОАО «Уралгидромаш», при этом погрешность полного давления с экспериментом составила для режима оптимальной подачи 1,2 %, нулевой подачи 2,2 % и для режима максимальной подачи насоса 1,6 %. Апробация расчета проведена на серийно выпускаемом для ЦБП центробежном насосе консольного типа.

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

Результаты работы следует использовать при проектировании ЦН на машиностроительных заводах выпускающих насосы специального назначения для ЦБП, а также при модернизации насосов на предприятиях ЦБП. Использование разработанных методик расчета динамики и гидродинамики ЦН подтверждены актом внедрения работы на ОАО «Уралгидромаш» и ОАО МБП «Сыктывкарский ЛПК». Результаты исследования используются в учебном процессе для подготовки инженеров-механиков по специальности 150504 в курсе «Теория и конструкция машин и оборудования отрасли», а также в разделе диагностика потокообразующих и потокопроводящих систем курса «Диагностика машин и оборудования».

Апробация работы. Основные положения диссертационной работы доложены на научно-технической конференции студентов и аспирантов Уральского государственного лесотехнического университета в апреле 2004г. (УГЛТУ, Екатеринбург); всероссийской научно-технической конференции, посвященной 90-летию Уральской государственной горно-геологической академии «Математическое моделирование механических явлений» 28-29 мая 2004г. (УГГГА, Екатеринбург); II всероссийской научно-технической конференции студентов и аспирантов 17-26 апреля 2006г. (УГЛТУ, Екатеринбург); евразийском симпозиуме «Деревообработка: технологии, оборудование, менеджмент XXI века» 20-21 сентября 2006г. (УГЛТУ, Екатеринбург).

Публикации. По теме диссертации опубликовано 9 печатных работ. Работа выполнялась в рамках госбюджетной темы по единому наряд-заказу Минобразования РФ [8].

Научные положения выносимые на защиту:

- результаты моделирования собственных и вынужденных колебаний ротора с учетом упругих свойств опор, под воздействием центробежных сил инерции неуравновешенных масс ротора, присоединенной массы волокнистой суспензии и гироскопического момента;

- результаты моделирования трехмерного течения жидкости в проточной части ЦН;

- методика определения параметров напряженно-деформированного состояния ротора и спирального отвода ЦН;

- результаты расчета колебаний спирального отвода на лопастной частоте.

Объём и структура диссертационной работы. Диссертационная работа состоит из введения, четырех разделов, заключения, списка литературы и двух приложений. Она включает 176 страниц машинописного текста, 95 ил., 35 табл., 134 наименования использованных источников, в том числе 18 иностранных.

В первом разделе проведен анализ эксплуатируемых ЦН в массоподво-дящих системах БМ, проанализированы отказы ЦН, выявлены факторы валяющие на динамику и прочность основных деталей ЦН, выполнен обзор работ по динамике, гидродинамике и прочности ЦН, сделаны выводы и произведена постановка задачи исследований.

Во втором разделе выявлены основные динамические воздействия на ротор, построена объемная балочная конечно-элементная модель ротора, приведен анализ статической прочности и определены параметры свободных и вынужденных колебаний ротора. Приведены частотные спектры вибрации опор ротора.

В третьем разделе проведено моделирование гидродинамических процессов, построена конечно-объемная модель сплошной жидкой среды, определены параметры потока и получена трехмерная картина течения жидкости в проточной части ЦН для режимов оптимальной, нулевой и максимальной подач ЦН, а также картина течения в боковых пазухах ЦН. Определена объемная картина НДС ротора и спирального отвода и проведен расчет колебаний спирального отвода на лопастной частоте.

В четвертом разделе приведены средства и методика определения основных параметров ЦН на испытательном стенде завода-изготовителя ОАО «Урал-гидромаш», а также методика измерения вибрации на месте эксплуатации ЦН в массоподводящей системе БМ №14 ОАО «МБП Сыктывкарский ЛПК». Результаты численного расчета подтверждены путем сравнения с экспериментальными данными.

Заключение диссертация на тему "Моделирование и исследование динамических и гидродинамических процессов в центробежных насосах массоподводящих систем бумагоделательных машин"

162 ВЫВОДЫ

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

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

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

3. Расчет динамических характеристик ротора на основе конечно-элементной модели показал, что с учетом упругих свойств опор, центробежных сил инерции неуравновешенных масс ротора, а также присоединенной массы водоволокнистой суспензии происходит снижение собственных частот изгибных колебаний вдоль оси X на 22,3 %; вдоль оси У на 19,2 %. При этом частота крутильных колебаний уменьшаются на 9,9 %, осевых колебаний ротора на 68,9 % относительно собственных частот неподвижного ротора на шарнирно-жестких опорах. Гироскопический момент на рабочем колесе приводит к увеличению собственной частоты изгибной формы колебаний вдоль оси X на 1,1 % и вдоль оси У на 4,2 %.

4. Вынужденные колебания определены с учетом коэффициента неупругих сопротивлений в виде частотной диаграммы ротора. Резонансные и критические режимы ротора в низкочастотной области получены путем построения частотной диаграммы, в которой основные частоты воздействия равны и кратны оборотной и лопастной. Установлено, что наиболее опасный резонанс ротора происходит по форме крутильных колебаний, частота которого экспериментально подтверждена путем идентификации гармоник в частотном спектре вибрации.

5. Гидродинамические процессы в насосах обуславливают нагруженность основных его деталей. На основе конечно-объемной модели сплошной жидкой среды получены параметры потока и установлены картины трехмерного течения чистой воды, реологическая характеристика которой мало отличается от полной реологической характеристики волокнистой суспензии с низкой концентрацией волокна до 2 % в диспергированном турбулентном режиме движения. Интегральные параметры потока на выходе и входе конечно-объемной модели подтверждены экспериментально путем сравнения с параметрами центробежного насоса полученных на испытательном стенде ОАО «Уралгидромаш». При этом погрешность полного давления с экспериментом составила для режима оптимальной подачи 1,2 %, для режима нулевой подачи 2,2 % и для режима максимальной подачи 1,6 %.

6. По результатам моделирования трехмерного течения в проточной части насоса определена амплитуда импульсов статического давления на лопастях рабочего колеса. Для режима оптимальной подачи насоса амплитуда равна 0,04 МПа, что составляет 12 % относительно его давления; для режима нулевой подачи 0,055 МПа, что составляет 15,4 % относительно давления насоса; для режима максимальной подачи 0,05 МПа, что составляет 17,6 % относительно давления насоса.

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

8. Разработанные методы расчета собственных и вынужденных колебаний пригодны при анализе вибрации и диагностике рабочих режимов центробежных насосов массопододящих систем, а также при проектировании с целью снижения их общей виброактивности. Методика расчета параметров напряженно-деформированного состояния применима в целях снижения механических напряжений в роторе и спиральном отводе. Методика моделирования трехмерного течения жидкости в проточной части и определения амплитуды импульсов статического давления пригодна для массных и смесительных насосов массоподводящих систем, а также напорных сортировок, дисковых мельниц и других гидравлических машин целлюлозно-бумажного производства.

Библиография Соколов, Евгений Вячеславович, диссертация по теме Технология и оборудование химической переработки биомассы дерева; химия древесины

1. Терентьев O.A. Массоподача и равномерность бумажного полотна / O.A. Терентьев. - М.: Лесная промышленность, 1986. - 264 с.

2. Терентьев O.A. Гидродинамика волокнистых суспензий в целлюлозно-бумажном производстве / O.A. Терентьев. М.: Лесная промышленность, 1980. - 680 с.

3. Климов В.И. Гидротранспорт волокнистых материалов в целлюлозно-бумажном производстве / В.И. Климов. М.: Лесная промышленность, 1971.-280 с.

4. Хенрикссон X. Новые системы перемещения и обработки волокон в целлюлозно-бумажной промышленности Текст. / X. Хенрикссон, Т. Нисканен //Бумажная промышленность. 1989. №9 - С. 29-31.

5. Львов А.И. Теоретическое обоснование условий работы смесительных насосов бумагоделательных машин при минимальном уровне пульсаций давления Текст. / А.И. Львов, B.C. Куров, O.A. Терентьев // Известия вузов. Лесной журнал 1997. №6 - С. 64-69.

6. Косторной A.C. Исследование нестационарных характеристик потока в проточной части центробежного насоса Текст. / A.C. Косторной // Насосы и оборудование. 2005. №1 - С. 18-19.

7. Зимницкий В.А. Лопастные насосы / В.А. Зимницкий, В.А. Умов. Справочник. Л.: Машиностроение, 1986. - 334 с.

8. ГОСТ 17398-72. Насосы. Термины и определения. М.: Изд-во стандартов, 1989. - 35 с.

9. Ломакин A.A. Центробежные и осевые насосы / A.A. Ломакин. М.: Машиностроение, 1966. - 364с.

10. Михайлов А.К. Конструкции и расчет центробежных насосов высокого давления / А.К. Михайлов, В.В. Малюшенко. М.: Машиностроение, 1971. - 3041с.

11. Руднев С.С. Лабораторный курс гидравлики, насосов и гидропередач / С.С. Руднев, Л.Г. Подвидз М.: Машиностроение, 1974. 304 с.

12. ГОСТ 6134-87. Насосы динамические. Методы испытания. М.: Изд-во стандартов. 1998. - 21 с.

13. Яременко О.В. Испытания насосов / О.В. Яременко. Справочник. -М.: Машиностроение. 1976. 225 с.

14. Справочник по расчетам гидравлических и вентиляционных систем / под ред. A.C. Юрьева. С.-Пб, AHO. НПО «Мир и семья» 2001. 1154'с.

15. Идельчик И.Е. Справочник по гидравлическим сопротивлениям / И.Е. Идельчик. М.: Машиностроение, 1990. - 672 с.

16. Климов В.И. Расчет массопроводов при транспортировании волокнистых суспензий Текст. / В.И. Климов // Бумажная промышленность 1986. №9 - С. 28-31'.

17. Киприанов А.И: Расчет систем гидротранспорта волокнистых суспензий высокой концентрации Текст. / А.И. Киприанов, H.H. Калинин и. др. //Бумажная промышленность 1988. №4 - С. 16-17.

18. Волков А., Панкратов С, Парыгин А. Анализ гидродинамических качеств сетевых насосов в энергоблоках мощностью 100 МВт Текст. / А. Волков, С. Панкратов, А. Парыгин, // Насосы и оборудование. 2006. №2-С. 25-27.

19. Абрамов И.В. Совершенствование структуры планово-предупредительных ремонтов оборудования Текст. / И.В; Абрамов, Ю.В. Турыгин, Б.М. Ценципер // Бумажная промышленность 1988. №9 - С. 27-28.I

20. Карелин В.Я. Насосные станции гидротехнических систем с осевыми и диагональными насосами / В.Я. Карелин, P.A. Новодережкин. М.: Энергия, 1980.-288 с.

21. Барков A.B. Мониторинг и диагностика роторных машин по вибрации / A.B. Барков, H.A. Баркова, А.Ю. Азовцев. СПб.: СПбГМТУ, 2000. - 169 с.

22. Розенберг Г.Ш. Вибродиагностика / Г.Ш. Розенберг и др. — СПб.: ПЭИПК, 2003.-284 с.

23. Гривнин Ю.А. Математическая модель динамического воздействия кавитации Текст./ЮА.Гривнин//Известия вузов. Энергетика- 1989.№10-С. 112-117.

24. Карелин В.Я. Кавитационные явления в центробежных и осевых насосах / В.Я. Карелин. М.: Машиностроение, 1975. - 432 с.

25. Покровский Б.В. Расчет уровней лопастных вибраций центробежных насосов Текст. / Б.В. Покровский, В.Я. Рубинов // Сб. научн. тр. ВНИИгидромаш. Совершенствование насосного оборудования. 1982. - С. 115-121.

26. Ширман А.Р. Практическая вибродиагностика и мониторинг состояния* механического оборудования / А.Р. Ширман, А.Б. Соловьев. М.: t 1996.-276 с.

27. Плотников И.А. К оценке реологических свойств волокнистых суспензий Текст. / И.А. Плотников, В.Н. Непенин, А.И. Киприанов // Известия вузов. Лесной журнал 1986. №1 - С. 78-80.

28. Терентьев O.A. Влияние турбулентности потока на реологические параметры бумажной массы Текст. / O.A. Терентьев, Э.А. Смирнова, B.C. Куров // Известия вузов. Лесной журнал 1990. №2 - С. 97-101.

29. Терентьев O.A. Реологические основы процесса массоподачи на бумагоделательную машину Текст. / O.A. Терентьев, B.C. Куров, Э.А. Смирнова//Известия вузов. Лесной журнал 1992. №1 - С. 85-89.

30. Евдокимов В.П. Демпфирование гидродинамических пульсаций давления напорных сортировок с целью снижения массоемкости бумаги.: Дисс. на соиск. учен, степени канд. техн. наук. Ленинград, 1989. - 207 с.

31. Морозов В.А. Исследование работы центробежных насосов на вязкопластичных суспензиях Текст. / В.А. Морозов, А.Н. Шерстюк // Известиячвузов. Энергетика 1987. №11 - С. 116-120.

32. Морозов В.А. Расчет характеристик центробежных насосов при работе на вязкопластичных жидкостях Текст. / В.А. Морозов, А.Н. Шерстюк // Известия вузов. Энергетика 1988. №1 - С. 123-125.

33. Хитрых Д. Проектирование турбомашин: обзор моделей турбулентностей Текст. / Д. Хитрых // ANSYS Solutions. Русская редакция. 2005. №1-С.9-11.

34. Андерсон Д. Вычислительная гидродинамика и теплообмен. В 2-х т. Т. 2: Пер. с англ. / Д. Андерсон, Дж. Таннехил, Р. Плетчер. М.: Мир, 1990. - 392 с.

35. Емцев Б.Т. Техническая гидромеханика. / Б.Т. Емцев. М.: Машиностроение, 1978.-463 с.

36. Ануфриенко О.С. К определению радиальных сил, действующих на рабочие колеса центробежных насосов и компрессоров Текст. / О.С. Ануфриенко, А.Н. Шерстюк // Известия вузов. Энергетика 1988. №2 - С. 106-108.

37. Ануфриенко О.С. К определению радиальных сил, действующих на рабочие колеса центробежных машин с учетом нагрузки на покрывающий'диск Текст. / О.С. Ануфриенко, А.Н. Шерстюк // Известия вузов. Энергетика- 1989. №6 С. 110-114.

38. Шнепп В.Б. Расчет осевых сил в центробежной ступени с учетом негерметичности уплотнений Текст. /В.Б. Шнепп // Химическое и нефтяное машиностроение 1983. №8-С. 21-23.

39. Евгеньев С.С. Экспериментальное определение осевого усилия действующего на вращающийся диск Текст. / С.С. Евгеньев, В.Б. Шнепп // Химическое и нефтяное машиностроение 1970. №8 - С. 10-12.

40. Рубинов В.Я. Влияние чисел лопаток рабочего колеса и направляющего аппарата на виброакустические характеристики центробежного насоса Текст. / В.Я. Рубинов, Б.В. Покровский // Сб. научн. тр. ВНИИгидромаш. Гидромашиностроение. 1975. №46- С. 71-89.

41. Покровский Б.В. О влиянии шероховатости на энергетические и виброакустические характеристики центробежных насосов Текст. / Б.В. Покровский, В.Н. Жемчугов // Сб. научн. тр. ВНИИгидромаш. Исследование гидромашин. 1967. №36- С. 141-147.

42. Рудис М.А. К расчету на прочность корпуса центробежного насоса с двухзавитковым спиральным отводом Текст. / М.А. Рудис, Ю.Л. Креймер // Сб. научн. тр. ВНИИгидромаш. «Повышение технического уровня центробежных насосов». 1980. - С. 77-82.

43. Духовный A.R Расчет на прочность двухзавитковых спиральных отводов с применением ЭВМ Текст. / АЛ Духовный, Н.С. Федяева // Сб. научн: тр. ВНИИгидромаш. «Эксплуатационная надежность насосного оборудования» -1986. С. 141-146.

44. Марцинковский В.А. Гидродинамика и прочность центробежных насосов /В.А. Марцинковский. М.: Машиностроение, 1970. -272 с.

45. Айзенштейн М.Д. Центробежные насосы для нефтяной промышленности / М.Д. Айзенштейн. Ml: Гостехиздат, 1967. - 363 с.

46. Териков Ю.П. Определение резонансных режимов рабочего колеса центробежного насоса Текст. / Ю.П. Териков, А.Н. Нацвин // Сб. научн. тр. ВНИИгидромаш. «Совершенствование насосного оборудования». -1982. С. 134-142.

47. Блюм Э.Э. Повышение стойкости деталей насосов Текст. / Э.Э. Блюм, Э.Э. Нугуманов // Бумажная промышленность. 1986. №10 - С. 22.

48. Раевская JI.T. Унификация центробежных насосов Текст. / JI.T. Раевская, М.Б. Зуншайн, Н.К. Джемилев // Бумажная промышленность. -1987. №4-С. 26.

49. Буренин В.В. Унификация современных насосов Текст. / В.В. Буренин, Н.И. Бражников, А.Я. Новиков // Химическое и нефтяное машиностроение. 1970. №9 - С. 28-29.

50. Витвинин А.М. Исследование колебаний быстроходных бумагоделательных машин: дисс. на. соиск. учен, степени канд. техн. наук. / A.M. Витвинин-Свердловск, 1974. 147 с.

51. Санников А.А. Решение проблемы виброзащиты и вибродиагностики бумагоделательного и лесопильного оборудования: дисс. на соиск. уч. степ. докт. техн. наук / Александр Александрович Санников. Екатеринбург, 2002. - 424 с.

52. Санников А.А. Колебания трубопроводов системы подачи массы в напорные ящики бумагоделательных машин Текст. / А.А. Санников, A.M. Витвинин // Бумажная промышленность. 1972. - №21. - С. 9-10.

53. Санников А.А. Колебания смесительных насосов и узлоловителей Текст. / А.А. Санников, А.М. Витвинин //Бумажная промышленность. -1972. №20. - С. 12-13.

54. ГОСТ 26563-85. Вибрация. Технологическое оборудование целлюлозно-бумажного производства. Методы и средства защиты. М.: Изд-во стандартов, 1985. - 13 с.

55. Sulzer Pumps delivers 100 pumps to Quqihar Text. // Paper Equipment & Materials International. 2001. December. - p. 20.

56. Samarasekera H. Pipeline Pumps: Efficient and Reliable Text. / H. Samarasekera // Sulzer technical review. 1998. №4 - p. 8-11.

57. Medium Consistency technology from Sulzer Pumps Text. // Pumps & Valves. 2003. №3 - p. 22-23.

58. Vesala R. Medium Consistency MC® Pumping Systems Text. / R. Vesala // Sulzer technical review. 2002. №3 - p. 7-9.

59. Dupont P. Numerical Prediction of Cavitation: Improving Pump Design Text. / P. Dupont // Sulzer technical review. 2001. №2 - p. 24-27.

60. Dupont P. Expertise Improves Products Text. / P. Dupont, M. Cugal // Sulzer technical review. 2006. №2 - p. 22-25.

61. Lienau W. Numeric Simulation Improves Operational Reliability Text. / W. Lienau, P. Meuter // Sulzer technical review. 2004. №1 - p. 8-11.

62. Salis D. Dynamic simulation of multiphase pumps / D. Salis, E. Heintze, Y. Charron // BHR Group. 1999. - p. 11- 44.

63. Védrine D. Improving Safety & Reliabilityin Refinery Operations Text. / D. Védrine, L.Donnat // Fluent NEWS. 2003. №2 - p. 30.

64. Blanco-Marigorta E. Numerical Simulationin a Centrifugal Pump with Impeller-Volute Interaction Text. / E. Blanco-Marigorta // Proceedings of ASME FEDSM. 2000. June — p. 11-15.

65. Carey C. QNET Building Quality and Trust in Industrial CFD Text. / C. Carey // Fluent NEWS. - 2002. №2 - p. 43.

66. Cropper M. Better Pumps Right,Away Text. / M. Cropper // Sulzer technical review. 1998. №1 - p. 10-13.

67. Lienau W. Behavior of Pumps Subject to Thermal Shock Text. / W. Lienau, R. Kaiser// Sulzer technical review. 1998. №2 - p. 24-27.

68. Lienau W. Early Optimization of Large Water Transport Pump Casing Text. / W. Lienau, T. Welschinger // Sulzer technical review. 2005. №2 - p. 4-7.

69. Анкудинов Д.Т. Шахтные и пневмоколесные самоходные машины / Д.Т. Анкудинов. М.: Недра. 1984. - 255 с.

70. Бендат Дж. Прикладной анализ случайных данных / Дж. Бендат, А. Пирсол: Пер. с англ. М.: Мир, 1989. - 540 с.

71. Светлицкий В.А. Случайные колебания механических систем / В.А. Светлицкий. М.: Машиностроение, 1976. 216 с.

72. Левитский Н.И. Колебания в механизмах / Н.И.' Левитский. М.: Наука. 1988. - 336 с.

73. Вейц В.Л. Динамика машинных агрегатов / В.Л. Вейц М.-Л.: Машиностроение, 1969. - 368 с.

74. Вибрации в технике: справочник, в 6 т. Т. 3 / Колебания машин, конструкций и их элементов / Под ред. Ф.М. Диментберга. М.: Машиностроение. 1980. - 544 с.

75. Уравновешивание роторов и механизмов / Под ред. В.А. Щепе-тильникова. — М.: Машиностроение. 1978. 320 с.

76. Прочность, устойчивость, колебания. Справочник в 3 т. Т. 1 / Под ред. И.А. Биргера. М.: Машиностроение. 1968. - 570 с.

77. Диментберг Ф.М. Колебания машин / Ф.М. Диментберг, К.Г. Шаталов, А.А. Гусаров. М.: Машиностроение. 1964. - 307 с.

78. Бидерман В.Jl. Прикладная теория механических колебаний / В.Л. Бидерман. М.: Высшая школа. 1972. - 416 с.

79. Пановко Я.Г. Основы прикладной теории упругих колебаний / Я.Г. Пановко. М.: Машиностроение. 1967. - 315 с.

80. Пановко Я.Г. Введение в теорию механических колебаний / Я.Г. Пановко. М.: Наука. 1971.-240 с.

81. Тимошенко С.П. Колебания в инженерном деле / С.П. Тимошенко. -М.: Наука. 1967.-444 с.

82. Зенкевич О. Метод конечных элементов в технике / О. Зенкевич. М.: Мир. 1975. - 540 с.

83. Зенкевич О. Метод конечных элементов в теории сооружений и в механике сплошных сред / Зенкевич О, Чанг И. Пер. с англ. М.: Недра. 1974. - 240 с.

84. Сегерлинд Л. Применение метода конечных элементов / Л. Сегерлинд. М.: Мир. 1979. - 392 с.

85. Норри Д. Введение в метод конечных элементов / Д. Норри. М.: Мир. 1981.- 155 с.

86. Галлагер Р. Метод конечных элементов / Р. Галлагер. Пер. с англ. -М.: Мир. 1984.-428 с.

87. Розин Л.А. Метод конечных элементов Текст. / Л.А. Розин // Соросовский образовательный журнал. 2000. Т. 6. №4 - С. 120-127.

88. Каплун А.Б. ANSYS в руках инженера / А.Б. Каплун, Е.В. Морозов, М.А. Олферьева. М.: Едиториал УРСС, 2003. - 272 с.

89. Чигарев A.B. ANSYS для инженеров / A.B. Чигарев, A.C. Кравчук, А.Ф. Смалюк. М.: Машиностроение, 2004. - 512 с.

90. Басов К.A. ANSYS в примерах и задачах / К.А. Басов; под ред. Д.Г. Красковского. М.: Компьютер пресс, 2002. - 224 с.

91. Программа ANSYS. Краткий курс: учеб. пособие / Представительство CAD-FEM GmbH в СНГ. Снежинск, 1996. - 28 с.

92. Биргер И.А. Расчет на прочность деталей машин / И.А. Биргер, Б.Ф. Шорр, Г.Б. Иосилевич. Справочник. М.: Машиностроение, 1979. - 702 с.

93. Гольдин A.C. Вибрация роторных машин / A.C. Гольдин. М.: Машиностроение, 1999. - 344 с.

94. Ден Г.Н. Механика потока в центробежных компрессорах / Г.Н. Ден. Л.: Машиностроение, 1972. - 272 с.

95. Серенсен C.B. Несущая способность и расчеты деталей машин на прочность / C.B. Серенсен, В.П. Когаев, P.M. Шнейдерович. -М.: Машиностроение, 1975.-488 с.

96. Феодосьев В.И. Сопротивление материалов / В.И. Феодосьев. М.: Физматгиз, 1962.- 536 с.

97. Прокофьев Ю.В. К определению потерь дискового трения Текст.'/ Ю.В. Прокофьев // Сб. научн. тр. ВНИИгидромаш. Исследование гидромашин. 1967. №36- С. 15-27.

98. Вибрации в технике: справочник в 6 т. Т. 1 / Колебания линейных систем / Под ред. В.В. Болотина. М.: Машиностроение. 1978. - 352 с.

99. Лихнов П.П. Динамика системы цилиндр-оболочка / П.П. Лихнов. -М.: Машиностроение. 1988. 152 с.

100. Kothnke Р. / ANSYS, Inc. Theory Reference. Release 8.0 / P. Kothnke. Houston, 1994.

101. Конюхов A.B. Основы анализа конструкции в ANSYS / A.B. Конюхов.-Казань.: КГУ. 2001.- 102 с.

102. Перель Л.Я. Подшипники качения: Расчет, проектирование и обслуживание опор / Л.Я. Перель. Справочник. М.: Машиностроение, 1992. - 606 с.

103. Гинесин Л.Ю. Применение MSC.NASTRAN для анализа динамики роторов / Л.Ю. Гинесин. M.: MSC Software. - 2000. - 28 с.

104. Чимде А.Г. Вибрационное проектирование и диагностирование дисковых мельниц: дисс. на соиск. уч. степ. кап. техн. наук / A.F. Чимде. Екатеринбург, 2004. - 164 с.

105. Мэтыоз Дж. Численные методы. Использование MatLab / Дж. Мэтьюз, Г. Финк. Пер. с англ. - Mi: Вильяме. 2001. - 720 с.

106. Кулинич С.И. Вибродиагностика неисправностей; машинного оборудования; http://www.vibration.ru/rotdin:shtml#top.l 1 .

107. Лойцянский ЛЖ. Механика; жидкости, и газа / Л-.Г. Лойцянский. — М.: Наука, 1973.-848 с.113; Коннор Дж. Метод конечных элементов в механике жидкости. Пер. с англ. / Дж. Коннор, К.Бреббия. - Л.: Судостроение, 1979. - 264 с.

108. FLUENT 6.1 Tutorial Guide / Fluent Inc. Release 6.0 / FLUENT, Belgium.-2003. p. 864.

109. ANSYS CFX Tutorial / ANSYS Inc. Release 10.0 / ANSYS CFX, ANSYS Canada Ltd. 2004. p. 368.

110. Воронин И.В. Газодинамический расчет первой ступени экспериментального двухступенчатого компрессора с использованием программного комплекса CFX Текст. /ИВ! Воронич, ЛФ. Ивчик, В.Н. Коныпин и др. // САПР и Графика. Апрель. 2005. С. 26-30.

111. Хитрых Д. Модуль Pre/Post ANSYS CFX Текст. / Д. Хитрых, М-. Плыкин // ANSYS Solutions. Русская редакция. 2006. №1-2-С.24-27.

112. Соколов Е.В. Повышение эффективности работы центробежных насосов массоподводящих систем Текст. / Е.В. Соколов // Целлюлоза. Бумага1. Картон. 2008. - №2 - С. 70-73.

113. Соколов Е.В. Численное моделирование рабочих процессов центробежных насосов как способ повышения эффективности' его работы Текст. / Е.В. Соколов, Д.Т. Анкудинов // Вестник Московск. государ, универ. леса. Лесной вестник.-2007.-№8-С. 165-169.

114. Слезкин H.A. Динамика вязкой несжимаемой жидкости / H.A. Слезкин. -М.: Физматгиз, 1955. 519 с.

115. Соколов Е.В. Моделирование гидродинамических процессов в гидротранспорте Текст. / Е.В. Соколов, Д.Т. Анкудинов // Известия вузов. Горный журнал. 2007. - №1. - С. 84-88.

116. ГОСТ 26493-85 Вибрация. Технологическое оборудование целлюлозно-бумажного производства. Нормы вибрации. Технические требования. М.: Изд-во стандартов, 1985. - 8 с.

117. ГОСТ ИСО 2954-97. Вибрация машин с возвратно-поступательным и вращательным движением. Требования к средствам измерений. М.: Изд-во стандартов, 1998. - 6 с.

118. ГОСТ 30296-95. Аппаратура общего назначения для определения основных параметров вибрационных процессов. Общие технические требования.- М.: Изд-во стандартов, 1966. 16 с.

119. ГОСТ Р ИСО 5348-99. Вибрация и удар. Механическое крепление акселерометров. М.: Изд-во стандартов, 1999. - 12 с.

120. ГОСТ ИСО 10816-1-97. Вибрация. Контроль состояния машин по результатам измерений вибрации на не вращающихся частях. 4.1. Общие требования. М.: Изд-во стандартов, 1998. - 12 с.