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

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

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

- инерции к силам вязкости:

 

Fинерц

 

 

d 2 V

2

 

V d

 

 

Re ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fвя зк

 

V d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т.е. получаем критерий Рейнольдса.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- пластической прочности структуры жидкости к силам вязкости:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fпл

 

 

 

0 d 2

 

0 d

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fвя зк

 

V d

V

Sen

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

получаем критерий Сен-Венана- Ильюшина;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

 

 

 

 

P l 2

 

 

 

 

P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

давл

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Eu ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V 2l 2

V

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fинер

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

получаем критерий Эйлера;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- инерции к силам тяжести:

Fинерцб

 

V 2l 2

 

 

V 2

 

Fr ,

 

 

 

 

 

 

 

 

 

 

 

 

 

Fтяж

 

r 3 g

 

r g

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

получаем критерий Фруда;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- сопротивления движению тела на поверхности жидкости:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fс

 

 

Fс

Ne ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fинер

V 2l 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

получаем критерий Ньютона и т. д.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Методом анализа размерностей можно получить и многие известные

формулы. Рассмотрим применение на ряде задач.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задача 1. Определение периода колебаний математического маятника

 

 

 

 

 

 

 

 

 

 

Запишем уравнения для сил инерции F а m и тяжести

F mg .

 

 

 

 

 

 

 

 

 

 

 

Переходя к масштабам подобия, имеем F m g

 

 

m

lс

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

c

c

c t 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда gс

lс

, и, решая это уравнение относительно tс, получаем

tс

 

 

lс

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g

с

 

 

 

 

 

 

 

 

 

 

 

 

 

с

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Возвращаясь от масштабов подобия к параметрам, будем иметь

t с

 

 

 

 

l

,

 

 

 

 

 

 

 

g

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где с- некоторая постоянная.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение дифференциального уравнения колебаний математического

маятника

 

d 2

 

g

sin ,

где φ –

угол отклонения

 

 

маятника

от

положения

 

 

 

 

 

 

 

 

 

 

 

 

 

d t 2

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

равновесия (от вертикали), дает значение постоянной

 

 

с = 2π.

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда для периода колебаний данного тела получим хорошо известную из курса элементарной физики формулу

56

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

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

Т 2

 

l

 

(2.70)

 

g

 

 

 

 

- закон математического маятника.

Задача 2. Определение периода колебаний груза, подвешенного на пружине Запишем закон Гука для пружины: Fp = С l, где С- жесткость,

l – деформация, и закон Ньютона Fu = m a.

Перейдя к масштабам подобия, получим Fp.c= Сclc ; Fи.с.= mc ac.

Приравнивая масштабы сил, получаем: Сclc mc ac mc lс . tс

Отсюда

tс

 

mс

 

.

 

 

 

 

 

 

 

 

 

 

 

Сс

 

 

 

 

 

 

 

 

Совершая обратный переход к параметрам, получаем t с

m

,

С

 

 

 

 

 

 

 

 

 

где с - некоторая постоянная.

Решение дифференциального уравнения колебаний груза с массой m на пружине с жесткостью С:

d 2 х С х m g , d t 2

где х-

отклонение груза от положения равновесия, дает значение постоянной

с = 2π.

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

 

 

 

 

 

 

 

 

T 2

 

 

m

,

(2.71)

 

 

 

 

 

 

С

 

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

Задача 3. Моделирование центральной силы тяготения. Сила притяжения между массой центрального тела М и вращающегося

вокруг него на радиусе r тела с массой

 

m

равна

 

F f

M m

,

 

 

r 2

 

 

 

 

 

 

 

 

 

 

 

 

где f – постоянная всемирного тяготения.

 

 

 

 

 

 

 

 

 

 

 

 

Переходя к уравнению масштабов, получаем

F

f

 

M C mC

.

C

 

 

 

 

 

 

 

 

 

C

 

 

r 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

Так как f и М – const, то

f

 

M

 

1

и

F

mC

.

 

C

C

 

 

 

 

 

 

 

C

 

l 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

С другой стороны:

57

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

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

F ma;

 

FC mC aC ;

1

 

aC ;

1

 

lC

.

lC2

 

lC2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tC2

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

 

 

 

 

 

 

 

 

 

 

 

 

 

lC3

1

или

 

 

l13

 

t12

1.

 

 

 

 

 

 

 

 

 

 

l23

 

t22

 

 

 

 

 

 

tC2

 

 

 

 

 

 

Отсюда получаем известный третий закон Кеплера

 

