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

УДК 629.015; 531.391.5

Гаврилов А.А., Власов Ю.Л., Морозов Н.А.

Оренбургский государственный университет E-mail: [email protected]

ДИНАМИЧЕСКАЯ МОДЕЛЬ ОТНОСИТЕЛЬНОГО ПОКОЯ ЖИДКОСТИ В АВТОЦИСТЕРНЕ ПРИ РАВНОМЕРНОМ ДВИЖЕНИИ ПО КРИВОЛИНЕЙНОМУ УЧАСТКУ ДОРОГИ

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

Ключевые слова: автоцистерна, уровень жидкости, центр масс, сила инерции, криволинейное движение.

Перевозка жидкости с помощью автомобильного транспорта часто осуществляется автоцистернами, частично заполненными жидкостью. В процессе перевозки изменяется положение жидкости в автоцистерне, что влияет на движение автоцистерны. Так, в работе [1] рассмотрено влияние давления жидкого груза частично заполненной железнодорожной цистерны при столкновении с препятствием на горизонтальном участке пути. В работе [2] предложена математическая модель автоцистерны с перемещающейся в ней жидкостью при торможении. Тем не менее, в вышеуказанных работах не учитывалось изменение давления жидкости на боковую поверхность цистерны в поворотах.

В работах [3], [4] была предложена динамическая модель равномерного движения автоцистерны с частичным заполнением на криволинейном участке дороги. Однако, в данных работах не учитывалось изменение уровня жидкости в цистерне при движении в поворотах. Изменение уровня жидкости в повороте связано с действием на жидкость сил инерции. Данные силы приводят к перераспределению объема жидкости в цистерне путем смещения жидкости к стенке, удаленной от центра кривизны поворота (рисунок 1).

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

1) получить уравнение свободной поверхности жидкости в частично заполненной автоцистерне;

2) разработать методику определения положения центра тяжести жидкости в повороте;

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

В качестве объекта исследования примем автоцистерну круглого сечения (рисунок 2).

Рисунок 1. Уровень жидкости в движущейся автоцистерне: а) на прямолинейном участке дороги; б) в правом повороте

Рисунок 2. Поперечное сечение автоцистерны в повороте

Гаврилов А.А. и др. Динамическая модель относительного покоя жидкости в автоцистерне...

Площадь сечения жидкости 5 определяется диаметром цистерны Б и долей заполнения. Доля заполнения может быть выражена через половину угла раствора сегмента сечения а, получаемого из выражения: D2

S =—(2п-2а + sin(2а)). (1)

8

Также долю заполнения цистерны можно охарактеризовать высотой уровня жидкости (рисунок 2):

h = D(1 + cos a) 2

(2)

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

Дифференциальное уравнение поверхности уровня жидкости будет иметь вид:

а^у + gdz = 0 , (3)

где а; - переносное нормальное ускорение точки жидкости, м/с2.

g - ускорение свободного падения, м/с2. х, у - декартовы координаты точки жидкости, м/с2.

Переносное нормальное ускорение точки жидкости:

а_ =-

(4)

где V - скорость автоцистерны, м/с,

г - радиус кривизны поворота дороги, м. Решением уравнения (2) будет являться выражение:

z = -^y + с,

(5)

где C - постоянная интегрирования, м.

Определим постоянную интегрирования путем подстановки в уравнение (6) граничных условий: при y=0 будем иметь _ D cosa 2 cos в '

где в - угол наклона жидкости к горизонту.

,2

tgß =

F u

mae

mg r • g

(6)

где F" - сила инерции точки жидкости, Н, G - сила тяжести точки жидкости, Н, m - масса точки жидкости, кг. Тогда:

^ _ D cosa (у)

_ 2 cos в

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

,,2

z = -

-y+:

D cosa

(8)

r •g 2 cos в Из полученного уравнения прямой (8) видно, что на одном и том же повороте разные цистерны, движущиеся с одинаковой скоростью, имеют одинаковый угол наклона уровней находящихся в них жидкостей к горизонту не зависимо от размеров цистерны и объемов залитых жидкостей (рисунок 3).

С другой стороны при изменении скорости движения цистерны или радиуса кривизны поворота дороги угол наклона уровня жидкости к горизонту будет меняться, что вызвано изменением переносного нормального ускорения точек жидкости, а, следовательно, изменением сил инерции (рисунок 4).

Координаты крайних точек уровня жидкости в соответствии с рисунком 2 принимают значения:

zi,2 = D cos(a±ß), y 1,2 = m D sin(a±ß)

(9)

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

Рисунок 3. Уровни жидкости в автоцистерне 0=2м при движении по правому

повороту (г=100 м) со скоростью v=20м/с (1 - заполнение цистерны 20%, 2 - 50%, 3 - 80%)

2

по объему цистерны в повороте может иметь три различные формы в зависимости от первоначального уровня залитой в цистерну жидкости (рисунок 5). Данные формы различаются наличием определенных зон: зоны 1, зоны 2 и зоны 3.

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

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

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

