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

9845

.pdf
Скачиваний:
0
Добавлен:
25.11.2023
Размер:
3.35 Mб
Скачать

Bp = Bmax cos αe cos ωt =

 

 

 

 

=

Bmax

cos(α

 

ωt )+

Bmax

cos(α

 

+ ωt ) = B

+ B

(7.26)

 

e

 

e

 

 

2

 

2

 

l +

l

 

 

 

 

 

 

 

 

 

т.е. волнами с половинной амплитудой и вращающимися в противоположных направлениях - прямом (Bl + ) и обратном (Bl ).

В выражении (7.26) αe = z pα – угловая координата точки на периоде пульсирующей волны Ta ; z p – число пар полюсов; α = x R − угловая коорди-

ната точки на окружности радиуса R , выраженная в дуговых единицах (радианах); x − координата точки в линейных единицах.

Период пульсирующей волны Ta соответствует полному циклу электро-

магнитных процессов во временной области, поэтому угловая мера αe = z pα

соответствует угловой мере фазового угла.

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

Ta 2 друг от друга (рис. 7.8), где период пульсирующей волны Ta пропорцио-

нален длине коаксиального дроссельного канала.

Рис. 7.8. Пульсирующая волна магнитной индукции в коаксиальном цилиндрическом

дроссельном канале 8 с МРЖ

110

Здесь и в выражении (7.26) неподвижный пульсирующий поток заменяется суммой идентичных круговых полей, вращающихся в противоположных направлениях и имеющих одинаковые частоты вращения: n1пр = n1обр = n1, где n1пр

– скорость поля прямой последовательности, n1обр – скорость поля обратной последовательности.

Электромагнитные моменты Мпр и Мобр , образуемые прямым и обрат-

ным магнитными полями внешнего однофазного индуктора МРТ, направлены в

противоположные стороны, а его результирующий момент М рез равен ал-

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

Скольжение МРЖ относительно прямого потока магнитного поля Фпр

sпр

= (n1пр n2 ) n1пр = (n1 n2 ) n1 = 1− n2 n1 ,

(7.27)

где n2

скорость вращения МРЖ в вихревом магнитном поле.

 

Скольжение МРЖ относительно обратного потока магнитного поля Фобр

sобр

= (n1обр + n2 ) n1обр = (n1

+ n2 ) n1 = 1 + n2 n1 .

(7.28)

Из (7.27) и (7.28) следует,

что

 

 

sобр = 1 + n2 n1 = 2 − sпр .

(7.29)

Скольжение МРЖ в коаксиальном цилиндрическом дроссельном канале МРТ определяет величину снижения частоты вращения МРЖ относительно частоты вращения вихревого магнитного поля внешнего однофазного индукто-

ра. У стенок прямолинейного проходного каналаVж = 0 , а в его центре она мак-

симальнаV = maks . С расстоянием

r

 

 

от дроссельного канала

скорость

ж

 

var

 

 

Vж изменяется по параболическому закону (рис. 7.9):

 

 

 

r 2

 

 

 

 

 

 

 

 

 

(7.30)

 

 

2

 

Vж (r) = Vж0 1−

R

,

 

 

 

 

 

 

 

где R – радиус проходного канала;

Vж0 скорость рабочей жидкости на оси

 

 

111

 

 

 

проходного канала, определенная формулой:

 

=

P − P

=

p

 

 

 

V

1 2

R 2

 

R 2

.

(7.31)

 

 

ж0

 

4η l

 

4η l

 

 

 

 

 

 

 

Здесь P1 − P2 = p − разность давлений на концах проходного канала; l - длина проходного канала; η − динамическая вязкость рабочей жидкости.

Рис.7.9. Движение отдельных слоев жидкости с различными скоростями при ламинарном течении рабочей жидкости в проходном дроссельном канале.

Поток рабочей жидкости

Q, т. е. объем рабочей жидкости, протекающий

через поперечное сечение проходного канала рабочей длиной

l = 500 мм за

единицу времени определяется формулой Пуазейля [16]:

 

Q =

(P1 − P2 )

π R 4 .

(7.32)

 

 

8η l

 

 

 

 

Зная скорость Vпотока вещества и объемную плотность ρ вещества,

можно определить массовый расход:

 

Q кг

= S

К

ρ V

(7.33)

 

 

 

СР ,

 

с

 

 

 

 

 

где Sк - площадь поперечного сечения потока на измерительном участке.

