О движении твёрдого тела с криогенной жидкостью Текст научной статьи по специальности «Физика»

НАУКА и ОБРАЗОВАН И Е

Эл № ФС77 - 48211. Государственная регистрация №0421200025. ISSN 1994-0408

О движении твёрдого тела с криогенной жидкостью # 12, декабрь 2013 DOI: 10.7463/1213.0627898 Ай Мин Вин , Темнов А. Н.

УДК 531.38

Россия, МГТУ им. Н.Э. Баумана ayeminwin84@gmail.com antt45@mail.ru

Введение. Движение твёрдого тела с полостью наполненной однородной несжимаемой жидкостью и движений самой жидкости в подвижном твёрдом теле исследовалось такими известными учёными как Жуковский Н.Е. [1], Слудский Ф.А. [2], Хаф [3], Пуанкаре [4]. Особый интерес к подобной задаче возник в 50-ом годы прошлого столетия, в связи с развитием ракетно-космической техники, о чем свидетельствуют работы [5], [6], [7], [8], с полостями наполненными однородной жидкостью.

Широкое исследование криогенных жидкостей в ракетно-космическом технике стимулировало появлении новых работ по движении твёрдых тел с полостями заполненными неоднородной несжимаемой жидкостью [9], [10]. В недавней работе [11] рассмотрены вопросы взаимодействия неоднородной (стратифицированной) жидкости с твёрдым телом, совершающим малые периодические движения. В данной статье приведены результаты исследование инерционных характеристик твёрдого тела с неоднородной жидкостью полностью или частично заполняющей цилиндрическую полость и совершающей квазипотенциальное движение.

Основная часть статьи состоит из 4-ох пунктов. В 1-ом пункте приведена постановка задачи о влиянии воздействия жидкости на твёрдое тело, в следующем показан физический смысл квазипотенциала скоростей, в 3-ем приведен анализ движения твёрдого тела, а в последнем пункте решена задача об определении квазипотенциала скоростей неоднородной жидкости с заданным законом изменения плотности и приведены результаты расчётов инерционных характеристики.

1. Постановка задачи. Рассмотрим движение твердого тела с полостью произвольной конфигурации, полностью заполненной неоднородной несжимаемой жидкостью.

Введем подвижную систему Ох1х2х3, жестко связанную с твердым телом, с началом в т.О, совпадающей с центром масс твердого тела. Определим движение твердого тела при помощи вектора скорости Уо точки О и вектора угловой скорости вращения со . Запишем уравнения движения твердого тела при наличии жидкости в полости [12]

dt

- + сх QT = KT + P, (1)

M + со х Мт = LT + N, (2)

dtTT

где QT и MT - количество движения и кинетический момент твердого тела; KT, P, LT и

N - главные векторы и главные моменты относительно т. О активных сил, приложенных к точкам твердого тела, и сил давлений жидкости.

Выражение для P и N получим, используя уравнения гидродинамики идеальной жидкости, записанные в подвижной системе координат

dV о - 1 ^ о

— + сх V =--Vp + f,

dt р0

где V - вектор абсолютной скорости частиц жидкости

V = VO + сх r + и,

и - вектор относительной скорости частиц жидкости, r - радиус вектор какой либо точки системы относительно точки О, р( r, t ), p ( r, t ) - плотность и давление частиц жидкости,

f - вектор интенсивности поля внешних сил.

P = -\p°^dG-\p0è х VdG, (3)

G d G

N = -^por х^^dG -\po? хё х VdG, (4)

G dt G

здесь G - область, занимаемая жидкостью.

Из уравнений (3), (4) видно, что для определения векторов P и N необходимо знать движение жидкости в полости твердого тела, т.е. решить гидродинамическую задачу

- найти поле скоростей жидкости - V (r, t) .

2. Квазипотенциал скоростей неоднородной жидкости.

Под квазипотенциальным движением будем понимать такое движение

неоднородной жидкости, при котором гогУ ^ 0, но выполняется условие

V ■ гогУ = 0, (5)

где V - скорость частиц жидкости.

Равенство (5) является необходимым и достаточным условием [13] существования скалярной функции Г (г, г), называемой в дальнейшем квазипотенциалом скоростей

неоднородной жидкости. Скорость частиц жидкости тогда может быть выражено формулой

у = —^ , (6) Ро ( г )

где ро (г ) - известное поле плотностей жидкости.

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

Поясним физический смысл квазипотенциала. С этой целью рассмотрим уравнение движения и уравнение неразрывности неоднородной жидкости

-V 1 - йр -

=--Ур + /, = -р ■ V;

ш р Ш

Пусть к жидкости приложены мгновенные силы с интенсивностью /', действующие в

течение очень малого времени Т и порождающие большие градиенты давления р' [14].

Уравнение движения жидкости тогда приобретёт вид

-V 1 „ 7 7* 1 гт '

=--Ур + / + /'--Ур',

ш р р

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

уравнение неразрывности от г = го до г = Т .

Учитывая, что за малый промежутки времени плотность частиц жидкости не успевает изменится, имеем

Ш р

I-^йг = ро(г),У ■У = 0;