2

 

 

3

 

 

 

 

 

 

 

 

 

 

 

t1

 

 

l1

 

 

 

 

 

 

 

 

(2.72)

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

t2

 

 

l2

 

 

 

 

 

 

 

 

 

 

 

Отношение квадратов времен обращения планет вокруг Солнца равно отношению кубов полуосей их эллиптических орбит.

Рассмотрим теперь третью теорему подобия. Эта теорема имеет несколько формулировок.

Первая формулировка

Необходимым и достаточным условием обеспечения подобия является пропорциональность сходственных параметров, входящих в условия однозначности, и равенство критериев подобия изучаемого явления.

Третья теорема подобия именуется также обратной теоремой подобия или теоремой Кирпичева – Гухмана.

Условия однозначности - это условия, определяющие индивидуальные особенности процесса или явления и выделяющие из общего класса конкретный процесс или явление.

К ним относятся следующие, не зависящие от механизма самого явления, факторы и условия:

1)геометрические свойства системы, в которой протекает процесс;

2)физические параметры среды и тел, образующих систему;

3)начальное состояние системы (начальные условия);

4)условия на границах системы (граничные условия);

5)взаимодействие объекта и внешней среды.

Нельзя математически сформулировать условия однозначности в общем виде. В каждом конкретном случае они могут быть различны в зависимости от рода решаемой задачи (и вида уравнения - если уравнение известно).

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

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

Положение 1

58

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

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

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

Положение 2 В созданной согласно положению 1 модели осуществление подобия

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

Положение 3 Моделирование по 1 и 2 положениям возможно в сколь угодно сложных,

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

В заключении этого параграфа приведем справочные сведения по теории подобия

Первая теорема подобия имеет следующий вид: явления, подобные в том или ином смысле, имеют определенные сочетания параметров, называемые критериями подобия, численно одинаковые для подобных явлений.

Размерности некоторых производных величин будут следующие:

1)Площадь S = l2 = > [L]2.

2)Объем Q = l3 => [L]3.

3)Момент силы M = F∙l => [M]∙ [L]2 ∙[T]-2.

4)Момент инерции I = m∙l2 => [M]∙[L]2.

5)Мощность N = A / t => [M]∙[L]2∙[T]-3.

6)Давление P = F / S =>[M]∙[L]-1[T]-2.

7)Плотность ρ = m / V => [M]∙[L]-3 .

8)

Динамическая вязкость [M ] [L] 1

[T ] 1.

9)

Кинематическая вязкость

 

L 2

T 1 .

 

 

 

 

 

Вторая теорема подобия (π – теорема). Эта теорема позволяет получить критерии подобия и в тех случаях, когда уравнение процесса можно представить лишь виде функциональной зависимости между параметрами системы и процесса:

y = f (x1, x2,…,xj,…,xn-1),или

F (y, xl, x2,....,xj,....,xn - l) = 0,

где f и F означают лишь символы зависимости.

59

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

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

π – теорема имеет несколько формулировок, одна из которых (основная)

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

Таблица 2.2

Размерности физических величин

 

Параметр

 

Показатель степени

 

 

 

 

 

 

 

 

 

 

[L]

 

[M]

 

[T]

 

 

 

 

 

 

 

m

 

α1= 0

 

β1 = 1

 

γ1 = 0

 

 

 

 

 

 

 

С

 

α2 = 0

 

β2 = 1

 

γ2 = – 2

x

 

α3 = 1

 

β3 = 0

 

γ3 = 0

t

 

a4 = 0

 

β4 = 0

 

γ4 = 1

F0

 

а5 = 1

 

β5 = 1

 

γ5 = – 2

ω

 

a6 = 0

 

β6 = 0

 

γ6 = – 1

k

 

а7 = 0

 

β7 = 1

 

γ7 = – 1

 

 

 

 

 

 

 

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

1

Критерий Архимеда

Ar

 

d 3 т в

g

.

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ne

 

 

 

 

Fc

 

 

2

Критерий Ньютона

 

 

 

 

 

 

 

.

 

 

V 2 d 2

 

 

3

Критерий Ньютона

Ne

F t 2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m l

 

 

4

Критерий Ньютона

Ne

 

M вр

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J 2

 

 

5

Критерий Коши Сa

 

V

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

 

 

 

 

 

 

 

 

 

6

Критерий Прандтля

Pr

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7

Критерий Фурье Fo t .

 

 

 

 

 

l 2

 

 

 

 

 

 

 

 

 

60

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

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