Научная статья на тему 'Моделювання деформованих профілів лопатей осьового компресора газоперекачувального агрегату'

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

CC BY
51
23
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
МОДЕЛИРОВАНИЕ / ПРОФИЛЬ ЛОПАСТЕЙ / ГАЗОПЕРЕКАЧИВАЮЩИЙ АГРЕГАТ / MODELING / BLADE PROFILE / GAS PUMP

Аннотация научной статьи по механике и машиностроению, автор научной работы — Штаєр Л. О.

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

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

GAS PUMP AXIAL COMPRESSOR STRAINED BLADE PROFILE MODELING

The longterm exploited gas pump units constructive elements strained profile modeling question are considered.

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

Таблица 1

Данные измерений и расчетов износа диаметра дисков

Таблица 2

Данные измерений и расчетов толщины лезвия дисков

Время t, ч Да, мм Оа, мм nj А- Е- 3/D(Aj) 5p(Ej)

Стендовые испытания

8 0,27 0,23 24 0,86 2,49 0,04 1,1

16 0,40 0,25 24 1,32 2,65 0,04 1,1

24 0,55 0,30 24 1,42 2,83 0,04 1,1

32 0,64 0,32 24 1,58 3,28 0,04 1,1

Эксплуатационные испытания

8,2 0,20 0,17 24 3,55 0,95 1,3 1,8

16,4 0,37 0,32 24 3,92 1,00 1,3 1,8

24,7 0,48 0,37 24 4,44 1,13 1,3 1,8

32,5 0,65 0,43 24 5,22 1,26 1,3 1,8

Время t, ч И, мм Ои, мм nj А- Е- 3]D(A;) ^D(Ej)

Стендовые испытания

8 0,20 0,07 24 0,008 0,82 0,04 0,9

16 0,25 0,09 24 0,023 0,76 0,04 0,9

24 0,29 0,13 24 0,127 0,53 0,04 0,9

32 0,30 0,18 24 0,175 0,36 0,04 0,9

Эксплуатационные испытания

8,2 0,07 0,03 24 0,25 0,37 1,2 2,1

16,4 0,15 0,06 24 0,42 0,39 1,2 2,1

24,7 0,19 0,10 24 0,71 0,43 1,2 2,1

32,5 0,24 0,12 24 0,75 0,51 1,2 2,1

5. Выводы

Полученные данные позволяют сделать следующие выводы.

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

2. На основании полученных данных можно реко-

мендовать линейную модель ускоренных испытаний для ----------------□ □------------------

Розглянуто питання моделювання деформованих профілів елементів конструкції газоперекачувальних агрегатів, що тривалий час перебувають в експлуатації.

Ключові слова: моделювання, профіль лопатей, газоперекачувальний агрегат.

□--------------------------------□

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

Ключевые слова: моделирование, профиль лопастей, газоперекачивающий агрегат.

□--------------------------------□

The long-term exploited gas pump units constructive elements strained profile modeling question are considered.

Keywords: modeling, blade profile, gas pump.

----------------□ □------------------

1. Вступ

Питання оцінки реального технічного стану елементів конструкцій газоперекачувальних агрегатів (ГПА) компресорних станцій, які зазнають дії аеродинамічних

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

Литература

1. Качинский Н. А. Физика почв [Текст] / Н. А. Качин-ский. — М. : Высшая школа, 1985. — 224 с.

2. Sach G. Zeitschrift ftr Metallkunde. — T. t9, 1987. — 376 p.

УДК 519.876.5

навантажень, залишається актуальним, оскільки на даний момент основним методом такої оцінки є візуальний контроль. В той же час існують методи розрахунку аеродинамічних параметрів крилових профілів (до яких можна віднести досліджувані елементи конструкції ГПА), які

моделювання деформованих профілів лопатей осьового компресора газоперекачувального

АГрЕГАТу

Л. О. Штаєр

Кандидат технічних наук, асистент Кафедра комп’ютерних технологій в системах управління і автоматики Івано-Франківський національний технічний університет нафти і газу вул. Карпатська, 15, м. Івано-Франківськ, Україна, 76019 Контактний тел.: 099-094-16-05 E-mail: lida.shtayer@gmail.com

Е

дозволяють встановити залежність між геометрією профілю та вказаними параметрами, до яких слід віднести коефіцієнт індуктивного опору, коефіцієнт підіймальної сили та коефіцієнт обертового моменту відносно передньої кромки [1].

2. моделювання деформованих профілів лопатей

