9845
.pdfBp = 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 |
|
|
|
|
|
Зная скорость VCР потока вещества и объемную плотность ρ вещества, |
||||||
можно определить массовый расход: |
|
|||||
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 |
|
|
8π |
|
|
|
|
|
¶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+ |
4π |
[H[VH]]- |
4π [HrotH.] |
||
|
2 |
2 |
2 |
|
1 |
|
νm |
|
здесь ν m = c2 / 4πγ - коэффициент магнитной вязкости, u - удельная внутренняя энергия среды, К - коэффициент теплопроводности, η и ζ - коэффициенты пер-
113
вой и второй вязкостей МРЖ, T − абсолютная температура. |
|
|
|||||||||
|
В уравнении (7.34) напряженность магнитного поля Н изменяется по гар- |
||||||||||
моническому |
|
закону |
и |
направлена по оси z: H = x |
H |
e−ikz , тогда |
|||||
|
|
|
|
|
|
|
|
|
0 |
0 |
|
|
¶ |
2 |
|
|
2 H |
e−ikz и |
|
|
= -k 2 H 2e−i 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) |
|
|||||||||||
|
|
4π |
4π |
||||||||
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) |
|
4π |
2π |
ω 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 |
. |
|
|
|
|
|
|
|
|
|
||
|
4ηl |
|
|
|
|
|
|
|
|
|
|
|
|||||||||
3 |
|
3 |
|
9η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=2p×25 рад/сек и ко-
эффициенте поглощения альфеновской волны β = 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
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