-г

'о

Далее проинтегрируем по малому промежутку Т уравнение движения, принимая во внимание малость импульса обычных сил по сравнению с импульсом мгновенных сил http://technomag.bmstu.ru/doc/627898.html 257

К' -V = J -— Уп, Ро

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

х г 1

J = ^ Р й1- импульс мгновенных сил, п = ^ р1' йт - импульс мгновенных давлений.

1о 1о

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

состоянии покоя V = 0 , то при мгновенном изменении этого состояния получим

V =-—Уп, Ро

причем кинематические состояние жидкости будет вихревым

- —' 1

Ух V = П =— Ур0 хУп.

Ро

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

Р (х, 1) = п( х, 1) + С (1),

приводящего к возникновению вихревого движения жидкости.

Сформулируем краевую задачу для нахождения функции Р (г, 1) в случае когда жидкость заполняет полость твёрдого тела, внезапно пришедшего в движение

У- — УР = о в т (7)

Ро

1 дР - - - „ _ „

--= Уо -у+с- (г ху) на а (8)

Ро ду

где v - орт внешней нормали к границе Б области G , У - оператор Гамильтона. Функции Р (х, 1) будем искать в виде [15^18]

Р (г, 1) = Ф-V + ¥ с (9)

где: Ф(Ф1,Ф2,Ф3); ¥(^1,¥2,¥3) - векторные функции, компоненты которых Фг. и ¥.

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

Используя представление (9) и уравнение (7), получим следующие краевые задачи Неймана для функций Фг. и

1 1 дФ ^ -

V-(—УФг) = 0 —= 1/ • /; (10)

Ро Ро д1

1 1 -

V• (—= 0 ——^ = (г XV)•/.; (11)

Ро Ро д1

Сделаем несколько предположений относительно характера величин, входящих в задачи (10), (11). Будем считать область односвязной, областью с достаточной гладкой границей Б. Из физических соображений функций ро (г ) подчиним условиям:

0 < т < Ро(г) <М <сю (12)

3. Анализ уравнений движения твердого тела с неоднородной жидкостью.

Считая функции Ф. и найденными, составим уравнения движения твердого

тела с неоднородной жидкостью, совершающей квазипотенциальное движение. Удобной формой уравнений движения будет являться такая форма, в которой инерционные характеристики системы могут быть рассмотрены в виде суммы инерционных характеристик твердого тела, «затвердевшей» жидкости и добавки, происходящей за счет подвижности жидкости относительно твердого тела [6]. Представим с этой целью вектор Ф в виде:

Ф = Ф"+р0 (г )г , (13)

и заменим в формулах (3),(4) вектор скорости V формулой (6), где функция ^ определена формулой (9), в которой функция Ф представляется формулой (13). Используя формулу Гаусса-Остроградского, граничные условия (10), (11) и формулу Грина

Г (иЬ \и ]- иЬ \иЫв = <$ — (и — - и —,

Г М 1 ] ? Ро д1 д1 '

( Ь - линейный оператор), после громоздких, но несложных преобразований:

^(ф •У )lт = fv(ф'•VQ У? + \ Г •У, Vpo йт + У0 |Р0 йт

(14)

т

£ Vo; £ъфФ;+ £У;£е1 ГX} Рйт + тУ0

]=1 к=1 ж Р0 д1 ] =1 к=1 т дХк

т

= V0 (JФ; — I ei дФ- + e2 + e3 ■дФз]ds + V0 (JФ] — | ei ^^ + e2 + e3 ^^Ads + 1 s P0 V дv дv дv J 2 s P0 V дv дv дv J

+V0 ф Ф; -1Í ¡ дФ- + ¡2 дф2 + ¡3 ds + V0 J xL f ¡ P + ¡2 P + ¡3

3 ; P0 V öv дv дv J L т V SxL дх2 Sx3

+V02 J

Л

р + ¡J p + - dp

2 i ^ 2 ^ 3

т V дх; дх2 дхз J

dT + Къ J X

'i P + e; P + Í P

i ^ 2 ^ 3

T V дх; дХ2 дх3 J

dT + mVr, =

i дф

Po Sv

i дФ2

Po Sv

i дФ3

i дФ, j г.. i дФ .

- ds ф Ф3--L ds

3

§ ф;

(J Ф;--- ds (J Ф]--- ds (J Ф

J r> Я1/ J Г) Я1/ J

¿ ^ i дФ3 . f .. i дФ3 , f _, Ф Ф;--3 ds ф Ф2--3 ds ф Ф

j r> Я1/ J Г) Я1/ J

Po Sv

i дФ2

Po Sv

i дФ3

Po Sv

i дФ2

Po Sv

i дФ3

+

Po Sv S Po Sv SPo f „ SPo f.. SPo

Sx

Po Sv

ds ds

f x^ dT f x2P^ dT f J Av* J Av* j

Sx, SP

SxL SPo

dT

J ^ dT J x22 ^ dT J

P dT

l

0

Sx J 2 Sx J 3 Sx,

Sx.

SP<

Sx2 SPo

J x^ dT J x2P^ dT J