При реалізації методу розрахунку плоскопаралельної течії навколо крилового профілю одержується інтегральне рівняння Фредгольма II роду відносно дотичної компоненти швидкості газового потоку до профілю:

Ve(6o) = 1 / Ve(6)K(6•,6o)d6 + 2 дф (бо),

(1)

прямокутник. В такому випадку можна використати відому апроксимацію прямокутника строго випуклими фігурами, які описуються залежностями [3]:

х2п + у2п _,

а2п + Ь2п _ і.

(8)

При п = 1 (8) задає еліпс, при п — прямокутник зі сторонами 2а та 2Ь, в якого точка перетину діагоналей співпадає з початком координат.

Проведено тестові розрахунки при п = 2:

або, з урахуванням (6),

(х - 0,5)4 0,54

У4 -1 (2х - 1)4

' Ь4 -1 ^ 1 '

.її -1 Ь4 1

(9)

(10)

де Ve — дотична компоненти швидкості, К(0;0О) — ядро інтегрального рівняння, яке обчислюється за формулою:

де Ь — звуження лопатки. В такому випадку залежність (6) набуває вигляду:

к(9.9 ) = (у(9)-У(9оЖ(9о)-(х(9)-х(9о))у/(9о) (2)

К(9’9(і) (у(9) - у(9о))2 + (х(9) - х(9о))2 ’ (2)

де х(0); х(0о)'у(0); у(0о) — відповідні координати точок профілю; Ф (0) — потенціал набігаю чого потоку:

Ф(0о) = (x(eo)cos а + y(eo)sin а)^*,, (3)

V* — швидкість незбуреного газового потоку на достатній відстані від профілю, а — кут атаки профілю.

При чисельній реалізації вказаної моделі рівняння (1) з використанням деякого методу чисельного інтегрування [2] зводиться до схеми лінійних алгебраїчних рівнянь відносно невідомих У0(0і) — значень дотичної компоненти швидкості в точках розрахункової сітки 0І по профілю. При цьому залежності (2) набувають вигляду:

К(6ь Є і) =

_(Уі - Уі)Хі- (хі - хі)Уі

(Уі - У і)2 + (Хі - Х)2

(4)

якщо точки 0І та 0j не співпадають, та

о • о ) — 1 х іу і -у іхі (5)

К(9);9)) 2(хі)2+(у))2’ (5)

якщо 9І = 0 j.

Проведено розрахунки для класичних профілів (еліпс, коло), з відомими аеродинамічними параметрами, при цьому для функцій х(0) та у(0) використовується подання у вигляді:

|х(0) = 2 + 2-СОЗ 0 ;

1у(Є) = f(0),

(6)

де f(0) — функція, що задає форму профілю. В більшості випадків координати задаються у вигляді співвідношень: у = ^х), тому при знаходженні відповідних похідних в (4) та (5) необхідно використовувати формулу для диференціювання складеної функції:

dy _ dy dx _ df dx

аё- аХае_ зхае •

(7)

Розглядаються наступні моделі процесу деформування профілю:

І) В процесі експлуатації плавні границі лопаток затуплюються, форма перерізу лопатки стає схожою на

(1 1

х = ^ + 2"Соз 0;

[у = ±Ц/1 - (2х -1)4 = ±Ь 4^—еОІ40.

(11)

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

Використання (11) дозволяє одержати прості явні формули для обчислення похідних в (4) та (5).

ІІ) В процесі експлуатації профіль приймає конфігурацію, близьку до крилової, з формуванням загострень на границі зходу потоку з профілю.

В такому випадку для задання функції ^х) або ^0) використовується апарат апроксимації або інтерполяції профілю. При цьому слід зазначити, що залежності ^0) та ^х) повинні бути достатньо гладкими, з неперервними першими та другими похідними. Крім того, необхідно враховувати вплив на результат апроксимації похибки при визначенні координат точок профілю. Для подолання вказаних проблем пропонується наступний підхід до апроксимації деформованих профілів (рис. 1).

рис. 1. Схема апроксимації деформованих профілів лопатей

В околі передньої кромки профілю (зони І та ІІІ) профіль наближається еліптичною конфігурацією:

1 1

х = ^+^соз 0;

у = —>/2ха - х2 ,

0 < х < а,

ІІІ:

1 1

х = ^ + 2Cos 0;

0 < х < К.

У = -^тл/2хК - х2 , ■’ К

(12)

(13)

З

Залежності (12) та (13) диференціюються з урахуванням (7).

