Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Повх И.Л. Аэродинамический эксперимент в машиностроении

.pdf
Скачиваний:
31
Добавлен:
25.10.2023
Размер:
17.67 Mб
Скачать

Легко видеть, что при со = О или п = О формула (6.22) дает решение задачи об обтекании круглого цилиндра однородным потоком.

На рис. 6.27 приведены кривые распределения скоростей по поверхности цилиндра при различных значениях степени неод­

нородности набегающего потока (я =

0; 0,5 и 1,0). Видно, что раз­

личия в обтекании весьма существенны. При я =

1,0 даже

 

кри­

тических

точек становится в два раза больше, чем

при я

— 0.

 

 

 

 

 

 

Несимметричность

обте-.

 

 

 

 

 

 

кания существенно искажает

 

 

 

 

 

 

показания

цилиндрических

 

 

 

 

 

 

зондов (как для величины

 

 

 

 

 

 

угла, так и для величины

 

 

 

 

 

 

скорости). Изменение вели­

 

 

 

 

 

 

чины поправки Дф для изме­

 

 

 

 

 

 

ренного зондом угла, подсчи­

 

 

 

 

 

 

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

путем

 

 

 

 

 

 

в зависимости

от степени не­

 

 

 

 

 

 

однородности,

показано

на

 

 

 

 

 

 

рис. 6.28, а. На этом же

ри­

 

 

 

 

 

 

сунке приведена кривая, по­

 

 

 

 

 

 

лученная

экспериментально.

 

 

 

 

 

 

Видно, что поправка,

опре­

 

 

 

 

 

 

деленная

экспериментально,

Р и с . 6 .2 7 .

Р а с п р е д е л е н и е

ск о р о ст ей

п о

меньше, чем

теоретическая,

п о в ер х н о ст и

ц и л и н д р а п р и

р а зл и ч н о й н е ­

и что величина поправки мо­

о д н о р о д н о ст и

н а б ега ю щ его

п оток а :

 

жет быть

легко

определена

/ — дг =

0 ; 2

п ~ 0 ,5 ;

3

п = 1 ,0

 

аналитически

по формуле

 

 

 

 

Дф =

80

я.

 

 

 

 

 

Для получения правильного значения угла поправка должна прибавляться к величине угла, измеренного зондом, так, чтобы вектор скорости отклонялся в направлении положительного гра­ диента скорости.

Неоднородность потока существенно скажется и на величине скорости, измеренной зондом. Если истинное значение скорости обозначить через V, а ее значение по показанию зонда через Ѵ3, то величина относительной погрешности

АѴ

Ѵ — Ѵ3

V ~

V

откуда

 

Зависимость этой погрешности от степени неоднородности, определенная по предыдущим формулам и из эксперимента, пока­ зана на рис. 6.28, б. Из приведенных кривых видно, что погреш-

198

ности при измерении угла и величины скорости довольно велики. Так, при п — 0,3 поправка на угол Дф достигает 13°, а ошибка при определении величины скорости превышает 20%.

Р и с . 6 .2 8 . П о п р а в к и на у г о л (а) и в ел и ч и н у ск о р о ст и

(б) п р и и зм ер ен и и ц и л и н д р и ч еск и м зо н д о м в н е о д н о р о д ­ ном п о то к е

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

делить

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

и известным диаметру зонда

лу

 

и величине скорости нахо­

 

 

дится

степень

неоднородно­

 

 

сти и

затем

по

графикам

 

 

(рис. 6.28) определяются со­

 

 

ответствующие

поправки.

 

 

Для шарового

зонда при

 

 

наличии градиента

скорости

 

 

в плоскости, нормальной оси

 

 

державки, величина попра­

 

 

вок меньше, чем для цилин­

Р и с . 6 .2 9 . П о п р а в к и п о к а за н и й ш а р о в о го

дрического зонда. На рис. 6.29

зо н д а д л я

вел и ч и н ы ск о р о ст и и у г л а в з а ­

приведены кривые изменения

ви си м ости

от п а р а м ет р а н е о д н о р о д н о ст и

значений поправок

для с к о ­

 

 

ростеи и угла в зависимости от степени неоднородности, получен­ ные экспериментально для шарового зонда в лаборатории аэро­ динамики ЛПИ им. Калинина. Видно, что при п = 0,3 поправка

199

для угла Дф составляет 10°, а для величины скорости — около 2%. Аналитически поправку для угла можно определить по фор­

