Проблемы машиностроения и надежности машин, 2019, № 7, стр. 54-63
ОЦЕНКА ДОПУСТИМОЙ ОСЕВОЙ СЖИМАЮЩЕЙ НАГРУЗКИ НА ЦИЛИНДРИЧЕСКУЮ ОБОЛОЧКУ ВЕРТИКАЛЬНОГО КРИОГЕННОГО БАКА ПРИ ЕГО ЗАПОЛНЕНИИ
В. С. Зарубин 1, *, В. Н. Зимин 1, Г. Н. Кувыркин 1
1 Московский государственный технический университет имени Н.Э. Баумана
г. Москва, Россия
* E-mail: fn2@bmstu.ru
Поступила в редакцию 11.07.2019
Принята к публикации 26.08.2019
Аннотация
При захолаживании или заполнении вертикально расположенного цилиндрического криогенного бака вдоль образующей его оболочки возникает неравномерное распределение температуры, перемещающееся вместе с подъемом уровня жидкости. Это приводит к искривлению образующей, возникновению подвижной области краевого эффекта и снижению сопротивления оболочки деформированию под действием осевой сжимающей силы. На основе разработанной математической модели термонапряженного состояния круговой цилиндрической оболочки такого бака установлена зависимость между сжимающей силой, геометрическими параметрами оболочки и характеристиками ее материала, скоростью подъема уровня криогенной жидкости и ее свойствами. Эта зависимость позволяет оценить допустимое в пределах линейно-упругого поведения материала оболочки значение действующей на оболочку сжимающей нагрузки. Такая оценка полезна при расчете и проектировании оболочек криогенных баков летательных аппаратов.
1. Введение. Крупногабаритные емкости для хранения и транспортировки сжиженных газов находят широкое распространение в различных областях техники. Сжиженный природный газ [1], основную долю которого составляет метан, является перспективным горючим для двигателей летательных аппаратов [2]. Для жидкостных ракетных двигателей в качестве окислителя и горючего используют соответственно жидкие кислород и водород [3, 4]. Процессы заполнения или захолаживания криогенной емкости низкотемпературной жидкостью, регламентированные перед заполнением емкости [5], связаны с возникновением в ее оболочке неравномерного распределения температуры. В силу существенного различия интенсивности теплообмена с жидкой и газообразной фазами наиболее резко изменяется температура оболочки в зоне уровня криогенной жидкости (как подвижного, так и неподвижного).
Одним из распространенных вариантов крупногабаритных криогенных емкостей являются топливные баки первых ступеней жидкостных ракет-носителей, основным конструктивным элементом которых является тонкостенная круговая цилиндрическая оболочка [4, 6]. Неравномерное распределение температуры в окрестности уровня криогенной жидкости вызывает искривление образующей такой оболочки и снижает ее сопротивление действию приложенной осевой сжимающей нагрузки. В данной статье с использованием результатов математического моделирования термонапряженного состояния цилиндрической оболочки криогенного бака при его заполнении [7–9] получены соотношения, позволяющие оценить допустимую осевую нагрузку на оболочку, при которой ее материал остается в пределах линейно-упругой области.
2. Температурное состояние оболочки. При заполнении вертикально расположенного цилиндрического криогенного бака на температурное состояние его оболочки влияют ряд факторов. Среди них скорость подъема уровня жидкости и интенсивность теплообмена на смоченном и несмоченном участках внутренней поверхности оболочки, определяемые соответственно коэффициентом теплообмена ${{\alpha }_{1}}$ с жидкостью, имеющую температуру $T_{1}^{*}$, и коэффициентом теплообмена ${{\alpha }_{2}}$ с имеющей температуру $T_{2}^{*}$ газовой средой в полости бака над уровнем жидкости.
Для уменьшения подвода теплоты от внешней среды и снижения потерь криогенной жидкости при заполнении бака внешнюю поверхность его оболочки можно покрыть слоем малотеплоемкой теплоизоляции с коэффициентом теплопроводности $\lambda {\kern 1pt} {\text{'}}$ и толщиной $h{\kern 1pt} '$. Например, оболочка топливного бака многоразовой космической системы Space-Shattle [6] была покрыта слоем пенополиуретана, что предотвращает образование инея и наледи на внешней поверхности бака в период его заполнения. При наличии внешнего слоя теплоизоляции эффективное значение коэффициента теплообмена оболочки с внешней средой определяет равенство ${{\alpha }_{3}} = 1{\text{/}}(1{\text{/}}\alpha _{3}^{'} + h{\kern 1pt} '{\text{/}}\lambda {\kern 1pt} ')$, где $\alpha _{3}^{'}$ – коэффициент теплообмена поверхности теплоизоляции с внешней средой, имеющей температуру $T_{3}^{*}$.
Если в качестве материала оболочки бака использовать достаточно теплопроводный материал (например, алюминиевый сплав АМг6 с коэффициентом теплопрводности $\lambda = 122$ Вт/(м ⋅ К) [10]), то изменением температуры по толщине $h$ оболочки можно пренебречь. Тогда при наличии осевой симметрии конструкции бака достаточно учитывать зависимость температуры оболочки лишь от времени $t$ и координаты $x$, отсчитываемой от уровня криогенной жидкости вверх вдоль образующей оболочки. В окрестности уровня жидкости распределение температуры в оболочке удовлетворяет условиям непрерывности [7]
(1)
${{T}_{1}}(0,t) = {{T}_{2}}(0,t),\quad {{\left. {\frac{{\partial {{T}_{1}}(x,t)}}{{\partial x}}} \right|}_{{x = 0}}} = {{\left. {\frac{{\partial {{T}_{2}}(x,t)}}{{\partial x}}} \right|}_{{x = 0}}},$При наличии внешнего слоя теплоизоляции с удельной объемной теплоемкостью $c{\kern 1pt} '$ справедливо неравенство $h{\text{'}}c{\text{'}} \ll hc$, где $c$ – удельная объемная теплоемкость материала оболочки. Поэтому влиянием теплоемкости этого слоя можно пренебречь. Тогда при движении уровня жидкости вверх со скоростью $v$ распределения температуры участков оболочки будут удовлетворять, согласно закону сохранения тепловой энергии, дифференциальным уравнениям [8] (аргументы функций ${{T}_{1}}(x,t)$ и ${{T}_{2}}(x,t)$ опущены)
(2)
$c{\kern 1pt} \frac{{\partial {{T}_{1}}}}{{\partial t}} = \lambda \frac{{{{\partial }^{2}}{{T}_{1}}}}{{\partial {{x}^{2}}}} + cv{\kern 1pt} \frac{{\partial {{T}_{1}}}}{{\partial x}} + {{\alpha }_{1}}(T_{1}^{*} - {{T}_{1}}){\text{/}}h + {{\alpha }_{3}}(T_{3}^{*} - {{T}_{1}}){\text{/}}h,$(3)
$c\frac{{\partial {{T}_{2}}}}{{\partial t}} = \lambda \frac{{{{\partial }^{2}}{{T}_{2}}}}{{\partial {{x}^{2}}}} + cv\frac{{\partial {{T}_{2}}}}{{\partial x}} + {{\alpha }_{2}}(T_{2}^{*} - {{T}_{2}}){\text{/}}h + {{\alpha }_{3}}(T_{3}^{*} - {{T}_{2}}){\text{/}}h.$По мере удаления от уровня жидкости неравномерность распределения температуры вдоль образующей оболочки будет убывать. Тогда помимо условий непрерывности (1) для уравнений (2) и (3) можно задать еще два граничных условия
Если условия теплообмена не зависят от времени, а уровень жидкости неподвижен $(v = 0)$, то распределение температуры в оболочке будет стремиться при $t \to \infty $ к некоторому установившемуся температурному состоянию, для которого из уравнений (2) и (3) при $t \to \infty $ следуют условия
(4)
$\mathop {lim}\limits_{x \to - \infty } {{T}_{1}}(x,t) = {{\tilde {T}}_{1}} = \frac{{{{\alpha }_{1}}T_{1}^{*} + {{\alpha }_{3}}T_{3}^{*}}}{{{{\alpha }_{1}} + {{\alpha }_{3}}}},\quad \mathop {lim}\limits_{x \to - \infty } {{T}_{2}}(x,t) = {{\tilde {T}}_{2}} = \frac{{{{\alpha }_{2}}T_{1}^{*} + {{\alpha }_{3}}T_{3}^{*}}}{{{{\alpha }_{2}} + {{\alpha }_{3}}}}.$Это состояние в момент начала движения уровня жидкости определит начальное условие для уравнений (2) и (3) при $v > 0$.
При независимых от $t$ и $x$ параметрах теплообмена и постоянной скорости движения уровня жидкости левые части соотношений (2), (3) при $t \to \infty $ стремятся к нулю. В этом случае соотношения (2), (3) перейдут в обыкновенные дифференциальные уравнения второго порядка [9]
(5)
$\frac{{{{T}_{1}}(x) - {{{\tilde {T}}}_{1}}}}{{{{{\tilde {T}}}_{2}} - {{{\tilde {T}}}_{1}}}} = \frac{{exp({{K}_{1}}x)}}{{1 + {{K}_{1}}{\text{/}}{{K}_{2}}}},\quad x < 0;\quad \frac{{{{T}_{2}}(x) - {{{\tilde {T}}}_{1}}}}{{{{{\tilde {T}}}_{2}} - {{{\tilde {T}}}_{1}}}} = 1 - \frac{{exp( - {{K}_{2}}x)}}{{1 + {{K}_{2}}{\text{/}}{{K}_{1}}}},\quad x > 0.$Из формул (5) в силу $exp( - 3) \approx 0.05$ следует, что неравномерность квазистационарного распределения температуры в первом приближении достаточно учитывать лишь при $x \in ( - 3{\text{/}}{{K}_{1}},3{\text{/}}{{K}_{2}})$. За пределами этого промежутка температуру смоченной и несмоченной частей оболочки можно принять равными ${{\tilde {T}}_{1}}$ и ${{\tilde {T}}_{2}}$ соответственно. При $v \approx 10$ мм/с переход от однородного начального распределения температуры к квазистационарному в оболочке толщиной не более 10 мм, выполненной из сплава АМг6 (при этом $2hV \approx 1$), можно считать практически завершенным за период времени, не превышающий 10 с [9].
При неподвижном уровне жидкости температуру всего смоченного участка оболочки в первом приближении можно положить равной температуре жидкости, т.е. принять ${{T}_{1}}(x) \approx {{\tilde {T}}_{1}} \approx T_{1}^{*}$ [7], т.к. отношение ${{\alpha }_{1}}{\text{/}}{{\alpha }_{2}}$ велико благодаря высокой интенсивности теплообмена на смоченной поверхности оболочки, а значения ${{\alpha }_{2}}$ и ${{\alpha }_{3}}$ обычно сопоставимы между собой [8].
На рис. 1 приведены распределения безразмерной температуры Θ(ξ) = $\Theta (\xi ) = (T(\xi ) - {{\tilde {T}}_{1}}){\text{/}}({{\tilde {T}}_{2}} - {{\tilde {T}}_{1}})$ в зависимости от безразмерной координаты $\xi = x{\text{/}}h$ при фиксированных значениях параметров ${{\tilde {H}}_{1}} = {{H}_{1}}h = 1$, ${{\tilde {H}}_{2}} = {{H}_{1}}h = 0.02$ и различных значениях параметра $\tilde {V} = Vh$, в том числе при $\tilde {V} = 0$, что соответствует неподвижному уровню жидкости (5).
Рис. 1.
Распределения безразмерной температуры $\Theta (\theta )$ в зависимости от безразмерной координаты $\theta $ при ${{\tilde {H}}_{1}} = 1$, ${{\tilde {H}}_{2}} = 0.02$ и различных значениях параметра $\tilde {V}$: 1 – 1; 2 – 0.5; 3 – 0.25 и 4 – 0.
![](/issues/pminm/2019/vol_2019/iss_7/PMiNM1907016Zarubin/PMiNM1907016Zarubin-F1.gif)
3. Радиальное перемещение оболочки. В силу температурного расширения материала оболочки неравномерное распределение температуры вдоль образующей криогенного бака (5) вызывает радиальное перемещение оболочки. Непосредственным влиянием внутреннего давления в баке, избыточного по отношению к давлению окружающей среды, можно пренебречь, поскольку при заполнении бака с открытым дренажным клапаном [5] давление ${{p}_{0}}$ в полости бака над уровнем жидкости может превышать давление окружающего бак воздуха не более, чем на 0.2 МПа [11]. Гидростатической составляющей давления на смоченной части оболочки водородного бака также можно пренебречь в силу малой плотности жидкого водорода. Для бака с жидкими кислородом или метаном, плотности которых соответственно 1136 кг/м3 и 426 кг/м3 [12] при стандартном атмосферном давлении 101.3 кПа и температуре кипения, это давление более существенно, но его влияние будет заметно лишь на глубине за пределами области краевого эффекта в окрестности уровня жидкости.
При анализе радиального перемещения оболочки помимо распределения температуры вдоль образующей бака достаточно учитывать лишь влияние осевой сжимающей силы $P$ и разгружающий эффект давления ${{p}_{0}}$, определяемый в первом приближении силой ${{P}_{1}} = \pi \,{{R}^{2}}{{p}_{0}}$, где $R$ – радиус средней поверхности оболочки.
При линейно-упругом поведении материала осесимметричной оболочки отнесенные к ее средней поверхности погонные усилия ${{N}_{\theta }}$ и ${{N}_{y}}$ соответственно в меридиональном и окружном направлениях представим в виде [7, 9]
(6)
${{N}_{\theta }} = Eh\frac{{{{\varepsilon }_{\theta }} + \nu {{\varepsilon }_{y}} - (1 + \nu ){{\varepsilon }_{T}}}}{{1 - {{\nu }^{2}}}},\quad {{N}_{y}} = Eh\frac{{{{\varepsilon }_{y}} + \nu {{\varepsilon }_{\theta }} - (1 + \nu ){{\varepsilon }_{T}}}}{{1 - {{\nu }^{2}}}},$Приложенной к верхнему торцу оболочки осевой сжимающей силе $P$ и разгружающей силе ${{P}_{1}}$ отвечает сжимающая погонная нагрузка $N = (P - {{P}_{1}}){\text{/}}(2\pi R)$. С учетом этой нагрузки условие равновесия оболочки в радиальном направлении в линейном приближении имеет вид [10]
(7)
$\frac{{{{d}^{2}}{{M}_{x}}}}{{d{{x}^{2}}}} + N{\kern 1pt} \frac{{{{d}^{2}}w}}{{d{{x}^{2}}}} + \frac{{{{N}_{y}}}}{R} = 0,$Исключив из соотношений (6) деформацию ${{\varepsilon }_{\theta }}$ получим
После подстановки формул для ${{M}_{x}}$ и ${{N}_{y}}$ в уравнение (7) с учетом приближенного равенства ${{N}_{\theta }} \approx - N$ запишем
(9)
$\frac{{{{d}^{4}}w}}{{d{{x}^{4}}}} + 2({{\beta }^{2}} - {{\gamma }^{2}})\frac{{{{d}^{2}}w}}{{d{{x}^{2}}}} + {{({{\beta }^{2}} + {{\gamma }^{2}})}^{2}}w = \frac{{\nu N}}{{DR}} + 12{{\alpha }_{T}}\frac{{T - {{{\tilde {T}}}_{1}}}}{{{{h}^{2}}R}}.$Характеристическое уравнение, соответствующее однородной части уравнения (9), имеет две пары комплексно сопряженных корней $\gamma \pm \beta i$ и $ - \gamma \pm \beta i.$ Вне зоны краевого эффекта функция $w\,(x)$ должна совпадать с частным решением $w_{{}}^{*}(x)$ уравнения (9), а общее решение соответствующего однородного уравнения при ${\text{|}}x{\text{|}} \to \infty $ должно убывать по абсолютному значению и стремиться к нулю. Этому условию при $x < 0$ удовлетворяет лишь первая пара корней и общим решением однородного уравнения будет соотношение
а при $x > 0$ удовлетворяет лишь вторая пара корней, так что общее решение примет видОбозначенные литерой $C$ постоянные коэффициенты в равенствах (10) и (11) определим с учетом частных решений уравнения (9) из условия непрерывности при $x = 0$ функции $w\,(x)$ и первых трех ее производных.
При квазистационарном распределении температуры, описываемым соотношениями (5), частными решениями уравнения (9) будут при $x < 0$
(12)
$w_{1}^{*}(x) = \frac{{\nu NR}}{{Eh}} + \frac{{{{w}_{T}}exp({{K}_{1}}x)}}{{{{h}^{2}}{{R}^{2}}B'(1 + {{K}_{1}}{\text{/}}{{K}_{2}})}},$(13)
$w_{2}^{*}(x) = \frac{{\nu NR}}{{Eh}} + \frac{{{{w}_{T}}}}{{1 - {{\nu }^{2}}}} - \frac{{{{w}_{T}}exp( - {{K}_{2}}x)}}{{{{h}^{2}}{{R}^{2}}B{\kern 1pt} '{\kern 1pt} '(1 + {{K}_{2}}{\text{/}}{{K}_{1}})}},$Общее решение ${{w}^{{^{ \circ }}}}(x) + {{w}^{*}}(x)$ уравнения (9) и первые три производные этого решения должны удовлетворять при $x = 0$ условиям непрерывности, что с учетом равенств (10)–(13) приводит к системе линейных алгебраических уравнений
(14)
$\left. {\begin{array}{*{20}{l}} {C_{2}^{'} - C_{1}^{'} = {{A}_{1}},} \\ {\gamma C_{2}^{'} + \beta C_{2}^{{''}} + \gamma C_{1}^{'} - \beta C_{1}^{{''}} = {{B}_{1}},} \\ {({{\gamma }^{2}} - {{\beta }^{2}})C_{2}^{'} + 2\beta \gamma C_{2}^{{''}} + ({{\beta }^{2}} - {{\gamma }^{2}})C_{1}^{'} + 2\beta \gamma C_{1}^{{''}} = {{A}_{2}},} \\ {({{\gamma }^{3}} - 3\gamma {{\beta }^{2}})C_{2}^{'} - ({{\beta }^{3}} - 3\beta {{\gamma }^{2}})C_{2}^{{''}}({{\gamma }^{3}} - 3\gamma {{\beta }^{2}})C_{1}^{'} + ({{\beta }^{3}} - 3\beta {{\gamma }^{2}})C_{1}^{{''}} = {{B}_{2}},} \end{array}} \right\}$Решением системы уравнений (14) будут равенства
В итоге для неоднородного уравнения (9) общее решение ${{w}_{1}}(x) = w_{1}^{{^{ \circ }}} + w_{1}^{*}(x)$ при x < 0 определяем равенствами (10), (12), а общее решение ${{w}_{2}}(x) = w_{2}^{{^{ \circ }}} + w_{2}^{*}(x)$ при x > 0 – равенствами (11), (13).
Полученные соотношения использованы для построения графиков зависимостей (рис. 2) от координаты $x$ температурного приращения $w_{T}^{'} = {{\varepsilon }_{T}}R$ радиуса оболочки, ее полного радиального перемещения $w$ и осевого погонного изгибающего момента $M$ при заполнении криогенного бака жидким кислородом. В расчетах принято: h = 0.005 м; R = 1 м; ${{\tilde {T}}_{2}} - {{\tilde {T}}_{1}} = 182$ К; $v = 0.01$ м/с; ${{H}_{1}}h = 0.5$; ${{H}_{2}}h = 0.01$, $N = 3 \cdot {{10}^{5}}$ Н/м; материал оболочки – алюминиевый сплав АМг6, для которого $\lambda = 122$ Вт/(м ⋅ К), $a = 5 \cdot {{10}^{{ - 5}}}$ м2/с, E = 71 ГПа, $\nu = 0.31$ и ${{\alpha }_{T}} = 24.7 \cdot {{10}^{{ - 6}}}$ 1/К.
Рис. 2.
Графики зависимости от осевой координаты $x$ температурного приращения $w_{T}^{'} = {{\varepsilon }_{T}}R$ радиуса оболочки толщиной h = 0.005 м, ее полного радиального перемещения $w$ и осевого погонного изгибающего момента Mx (кривые 1, 2 и 3 соответственно) при скорости $v$ = 0.01 м/с подъема уровня жидкости.
![](/issues/pminm/2019/vol_2019/iss_7/PMiNM1907016Zarubin/PMiNM1907016Zarubin-F2.gif)
4. Напряженное состояние оболочки. Определитель матрицы системы уравнений (14) обращается в нуль при значении ${{N}_{{cr}}} = (2{\text{/}}R)\sqrt {DEh} $, т.е. при верхнем значении критического погонного усилия, соответствующего потере устойчивости достаточно длинной цилиндрической линейно-упругой оболочки, сохраняющей в докритическом напряженном состоянии прямолинейность образующей [13, 14]. Значению ${{N}_{{cr}}}$ отвечает критическое значение осевого сжимающего напряжения ${{\sigma }_{{cr}}} \approx 0.605{\kern 1pt} Eh{\text{/}}R$ [7], т.е. эта формула применима лишь при условии $h{\text{/}}R < 1.65{\kern 1pt} {{\sigma }_{Y}}{\text{/}}E$, где ${{\sigma }_{Y}}$ – предел текучести материала оболочки. Для сплава АМг6 при температуре 78 К ${{\sigma }_{Y}} = 185$ МПа и $E = 71$ ГПа [10], т.е. отношение $h{\text{/}}R$ не должно превышать 0.0043, что соответствует весьма тонкостенной оболочке.
Искривление образующей оболочки в зоне подвижного краевого эффекта приводит к возрастанию локальных абсолютных значений напряжений, которые могут превысить предел текучести материала оболочки. При этом известные подходы анализа упругой устойчивости оболочки станут неприменимы. Пластическое деформирование материала оболочки в зоне краевого эффекта, возникающего на стадии заполнения криогенного бака, вызовет остаточное формоизменение, влияющее на несущую способность конструкции на последующих этапах ее эксплуатации. В связи с этим актуальным является анализ напряженного состояния оболочки бака в зоне краевого эффекта с точки зрения возможности выхода материала за пределы области линейной упругости.
В рамках допущения о плоском напряженном состоянии оболочки примем, что линейно-упругое поведение материала согласно критерию Мизеса [7] имеет место при условии
(15)
${{\sigma }_{Y}} > {{\sigma }_{i}} = \sqrt {\sigma _{\theta }^{2} + \sigma _{y}^{2} - {{\sigma }_{\theta }}{{\sigma }_{y}}} ,$Изменение радиального перемещения средней поверхности оболочки в зоне подвижного краевого эффекта (рис. 2) позволяет заключить, что при фиксированном сжимающем погонном усилии $N$ и примерно равных ближайших к уровню жидкости экстремальных значениях осевого погонного изгибающего момента наибольшие значения ${{\sigma }_{i}}$ будут соответствовать окрестностям точек с абсциссами этих экстремумов на сжатой наружной поверхности смоченной части оболочки и на ее внутренней поверхности несмоченного участка. Осевое сжимающее напряжение в окрестностях этих точек, вызванное осевым погонным изгибающим моментом ${{M}_{x}}$, равно ${{\sigma }_{M}} = - 6\left| {{{M}_{x}}} \right|{\text{/}}{{h}^{2}}$, а осевой сжимающей силой – ${{\sigma }_{N}} = - N{\text{/}}h$. Следовательно, совпадающее в данном случае с меридиональным осевое напряжение, наибольшее по абсолютному значению, определяет равенство $\sigma _{x}^{*} = - \left| {6{{M}_{x}}} \right|{\text{/}}{{h}^{2}} - N{\text{/}}h$. C погонным изгибающим моментом ${{M}_{y}} = \nu {{M}_{x}}$ в окружном направлении связано сжимающее напряжение $\sigma _{y}^{'} = - 6\nu {\text{|}}{{M}_{x}}{\text{|/}}{{h}^{2}}$, а с радиальным перемещением $w$ – напряжение $\sigma _{y}^{{''}} = {{N}_{y}}{\text{/}}h$ (в формуле (8) для ${{N}_{y}}$ следует ${{N}_{\theta }}$ заменить на $ - N$). Значение ${{N}_{y}}$ может быть как положительным, так и отрицательным.
В линейном приближении, пренебрегая величиной ${{(dw{\text{/}}dx)}^{2}}{\text{/}}2$ по сравнению с единицей, в формуле (15) можно заменить напряжения ${{\sigma }_{\theta }}$ и ${{\sigma }_{y}}$ соответственно на $\sigma _{x}^{*}$ и $\sigma _{y}^{*} = \sigma _{y}^{'} + \sigma _{y}^{{''}}$ и использовать при вычислениях формулу
(16)
${{\sigma }_{i}} = \sqrt {{{{(\sigma _{x}^{*})}}^{2}} + {{{(\sigma _{y}^{*})}}^{2}} - \sigma _{x}^{*}\sigma _{y}^{*}} .$Аналогичные подсчеты напряжений в окрестности точки ближайщего к уровню жидкости экстремума осевого погонного изгибающего момента на несмоченном участке оболочки приводят к несколько меньшему значению ${{\sigma }_{i}} = 176.5$ МПа. При уменьшении скорости движения уровня жидкости и фиксированном значении погонной силы $N$ интенсивность напряжения ${{\sigma }_{i}}$ возрастает. Основной причиной является более пологое изменение температуры несмоченной части оболочки в окрестности уровня жидкости и, наоборот, более резкое изменение в этой окрестности температуры смоченной части оболочки. Это приводит к увеличению положительного значения ${{N}_{y}}$ и возрастанию ${{\sigma }_{i}}$.
На рис. 3 представлены графики при неподвижном уровне жидкости и прежнем значении $N = 3 \cdot {{10}^{5}}$ Н/м. В этом случае на несмоченной части оболочки наибольшее значение 226.3 МПа интенсивности напряжения превышает предел текучести, но и на смоченной части соответствующее значение 186.1 МПа также несколько выше допустимого. Условие ${{\sigma }_{Y}} > {{\sigma }_{i}}$ на несмоченной части оболочки может быть выполнено лишь при значении $N < 1.74 \cdot {{10}^{5}}$ Н/м.
Рис. 3.
Графики зависимости от осевой координаты $x$ температурного приращения $w_{T}^{'} = {{\varepsilon }_{T}}R$ радиуса оболочки толщиной h = 0.005 м, ее полного радиального перемещения $w$ и осевого погонного изгибающего момента Mx (кривые 1, 2 и 3 соответственно) при неподвижном уровне жидкости.
![](/issues/pminm/2019/vol_2019/iss_7/PMiNM1907016Zarubin/PMiNM1907016Zarubin-F3.gif)
Увеличение толщины оболочки позволяет даже в случае неподвижного уровня жидкости повысить значение $N$ при выполнении условия, определяемого соотношением (15). На рис. 4 приведены графики при толщине h = 0.01 м и прежних остальных исходных данных. Для несмоченной части оболочки выполнению условия ${{\sigma }_{i}} \leqslant {{\sigma }_{Y}} = 185$ МПа соответствует значение $N = 4.27 \cdot {{10}^{5}}$ Н/м. При этом на смоченной части наибольшее значение интенсивности напряжения не превышает 160 МПа. При заполнении бака допустимый уровень осевой сжимающей нагрузки возрастает, но более напряженным становится смоченный участок оболочки вблизи уровня жидкости. При скорости подъема уровня жидкости $v = 0.005$ м/с условие (15) выполнимо, если $N < 6.53 \cdot {{10}^{5}}$ Н/м, а при $v = 0.01$ м/с, если $N < 7.24 \cdot {{10}^{5}}$ Н/м.
Рис. 4.
Графики зависимости от осевой координаты $x$ температурного приращения $w_{T}^{'} = {{\varepsilon }_{T}}R$ радиуса оболочки толщиной h = 0.01 м, ее полного радиального перемещения $w$ и осевого погонного изгибающего момента Mx (кривые 1, 2 и 3 соответственно) при неподвижном уровне жидкости.
![](/issues/pminm/2019/vol_2019/iss_7/PMiNM1907016Zarubin/PMiNM1907016Zarubin-F4.gif)
Заключение. Количественный анализ комплексной математической модели, описывающей термодеформированное состояние подкрепленной стрингерами оболочки вертикально расположенного цилиндрического криогенного бака при его заполнении, позволяет оценить неравномерность распределения температуры в оболочке и ее радиального перемещения в зоне движущегося уровня криогенной жидкости, а также найти экстремальные значения погонного изгибающего момента. Это дает возможность уточнить постановку и решение задачи об устойчивости оболочки бака при действии осевой сжимающей нагрузки.
Полученные соотношения можно применить и для оценки допустимой нагрузки в случае оболочки бака, подкрепленной силовым набором (достаточно часто расположенными стрингерами и кольцевыми ребрами или вафельным подкреплением [15]), если привести характеристики жесткости такой оболочки к ее некоторой условной средней поверхности [16].
Конфликт интересов. Авторы заявляют, что у них нет конфликта интересов.
Список литературы
Альтернативные топлива для двигателей внутреннего сгорания / Под ред. А.А. Александрова, В.А. Маркова. М.: ООО НИЦ “Инженер”, ООО “Онико-М”, 2012. С. 791.
Рачук В.С., Коротеев А.С., Каторгин Б.И. и др. Обоснование выбора компонентов ракетного топлива для двигательных установок первой ступени многоразовой ракетно-космической системы // Ж. Авиакосмическая техника и технология. 2010. № 1. С. 19.
Феодосьев В.И. Основы техники ракетного полета. М.: Наука, 1979. С. 496.
Ковалев Б.К. Развитие ракетно-космических систем выведения. М.: Изд-во МГТУ им. Н.Э. Баумана, 2014. С. 398.
Архаров А.М., Кунис И.Д. Криогенные заправочнные системы стартовых ракетно-космических комплексов / Под ред. И.В. Бармина. М.: Изд-во МГТУ им. Н.Э. Баумана, 2006. С. 252.
External TankThermal Protection System. https://ru.scribd.com/document/349711/NASA-114022main-TPS-FS (26.01.2019).
Балабух Л.И., Колесников К.С., Зарубин В.С., Алфутов Н.А., Усюкин В.И., Чижов В.Ф. Основы строительной механики ракет. Москва, Высшая школа, 1969. С. 496.
Зарубин В.С., Зимин В.Н., Кувыркин Г.Н. Математическое моделирование температурного состояния оболочки цилиндрической криогенной емкости при заполнении и опорожнении // Математика и математическое моделирование. 2015. № 6. https://doi.org/10.7463/mathm0615.08829350
Александров А.А., Зарубин В.С., Зимин В.Н. Термонапряженное состояние оболочки цилиндрического криогенного бака при его заполнении // Изв. вузов. Машиностроение. 2016. № 4. С. 5.
Конструкционные материалы: Справочник / Под общ. ред. Б.Н. Арзамасова. М.: Машиностроение, 1990. С. 688.
Александров A.A., Бармин И.В., Кунис И.Д., Чугунков В.В. Особенности создания и развития криогенных систем ракетно-космических стартовых комплексов “Союз” // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. 2016. № 2. С. 7. https://doi.org/10.18698/0236-3941-2016-2-7-27
Справочник по физико-техническим основам криогеники / Под ред. М.П. Малкова. М.: Энергоатомиздат, 1985. С. 432.
Алфутов Н.А. Основы расчета на устойчивость упругих систем. 2-е изд. М.: Машиностроение, 1991. С. 336.
Алфутов Н.А., Колесников К.С. Устойчивость движения и равновесия. М.: Изд-во МГТУ им. Н.Э. Баумана, 2003. С. 252.
Лизин В.Т., Пяткин В.А. Проектирование тонкостенных конструкций. 3-е изд. М.: Машиностроение, 1994. С. 384.
Кан С.Н., Бырсан К.Е., Алифанова О.А., Бутенко Ю.И., Ингульцев В.Л. Устойчивость оболочек. Харьков: Изд-во Харьковского ун-та, 1970. С. 154.
Дополнительные материалы отсутствуют.
Инструменты
Проблемы машиностроения и надежности машин