Sx2 SPo

dT

V0i V02 V

+ mV0 =

3 3

= ZZ e, e;ßß V + mVo = M • Vo + mVo ;

j=i k=i

3 3

гДе M = ZZ ek ej Vjk, m = J PodT;

j=i k=i

J r X v(Ф • V0 ) dT = J r X v(^ • V0 ) dT + J r xVp0 (r -V0 ) dT + J r xV0p0 dT

Z Vo j Z ek fy1- ^тФ] ds + Z Voj Z ek J r, ( r x Vpo )k dT + L x V

j=i k=i

Po Sv

j=i k=i

= V0 JФ] L f ¡S^ + e2^ + ¡зds + V0 JФ2 -1 f ¡S^ + e2^ + ¡3 +

L S P0 V Sv Sv Sv J 2 S P0 V Sv Sv Sv J

+V03 J Ф]

, i Í — S^ — дщ2 — S^31 — e —— + e —— + e —— \ds +

i -л 2 ^ 3 -

p0 V Sv Sv Sv

+V0i J I x2 SP ~ x3 SP l( xi ei + x2 e2 + x3 e3 )dT+ V02 J

x Spo - x Spo ]( xi e; + x2 e2 + x3 e3 )dT +

SxL 1 Sx3

+V03 J! xi tr " x2 "P^](xi eei + x2 e2 + x3 e3 )dT+ L0 x V0

I

I

ф Ф* — -d^L ds ф Ф* — -d^L ds ф Фз*

* /л Я1/ J г\ г) л/ *

1 d^1

Po dv

1 d^2

Po dv

1 d^3

Ф Ф1вТГ ^ ds í

ф ф; -1 ËlL ds ф J Po dv J

Ф*

Ф*

Po dv

1 d¥2

Po dv

1 d^3

Po dv

ds i

S

ds i

Ф*

*

Ф*

Po dv

1 d¥2

Po dv

1 d^3

Po dv

ds ds ds

+

i

x P - x P

dx3 dx.

dPo „ dP

dx,

2 J Л

- x,

dx,

dPo „ dP

dx„

2 J Л

- x,

dx.