муле Дф

34,8 п.

 

 

 

 

 

 

 

 

 

 

 

 

 

Существует ряд теоретических работ по вопросу об обтекании

цилиндра и шара вихревым потоком. На рис.

6.30 приведены кри­

 

 

 

 

 

 

 

 

 

вые величины эффективного сме­

 

 

 

 

 

 

 

 

 

щения для цилиндра [205] и

 

 

 

 

 

 

 

 

 

сферы

[202 ]

при

постоянном

 

 

 

 

 

 

 

 

 

значении вихря, т. е. с линей­

 

 

 

 

 

 

 

 

 

ным распределением

 

скоростей,

 

 

 

 

 

 

 

 

 

и

цилиндра,

расположенного

 

 

 

 

 

 

 

 

 

в

турбулентном

пограничном

 

 

 

 

 

 

 

 

 

слое

[87 ]

(диаметр

 

цилиндра

 

 

 

 

 

 

 

 

 

при этом должен быть во много

 

 

 

 

 

 

 

 

 

раз меньше толщины слоя).

Там

 

 

 

 

 

 

 

 

 

же приведены результаты, по­

Р и с .

6 .3 0 .

Э ф ф ек т и в н о е см ещ ен и е в з а ­

лученные для

насадка по фор­

в и си м ости

от п а р а м ет р а

н е о д н о р о д н о ­

муле (6.21).

 

 

 

 

 

 

сти

д л я

ц и л и н д р а в п о гр а н и ч н о м с л о е

 

 

 

 

 

 

 

Из кривых

видно,

что вели­

( /) с п остоян н ы м

в и х р ем

(2)

и (3 ) — по

 

чина

Д для насадка

[формула

ф о р м у л е

 

(6 .2 1 ),

д л я ш а р а

с

п о с т о я н ­

 

 

 

ным

в и х р ем

(4)

 

 

(6.21) 1значительно больше, чем

 

 

 

 

 

 

 

 

 

для цилиндра,

обтекаемого

по­

перечным

потоком,

и

 

близка

к соответствующей

величине для

сферы.

 

Последнее позволяет утверждать, что величина

эффектив­

ного смещения существенно зависит от того,

плоский поток или

пространственный.

 

 

зонд

Kg

 

 

 

 

 

 

 

 

 

Цилиндрический

 

 

 

 

 

 

 

 

 

 

предназначен для измерений

 

 

 

 

 

 

 

 

 

в плоских потоках. Если по­

 

 

 

 

 

 

 

 

 

ток не плоский, показания ци­

 

 

 

 

 

 

 

 

 

линдрического зонда не соот­

 

 

 

 

 

 

 

 

 

ветствуют

истинным

значе­

 

 

 

 

 

 

 

 

 

ниям скоростей, и измерение

 

 

 

 

 

 

 

 

 

таких полей следует произво­

 

 

 

 

 

 

 

 

 

дить шаровым зондом. При­

 

 

 

 

 

 

 

 

 

чиной

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

 

пока­

 

 

 

 

 

 

 

 

 

заний

цилиндрического

зон­

 

 

 

 

 

 

 

 

 

да

при

наличии

градиента

 

 

 

 

 

 

 

 

 

скоростей

вдоль

державки

 

 

 

 

 

 

 

 

 

является

известный в гидро­

 

 

 

 

 

 

 

 

 

динамике эффект «устоя мос-.

 

 

 

 

 

 

 

 

 

та». Смысл этого явления за­

 

 

 

 

 

 

 

 

 

ключается в том, что гра­

 

 

 

 

 

 

 

 

 

диент скоростей вызывает гра­

 

 

 

 

 

 

 

 

 

диент давлений, под дейст­

 

 

 

 

 

 

 

 

 

вием которого жидкость на­

 

 

 

 

 

 

 

 

 

чинает перетекать ВДОЛЬ

дер-

р ис.

б ,з ц

К о эф ф и ц и ен т ,

уч и ты в а ю щ и й

жавки.

 

Из-за

наличия

этих

с к о с п о то к а

в до л ь

оси

ц и л и н д р а

200

перетеканий появляются погрешности при проведении замеров цилиндрическим зондом.

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

потоке можно вычислить

по

формуле

у __ р

| /

2 (h 2 h a) 6

где k6 — коэффициент цилиндрического зонда, определенный при

некотором

угле

б (рис. 6.28); (/і2 — h3) б — разность давлений

в центральном

и боковом отверстиях зонда при том же угле.

На рис.

6.31

приведено изменение k6 в зависимости от б в пре­

делах от —35 до +35°.

Очевидно, что составляющую скорости, перпендикулярную

оси державки, можно определить по

формуле

=

V cos б.

 

Относительная погрешность

при

определении скорости V

с изменением углов б в диапазоне ±20° не превышает 2,5%, а при изменении углов от 5 до —5° пренебрежимо мала. Следова­ тельно, отклонения от нормали до ±5° при установке зонда в по­ ток не вносят каких-либо существенных погрешностей.

37. Измерения насадками нестационарных потоков

При исследовании потока через рабочие колеса турбомашин измерение полного давления за работающим колесом обычно про­ изводится неподвижно установленным насадком. При вращении колеса на неподвижный насадок набегает поток с изменяющимся направлением скорости и величиной полного напора.

Считая, что величина давления за колесом меняется слабо, можно задачу о нестационарном обтекании насадка разделить на две более простые: первая — насадок обтекается потоком с по­ стоянными по шагу давлением р и углом набегания ß2 и с меняю­ щейся лишь величиной скорости потока; вторая — давление и ве­ личина скорости постоянны, изменяется лишь направление ско­ рости по шагу.

Решение первой задачи выполнено для двух законов распре­ деления скоростей: парциального и параболического [146]. Если V 2— скорость потока, набегающего на насадок, т. е. скорость за колесом, а Ѵ2Х— ее расходная составляющая, то среднее значение

полного напора за рабочим колесом будет равно

 

t

v\v2x dx

 

I

 

Ро2Cp =~р2 "~f~ ~~

t--------- •

(6.23)

 

j v.,x dx

 

 

О

 

201

Интегрирование производится по шагу от 0 до t — значения х

равного шагу t. Определим Ѵ2 из треугольника скоростей пока­ занного на рис. 6.32,

 

 

У2Ѵ2х +

( и ^2и )2 = Ѵ\х +

—•Ѵ2х Ctg ß2)2 =

 

 

 

 

 

