Изменение формы вертикально стоящего упругого резервуара с жидкостью
- Авторы: Волобуев А.Н.1, Краснов С.В.1, Адыширин-Заде К.А.1, Антипова Т.А.1, Александрова Н.Н.1
-
Учреждения:
- Самарский государственный медицинский университет
- Выпуск: Том 25, № 1 (2022)
- Страницы: 80-86
- Раздел: Статьи
- URL: https://journals.ssau.ru/pwp/article/view/10152
- DOI: https://doi.org/10.18469/1810-3189.2022.25.1.80-86
- ID: 10152
Цитировать
Полный текст
Аннотация
Обоснован вид уравнения гидростатики для упругого, вертикально стоящего резервуара, например топливного бака ракеты на стартовом столе. Уравнение гидростатики получено на основе уточненного уравнения Бернулли. Обоснование проведено с помощью формулы Лапласа для давления под упругой поверхностью жидкости, которая может возникнуть как за счет сил поверхностного натяжения, так и за счет упругой тонкостенной оболочки, как в настоящей задаче. Найдена форма вертикально расположенного упругого резервуара с жестким дном и жестким верхним обручем, заполненного покоящейся жидкостью. Показано, что необходимо использовать особую запись закона Гука для получения формы резервуара. Проведен анализ этой формы. Показано распределение по высоте резервуара гидростатического давления, объемной плотности энергии растянутой упругой стенки, а также суммы этих величин. Найдено гидростатическое давление, на уровне которого возникает максимальное увеличение площади упругого резервуара.
Полный текст
Введение
В процессе исследования покоящейся жидкости в вертикально стоящих упругих резервуарах, часто возникает вопрос о форме, которую принимает резервуар. Например, как изменяется форма заправленной топливом ракеты на стартовом столе, какова роль гидростатического давления в изменении формы грудного отдела аорты человека и т. д.
Решение этой задачи в классическом труде [1] основывается на т. н. «мембранной теории». В этой теории рассматриваются силы и моменты сил в элементе оболочки вертикально стоящего резервуара, заполненного покоящейся жидкостью. Составляется дифференциальное уравнение 4-го порядка для радиальной деформации резервуара. Несмотря на достаточно большую сложность и громоздкость преобразований, результат частного решения дифференциального уравнения оказывается достаточно тривиальным: относительная деформация площади сечения резервуара на определенном уровне пропорциональна гидростатическому давлению на этом уровне. В дальнейшем будет показано, что этот результат слишком грубый для правильной оценки формы резервуара.
Используется также другой подход, основанный на т. н. «безмоментной теории» [2]. В этой теории используется уравнение Лапласа, которое позволяет приближенно оценить напряженное состояние в стенках резервуара под действием внутреннего давления.
Достаточно разнообразно данная задача исследуется в биомеханике, например, [3–5]. Задача исследуется с разной степенью сложности, но основное внимание часто уделяется биомеханическим особенностям стенок аорты и других кровеносных сосудов.
Целью настоящей работы является исследование способа нахождения формы вертикально расположенного упругого цилиндрического резервуара, заполненного покоящейся жидкостью, определение условий, необходимых для достижения корректного результата решения поставленной задачи.
Уравнение Бернулли для упругого трубопровода
Рассмотрим модельную задачу нахождения формы открытого вертикального цилиндрического резервуара с упругими стенками и жестким дном, рис. 1, в который налита жидкость (вода).
Рис. 1. Изменение формы упругого резервуара
Fig. 1. Changing the shape of an elastic reservoir
Пусть в сечении 1, у дна резервуара координата а гидростатическое давление На поверхности жидкости в сечении 2 находится начало координат где давление Используем давление, избыточное над атмосферным. Так как жидкость находится в состоянии покоя, то ее вязкость роли не играет.
На дне резервуара гидростатическое давление равно:
(1)
где - ускорение свободного падения, – плотность жидкости.
Распределение гидростатического давления по высоте резервуара имеет вид:
(2)
Распределение гидростатического давления по высоте упругого резервуара, согласно формуле (2) носит линейный характер. Оно, очевидно, не должно зависеть от того жесткие у резервуара стенки или упругие.
Для учета упругости стенок резервуара найдем уравнение Бернулли в упругом трубопроводе с движущейся в нем жидкостью.
В [4; 6; 7] рекомендуется для волнового процесса в упругом трубопроводе использовать уравнение импульса в виде:
(3)
где V и W – продольная и поперечная составляющие скорости жидкости, r – радиальная координата, Х – продольная координата, t – время, S – поперечное сечение тонкостенного упругого трубопровода. Вязкость жидкости не учитываем.
Преобразуем уравнение (3) к виду:
(4)
При записи (4) использован закон Гука для упругого трубопровода в виде [3]:
(5)
где D – упругость стенки трубопровода. Формула (5) характеризует связь распределенной реакции упругой стенки трубопровода Р и его площади сечения S, поэтому знаки дифференциалов разные. Эта форма записи закона Гука фактически получена в [1] на основе «мембранной теории» при решении уравнения для радиальной деформации тонкостенного упругого резервуара, если принять [8], где – толщина стенки резервуара, d – его диаметр, Е – модуль упругости материала оболочки резервуара.
Принято также, что вихрей в потоке не образуется, течение потенциальное [9], т. е. где V – вектор скорости, следовательно,
Преобразуем уравнение (4), считая жидкость несжимаемой, т. е.
(6)
Если течение жидкости стационарное, т. е. то уравнение (6) можно проинтегрировать:
(7)
Проанализируем более подробно причину возникновения последнего слагаемого в (7). Покажем, что оно является следствием закона Гука в форме (5), который для удобства преобразований запишем в виде где – начало отсчета давления.
Рис. 2. Участок оболочки вертикального упругого резервуара
Fig. 2. Section of the shell of a vertical elastic tank
Исследуем более детально связь между упругостью стенок резервуара D и модулем упругости вещества стенок Е, а также геометрическими размерами резервуара. Разрежем мысленно вдоль небольшой участок резервуара длиной l и средним диаметром d, рис. 2.
Сила давления, избыточного над растягивающая резервуар, равна силе сопротивления его стенок:
где σ - механические напряжения в стенке резервуара толщиной δ, – механические напряжения в стенке резервуара при давлении
Следовательно, давление в резервуаре равно Полученное соотношение носит название формулы Лапласа. Однако, формула Лапласа принципиально неточная. Термодинамический анализ показывает, что более точная формула [10]. Разлагая экспоненту в ряд, получаем
Из первого слагаемого в правой части, как первого приближения (формула Лапласа), следует
Учитывая закон Гука где - относительная деформация длины окружности L стенки резервуара d – изменение диаметра резервуара при его растяжении, получаем:
Найдем связь между относительным изменением площади сечения резервуара и относительной деформацией его диаметра Учитывая связь между площадью сечения и диаметром резервуара находим производную следовательно,
Поэтому для давления, избыточного над получаем:
или
Полученный результат, при показывает, что если использовать закон Гука в виде (5), то более правильно вместо давления Р использовать величину как это принято в формуле (7). Данный вывод не связан с течением жидкости.
В случае, если на текущую по упругому трубопроводу жидкость действует гравитационная сила, в левую часть уравнения (7) необходимо добавить объемную плотность потенциальной энергии жидкости в гравитационном поле. В этом случае уравнение (7) приобретет вид:
(8)
где – высота рассматриваемого элемента жидкости над уровнем отсчета.
Величина представляет собой объемную плотность энергии растянутой упругой стенки. Величину Рс - можно отождествить с некоторым давлением.
Полагая статическое давление:
(9)
получаем стандартную форму уравнения Бернулли. Будем использовать термин давление для величины Р, учитывая, что оно не тождественно измеряемому статическому давлению.
Таким образом, уравнение Бернулли для упругого трубопровода ничем не отличается от такового для жесткого трубопровода. Так же как и в жестком трубопроводе, сумма статического гидростатического и динамического давлений в каждом поперечном сечении упругого трубопровода остается постоянной.
Однако в уравнение Бернулли (8) входит не статическое давление а давление которое равно сумме статического давления и некоторого давления, которое равно объемной плотности энергии растянутой упругой стенки. Отметим, что в жестком трубопроводе при давление
Найдем распределение всех давлений по высоте упругого резервуара.
Используя формулы (2) и (9), с учетом получим:
(10)
В отличие от уравнения (8), координата h направлена сверху вниз, рис. 1.
Решая квадратное уравнение (10) относительно давления Р, найдем:
(11)
Знак плюс перед корнем неприемлем, т. к. при (жесткие стенки резервуара) величина Р должна стремиться к гидростатическому давлению
Давление, определяемое объемной плотностью энергии растянутой упругой стенки резервуара, найдем по формуле (9):
(12)
На рис. 3 показано распределение давлений по высоте резервуара с упругими стенками. Для расчета приняты следующие параметры: плотность жидкости (воды) упругость стенки резервуара (для наглядности, как у резины) Эта величина зависит от диаметра резервуара, толщины его стенки, модуля Юнга вещества стенки [8].
Рис. 3. Распределение давлений по высоте вертикального резервуара с упругими стенками: прямая 1 – гидростатическое давление 2 – давление 3 – объемная плотность энергии растянутой упругой стенки резервуара –
Fig. 3. Pressure distribution along the height of a vertical tank with elastic walls: straight line 1 – hydrostatic pressure 2 – pressure 3 – volumetric energy density of the stretched elastic tank wall –
Все давления равны нулю на поверхности жидкости в упругом резервуаре и достигают максимального значения на дне резервуара.
Найдем форму стенки вертикального упругого резервуара, в который налита вода. Пусть высота жидкости в резервуаре H, а площадь жесткого дна резервуара SH.
Заметим, что точность формул (10), (11) и (12) определяется справедливостью использованного закона Гука в виде (5).
Нахождение формы вертикального резервуара с упругими стенками
Применение закона Гука для упругой стенки резервуара в виде (5) при решении поставленной задачи невозможно, т. к. эта формула слишком грубо описывает зависимость давления и площади поперечного сечения резервуара. Действительно, если рассматривать формулу (5) как дифференциальное уравнение первого порядка, то в результате интегрирования возникает только одна постоянная интегрирования. Поэтому невозможно удовлетворить сразу двум граничным условиям: площадь поверхности жидкости в области жесткого верхнего обруча и площадь жесткого дна упругого резервуара
Применим закон Гука в виде, аналогичном [2], где подобная запись используется при анализе относительного удлинения стержня:
(13)
где - изменение дифференциала площади сечения резервуара. Знак минус определяется тем же, что и в случае формулы (5), – это реакция упругой стенки резервуара на жидкость.
Учитывая зависимость площади сечения резервуара от высоты перейдем от приращений к дифференциалам:
(14)
Преобразуем формулу (14):
(15)
Уравнение (15) можно один раз проинтегрировать:
(16)
где - постоянная интегрирования.
Следовательно:
(17)
Используем связь между давлением Р и высотой h в виде (10). Нужно отметить, что формула (10) найдена из более простой формы закона Гука (5). Поэтому на данном этапе в проводимый анализ вносится некоторое приближение. Это приближение связано с тем, что решение уравнения импульса (3) совместно с законом Гука (13) в аналитическом виде затруднительно. Изменением высоты жидкости при изменении формы резервуара пренебрегаем, что предполагает небольшую деформацию формы резервуара.
Дифференцируя (10), имеем:
(18)
Умножив (17) на (18), получим:
(19)
Интегрируя (19), находим:
(20)
где - вторая постоянная интегрирования.
Если площадь поверхности жидкости в упругом резервуаре задана за счет жесткого обруча то, учитывая давление на поверхности жидкости найдем Следовательно, формулу (20) можно записать в виде:
(21)
Кроме того, учитывая площадь жесткого дна резервуара найдем:
(22)
где РН - давление Р на дне резервуара.
Поделив (22) на (21), получим:
(23)
В технике часто возникает задача определения деформации высокого вертикально стоящего резервуара с периодическими или непериодическими подкреплениями оболочки горизонтальными круговыми внутренними обручами с заданной площадью Формула (23) позволяет оценить изменение формы оболочки на каждом участке такого резервуара. В этом случае постоянные и будет изменяться от участка к участку. Эти постоянные определяются по формуле (20) последовательно, начиная с верхнего участка, в соответствии с давлениями на верхней и нижней границах участков.
Проведем анализ полученной зависимости (23). Найдем, при каком давлении Р площадь сечения резервуара будет максимальной. Находя производную от (23) и приравнивая ее к нулю, имеем, что при давлении площадь сечения резервуара будет максимальной. Из формулы (23) находим:
(24)
Формула (24) фактически характеризует изгибные свойства стенок упругого вертикального резервуара и позволяет найти в расчете связь величин площадей и
Рис. 4. Зависимость относительной площади поперечного сечения резервуара от безразмерной высоты жидкости в нем
Fig. 4. Dependence of the relative cross-sectional area of the tank on the dimensionless height of the liquid in it
Приближенно, при достаточно жесткой стенке резервуара, т. е. при достаточно большой величине D, на данном заключительном этапе анализа, чтобы не увеличивать громоздкость формул, примем самую простую зависимость а именно (10) при условии Считаем, что давление Р пропорционально высоте жидкости h, отсчитанной от поверхности ко дну резервуара, Фактически, мы приравняли давление Р к гидростатическому давлению
Для расчетного примера положим что где - положение максимальной площади поперечного сечения упругого резервуара. Таким образом, мы приняли, что максимальная площадь поперечного сечения упругого резервуара возникает на высоте 1 / 3 от его дна.
Следовательно:
(25)
Учитывая, что для данного случая найдем связь высоты жидкости в резервуаре и упругости его стенок
При практических расчетах порядок анализа обычно иной. По заданным параметрам: высоте жидкости в резервуаре Н и упругости его стенок D находится положение максимальной площади сечения резервуара. Кроме того, в практических расчетах предположение вряд ли допустимо. По-видимому, для зависимости допустимо использовать также приближенную формулу (11), однако это делает выкладки значительно более громоздкими.
Учитывая и записываем формулу (25) в виде:
(26)
На рис. 4 показан график зависимости относительной площади сечения упругого резервуара от величины построенный по формуле (26).
Заключение
На основе простейшей формы закона Гука, связывающей давление в упругом трубопроводе по которому течет жидкость и относительную деформацию площади поперечного сечения упругого трубопровода найдено уравнение Бернулли для этого трубопровода. В это уравнение входит объемная плотность энергии растянутой упругой стенки, которую можно отождествить с некоторым давлением.
Показано, что используемый вид закона Гука не позволяет найти форму модельного вертикально стоящего упругого резервуара с жестким дном и жестким верхним обручем, заполненного жидкостью.
Переход к более точной записи закона Гука позволил найти форму такого упругого резервуара. Проведен анализ этой формы. Найдено давление, на уровне которого возникает максимальное увеличение площади упругого резервуара.
Об авторах
Андрей Николаевич Волобуев
Самарский государственный медицинский университет
Email: volobuev47@yandex.ru
ORCID iD: 0000-0001-8624-6981
доктор технических наук, профессор
Россия, СамараСергей Викторович Краснов
Самарский государственный медицинский университет
Email: krasnovtlt@mail.ru
доктор технических наук, профессор, заведующий кафедрой медицинской физики, математики и информатики
Россия, СамараКаира Алимовна Адыширин-Заде
Самарский государственный медицинский университет
Email: adysirinzade67@gmail.com
ORCID iD: 0000-0003-3641-3678
кандидат педагогических наук, доцент кафедры медицинской физики, математики и информатики
Россия, СамараТатьяна Александровна Антипова
Самарский государственный медицинский университет
Email: antipovata81@gmail.com
ORCID iD: 0000-0001-5499-2170
кандидат физико-математических наук, доцент кафедры медицинской физики, математики и информатики
Россия, СамараНаталья Николаевна Александрова
Самарский государственный медицинский университет
Автор, ответственный за переписку.
Email: grecova81@mail.ru
ассистент кафедры медицинской физики, математики и информатики
Россия, СамараСписок литературы
- Тимошенко С.П., Войновский-Кригер С. Пластинки и оболочки. М.: Наука, 1966. С. 535.
- Феодосьев В.И. Сопротивление материалов. М.: Наука, 1986. С. 37.
- Педли T. Гидродинамика крупных кровеносных сосудов. М.: Мир, 1983. 400 с.
- Биофизика для инженеров. Т. 2 / Е.В. Бигдай [и др.]. М.: Горячая линия – Телеком, 2008. С. 126.
- Бранков Г. Основы биомеханики / пер. с болг. М.: Мир, 1981. 256 с.
- Будак Б.М., Самарский А.А., Тихонов А.Н. Сборник задач по математической физике. М.: Наука, 1980. С. 161.
- Волобуев А.Н. Нелинейные особенности течения жидкости в упругом трубопроводе // Математическое моделирование. 2019. Т. 31, № 6. С. 43–54. DOI: https://doi.org/10.1134/S0234087919060030
- Механика кровообращения / К. Каро [и др.]. М.: Мир, 1981. С. 121.
- Ландау Л.Д., Лифшиц Е.М. Гидродинамика. T. 6. М.: Наука, 1986. 736 с.
- Левич В.Г. Курс теоретической физики. Т. 1. М.: Физматгиз, 1962. С. 583.
Дополнительные файлы
![](/img/style/loading.gif)