Для апроксимації зон ІІ та IV використовується техніка многочленів Ерміта [2], при цьому в точках А та В відповідні граничні умови задаються умовами (12) та (13), а в точці С передбачається, що у(С) = 0 . З урахуванням вказаного, одержується:

ІІ:

ІІ

x = 2 + 2 cos 0;

y = £i(x- 1)(aix2 + bx + c)-11 (X 2(;)2-(i) 1) -

де аІ =

IV:

a < x < 1, І

; b =

І - За

(a -1)3 (a - І)3’

1 1

x = 2 + 2 cos 0;

3a2 - 3a + І (a - І)3 '

y = _£2(x _ 1)(aix2 + bx + c)+K2 (x 2K)2(1) 1) -

де ai = -

K < x < 1, І

b=

1 - 3K

3K2 - 3K +1

коефіцієнт аеродинамічного опору зростає за величиною на порядок незалежно від величини малої піввісі профілю та відповідного кута атаки, причому такий ріст характерний лише для коефіцієнту опору — коефіцієнт підіймальної сили та обертовий момент змінюються мало (табл. 1);

їаблиця 1

Коефіцієнт опору для різних моделей профілю

(14)

(15)

Єї = £2 Модель (11) Еліпс; а = 0 Еліпс; а = 8

0,01 -2,63 0,00128 -0,00604

0,014 -3,68 0,00183 0,00119

0,018 -4,73 0,00234 -0,50527

0,022 -5,79 0,00288 0,00149

0,026 -6,84 0,00340 -0,06419

0,03 -7,89 0,00393 -0,52543

0,034 -8,94 0,00445 -0,20294

(К -1)3’“ (К -1)3’ (К -1)3 •

Перевагою такого подання є можливість одержання аналітичних залежностей для визначення похідних в (14)—(15) по змінній 0 за формулою (7).

Вказані два підходи до побудови геометрії профілю були реалізовані у вигляді відповідних програмних продуктів, проведено розрахунки для різних конфігурацій профілів елементів конструкції ГПА.

3. Висновки

За результатами розрахунків можна зробити наступні висновки:

■ при використанні подання (11) для моделювання деформованого профілю встановлюється, що профіль з затупленими границями має суттєво інші аеродинамічні характеристики, ніж відповідний еліптичний профіль —

■ використання моделі деформування (12)—(15) дозволяє одержати значення аеродинамічних характеристик профілю, які співпадають з відповідними теоретичними результатами, крім того, вказана модель процесу деформування дозволяє уникнути впливу похибок визначення координат точок профілю;

■ використання (12)—(15) фактично задає процедуру згладжування даних, причому розрахунки показують, що результати суттєво залежать від вибору параметрів а; К; Єї та є2 в цих залежностях, які, фактично, можуть виступати в ролі параметрів оптимізації моделі.

Література

1. Флетчер К. Численные методы на основе метода Галёр-кина [Текст] / К. Флетчер. — М. : Мир, 1988. — 352 с.

2. Самарский А. А. Численные методы [Текст] / А. А. Самарский, А. В. Гулин. — М. : Наука, 1989. — 432 с.

3. Федоренко Р. П. Приближенное решение задач оптимального управления [Текст] / Р. П. Федоренко. — М. : Наука, 1978. — 488 с.

Розглядаються недоліки сучасної теорії повітря-реактивних двигунів, які приводять до заглушання двигунів при зльоті, польоті та посадці, до знижування безпеки польотів.

Ключові слова: кінематичний аналіз, тяга турбореактивного двигуна.

□----------------------------------------□

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

Ключевые слова: кинематический анализ, тяга турбореактивного двигателя.

□----------------------------------------□

Are looked through the drawbacks of modern theory of air-jet engines, which come to engine stoping while takeoff, cruise and takedown, to the lowering of the fly security.

Keywords: kinematical analyse, turbojet engine thrust.

-------------------□ □------------------------

УДК 629.7.036.001

ОСНОВЫ ЕДИНОЙ ТЕОРИИ ДВИЖИТЕЛЕЙ НА НЕПРЕРЫВНЫХ ПОТОКАХ. ВЫВОД ФОРМУЛЫ ТЯГИ, ПОЛЕТНОГО (ТЯГОВОГО) КПД ТУРБОРЕАКТИВНЫХ

ДВИГАТЕЛЕЙ

Б. Ш. Мамедов

Кандидат технических наук, доцент Запорожский национальный технический университет ул. Жуковского, 64, г. Запорожье, Украина, 69063

Е

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