=

v l

(1 + Ctg2 ß2) +

и -

2uV2x ctg ß2,

 

 

где

и — окружная

скорость

рабочего

колеса

на том

радиусе

где

расположен

насадок;

 

W 2и — окружная

составляющая ско­

рости в

относительном движении.

 

 

 

Ѵ2

в

 

 

 

 

Подставив

найденное

выражение для

формулу (6 23)

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

(ctg2ß2 +

l ) \ v l x dx.

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

— 2« ctg ß2 j Vlx dx

 

 

 

 

 

(6.24)

Считая,

что

насадок

 

осредняет полный

напор по импульсам

получим

величину давления,

измеренную

насадком,

 

 

 

 

 

 

 

 

 

 

 

 

і

 

 

 

 

 

 

 

 

 

 

 

 

 

Р02ср = № +

- J Vlxdx = р2+

P f +

 

 

 

+

 

Р

( 1 +

ctg2

 

 

Г

v \xdx --

 

и ctg

 

1

Vxdx

 

(6.25)

 

 

2/

ß2) J

2

ß2 J

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Разность

Ар02 =

р02 кр —

 

 

 

 

 

 

 

 

 

 

 

—Р»2ср характеризует ошибку

 

 

 

 

 

 

 

 

 

 

 

показания прибора из-за на­

 

 

 

 

 

 

 

 

 

 

 

личия нестационарности. Вве­

 

 

 

 

 

 

 

 

 

 

 

дем безразмерную

величину

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А

 

Р“к

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

Рис. 6.32.

Треугольник скоростей

 

где ик — окружная

скорость

Вычитая из (6.25)

уравнение (6.

периферийной точки

колеса,

t) и подставляя безразмерные

значения х = ~ ;

и =

 

и V.

 

Сѵ

получим

 

 

 

= ~ ~ ,

 

 

 

 

 

 

t

 

 

U„

 

 

 

ик

 

 

 

 

 

 

202

1

\

 

Vtx dx

оо

 

I

l v**d~x

\

 

 

-— // Ctgß2|

1

Ѵ2хсіх— - \ ------- .

\

J

1 г » “'*

,

о

/

Если Ѵ2х представить как произведение максимального значе­ ния этой скорости Ѵ2хтах на некоторую функцию / (я), т. е.

V2x = V2xf(x),

и подставить полученное выражение в последнее уравнение, то будем иметь

 

 