xxdr i| 2dT i( 3dT il

dPo „ dP

Л

x, _ _ xi i I

1

dx1 dx3 j dPo „ dPo Л

■-x

3 dx1 1 dx3 j

dPo ^ dPo Л

- x

2

dx 1 dx,

2dr x1 3 dr i| x1

3 J

dPo x dPo 1

dx2 dx1 J

dPo x dPo Л

dx2 dx1 J

dPo x dPo Л

dx2 dx1 J

x1 dr

x2 dr

x3dr

+ Lo x Vo

3 3

= ZZe.ekljk VV + Lo x Vo = L- Vo + Lo x Vo;

J=1 k=1

3 3

где L = ZZjkhk, Lo =\p rdr

= mrc ;

J=1 k=1

iv(v^)dr = i® -y/vds = i>-co —Pds =

т

3 3

P dv

Z® Z ï fa P v=Z ® Z e iф P dv=

j=1 k=1 s Po dv J=1 k=1 s P dv

3 3 33

= Z® ZÏiФ (r x v) ,ds = Z® Zïi(r x v) Фdr =

J=1 k=1 s J J=1 k=1 r J

3 3

= Z ® Z ek i (r xv), (Ф* +Po г, ) dr =

J=1 k=1 r J

3 3 3 _ 3 3

Z® Z ek i (r x v) +Z c Z e i rk (r x vPo )dr+Z ® Z e i Po (r x v), rkdr

J=1 k=1

J=1 k=1

J=1 k=1

:Ь|

л

2 дх3 3 дх2 J

г ( д д^

+1 03 х--х2 —

Т V дх2 дх1 J

д Ф е + Ф2 ^ + Ф; ¿3 + ! Ф е + Ф* ^ + Ф* ёз) dт +

0

' д д^

2

Т V

дх1 1 дх,

г)(фГ +ф£+Фй Ут+

3 J

Т V

(

{ . дРо „ дРо Л

+ 0.

3

Т V

I

Т

+Ро 101

Т

+Ро 103

дх3 л3 дх2 J

дРо х дРо Л

дх2 дх1 J

х д д

2 дх3 - х3 — дх2

д д

х1- дх2 - х2 - дх1

(х1 е1 (х1 е1

+ х2 е2 + х3 е3

+ х2 е2 + х3 е3

Т V ^

( х1 е1

+ х2 е2 + х3 е3

К

( х1 е1

дФ* дФ* Л

+ х2 е2 + х3 е3

Ут + \

) ёт +

) ёт+ ро I ) ёт

{ , дРо „ дРо Л

- х

дх 1 дх.

" |(х1

3 J

+ х2е2 + х3е3 щт +

а2

т V

д д —х —

дх1 дх,

- (х1 е1

3 J

+ х2 е2 + х3 е3

) ёт +

дх,

дх,,

2 J

дФ* дФ*Л

- х„

дх0

дх3 ^2 J

дФ* дФ*Л

- х,

2 дх3 3 дх.

дРо „ дР

2 J

Л

ёт Ц

ёт к

ёт II

- х„

дх у

.. дРо „ дРо

- х

3 дх1 1 дх3

1Ь

дРо „ дРо

дх,,

-- х„

дх.

х1ёт 11 х1ёт I ( х1ёт 11

1

дФ; дх1 - х1 дФ1* Л дх3 J

дФ* дФ* Л

дх1 1 дх3 J

дФ; дФ3* Л

дх1 1 дх3 J

дРо дРоЛ

дх3 3 дх2 J

дРо дРо Л

дх1 1 дх3 J

дРо дРо Л

дх2 2 дх1 J

ёт ^ х1 ёт !(х1 ёт ^ х1

х2ёт ^ х2ёт ^ х2ёт II

дФ* дФ* Л

дх,,

дх

ёт

1

дФ* дФ*Л

дх0

дх,

ёт

1 J

дФ* дФ*Л

дх,,

-- х„

дх

ёт

1

0

0

0

дРо „ дРо

- х„

дх д^2 ^ л

х3ёт

дРо „ дР

- х

3 дх1 1 дх3 J

дРо „ дРо

х3ёт

дх.

-- х„

дх.

х3йт

1

0

0

0

4

Ро

х2&1 е3 хе2 + х32 е1 е3 + 0^3 е2 0^3 е1

ёт =

= Ь* -0-Ь хш;

где Ь - тензор сопряжённой к тензору Ь ;

IгхУ(^-а)= Xек^гхУ) уоёт =

т 1=1 к=1 т

=^0] х ^ <£(г )к =£0 х ? I-1 д¥к

1=1 к=1

1=1 к=1

Ро

у, ёя =

1 д¥1

Ро ду

1 д¥

Ро ду

1 д¥

Ро ду

¥1 ds </ — ¥1 ds (/ Ро ду

s ~о s

¥2 ds (/—д¥2- у/2 ds </

* п г) 1/ " *

s

I

s

(/—■д¥ \ffjds <£—<д¥ \ffjds <£ У р д V л р0 ду

Ро ду

1 д¥2

Ро ду

1 д ¥2

1 д¥3

Ро ду

1 д¥2

Ро ду

1 д¥3

s

3 3

Ро ду

¥1 ds ¥2 ds ¥3 ds

s го

= Е-X с =7

у=1 ¿=1

уравнения движения (1), (2) запишутся в следующем виде:

d

dt

(д + д) + Сх (д + д) = кт,

(15)

(16)

(Л/ + М') + сС х (ММ + М') + V х 0' = 4,

dt

д = Ст + С?,.«., б' = МV -1-С - 1о хс; М = Мт + М' = IV +10 XV0 + 7-с;

где С,« = dG, М= [Рг х V dG - количество и момент количества

G G

«затвердевшей» неоднородной жидкости;

М, 7,7 - тензоры, характеризующие подвижность жидкой массы, причем М,7 -тензоры присоединенной массы и присоединенных моментов инерции жидкости [10].

3 3

М=ЕЕе—; 1-;

*=1 У=1

к=1 У=1

3 3

7=ЕЕ еке/к- ; 4=$ро^т

=;

к=1 ]=1

*=/ ^+е ду ф ^;

1 дФ,

dxl

G к ^ У о

Ро ду

Г Г 1

= /х-(г хУро+ <|--^dS;

^ =

/I

Ро ду 1

— ^ Ж;

Ро ду

(17)

(18)

(19)

(20) (21)

Ьо - статический момент «затвердевшей» неоднородной жидкости массой т, Гс - радиус-вектор центра масс «затвердевшей» неоднородной жидкой массы.

Формулы (17) - (21) позволяют рассмотреть и частный случай, когда плотность частиц постоянна и движение жидкости безвихревое. При этом краевые задачи для функций Yi переходят при ро (r) = const в задачи Неймана для гармонических функций, называемых потенциалами Жуковского для данной полости. Как видно из выражений (17) - (21) компоненты тензоров M, L, состоят из двух, качественно разных, слагаемых. Первое слагаемое учитывает влияние изменения плотности жидкости на инерционные характеристики системы, последние же слагаемые, как и в случае однородной жидкости, характеризуют воздействие вихревого движения жидкости, возникающего вследствие движения твердого тела.

Из уравнений движения (15) - (16) также следует, что при приложении в центре масс системы к твердому телу, система не будет совершать поступательное движение коллинеарное вектору действия силы. Допустим, что к твердому телу с неоднородной жидкостью приложена такая система сил, при действии которой рассматриваемая механическая система будет совершать только поступательное движение. Тогда отличие от соответствующего движения твердого тела с однородной жидкостью будет состоять в том, что инерционные характеристики системы при приложении совокупности сил, компланарной одной из плоскостей Oxixj (i ф j i, j = 1,2,3 ) оказываются разными.

Указания анизотропность инерционных характеристик системы тело -неоднородная жидкость, а также существование тензоров M и L, можно объяснить возникновением вихревого движения жидкости. Так для относительной скорости частиц жидкости при наличии только поступательного движения полости имеем:

u = V-К =— У(Ф-К)-Va; (22)

Ро

Применим к обеим частям равенства (22) операцию rot

Qr = rotU = --РVpo хУ(Ф •V0) = Пу; (23)

ро

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

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

Ог = Ши = —\т Уро хУ^Р • сс) - 2сс Р2

или

Ог -2с; * о;

(25)

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

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

Если форма полости твердого тела, изменение плотности жидкости ро (х1, х2, х3) и

направление движения тела таковы, что векторное произведение в выражениях (23), (24) равно нулю, то в этих случаях скорости жидкости при поступательном движении твердого тела будут иметь однозначную потенциальную функцию, и движение тела не будет оказывать влияния на ее внутреннее движение.

4. Инерционные характеристики твердого тела с цилиндрической полостью, наполненной жидкостью, совершающей квазипотенциальное движение

В качестве примера рассмотрим движение твердого тела, имеющего цилиндрическую полость радиуса го, заполненную криогенной жидкостью глубиной Н с экспоненциальным распределением плотности вдоль оси ОХ3(р00 = Р0(х3)). Ограничимся рассмотрением движения твердого тела в плоскости Х10Х3 так что Ф1(Ф1,о,о) и

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

Для функций Ф1, имеем краевые задачи, записанные в цилиндрической системе координат (г, п, Х3)

1. Для твердого тела с полностью заполненной полостью

V- (— УФ,) = 0, V- (— VT2) = 0, Ро Ро

1 ЗФ,

Ро дг 1 дФ,

1 дТ

г = sinn ,--

Го Ро дг

Ро дХ3

= sin Л .

1 дТ

Ро дХ3

= X3 sin n

о,я =-г sinn-

г=г

х, =о,я

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

У- (— УФ1) = о, У- (— У^2) = о,

Ро Ро

1 дФ1

Ро дг 1 дФ1

1 дТ

г = sinn ,--

го Ро дг

= X3 sin n

(27)

Ро дх3

X 3 = Я = sin^

Ф,

1 дТ

X 3 =о

X 3 = я =-г sinn

= о.

X 3 =о

Ро дX3

= о, Т2

Решение обеих задач ищем в виде:

Ф1 = Ф* + Ро(X3) - г - sinn; Т2 = Т2* + Ро(X3) - г - X3 - sinn .

Функции Ф* и Т2* находим из краевых задач (26)-(27), используя метод разделения переменных. Найденные решения запишем как

Ф* = 2го sinnX

( г3

V го У

1

znw)( X3)

(X3):

(28)

J1

Т2 = 2го sin ЛХ

L

'о У

1 (L2 - 1)J! Ln )

znS( X3)

■(2)( nS *

Л*(X3)

(29)

где (£пг/го)/(%п) - нормированные решения уравнения Бесселя 1-го порядка [11], причем 1;п - корень уравнения ) = о;

KS)

Z

Z (S)

nw > nS '

ZS - есть сопутствующие функции, полученные при решении

краевых задач (26)-(27) методом разделения переменных (значок (2) означает, что функция получена при наличии свободной поверхности жидкости).

т = т.

г

X 3) = -

во якЛпН + /лпск^пН

во • го2

в *кцп (н - г )+лпскцп (н - г)

+ Сяклп2

+ е

вг Но

во • Го2

(3о)

е

воХ 3/

О X3) = -

в) ЯкЛпН + ЛпРкЛпН

в (Н - Х3) + ¡псклп (Н - Х3) + во'

РоН/

1**кМпХ3

(31)

е

г.а =

Ро X 3 /

е /2 • Го2

&кцпН

РвзИМп(Н -X3)+ ¡псклп(Н -X3)

+ С

в X3 -л„склпX3

+ е

во • Го2

во X 3 Но 'о

=во

во-^3 /

е /2 • Го2

^ п^ ^о £> 2

£ *к»я Н

в / ч / ч воН/

^¡(Н -Xi) + ^cк^{H -X3)+е /2

у X3 - ЛскЛX3

(32)

(33)

Теперь воспользуемся полученными выше выражениями для величин 2, ¡, т, 7, к0 и оценим влияние подвижности жидкости на "сопротивление" движению тела в плоскости

X!0X3

Введем безразмерные тензоры:

Т 1 I ~ 1 Л ~т 1 7

Ь = 1 + — ; л = 1 + —; 7 = 1 +-^

к

т

того

(34)

где г00 - характерный размер полости.

При движении твердого тела с цилиндрической полостью только в плоскости X¡0X3 вместо тензоров (34) будем иметь.

/12 =1 +

1л = 1+—; 722 = 1+

7.

22

т г

22 2 т т • гп

(35)

Вычисляя входящие в (35) величины по формулам (19) - (21), получим:

2во

Л11 1 + ево Н 2во

то 1

-У 1 1 -1

Л

^ 1 + во Н

евоН -1Е гс(£ -1)

+ А - воНевоН 1 во2 (евоН -1)-4

п = 1 Т> п

^ = НевонН - 2/во , 2 воНевоН

ево Н -

+

2во

ево Н -

7 Е /о

1 п=1

(36)

(37)

(38)

где:

/1 = /2т(Xз)dXз, /2 = /X32пИ,(Xз)dX3 ,

о

/ о =

£ -1

Ьп

Н 1

/ Xз - +£ \гп9(Н) - гп&(о)]

п

1

Для твердого тела, имеющего полость, частично заполненную жидкостью, интегралы 11, 12, 10 вычисляются по тем же формулам, с заменой 2пк, 2п9 на функции ,

Зависимости величин /11(в0,Н), /12(в0,Н), 322(в0,Н) даны на рис. 1-5 как для случая полости, целиком заполненной жидкостью, так и при частичном заполнении сосуда.

Из рис. 1 для величины /11 видно, что поступательное движение твердого тела оказывает заметное влияние на внутреннее движение жидкости при экспоненциальном распределении плотности и значении параметра в0 - 0.3 , (в0 = в0 • г0). До этого значения инерционность системы при движении твердого тела вдоль осей 0Х3 и 0Х1 практически оказывается одинаковой как и для твердых тел, имеющих полости, целиком заполненные

однородной жидкостью. При постоянной высоте сосуда увеличение числа в0 равноценно увеличению угловых скоростей частиц жидкости и, как следствие, приводит к уменьшению величины /11.

Рис. 1: Изменение ¿¿ц от глубины заполнения бака:

полное заполнение,

свободной поверхностью; 1 - /?0 = 0, 2 - = 0.1, 3 - /?0 = 0.3, 4 - = 0.9;

со

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

заполненной жидкостью, здесь сказывается на уменьшении величины /и11 с ростом параметра в 0.

На рис. 2 дан график величины /12 для обоих случаев заполнения полости.

Величина /12 не имеет такого явного физического смысла, как величины /и11 и /22 и служит в основном для определения положения центра инерции жидкости [18].

Рис. 2: Изменение 112 от глубины заполнения бака:---полное заполнение, —

свободной поверхностью; 1 - /?0 = 0, 2 - = 0.1, 3 - /?0 = 0.5, 4 - = 0.9;

со

Под центром инерции здесь понимается такая точка системы тело - жидкость, при приложении в которой некоторой силы Г твердое тело будет совершать поступательное движение в направлении действия приложенной силы, а при действии момента твердое тело будет совершать вращательное движение вокруг оси, проходящей через эту точку. Из характера распределения /12 для полости целиком заполненной жидкостью, следует, что центр инерции жидкости при любой высоте сосуда и при любом значении параметра в0 находится выше положения центра масс "затвердевшей" жидкости. График величины ё = Ьц - 2с (ё = ё/г0, Ьц - расстояние от оси вращения до места положения центра

инерции жидкости) подтверждает это предположение (рис. 3.) Такое поведение положения центра инерции можно объяснить тем, что в невозмущенном состояния центр масс жидкости находится ниже середины высоты заполнения.

0.3

0.2 0.1 0.0 -0.1 -0.2 -0.3 -0.4 -0.5 -0.6

Рис. 3: Изменение (1 от глубины заполнения бака:---полное заполнение,-со

свободной поверхностью; 1 - /?0 = 0, 2 - /50 = 0.1, 3 - /?0 = 0.5, 4 - /50 = 0.9;

Наличие свободной поверхности и изменение параметра в0 довольно сильно

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

параметра в0, характеризующего степень неоднородности жидкости, влияние свободной

поверхности становится меньше и центр масс "затвердевшей" жидкости остается ниже положения центра инерции при любой глубине.

Зависимости /22 (в0, Н) для полостей, целиком или частично заполненных

жидкостью, располагаются довольно близко (рис. 4). Так как ось вращения в рассматриваемом примере, совпадает с осью 0Х1, то такое поведение зависимостей

322 (в0, Н ) не требует дополнительных пояснений.

Рис. 4: Изменение

от глубины заполнения бака:

полное заполнение,

со свободной поверхностью; 1 - /?0 = 0.1, 2 - /50 = 0.5, 3 - /?0 = 0.9;

На рис. 5 представлен график относительного момента инерции жидкости при вращении вокруг оси, проходящей через центр инерции жидкости, отнесённый к моменту инерции затвердевшей жидкости, при вращении вокруг оси 0Х1. Роль свободной поверхности здесь также невелика и не вносит качественных изменений в поведении момента инерции У0 при увеличении глубины заполнения. Из этого же графика видно:

изменение степени неоднородности жидкости меняет поведение момента инерции J0 и приводит к его дальнейшему уменьшению с ростом глубины жидкости.

О 2 4 6 8 10

Рис. 5: Изменение |о от глубины заполнения бака:---полное заполнение,-со

свободной поверхностью; 1 - /50 = 0.1, 2 - /50 = 0.7, 3 - /?0 = 0.9;

Заключение. Из изложенного видно, что движение твердого тела с полостью, целиком заполненной неоднородной жидкостью, отличается от случая, рассмотренного Жуковским Н.Е. в работе [1] о движении твердого тела с однородной жидкостью. Механический эффект неоднородной жидкой массы, совершающей квазипотенциальное движение в полностью наполненном сосуде, будет эквивалентен действию «затвердевшей» жидкой массы и твердого тела, обладающего количеством движения 0 и

моментом количества движения M ', присоединенных к несущему твердому телу и образующих в совокупности систему твердых тел.

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

Список литературы

1. Жуковский Н.Е. О движении твердого тела, имеющего полости, наполненные однородной капельной жидкостью. Собр. соч. Т. 2. М.-Л.: Гостехиздат, 1948.

2. Слудский Ф. А. De la rotation de la terre suppose fluide a son interieus // Bulletin de la Societe des naturalists de Moscow. 1895. Vol. IX. P. 285-318.

3. Houqh S.S. The Oscillations of a Rotating Ellipsoidal Shell containing Fluid // Philosophical Transactions of the Royal Soc. of London (A). 1895. Vol. 186, pt. 1. P. 469-506.

4. Poincare H. Sur la precession des corps deformables // Bulletin Astronomique. 1910. Vol. XXVII. P. 321-356.

5. Черноусько Ф.Л. Движение твердого тела с полостями, содержащими вязкую жидкость. М.: ВЦ АН СССР, 1968. 232 с.

6. Колесников К.С. Динамика ракет. М.: Машиностроение, 2003. 520 с.

7. Рабинович Б.И. Введение в динамику ракет-носителей космических аппаратов. М.: Машиностроение, 1975. 416 с.

8. Луковский И.А. Введение в нелинейную динамику твердого тела с полостями, содержащими жидкость /АН УССР, Ин-т математики. Киев: Наукова думка, 1990. 296 с.

9. Ганичев А.И., Качура В.П., Темнов А.Н. Малые колебания двух несмешивающихся жидкостей в подвижном сосуде // В кн.: Колебания упругих конструкций с жидкостью. Новосибирск: НЭТИ, 1974. С. 82-88.

10. Темнов А.Н. Уравнения движения твердого тела с неоднородной жидкостью, совершающей квазипотенциальное движение // Труды МВТУ им. Н.Э. Баумана № 442. Динамика гидроупругих систем / под ред. К.С. Колесникова. М.: Изд-во МГТУ им. Н.Э. Баумана, 1985.

11. Ай Мин Вин, Темнов А.Н. О движении стратифицированной жидкости в полости подвижного твёрдого тела // Вестник МГТУ им. Баумана. Сер. Естественные науки. 2012. Спец. вып. С. 86-101.

12. Моисеев Н.К., Румянцев В.В. Динамика тела с полостями, содержащими жидкость. М.: Наука,1966. 440 с.

13. Смирнов В.И. Курс высшей математики. В 5 т. Т. 2. М.: Наука, 1974. 656 с.

14. Кочин Н.Е., Кибель И.А., Розе Н.В. Теоретическая гидромеханика. В 2 ч. Ч. 1. М.: Физматгиз, 1963. 584 с.

15. Нариманов Г.С., Докучаев Л.В., Луковский И.А. Нелинейная динамика летательного аппарата с жидкостью. М.: Машиностроение, 1977. 203 с.

16. Моисеев Н.Н. Вариационные задачи теории колебаний жидкости и тела с жидкостью // Вариационные методы в задачах о колебании жидкости и тела с жидкостью: сб. статей. М.: ВЦ АН СССР, 1962. С. 7-118.

17. Петров А.Г. Аналитическая гидродинамика: учеб. пособие для вузов. М.: Физимтлит, 2010. 520 с.

18. Охоцимский Д.Е. К теории движения тела с полостями, частично заполненными жидкостью // Прикладная математика и механика. 1956. Т. 20, вып. 1. С. 3-20. 6, т. 20, вып. I, Page 3-20.

SCIENTIFIC PERIODICAL OF THE BAUMAN MSTU

SCIENCE and EDUCATION

EL № FS77 - 48211. №0421200025. ISSN 1994-040S

electronic scientific and technical journal

On motion of a solid body with a cryogenic liquid # 12, December 2013 DOI: 10.7463/1213.0627898 Ai Min Vin, Temnov A.N.

Bauman Moscow State Technical University, 105005, Moscow, Russian Federation

ayeminwin84@gmail.com antt45@mail.ru

Relevance of this work is associated with the problem caused by ever-increasing usage of cryogenic fluids in the space rocket technology. Further exploration of space is impossible without creation of orbital cryogenic gas stations. A distinctive feature of all cryogenic liquids is a non-uniform distribution of density and temperature observed during all modes of operation and storage. A suitable model for studying motion of a mechanical system is a heterogeneous incompressible fluid. In this work interaction between a heterogeneous incompressible fluid and a movable rigid body's cavity was discussed. Cases in which motion of a rigid body is caused by the action of momentary forces or cases of abrupt movement were considered. The authors proposed a concept of quasi-potential of fluid velocities was introduced; solution to the problem of motion of a rigid body with a cylindrical cavity partially or completely filled with a cryogenic liquid.

Publications with keywords: moment number of movement, center of masses, moments of inertia, cryogenic liquid, quasipotential speeds, number of movement Publications with words: moment number of movement, center of masses, moments of inertia, cryogenic liquid, quasipotential speeds, number of movement

References

1. Zhukovskiy N.E. O dvizhenii tverdogo tela, imeyushchego polosti, napolnennye odnorodnoy kapel'noy zhidkost'yu. Sobr. soch. T. 2 [About motion of the rigid body having cavities, filled with uniform drop liquid. Collected works. Vol. 2]. Moscow - Leningrad, Gostekhizdat, 1948.

2. Sludskiy F.A. De la rotation de la terre suppose fluide a son interieus. Bulletin de la Societe des naturalists de Moscow, 1895, vol. 9, pp. 285-318.

3. Houqh S.S. The Oscillations of a Rotating Ellipsoidal Shell containing Fluid. Philosophical Transactions of the Royal Soc. of London (A), 1895, vol. 186, pt. 1, pp. 469-506.

4. Poincare H. Sur la precession des corps deformables. Bulletin Astronomique, 1910, vol. 27, pp. 321-356.

5. Chernous'ko F.L. Dvizhenie tverdogo tela spolostyami, soderzhashchimi vyazkuyu zhidkost' [Motion of rigid body with cavities containing viscous liquid]. Moscow, VTs AN SSSR Publ., 1968. 232 p.

6. Kolesnikov K.S. Dinamika raket [Dynamics of rockets]. Moscow, Mashinostroenie, 2003. 520 p.

7. Rabinovich B.I. Vvedenie v dinamiku raket-nositeley kosmicheskikh apparatov [Introduction to dynamics of carrier rockets of spacecrafts]. Moscow, Mashinostroenie, 1975. 416 p.

8. Lukovskiy I.A. Vvedenie v nelineynuyu dinamiku tverdogo tela s polostyami, soderzhashchimi zhidkost' [Introduction to nonlinear rigid body dynamics with the cavities containing liquid]. Kiev, Naukova dumka, 1990. 296 p.

9. Ganichev A.I., Kachura V.P., Temnov A.N. Malye kolebaniya dvukh nesmeshivayushchikhsya zhidkostey v podvizhnom sosude [Small oscillations of two immiscible liquids in a mobile vessel]. In: Kolebaniya uprugikh konstruktsiy s zhidkost'yu [Fluctuations of elastic designs with liquid]. Novosibirsk, NETI Publ., 1974, pp. 82-88.

10. Temnov A.N. Uravneniya dvizheniya tverdogo tela s neodnorodnoy zhidkost'yu, sovershayushchey kvazipotentsial'noe dvizhenie [The equations of movement of a rigid body about the non-uniform liquid making quasipotential movement]. In: Kolesnikov K.S., ed. Trudy MVTU im. N.E. Baumana No. 442. Dinamika gidrouprugikh sistem [Proc. of the Bauman MSTU No. 442. Dynamics of hydroelastic systems]. Moscow, Bauman MSTU Publ., 1985.

11. Ay Min Vin, Temnov A.N. O dvizhenii stratifitsirovannoy zhidkosti v polosti podvizhnogo tverdogo tela [On the motion of stratified liquid in the cavity of movable solid]. VestnikMGTU im. Baumana. Ser. Estestvennye nauki [Herald of the Bauman MSTU. Ser. Natural science], 2012, spec. iss., pp. 86-101.

12. Moiseev N.K., Rumyantsev V.V. Dinamika tela spolostyami soderzhashchimi zhidkost' [Dynamics of body with cavities containing liquid]. Moscow, Nauka, 1966. 440 p.

13. Smirnov V.I. Kurs vysshey matematiki. V5 t. T. 2 [Course of Higher Mathematics. In 5 vols. Vol. 2]. Moscow, Nauka, 1974. 656 p.

14. Kochin N.E., Kibel' I.A., Roze N.V. Teoreticheskaya gidromekhanika. V2 ch. Ch. 1 [Theoretical hydromechanics. In 2 parts. Pt. 1]. Moscow, Fizmatgiz, 1963. 584 p.

15. Narimanov G.S., Dokuchaev L.V., Lukovskiy I.A. Nelineynaya dinamika letatel'nogo apparata s zhidkost'yu [Nonlinear dynamics of the aircraft with liquid]. Moscow, Mashinostroenie, 1977. 203 p.

16. Moiseev N.N. Variatsionnye zadachi teorii kolebaniy zhidkosti i tela s zhidkost'yu [Variational problems of the theory of oscillations of liquid and liquid-containing body]. In: Variatsionnye metody v zadachakh o kolebanii zhidkosti i tela s zhidkost'yu: sb. statey [Variation methods in problems of oscillations of liquid and liquid-containing body: Collected articles]. Moscow, VTs AN SSSR Publ., 1962, pp. 7-118.

17. Petrov A.G. Analiticheskayagidrodinamika [Analytical hydrodynamics]. Moscow, Fizimtlit, 2010. 520 p.

18. Okhotsimskiy D.E. K teorii dvizheniya tela s polostyami, chastichno zapolnennymi zhidkost'yu [To the theory of motion of body with cavities partially filled with liquid]. Prikladnaya matematika i mekhanika, 1956, vol. 20, no. 1, pp. 3-20.