CC BY
149
НАУКА и ОБРАЗОВАН И Е
Эл № ФС77 - 48211. Государственная регистрация №0421200025. ISSN 1994-0408
Колебания жидкого топлива в цилиндрических и конических
ёмкостях
# 11, ноябрь 2013
DOI: 10.7463/1113.0623923
Дьяченко М. И., Орлов В. В., Темнов А. Н.
УДК 621.3
Россия, МГТУ им. Н.Э. Баумана antt45@mail. ru s [email protected]
Введение.
Исследованию задач о колебаниях жидкости в баках ракет-носителей посвящено много работ [1-11]. В предыдущих работах авторов [12, 13] приведена постановка модельной задачи о малых движениях несжимаемой жидкости, вытекающей из топливного бака с заборными устройствами (ЗУ). Исследование рассматриваемой задачи в этих работах, проведенные методами функционального анализа и спектральной теории
операторных пучков, показало, что спектр нормальных движений
Ät
несжимаемой жидкости (движений, подчиняющихся закону е ) обладает двумя ветвями собственных значений: дискретным множеством вещественных чисел, расположенных на положительной части вещественной оси и дискретным множеством комплексно сопряженных чисел, расположенных вблизи мнимой оси. Действительным собственным значениям отвечают апериодические режимы движения жидкости, происходящие преимущественно вблизи дна топливного отсека, а комплексно сопряженным числам отвечает режимы затухающих колебаний, преобладающих в основном на свободной поверхности жидкости.
В предлагаемой статье приведены примеры рассматриваемых задач о движении несжимаемой жидкости, вытекающей из наиболее распространенных на практике топливных отсеков в виде цилиндра и конуса.
1. Постановка задачи.
Пусть несжимаемая жидкость, частично заполняющая неподвижный бак произвольной формы вытекает через ЗУ и может совершать малые движения. Рассматриваемая проблема малых движений жидкости может быть описана уравнениями гидродинамики, линеаризованными вблизи невозмущённого состояния. Подробная постановка задачи приведена в работах [12-14].
Предполагая возмущенное движение жидкости потенциальным, сформулируем краевую задачу о колебаниях вытекающей жидкости для потенциала скоростей Ф( х, г, п, х) в виде
АФ = 0 в области, заполненной жидкостью, 0,
дФ
— = 0 на смачиваемой поверхности £,
дп
— _ — -У0(Н) + gwГ = / на свободной поверхности Г0, (1)
дХ дХ^
дФ дФ
----¥0(0)- у^Е = /2 на поверхности сливаЕ
дх дХз
Ф ( х1. х2. х3.0 ) = Ф(0), при Х = 0 ,
д сдФ„
где — - производная по внешней нормали к поверхности £; = 1 —Л,
дп 1 д х3
дФ г =--малые смещения и скорость частиц жидкости на поверхностях Г0 и
д х3
Е соответственно, у = дУ^ = д¥0(0)- обобщенный коэффициент сопротивления поверхности слива, д - коэффициент гидравлического сопротивления поверхности слива, g - величина интенсивности внешнего однородного поля массовых сил, /1, /2 - заданные поля внешних воздействий соответственно на
поверхностях Г0 и Е.
Проинтегрировав уравнение неразрывности по объёму, занимаемому жидкостью, для любого момента времени ^ получим
дополнительное интегральное условие | V ■ пГ ёГ = | V ■ пЕ dS, которому
Г I
должны подчиняться поле скоростей и поле смещений в рассматриваемой задаче.
Впервые исследование колебаний жидкости с учетом вытекания было предложено Кирилловым В.В. [15] и продолжено в работах [5], [11], [16]. В упомянутых работах рассматривались задачи для жидкости, занимающей часть цилиндрического бака, на дне которого ставилось кинематическое условие вытекания.
2. Колебания жидкого топлива в цилиндрической ёмкости
Рассмотрим задачу о собственных движениях жидкости, вытекающей через заборное устройство из цилиндрической ёмкости при наличии свободной поверхности. Используя цилиндрические координаты г, П, х с началом на поверхности слива, получим задачу для определения потенциала скоростей Ф( х, г ,п, t):
д2Ф 1 дФ т2 _ дФ Л _
—г +----гФ + —г = 0 в Q , т = 0,±1,±2,...,
дг г дг г д х
дФ дФ сдФ . Л ^
^(Н) + 8= 0 на Го, (2)
дг дх •> д х
дФ дФ ( дФ
---V0(0)-у-= 0 на I,
дг дх дх
дФ
— = 0 на S,
дг
Ф( х, г ,п, t )= Ф(0) при г = 0.
Функцию Ф( х, г ,п, t) представим в виде суммы двух функций Ф(х, г ,п, t) = Ф1(х, г ,п, t) + Ф2(х, г ,п, t), и будем искать каждый потенциал Ф1( х, г ,п, t) и Ф2( х, г,п, t) в виде
произведения четырёх функций:
Фх(х,г,",t) = Z(1)(х)R (г)Н(")Я ($), Ф2(х,г,",t) = Z(2)(х)R (г)Н(")Д (Г), где Н(") - функция периода 2п, функция Я (г) удовлетворяют уравнению
Бесселя, а функции 2 (1)( х), 2 (2)( х) есть решения задач
—2 7(1)
^^^ - к27(1) = 0, (0 < х < Н),
— х (3)
— ,7(1) = 0, х = 0, —7(1) = 1, х = Н, —х —х
— 27(2)
—- к27(2) = 0, (0 < х < Н),
— х (4)
7(2) = 1, х = 0, —7(2) = 0, х = Н .
—х —х
Используя метод разделения переменных, находим выражения для потенциалов скоростей
ф (х, г, п t) = £ £ Зт (ктп г )в1т11 ^ ^), (5)
т=0 п=1 (ктпН )
Ф2( х, г,,, t) = ££ Л, (ктп г У" Ск(ктП;Х, )) >1тп С), (6)
т=0 п=1 ктпН )
где Jm(ктпг) - функция Бесселя 1-го рода, т-го порядка, ктп = %тп - корни уравнения = а функции Ятп(*), ¿тп ) удовлетворяют системе
уравнений
^тп + атп^тп + ^0£тп^тп + ^^тп^тп + 0 тп^ тп = 0,
^тп + атп^тп + ^^тп^тп + (^0£тп + Утп )^тп = 0, (7)
т = 0,1,2..., п = 1,2,3..., где атп, £тп - коэффициенты инерционных и диссипативных связей, 0)2тп, Утп, $тп, Б'тп - динамические характеристики парциальных подсистем
amn = , / 7 , Smn = kmnSh l(kmnИ) , Уmn = Уkmnth(kmnИ),
СК КпИ )
Б = k ^И), 8 = 2k /sh(2k И), О)2 = ?к гh(k И).
^п mn \ mn /' mn mn V mn /' mn о mn V mn /
Для определения собственных частот рассматриваемой механической системы, положим £ = А еО г, Л = А еО г. Из уравнений
5 mn ^п ? mn 2mn ^ ±
(7) получаем характеристической уравнение
ОЩ^пИ) + О2^ (V + г) + О^ (Vо8mnУ + 8) + ^п + Гтп) = 0. (8) Здесь О - комплексный коэффициент затухания волновых движений жидкости m = 0,1,2..., п = 1,2,3...,.
Уравнение (8) имеет две ветви решений: действительные корни и ветвь комплексно сопряженных корней. Для исследования уравнения введём безразмерные параметры
о=О /481^0, 7=7/480. р=^г>, И=И/
и перепишем уравнение (8) в безразмерном виде
'Г
- - - 2 V Р
ОЩ^пИ) + ор V + У) + ОР г +1) +
(9)
+4 (КСЧ^ПИ ) + )) = 0.
Результаты вычислений корней кубического уравнения (9) приведены в таблице (1).
В случае глубокой жидкости (И > 1, ^(ртпИ) ^ 1, sh(2pmnИ) ^ да) уравнение (19) принимает вид
О3 + о24п (V; + 7) + Ор + р2П V + 7) = 0
и имеет корни
О01,2)=±^ , О03)=-(^ + 74. (10)
В таблице № 1 представлены результаты вычислений безразмерных собственных частот и коэффициентов распределения х( k) = А1 / А2 при И = 1, m=1, п=1,о11 = 1,323 для различных значений V*, в = 7~1, а на рис. 1-4 показаны формы колебаний.
Таблица № 1 - Собственные частоты колебаний жидкости в цилиндрической емкости при перераспределении жидкости из бака
V* п 1 ^2,3 /(1) 1
0 0,05 -12,348 -0,00742±1,324i -0,306 -0,00408±0,036i
0 0,15 -4,080 -0,021±1,330i -0,280 -0,0338±0,0990i
0,04 0,01 -61,834 -0,00268±1,323Ш -0,309 -0,0003±0,0073i
0,04 0,05 -12,37 -0,00857±1,324Ш -0,305 -0,0047±0,0360i
0,04 0,1 -6,175 -0,01544±1,3266i -0,295 -0,0170±0,0690i
0,04 0,15 -4,103 -0,02143±1,330i -0,279 -0,0340±0,0980i
0,07 0,01 -61,85 -0,00357±1,32322i -0,309 -0,0004±0,0073i
0,07 0,05 -12,387 -0,00942±1,32432i -0,305 -0,0052±0,0360i
0,07 0,1 -6,192 -0,0162±1,326946i -0,294 -0,0178±0,0690i
0,07 0,15 -4,12 -0,022±1,3307167i -0,278 -0,036±0,097i
0,1 0 -6,18^-9 -0,003±1,323i -0,309 -3,2бс(-12) ±7,33с(-Щ
0,1 0,01 -61,867 -0,00445±1,3232i -0,309 -0,0005±0,0073i
0,1 0,05 -12,404 -0,0103±1,3245i -0,305 -0,0056±0,036i
0,1 0,1 -6,209 -0,0169±1,3272i -0,293 -0,019±0,069i
0,1 0,15 -4,137 -0,0227±1,331Н -0,276 -0,037±0,096i
л,
Н=0.5 Н =1
Рис. 1. Формы колебаний жидкости в цилиндрическом сосуде при 0*=-61.834, , m=0, п=1, £01 =3. 83, А1/А2=-0.309.
Н=0.5 Н=1
Рис. 2. Формы колебаний жидкости в цилиндрическом сосуде при различных Н при 0*=-0.0027+1.321, ш=1, п=1, Р11=1.841, А2/А1=0.007.
Н=0.5 Н=1
Рис. 3. Формы колебаний жидкости в цилиндрическом сосуде при различных Н при 0*=-122.45, ш=1, п=1, Р11=1.841, А1/А2=-0.043.
Н=0.5 Н=1
Рис. 4. Формы колебаний жидкости в цилиндрическом сосуде при различных Н при 0*=-0.000076+1.961, ш=0, п=1, Р01 =3. 83, А2/А1=0.0007.
3. Колебания жидкого топлива в конической ёмкости
Рассмотрим круговой конус, частично заполненный жидкостью.
Введем сферическую систему координат R,3,r с центром в вершине конуса (рис. 5). Угол 3 отсчитывается от положительного направления оси OZ, угол П измеряется в плоскости OXY от оси OX, сторону оси OY . Из рис. 5 следует. что Z = R cos 3, y = R sin 3 sin r, x = R sin 3- cos r. Если ограничиться малым углом конусности. то оказывается правомерным граничное условие для плоской невозмущённой свободной поверхности рассматривать на части сферической поверхности R=R1. Очевидно, что уравнение для конической полости R=R2, боковой поверхности - 3 = 30. Условием применимости полученных формул для плоской свободной поверхности является следующее ограничение, накладываемое на угол конусности: sin 3 0 >> 1 - cos 30. Наиболее полное исследование колебаний жидкости в усечённых конических баках без учёта эффекта вытекания приведено в работе [18].
Рис. 5. Геометрические характеристики обратного и прямого конусов
Примем за характерный размер радиус дна конической полости R=R2 и введем безразмерный радиус г = R/R2. Обозначим q = R1/R2. Допустим далее, что жидкость либо достигает вершины конуса, либо конус усеченный и имеет дно R=R2 (рис. 5 а, б). Для сокращения будем называть полость на рис. 5 а - обратным усеченным (не усеченным) конусом ^ > 1).
а) с] > 1
б) q < 1
Краевая задача о свободных колебаниях жидкости в конической полости имеет вид
дФ
ДФ = 0.
=0
З=З
дЗ
дФ 1 дФ . g гдФ
1 я Г дФ
-У0(ч) + ^ ¡-¡г* = 0 на Г
дг я2 дг к,* дг 0'
дФ 1 дФ т у дФ „ пп ----К(1) ----= 0 на I, (11)
дг к2 дг 0 Я2 дг
где Д - лапласиан, записываемый в сферических координатах [17].
Потенциал скоростей Ф( г ,З,п, t) для конической полости представим в виде суммы потенциалов
Ф(г ,З,п, t) = Ф1(г ,З,п, t) + Ф2(г ,З,п, t), (12)
каждый из которых есть произведение функций
Ф1(г ,З,п, t) = У (З) R 5 (г) Н (п) 5 (t), Ф2{ г ,З,п, t) = У (З) RЛ(г) Н (n)Л(t),
где У(З)-присоединённая функция Лежандра, а функции К(5)(г), К(Л)(г), есть решения краевых задач
* (г2 * (К5 (г))) -у(у +1) К5 (г) = 0, аг аг
* (Я5 (г)) = 0, г = 1, * (К5 (г)) = 1, г = д, аг аг
* (г2 * (КЛ(г)))-у(у +1) КЛ(г) = 0, аг
* (КЛ( г)) = 1, г = 1, * (КЛЛ г)) = 0, г = д.
(13)
(14)
Используя метод разделения переменных, находим представления решений поставленных краевых задач
Фх{т Лп, t) = SSYm» ЛУ^ (r) Sm (t), (15)
m=0 n=1
Ф2(Г,Л, П, t) = W^Rl (r)ln (t)
m=0 n=1
PL (cosW)
где Ymn (Л) = —--, - присоединенные функции Лежандра первого рода,
Pvm (cOsW0)
mn ^
v -го порядка, vmn - n -ый корень уравнения
d
(Pvmmn (cos Л)) = 0 .
dW
Функции RSmn (r). Rl (r) соответственно равны
RL(r) = v +i , ^ ■■—, (16)
qvmn +2 (v + 1) rV<mn + V Г (Vmn +1)
_ 4_\Vmn "*" 4r "*" Vmn '
mn V 'J 2v +1 1 / 1 \
q Vmn 1 -1 v (v +1)
mn mn
Rl (r) = q2 (vm„ + 1)rvmn + vm„q^+1)r"^ +1)
mn
а (2^+1) -1 V (V +1)
1 тп\ тп /
Подставив полученные решения (16) в граничные условия (11), и предполагая, что 8тп = 80вп 1, Лтп = Л0в°1, получим уравнение для определения собственных чисел:
а0П3 +а1П2 + а2П + а3 = 0, (17)
где коэффициенты вычисляются по формулам
а0 = а11а22 ~а12а2Ъ
а1 = а11@22 + а22@11 - а12@21 - а21@21,
а2 = Р11Р22 +а22®2 -Р12Р21, (18)
а3 =Р22(®тп )2.
Коэффициенты инерционных и диссипативных связей и квадрат парциальной частоты а>2тп = (а*тп)2/Rsmn(а) колебаний свободной поверхности жидкости в коническом баке представлены в таблице 2.
Таблица № 2 - Коэффициенты а, , (а>*тп )2 при обобщенных координатах
для вычисления ак при /, ] = 1,2, к = 0,1,2,3
а11 Кп (я) я( Я 2"тп +1 • (V + 1) + У ) Ч. хЧ. \ тп / тп / Утп V + 1)(Я2Утп +1 - 1)
«12 ^п (я) Я'тп+2(2Утп + 1) Vтп (Vтп + 1)(Я^ " - 1)
«21 К. (1) Я'тп+2(2Утп + 1) V (V + 1)( Я2"тп- 1) т т
а22 *!пп (1) Я2 Г (V + 1) + V Я2Утп 1 ^ V т п / т пЧ. V (V + 1)( я'2Vmn +1 - 1) т т
в„ Кп (Я) - Я Г (1 - V 2)я^^^^+2 - 2^ 2 + V +1)я2у™+1 - V (V + 2) 1 Я тп /Ч. \ тп тп /Ч. тп\ тп / 1
Г V V + 1)( я2Ут"+1 -1) ]
в\2 я,! (я) 2Vmn + 1 Я+ 2(1 Vтп ) - Я"" +1 (Vmn + 2) V™ Кп +1) (Я2Утп+ - 1)2
021 Ят! (1) я 2Vmn + 1 ЯЗ'тп + 2(1 Vтп ) - Я"" +1 (Vmn + 2) V (V + 1) (Я2"тп- 1)2 тп\ тп / /
в22 кп (1) - я/ Я2 2v я4гтп+3 + Я2'тп+2 Г1 - 4v - 4v 21 - 2я(v + 1) у т п± ^ т п т п т п / . 2/
Vтп V + 1)(Я2^ + - 1)2 ^ ^2
* 2 (^тп ) £ Я2 £ Я2
В таблице 3 представлены результаты расчетов собственных чисел в зависимости от скорости Vо и параметра в .
Таблица № 3 - Собственные частоты колебаний жидкости в конической
емкости при 3о = 45°, V = 2, д = 2 , т=1, п=1.
V 0 в о1 ^2,3
0 0 0 0,000 ± 0,9741
0 0,2 -15,798 -0,004 ± 0,9741
0,05 0,2 -15,704 -0,008 ± 0,9751
0,1 0,2 -15,610 -0,012 ± 0,9751
0 0,5 -6,302 -0,010 ± 0,9761
0,05 0,5 -6,208 -0,014 ± 0,9761
0,1 0,5 -6,114 -0,018 ± 0,9771
0 1 -3,123 -0,019 ± 0,9801
0,05 1 -3,028 -0,023 ± 0,9811
0,1 1 -2,934 -0,027 ± 0,9821
Также рассматриваемая модельная задача для конической области была решена вариантом метода конечных элементов (МКЭ), численно реализующим нахождение обобщённого решения некоторого интегрального тождества, полученного методом Галёркина. В следующих таблицах приведены результаты расчетов для конусов, аналогично изображенных на рисунке 5, но с плоским дном (таблица № 4) и с плоской свободной поверхностью (таблица № 5).
Таблица № 4 - Собственные частоты колебаний жидкости в конической емкости с плоским дном при За = 45°, V = 2, д = 2 , т=1, п=1.
V 0 в П1 02,3
0 0 0 0,00 ± 0,99091
0 0,2 -17,569 -0,001861 ± 0,99091
0 0,5 -7,02015 -0,004585 ± 0,991421
0 1 -3,4977 -0,008707 ± 0,993161
Таблица № 5 - Собственные частоты колебаний жидкости в конической емкости с плоской свободной поверхностью при Зо = 45°, V = 2, д = 0.5 , т=1,
n=1.
V 0 ß Q1 Q23
0 0 0 0,00 ± 2,5746i
0 0,2 -10,47867 -0,01931 ± 2,57469i
0 0,5 -4,11859 -0,03727 ± 2,59303i
0 1 -2,011915 -0,0401 ± 2,62067i
В заключение необходимо отметить, что полученные результаты подтверждают выводы, сформулированные ранее в работах авторов [12, 13], то есть спектр нормальных движений несжимаемой жидкости, вытекающей из ограниченной области, обладает двумя ветвями собственных значений: дискретного множества вещественных чисел, и дискретного множества комплексно-сопряженных чисел, расположенных вблизи мнимой оси.
Список литературы
1. Колесников К.С. Динамика ракет. М.: Машиностроение, 2003. 500 с.
2. Луковский И.А. Введение в нелинейную динамику твердого тела. Киев: Наукова думка, 1990. 296 с.
3. Рабинович Б.И. Введение в динамику ракет-носителей космических аппаратов. М.: Машиностроение, 1975. 416 с.
4. Микишев Г.Н.. Рабинович Б.И. Динамика тонкостенных конструкций с отсеками, содержащими жидкость. М.: Машиностроение, 1971. 563 с.
5. Моисеев Н.Н., Петров А.А. Численные методы расчёта собственных частот колебаний ограниченного объёма жидкости. М.: ВЦ АН СССР, 1966. 270 с.
6. Моисеев Н.Н. Румянцев В.В. Динамика тел с полостями, содержащими жидкость. М.: Наука, 1965. 440 с.
7. Черноусько Ф.Л. Движение твердого тела с полостями, содержащими вязкую жидкость. М.: ВЦ АН СССР, 1968. 232 с.
8. Фещенко С.Ф., Луковский И.А., Рабинович Б.И., Докучаев Л.В. Методы определения присоединенных масс жидкости в подвижных полостях. Киев: Наукова думка, 1969. 250 с.
9. Нариманов Г.С., Докучаев Л.В., Луковский И.А. Нелинейная динамика летательного аппарата с жидкостью. М.: Машиностроение, 1977. 203 с.
10. Колесников К.С., Шкапов П.М., Пожалостин А.А. Задачи динамики гидромеханических систем в трудах кафедры теоретической механики имени профессора Н.Е. Жуковского // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2012. Спец. вып. № 8. С. 15-30.
11. Лимарченко О.С., Матараццо Д., Ясинский В.В. Динамика вращающихся конструкций с жидкостью. Киев: ГНОЗИС, 2002. 304 с.
12. Степанова М.И., Темнов А.Н. Малые движения жидкости с поверхностной диссипацией энергии // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2011. № 4. С. 99-110.
13. Дьяченко М.И., Темнов А.Н. Собственные колебания жидкого топлива в условиях перераспределения // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. 2012. № 3. С. 31-38.
14. Дьяченко М.И., Темнов А.Н. Проблемы динамики перераспределения топлива в крупногабаритных ракетно-космических объектах // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. 2012. Спец. вып. С. 164-174.
15. Кириллов В.В. Исследование колебаний жидкости в неподвижном сосуде с учётом её вытекания // Труды МФТИ. 1960. Вып. 5. С. 19-25.
16. Орлов В.В., Темнов А.Н. Малые движения жидкости, вытекающей из бака // Современные методы теории функций и смежные проблемы : тез. докладов. Воронеж: Изд-во ВГУ, 1997. С. 124.
17. Кошляков Н.С., Глинер Э.Б., Смирнов М.М. Основные дифференциальные уравнения математикой физики. М.: Физматлит, 1962. 767 с.
18. Луковский И.А., Солодун А.В., Тимоха А.Н. Собственные колебания жидкости в усеченных конических баках // Акустический вестник. 2006. Т. 9, № 3. С. 42-61.
SCIENTIFIC PERIODICAL OF THE BAUMAN MSTU
SCIENCE and EDUCATION
EL № FS77 - 48211. №0421200025. ISSN 1994-040S
electronic scientific and technical journal
A problem of propellant oscillations in cylindrical and conical tanks
# 11, November 2013
DOI: 10.7463/1113.0623923
D'yachenko M.I., Orlov V.V., Temnov A.N.
Bauman Moscow State Technical University, 105005, Moscow, Russian Federation
antt45@mail. ru s [email protected]
A problem of fluid oscillations in tanks of space launch vehicles is quite well known and has been investigated extensively during the last decades. This is caused by the fact that the major mass of a liquid-propellant rocket is fluid and its forced motion affects vehicle controllability significantly. In their recent works the authors considered small oscillations in tanks with a propellant intake and with fuel reallocation inside a vehicle. Such problems have arisen due to the purpose of increasing efficiency and output performance of multistage cluster rockets by central stage fueling from side stages of a cluster. The problem of propellant free oscillations in cylindrical and conical tanks was considered, and solutions were presented with boundary conditions on a free surface and the surface of intake that, figuratively speaking, provides resistance to the fluid descent. Eigen-frequencies and Eigen-values of perturbed motion equations with dissipation on boundary surfaces were found. It was shown that slow fluid decent and intake surface boundary conditions may affect both oscillating and damped terms of the solution.
Publications with keywords: the cylindrical and conical tanks, rearrangement, cluster configuration, eigen values and eigenvibrations
Publications with words: the cylindrical and conical tanks, rearrangement, cluster configuration, eigen values and eigenvibrations
References
1. Kolesnikov K.S. Dinamika raket [Dynamics of missiles]. Moscow, Mashinostroenie, 2003.500 p.
2. Lukovskiy I.A. Vvedenie v nelineynuyu dinamiku tverdogo tela [Introduction to nonlinear rigid body dynamics]. Kiev, Naukova dumka, 1990. 296 p.
3. Rabinovich B.I. Vvedenie v dinamiku raket-nositeley kosmicheskikh apparatov [Introduction to the dynamics of carrier rockets of space vehicles]. Moscow, Mashinostroenie, 1975. 416 p.
4. Mikishev G.N.. Rabinovich B.I. Dinamika tonkostennykh konstruktsiy s otsekami, soderzhashchimi zhidkost' [Dynamics of thin-walled structures with compartments containing liquid]. Moscow, Mashinostroenie, 1971. 563 p.
5. Moiseev N.N., Petrov A.A. Chislennye metody rascheta sobstvennykh chastot kolebaniy ogranichennogo ob'ema zhidkosti [Numerical methods of calculation of natural frequencies of oscillations of limited volume of fluid]. Moscow, Publ. of VTs AN SSSR, 1966. 270 p.
6. Moiseev N.N. Rumyantsev V.V. Dinamika tel spolostyami, soderzhashchimi zhidkost' [Dynamics of bodies with cavities containing liquid]. Moscow, Nauka, 1965. 440 p.
7. Chernous'ko F.L. Dvizhenie tverdogo tela spolostyami, soderzhashchimi vyazkuyu zhidkost' [Motion of rigid body with cavities containing viscous fluid]. Moscow, Publ. of VTs AN SSSR, 1968. 232 p.
8. Feshchenko S.F., Lukovskiy I.A., Rabinovich B.I., Dokuchaev L.V. Metody opredeleniya prisoedinennykh mass zhidkosti vpodvizhnykhpolostyakh [Methods of determination of added mass of liquid in moving cavities]. Kiev, Naukova dumka, 1969. 250 p.
9. Narimanov G.S., Dokuchaev L.V., Lukovskiy I.A. Nelineynaya dinamika letatel'nogo apparata s zhidkost'yu [Nonlinear dynamics of aircraft with liquid]. Moscow, Mashinostroenie, 1977. 203 p.
10. Kolesnikov K.S., Shkapov P.M., Pozhalostin A.A. Zadachi dinamiki gidromekhanicheskikh sistem v trudakh kafedry teoreticheskoy mekhaniki imeni professora N.E.Zhukovskogo [Problems of hydromechanical systems dynamics in proceedings of the zhukovsky theoretical mechanics department]. VestnikMGTUim. N.E. Baumana. Ser. Estestvennye nauki [Herald of the Bauman MSTU. Ser. Natural science], 2012, spec. iss. no. 8, pp. 15-30.
11. Limarchenko O.S., Mataratstso D., Yasinskiy V.V. Dinamika vrashchayushchikhsya konstruktsiy s zhidkost'yu [Dynamics of rotating structures with liquid]. Kiev, GNOZIS, 2002. 304 p.
12. Stepanova M.I., Temnov A.N. Malye dvizheniya zhidkosti s poverkhnostnoy dissipatsiey energii [Small motions of liquid with surface energy dissipation]. Vestnik MGTU im. N.E. Baumana. Ser. Estestvennye nauki [Herald of the Bauman MSTU. Ser. Natural science], 2011, no. 4, pp. 99-110.
13. D'yachenko M.I., Temnov A.N. Sobstvennye kolebaniya zhidkogo topliva v usloviyakh pereraspredeleniya [Natural oscillations of liquid propellant under redistribution conditions]. Vestnik MGTU im. N.E. Baumana. Ser. Mashinostroenie [Herald of the Bauman MSTU. Ser. Mechanical Engineering], 2012, no. 3, pp. 31-38.
14. D'yachenko M.I., Temnov A.N. Problemy dinamiki pereraspredeleniya topliva v krupnogabaritnykh raketno-kosmicheskikh ob'ektakh [Problems of fuel redistribution dynamics in large-sized rocket and space objects]. Vestnik MGTU im. N.E. Baumana. Ser. Mashinostroenie [Herald of the Bauman MSTU. Ser. Mechanical Engineering], 2012, spec. iss., pp. 164-174.
15. Kirillov V.V. Issledovanie kolebaniy zhidkosti v nepodvizhnom sosude s uchetom ee vytekaniya [Study of oscillations of liquid in stationary vessel taking into account its leakage]. Trudy MFTI, 1960, iss. 5, pp. 19-25.
16. Orlov V.V., Temnov A.N. Malye dvizheniya zhidkosti, vytekayushchey iz baka [Small movement of the fluid flowing out of the tank]. Sovremennye Metody Teorii Funktsiy I Smezhnye Problemy : tez. dokladov [Modern Methods of the Theory of Functions and Related Problems: abstracts]. Voronezh, VSU Publ., 1997, p. 124.
17. Koshlyakov N.S., Gliner E.B., Smirnov M.M. Osnovnye differentsial'nye uravneniya matematikoyfiziki [Basic differential equations of mathematics physics]. Moscow, Fizmatlit, 1962. 767 p.
18. Lukovskiy I.A., Solodun A.V., Timokha A.N. Sobstvennye kolebaniya zhidkosti v usechennykh konicheskikh bakakh [Natural vibrations of a liquid in truncated conical tanks]. Akusticheskiy vestnik, 2006, vol. 9, no. 3, pp. 42-61.