(6.26)

где

 

 

j p(x)dx\

/ 2=

j fs(x)dx',

о

 

о

коэффициент

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

о

 

 

расходных составляющих скоростей.

 

В зависимости от вида функции / (х) общая формула (6.26) будет иметь тот или иной вид. Рассмотрим вид формулы (6.26) для парциального поля. Пусть на участке от 0 до а скорость по­ стоянна и равна Ѵ2х, а на участке от а до t скорость равна нулю, тогда на первом участке

f (x) = 1; Ѵ2хср = Ѵ2х.

Все интегралы, входящие в формулу (6.26) равны а и, следо­ вательно,

2

203

или

 

 

Л р о 2 =

 

 

 

=

 

 

 

При а =

1 насадок

показывает действительное полное давле­

ние. При а <

1 и И2Х<

и насадок показывает величину, большую

действительного полного давления. Чем меньше

величина а,

тем больше искажаются показания насадка.

за

колесом

Для

параболического распределения

скоростей

 

 

 

 

/ (х) =

1 — хп' ,

 

 

 

где /4 =

- ^ .

 

 

 

 

 

 

 

 

В этом случае

 

 

 

 

 

 

 

 

 

 

h

п (п + 1)

И

/ о

п ( п +

1 )(я + 2)

 

 

 

 

 

 

 

 

и формула (6.26)

принимает вид

 

 

 

 

Дрог = V.2* ср

ft-- 1

Hctgß2 — 2П2Х~ 2 - ( 1 + Ctg2 ß2)

(6.27)

ft -f“ 1

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

Аналогичные решения получены для синусоидального распре­ деления скоростей в несжимаемой жидкости и парциального рас­ пределения скоростей в сжимаемом газе [189]. Во всех случаях погрешность измерений обращается в нуль только при определен­

ных соотношениях между окружной скоростью и, расходной со­ ставляющей V 2х и углом выхода потока ß2, т. е. только на одном режиме движения колеса. На других режимах измерения непод­ вижным насадком дают значения полного напора, отличные от действительных.

Для экспериментальной проверки погрешностей измерения полного напора были поставлены опыты, при которых в устано­ вившемся однородном потоке аэродинамической трубы произво­ дилось измерение полного напора насадком, совершающим коле­ бательное движение. В этом случае при постоянных значениях давления р 2 и величины скорости, а следовательно, и полного напора изменялась величина угла ß2. Насадок, будучи чувстви­ тельным к углу поворота, отклонялся относительно средних углов на ±12,5 и 21°. Частота колебания насадка характеризо­ валась числом Струхала, составленным следующим образом:

где d — диаметр приемного канала насадка; V 2 — средняя вели­ чина скорости потока; f — частота колебаний.

204

Рис. 6.33. Зависимость показания на­ садка h от числа Sh насадка:
д — показания насадка;------— средние по времени значения мгновенных скоро­ стей

На рис. 6.33 показана зависимость результатов измерения на­ садком h и средних значений мгновенных давлений в нестационар­ ном потоке от числа Sh.

Из рисунка видно, что при числах Sh > 2 0 показания насадка совпадают со средним значением мгновенных давлений, в то время как при числах Sh < 20 они

расходятся.

Эти разультаты получены при малых числах М и для оп­ ределенного типа насадка. Для других типов насадков и боль­ ших чисел М область чисел Sh, в которой имеет место указан­ ное расхождение, может изме­ ниться.

Таким образом, полное дав­ ление, измеренное за рабочими колесами неподвижными насад­ ками, отличается отдействитель­

ного. Измеренная величина полного давления может быть как мень­ ше, так и больше действительной. Если учесть, что при измерении давлений и скоростей неподвижными насадками за рабочим колесом изменяется не только угол, но и величина скорости, то погрешно­ сти будут иметь еще большие значения. Отсюда следует, что изме­ рение параметров потока должно производиться насадками, вра­ щающимися вместе с рабочим колесом [189]. Методика измерений на вращающихся элементах машин изложена в гл. VIII.

38. Инерционные свойства насадков

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

зафиксировано прибором.

В неудачно спроектированных насад­

а)

 

ках это время может оказаться очень

 

 

большим.

 

 

 

может

быть

 

 

Математически задача

 

 

поставлена следующим образом: про­

5)

 

странство,

в котором

производится из­

 

мерение давления, имитируемое беско­

 

 

------- \ —

нечным резервуаром,

