Материал: МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ В ЗАДАЧАХ НЕФТЕГАЗОВОЙ ОТРАСЛИ

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

СПБГУАП| Институт 4 группа 4736

Соответственно, запись u r в предыдущей формуле означает, что местные скорости в сечении трубы изменяются по радиусу, а u r

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

скоростей

равен

объему старой. Таким образом, реализуется условие:

u r V const и, соответственно, получаем:

Q VdS V dS VS u r dS

S

S

S

V QS

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

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

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

131

Контакты | https://new.guap.ru/i03/contacts

uy , uz

СПБГУАП| Институт 4 группа 4736

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

дифференциальных уравнений Навье-Стокса принимает вид:

 

1 p

X

 

 

 

x

 

 

 

 

 

1 p

Y

 

 

 

 

 

y

 

 

 

 

 

1 p

Z

 

 

 

z

 

ux uxx

0

0

Последние два уравнения полностью совпадают с уравнениями гидростатики, а это означает, что в пределах живого сечения движущейся жидкости давление распределено по гидростатическому закону. В частности, если объемной силой является сила тяжести, то для произвольной точки, находящейся в живом сечении, имеем Y=0, Z= - g cos α, где α – угол между направлением вертикали и нормалью к вектору скорости. Следовательно, в пределах живого сечения (для Х= const) гидростатический напор постоянен:

z

p

const , либо gz

p

const

g

 

 

 

 

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

132

Контакты | https://new.guap.ru/i03/contacts

СПБГУАП| Институт 4 группа 4736

4. ВЫВОД ОСНОВНЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ КОЛЕБАНИЙ МЕХАНИЧЕСКИХ СИСТЕМ И ТЕПЛОВЫХ ПРОЦЕССОВ

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

Рассмотрим стержень (Рис.4.1) длиной l, площадью поперечного сечения F(x), модулем упругости E, плотностью ρ(x). Пусть на стержень вдоль оси Ox действует продольная сила с интенсивностью q(x,t). Ось Ox направлена вдоль оси стержня. Функция F(x), ρ(x) и q(x) непрерывны, или кусочно-непрерывны

[3,5,21].

 

Из стержня

выделим малый

 

участок x сечениями I-I и II-II, на

 

который

действуют

следующие

 

силы:

 

 

 

 

 

а)

продольная

внешняя

 

нагрузка,

равная

q(x,t) x

 

направленная вдоль оси Ox;

 

б)

равнодействующая

 

натяжения со

стороны

левого

 

участка

 

I-I,

 

равная

Рис 4.1

S(x;t)=F(x)Δ(x;t)

и

направленная

противоположно оси Ox . Здесь Δ(x;t) - напряжение в сечении F(x), согласно закону Гука равное

,

где u – упругое перемещение;

в) равнодействующая натяжения со стороны правого участка II-II, равная S(x+Δx,t)=F(x)Δ(x+ x,t) и направлено по оси Ox . Здесь (x+ x,t) напряжение в сечении II-II;

г) сила сопротивление среды. Если её интенсивность по длине стержня обозначить через Ф (x,t,υ), то на участок x придётся сила Ф (x,t,υ) x.

Так как сила сопротивления направлена противоположно скорости, то её знак определяется характером её зависимости. Например, для кулонова трения

сила сопротивления запишется в виде

 

 

 

 

 

 

 

 

 

 

Ф sign

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

,

 

 

 

 

 

 

 

 

 

Для случая, когда сила сопротивления пропорциональна первой степени

скорости, знак минус перед Ф сохраняется и опускается sign

 

.

 

Сила инерции малого элемента стержня определяется выражением

133

Контакты | https://new.guap.ru/i03/contacts

СПБГУАП| Институт 4 группа 4736

Тогда по принципу Даламбера

–

(4.1)

Разлагая S(x+ x,t) в ряд Тейлора поставляя в (6.1.1), после деления на x и перехода к приделу при x 0 получим

(4.2)

Аналогичный результат может быть получен, если в результате разложения S(x+ x;t) пренебречь членами, содержащими x степени выше второй. В дальнейшем, учитывая это обстоятельство при выводе других дифференциальных уравнений, в разложениях будем оставлять два члена разложения, подразумевая, что имеет место такой предельный переход.

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

Например, если принять ρ(x)=const, F(x)=const, а также принять пропорциональность силы сопротивления скорости, то получится дифференциальное уравнение вынужденных продольных колебаний однородного стержня

 

 

 

 

 

 

 

 

–

 

g ( x , t ) (4.3)

 

 

 

 

 

 

где

 

; g ( x , t ) =

 

 

; h – коэффициент сопротивления.

 

 

Приняв g ( x ; t ) = 0 , получим дифференциальное уравнение свободных колебаний однородного стержня

–

(4.4)

При h=0 из (4.4) получится уравнение свободных колебаний стержня в идеальной среде.

= a2

 

(4.5)

 

Для решения дифференциальных уравнений в частных производных необходимо задание начальных и граничных условий.

В качестве начальных условий обычно задается значения функции и ее первой производной по аргументу t в момент времени . При задании только начальных условий совместно с дифференциальным уравнением имеет место задача Коши. Начальные условия для уравнения типа (1.5) могут быть записаны в виде

134

Контакты | https://new.guap.ru/i03/contacts

СПБГУАП| Институт 4 группа 4736

 

 

 

(4.6)

где

заданные функции.

В качестве граничных условий задаются значения функции и ее производных в некоторых фиксированных точках. Так, например если известен закон движения обоих концов стержня, то эти условия задаются в виде

( 4.7)

Если на одном из концов известна внешняя действующая сила, то согласно закону Гука

или

 

 

 

 

 

, где

 

 

.

 

 

 

 

Следовательно

 

 

Тогда граничное условие имеет вид

 

 

 

(4.8)

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

(4.9)

где с – жесткость пружины.

Если точки закрепления являются подвижными и их перемещения задаются функциями типа S(t), то граничные условия (4.9) примут вид

(4.10)

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

(4.11)

Для неограниченного стержня вместо полной задачи можно рассматривать задачу только с начальными условиями – задачу Коши:

135

Контакты | https://new.guap.ru/i03/contacts

Источник: https://studfile.net/preview/16405273/