При составлении уравнений движения магнитореологических сред в магнитных полях гидроопоры предполагаются следующие условия [5]. Магнитная проницаемость магнитной среды однородна и изотропна во всем объеме действия и не зависит от напряженности магнитного поля Н. Это условие имеет место при ω0τ <<1, где ω0 ларморова частота прецессии для полярных частиц

112

МРЖ, τ − среднее время свободного пробега частицы, электропроводность -γ

достаточно мала, так как нет свободных носителей заряда.

Система уравнений, описывающих движение МРЖ в магнитном поле дроссельных каналов, включает в себя уравнения Максвелла:

 

 

 

H = rot[VH] +ν m DH, divH = 0,

 

 

 

t

где ν m

=

c2

- коэффициент магнитной вязкости, который тем больше, чем ниже

4πγ

 

 

 

электропроводность среды.

Вгидродинамические уравнения движения среды входят:

-обобщенное уравнение Навье-Стокса:

V

+ (VÑ)V = -

1

 

1

[HrotH]+

η

 

1

 

ξ +

η

 

 

gradp -

 

 

DV +

 

 

graddivV ;

t

ρ

4πρ

ρ

 

 

 

 

 

ρ

 

3

- уравнение неразрывности: ρ + div( ρV) = 0;

t

- уравнение состояния среды: p = p(ρ, T) ;

 

 

ρV

2

+ ρu +

H

2

 

 

 

- уравнение закона сохранения энергии:

 

 

 

 

= -divW , где

ρ − плот-

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

t

 

 

 

 

 

ность среды, V − скорость движения среды в дроссельных каналах, u внутренняя энергия, H напряженность внешнего магнитного поля, W - плотность потока энергии (вектор Умова-Пойнтинга).

Диссипация плотности потока энергии W резко возрастает, если в потоке МРЖ нарушается химическое равновесие. С учетом параметров реологического заполнителя и внешнего магнитного поля плотность потока энергии можно представить в виде [6]:

W= ρ×V u

 

 

 

2

 

 

 

 

2

 

 

2

 

2

 

 

 

 

Vy

2

 

 

2

 

 

p

 

V

 

 

V

 

Vy

 

V

 

V

 

V

V

 

+

 

+

 

 

-KÑT+η 2

x

 

 

+

 

 

+

z

 

+

 

x

+

 

 

+

x

+

z

 

+

ρ

2

 

 

y

 

 

x

 

 

 

 

 

 

 

x

 

 

 

 

z

 

 

y

 

 

z

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(7.34)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Vy

 

V 2

 

 

+

z

 

 

z

 

 

 

 

y

]-

3(divV)

++ζ (divV)

}V+

[H[VH]]-

[HrotH.]

 

2

2

2

 

1

 

νm

 

здесь ν m = c2 / 4πγ - коэффициент магнитной вязкости, u - удельная внутренняя энергия среды, К - коэффициент теплопроводности, η и ζ - коэффициенты пер-

113

вой и второй вязкостей МРЖ, T − абсолютная температура.

 

 

 

В уравнении (7.34) напряженность магнитного поля Н изменяется по гар-

моническому

 

закону

и

направлена по оси z: H = x

H

eikz , тогда

 

 

 

 

 

 

 

 

 

0

0

 

 

2

 

 

2 H

eikz и

 

 

= -k 2 H 2ei 2kz .

 

 

rotH =

 

H = -k

HrotH

 

 

 

z

 

 

0

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

Будем рассматривать движение МРЖ в коаксиальном цилиндрическом зазоре МРТ при div = 0 , так как объём рабочей среды в коаксиальном дроссельном канале не изменяется.

Пусть напряженность магнитного поля ортогональна направлению скоро-

сти в коаксиальном зазоре H ^ v , тогда [vH]= v0 H .

С учетом этих допущений перепишем уравнение (7.34)

W

 

=ρ × v (u +

Dp +

v02

) - kÑT +

1

H 2 v

+

ν m

H 2 .

(7.35)

 

 

 

x=0

 

0

ρ 2

 

0 0

 

0

 

Коэффициент поглощения β альфеновской волны (т.е. электромагнитной волны, направление распространения которой, понимаемое как направление ее групповой скорости, совпадает с направлением поля H ) можно найти из формулы [8]:

β =

ω 2

(

η

+ν

 

)

(7.36)

2U 3

ρ

m

 

 

 

 

 

 

A

 

 

 

 

 

 

где ω - частота магнитного поля, UA – скорость распространения альфеновской волны,

UA =

 

H

 

(7.37)

 

 

 

4πρ

 

 

 

 

из уравнений (7.36) и (7.37) выразим магнитную вязкость νм :

ν

 

=

β 2H 3