соединено

с из­

мерителем

давления

тонкой

трубкой

4

J

(рис. 6.34).

В начальный

момент вре­

Рис. 6.34.

Идеализированная

мени £ = 0 при входе в трубку создается

(а) и типичная (б) схемы из­

конечное давление,

отличное

от посто­

мерения давлений

 

 

янного давления в

трубке и в измери­

теле. Требуется определить изменение давления в трубке и в измери­ теле. Требуется определить изменение давления по длине трубки в зависимости от времени р (t, х) и время запаздывания, т. е. время,

205

в течение которого в измерителе установится давление, возникшее при t — 0 в начале трубки. В зависимости от среднего давления

в

трубке в ней может установиться три типа движения

газа:

1)

непрерывное движение или течение континиума; 2) течение со

скольжением и 3) свободное молекулярное течение.

 

 

В случае течения со

скольжением время запаздывания

[58 ]

будет существенно лишь

тогда, когда величина

 

0,02 < Я = ——- < 1,

где рср — среднее давление в мм рт. ст.; г — радиус трубки в см. Следовательно, при течении со скольжением среднее давление

имеет значения, ограниченные неравенствами

4,1 .

 

410

— <

Рср <

2л -

 

 

Такие давления могут встретиться при аэродинамическом экс­ перименте.

При значениях среднего давления в мм рт. ст., меньших об­

ратной величины диаметра трубки в мм, т. е. рср < — ■, будет

иметь место свободномолекулярный поток. В обычном аэродина­ мическом эксперименте свободомолекулярных течений не бывает.

В случае, когда длина свободного пробега молекул имеет по­ рядок внутренних размеров насадка, т. е. имеет место течение со скольжением, время установления давления в измерителе (или время запаздывания) определяется по формуле [113]

 

Q +

ял2/

4/2

(6.28)

 

 

- +

ял

где С = ~ ~ - у ~ = г (для

воздуха при t = 20° С С = 10"5

с/см);

Q — объем измерителя;

I — длина

трубки;

г0 — радиус

отвер­

стия.

 

 

 

 

Как видно из формулы (6.28), увеличение радиуса трубки г вызывает одновременно и падение, и рост времени запаздывания; первое происходит из-за уменьшения сопротивления трубки, а второе — из-за увеличения объема воздуха. Таким образом, для каждой установки существует радиус трубки, при котором время запаздывания будет минимальным.

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

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

206

нений. Считая, что движение ламинарное, можно определить мас­ совый расход газа по формуле

т = п

рг4

др

(6.29)

8

дх ’

 

где X — координата вдоль трубки

(рис.

6.34, а).

Уравнение неразрывности для одномерного нестационарного движения имеет вид

Ф I д(ри) _ п dt ^ дх

Допуская, что изменение состояния происходит по изотерми­ ческому закону, т. е. р = kp, и используя формулу (6.29),полу­ чим уравнение одномерного движения квазистационарного тече­ ния вязкого сжимаемого газа в виде

<6 ' 3 0 >

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

Граничные условия определяются из следующих физических условий: в точке, в которой производится измерение, давление постоянно, т. е. при х = I р = const.

В конце трубки, т. е. при л: = О, граничное условие полу­ чается из соображения равенства изменения массы газа в камере

измерительного прибора массовому

расходу через трубку, т. е.

др

я г4

др

d t

8 Q p Р

д х

Начальные условия будут

р (0, х) = р х и р (О, I) = р 2.

Уравнение (6.30) было проинтегрировано [175] численно для трех случаев. Во всех случаях начальное давление в измеритель­ ном приборе и в трубке бралось равным р = 10107 кгс/м2 (101070 Па) и длина трубки I = 61 см, давления в баке, соответ­ ствующие давлению в отверстии насадка, задавались равными

р2 = 272 кгс/м2 (2720 Па) и р 2 = 6904 кгс/м2 (69040 Па).

Результаты теоретического расчета и результаты эксперимента

сидеализированной схемой (рис. 6.34) для трубок с диаметрами 0,64 и 1,60 мм при р 2 — 272 кгс/м2 (2720 Па) и одинаковыми осталь­ ными параметрами показаны на рис. 6.35. Видно, во-первых, что теоретические данные блестяще совпадают с эксперименталь­ ными и, во-вторых, что время для создания измеряемого давления

визмерительной камере устанавливается асимптотически.

207

Соседние файлы в папке книги из ГПНТБ