Координаты центров тяжести определяются выражениями:

•'Ъл <>--I Ъ-СОБз--СОБа <>--<> 1 П •'Ъ? <>

•'ъл -1—-—I ъ-собв--соба

1 бшр! 2

П К

ъёуёъ

Г2 п . уйуйъ+81 П, уауёъ+§2 Г2Ауёуёъ

•*ъ1 •>——i ъ-собз--соба <>--<> 1 п г -> * 1 п г о

П - « ПТ-ъ2

•'ъ] -1 ъ-собз--соба

1 ¡щз[ 2

ус =-

П Р

(10)

IП Т-Ъ2

2

2

где 8Р 82 - коэффициенты, характеризующие наличие или отсутствие зон 1 и 3 в формах заполнения (зональные коэффициенты).

Зональные коэффициенты принимают значение 1 при наличие данной зоны и значение 0 -при ее отсутствии. Так при наличии зоны 1 81=1 (рисунки 5 б, 5 в), в противном случае 81=0 (рисунок 5 а); при наличии зоны 3 82=1 (рисунок 5 в), в противном случае 82=0 (рисунки 5 а, 5 б).

Рассмотрим цистерну диаметром Э=2 м. Проведем расчеты по формулам (10) для слу-

чаев заполнения цистерны на 20% (а=2,085 рад, Ь=0,508 м), 50% (а=1,571 рад, Ь=1 м) и 80% (а=1,057 рад, Ь=1,492 м). Графики зависимости положения центра масс от скорости и радиуса закругления поворота приведены на рисунке 6.

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

/ (г )=РЬ-

при z < z1 81.11 — I - г2

при z1 < z < z2 —1— I г - cos в -—соэа этуЗ [ 2

при z < z2

Рисунок 4. Уровнн жидкости в автоцистерне 0=2 м, заполненной на 50%, при движении по правому повороту (г= 100 м) (1 - V =10 м/с, 2 - V =20 м/с, 3 - V =30 м/с)

(11)

Рисунок 5. Зоны заполнения цистерны при различных формах заполнения

Гаврилов А.А. и др. Динамическая модель относительного покоя жидкости в автоцистерне..

-0.4

V , м/с

г , м

V , м/с

г , м

Рисунок 6. Зависимость коордннат центра тяжести жидкости от скорости и радиуса закругления дороги

(1 - заполнение цистерны 20%, 2 - 50%, 3 - 80%)

Рисунок 7. Эпюра распределения сил инерции

где р - плотность жидкости, кг/м3;

Ь - длина цистерны, м.

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

Таким образом, полученная динамическая модель относительного покоя жидкости в частично заполненной автоцистерне при равномерном движении на криволинейном участке дороги позволяет получить уравнение свободной поверхности жидкости. Данное уравнение

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

16.06.2014

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

1. Кельрих, М.Б. Определение давления жидкого груза на днище котла железнодорожной цистерны при столкновении ее с препятствием / М.Б Кельрих, С.И. Криль // Прикладна гщромехатка. - 2011. Том 13. - №4. - С. 30-36.

2. Шимановский, А.О. Динамическая нагруженность конструкции цистерны при торможении / А.О. Шимановский, Г.М. Куземкина // Автомобильный транспорт. - 2006. - №19. - С. 32-34.

3. Власов, Ю. Л. Динамическая модель равномерного движения автоцистерны без перегородок с частичным заполнением / Ю. Л. Власов, Н. А. Морозов, А. А. Гаврилов // Вестник Оренбургского государственного университета. - 2013. - №9. -С. 227-231.

4. Власов, Ю.Л. Расчет инерциальных характеристик автоцистерны с частичным заполнением [Электронный ресурс] / Ю.Л. Власов, А.А. Гаврилов. - Прикладная программа (0,36 Мб). Зарегистрирована в УФАП ОГУ №833 от 20 июня 2013 г. - Оренбург: ОГУ, 2013.

5. Власов, Ю.Л. Влияние расположения груза в кузове грузового автомобиля на устойчивость автомобиля в поворотах / Ю. Л. Власов [и др.] // Вестник Оренбургского государственного университета. - 2011. - №4(123) - С. 152-155.

Сведения об авторах:

Гаврилов Александр Александрович, ассистент кафедры сопротивления материалов и теоретической механики Оренбургского государственного университета, e-mail: [email protected]

Власов Юрий Леонидович, доцент кафедры сопротивления материалов и теоретической механики Оренбургского государственного университета, кандидат технических наук, e-mail: [email protected]

Морозов Николай Анатольевич, доцент кафедры сопротивления материалов и теоретической механики Оренбургского государственного университета, кандидат технических наук,

e-mail: [email protected]

460018, г. Оренбург, пр-т Победы, 13, ауд. 20404, тел (3532) 372513