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

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

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

 

Рис.2.9.

Пример.

Вывод уравнения колебаний натянутой упругой струны на основе

закона сохранения энергии.

Колебания принимаются малыми, что позволяет пренебречь

увеличением

длины струны и приводит к постоянству натяжения в любом

сечении и в любой момент времени. Пренебрегаем трением, т. е. диссипация (рассеяние) энергии принимается равной нулю (так как приток энергии извне принимается равным нулю). Предполагаем, что начальный момент времени соответствует состоянию равновесия. Например, в начальный момент времени струна натянута, находится в горизонтальном положении и концы ее закреплены. Для того чтобы начался процесс колебания струны, необходимо ее вывести из положения равновесия, положим, начальным импульсом. Обозначим отклонение струны от ее положения равновесия через x,t .

На рис. 8 цифрами 1—1 и 2—2 обозначены сечения, в которых закреплены концы струны. Выделим на расстоянии х от начала оси х элементарный отрезок струны x между сечениями, обозначенными цифрами

3—4.

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

точки 3 будет x,t а через интервал времени t оно составит

x,t t .

Средняя скорость этой точки равна

 

x,t t x,t

(2.136)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Мгновенная скорость в точке x определяется при t 0 соотношением

 

lim

x,t t x,t

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t 0

 

 

 

 

 

 

t

 

 

 

 

 

 

t

 

 

 

 

Учитывая, что масса отрезка струны l x ,

где l

линейная плотность,

получим изменение кинетической энергии элемента k имеет вид

 

 

 

 

1

 

 

 

 

2

 

 

 

 

1

 

2

k

 

 

l x

 

 

или в дифференциалах dk

 

 

ldx

. Суммируя по

 

 

 

 

 

 

 

 

 

2

 

 

 

t

 

 

 

 

 

2

 

t

длине получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 l

 

 

2

 

 

 

 

 

 

 

 

 

 

k

 

 

0

l

 

 

dx

(2.137)

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

t

 

 

 

 

 

 

 

 

 

 

 

Получим аналогичным образом соотношение для потенциальной энергии

d P dl dx

 

 

(2.138 )

 

 

 

 

Применяя известную из курса математики формулу длины дуги кривой ,

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

dl

 

 

 

 

 

dx .

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

91

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

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

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

2

d P

 

 

 

1 dx ,

1

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

2

P

 

 

 

1 dx .

1

 

 

0

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

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

 

2

 

1

2

1

 

1

 

 

 

 

 

x

 

2

x

и выражение для потенциальной энергии для всей струны

 

l

1

 

 

2

 

 

 

 

 

 

 

 

 

P

 

 

 

 

dx .

 

 

 

 

2

 

 

 

0

 

 

x

 

 

Поскольку сумма потенциальной и кинетической энергии составляет постоянную величину равную полной энергии получим k .

 

1 l

 

2

 

 

2

 

 

 

 

 

 

 

l

 

 

 

P

 

 

 

dx .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 0

 

t

 

 

 

x

 

 

 

 

 

Продифференцируем по t

 

 

 

 

1 l

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

l

 

t

 

P

 

x

 

dx

0

 

 

 

 

 

 

t

2

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

или

l

 

 

 

 

 

2

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

t

 

 

 

t

0

 

 

 

t

 

 

 

2

P

 

dx 0

 

x

 

 

 

 

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

l

 

 

2

 

2

 

 

 

l

 

2

P

 

 

dx 0

(2.139)

 

 

 

0

 

 

t t

 

 

x x t

 

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

 

 

 

 

 

 

 

 

 

0,

 

 

 

0

x

 

x 0

 

x

 

x l

 

 

 

 

 

92

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

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

получим

l

 

l

2

P

2

(2.140)

 

t t

2

x x

2

dx 0

0

 

 

 

 

 

 

 

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

l

2

P

2

0 (2.141)

 

t t2

 

x x2

 

Отсюда, поскольку 0 , можно разделить левую и правую части

x

уравнения на x . Тогда получим

l 2 P 2 . (2.142)

t2 x2

Это равенство представляет собой дифференциальное уравнение колебаний струны. Отметим, что здесь одно уравнение и одна неизвестная функция – поперечное перемещение струны. Натяжение струны было принято постоянным, поэтому оно исключено как неизвестное.

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

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

Дифференциальное уравнение для этого случая можно получить из уравнения движения.

В данном продольном стержне смещение происходит вдоль оси. Обозначим через x,t смещение вдоль оси ох (см. рис.2.10).

Скорость вдоль оси Ox

93

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

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

lim

x,t t x,t

 

 

u x,t

t

t

t0

 

 

Тогда ускорение вдоль оси Ox имеет вид

lim

u x,t t u x,t

 

u

 

2u

t

t

t2

t 0

 

 

Рис. 2.10

Упругое напряжение в сечении 1—1 в момент времени t будет F (x, t), а в сечении 2— 2 — будет

F x x,t .

Согласно второму закону Ньютона можно записать

m

2

F x x,t F x,t

1

t2

 

или

 

 

m

2

xF (2.142)

1

t2

,

где

— плотность; F — площадь поперечного сечения стержня.

Разделив уравнение (1.78) на х и перейдя к пределу при x 0 , получим

F

2

lim

F x x,t F x,t

(2.143)

t2

x

 

 

x0

 

F

2

 

F

(2.144)

 

 

t2

 

x

 

 

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

сохранения энергии.

В качестве примера возьмем вывод уравнения колебания груза на упругой нити на основе закона сохранения энергии.

Рассмотрим колебание груза, подвешенного к упругой нити. Примем, что сила упругости пропорциональна удлинению х (рис. 2.10).

На рис. 2.10 показаны: 1 — положение равновесия, 2 — отклонение от положения равновесия, которое изменяется в зависимости от времени. Если

94

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

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

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

0.

 

 

Кинетическая энергия в положении 2 равна

m x ' 2

, где т —масса;

 

 

2

 

 

 

 

 

x '

dx

-

скорость груза.

 

 

dt

 

 

 

 

 

 

 

Потенциальная энергия — это энергия упругости нити. Силы упругости приняты пропорциональными х, т. е. кх совершает элементарную работу kxdx, соответствующую элементарной потенциальной энергии. В положении 2 потенциальная энергия определяется как

 

x

 

 

kx

2

 

 

 

 

kxdx

 

 

.

 

 

2

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда

 

 

 

 

 

dx 2

 

 

 

 

 

 

 

m

 

 

 

 

kx

2

 

 

 

 

 

 

dt

 

 

 

const , (2.145)

2

 

 

2

 

где Ф — полная энергия в положении 2. Дифференцируя (2.145) по t, получаем

m

d 2 x

 

dx

kx

dx

0 .

dt2 dt

dt

 

 

 

При x ' 0

 

тх" + кх = 0.

(2.146)

Уравнение (2.146) и есть выражение второго закона Ньютона.

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

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

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

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

моделирования

широко

применяется

в

различных

отраслях

 

 

 

 

 

95

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

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