-

η .

 

 

 

m

 

(4πρ )

3 2

 

ρ

Подставляя (7.38) в (7.35) получим:

 

= ρ ×v (u +

Dp

 

v2

W

 

+

0

) - kÑT +

ρ

2

z =0

0

 

 

(7.38)

1

H 2v +

1

β

H05

-

η

H 2

(7.39)

ω 2 (4πρ )3 2

4πρ

0 0

 

 

0

 

Из выражения (7.39) следует, что плотность потока энергии зависит от Н

114

как функция пятой степени [8]. Рассмотрим конкретный пример.

Пусть: ν 0 = 0,5м/с; u = 30 дж;

k =

 

 

T2

 

 

o

 

а = 0,25×10

-3

7,0;

kÑ=

t

=2TdT;

0 £T £

70

;

м;

H0 = 200×103sin(wt+φ(t)); β = 0,1.

Чтобы упростить выражение (7.39), будем считать, что температура жидкости в канале постоянна (т.е. ÑT = 0 ), так же примем, что в первом приближении вязкость жидкости не зависит ни от напряженности магнитного поля ни от скорости течения (η = const ).

В качестве рабочей жидкости возьмем магнитореологический заполни-

тель на основе глицерина при 200С. Тогда

 

r=1,26×103

кг/м3,

h=1480×10-3

кг/м×сек. Пусть r0 =

1

a ; тогда v0

=

a 2 Dp

(1 -

1

)2

=

a 2 Dp

.

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

3

 

3

 

l

 

 

 

 

 

 

 

 

 

 

 

Перепад давления в гидроопоре Dp = p - p

 

 

(рис. 7.8)

p

 

=

F

а p

 

=

F

,

2

 

1

 

2

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

S1

 

 

Sкан

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где F - сила, действующая на гидроопору, Sкан =12,5×10-6м2

площадь попереч-

ного сечения коаксиального канала, S1 =36,7×10-3м2 – площадь нижнего сечения рабочей камеры гидроопоры. Пусть сила воздействия на гидроопору изменяется по гармоническому закону, с частотой, равной частоте магнитного поля F = F0 sin(ωt) . Таким образом:

Dp =

F0

sin(ωt) -

F0

sin(ωt) .

(7.40)

 

 

 

S1

 

Sкан

 

Сила воздействия F0=1000 Н. Учитывая все сказанное выше, построим зависимость плотности потока энергии в коаксиальном дроссельном канале гид-

роопоры с МРЖ W(H, t) x=0 при частоте магнитного поля w=225 рад/сек и ко-

эффициенте поглощения альфеновской волны β = 0,1 (рис. 7.10).

115

Рис. 7.10. Зависимость плотности потока энергии магнитного поля от времени и давлении при f1 = 30 Гц и Т1 = 20° С .

Рис. 7.11. Зависимость плотности потока энергии магнитного поля от времени и давлении при f1 = 30 Гц и Т1 = 50° С.

Вывод.

Из рис. 7.10 и 7.11 следует, что плотность потока энергии зависит от напряжённости магнитного поля как функция пятой степени, что так же следует и из уравнения (7.39).

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

116

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

117

ГЛАВА 8

ВЕРОЯТНОСТНЫЕ И СТАТИСТИЧЕСКИЕ МОДЕЛИ ВИБРАЦИОННЫХ ПРОЦЕССОВ

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

8.1. Случайные события и величины

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

118

£ a < b £ +¥ ,

pi ≈ νi ;

νi = ni / n;

n = ni ,

 

 

i

где ni - число i - го события в серии из n наблюдений. Множество всех возможных событий в опыте объединяется понятием случайной величины. Слу-

чайная величина Х есть числовая величина, которая при наблюдении прини-

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

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

X = {WX ; PX (x)} .

Здесь WX = { x} - множество возможных значений х, а PX (x) - закон рас-

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

Случайные величины могут быть дискретными или непрерывными, но колебательные и вибрационные явления, рассматриваемые ниже, являются непрерывными, и поэтому рассматриваем только непрерывные величины. Для их описания необходимо задать: WX = (a, b),

FX (x) = P( X < x) и f X (x) = dFX ³ 0 . dx

Здесь а, b - границы односвязной конечной или бесконечной непрерывной области, FX(x) - функция распределения случайной величины Х, определяющая вероятность того, что случайная величина принимает значения меньшие, чем ее аргумент х, а fX(x) - функция плотности распределения вероятности случайной величины. Каждая из этих функций и задает закон распределения случайной величины, для них выполняются следующие формулы:

119

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]