Московский Государственный Технический Университет имени Н.Э. Баумана.

Курсовая работа

По предмету

“Дифференциальные уравнения.”

Тема: Математическая модель всплытия подводной лодки

Выполнила:

студентка группы

ФН 2-31, Иванова А.

Научный руководитель:

профессор В.И. Ванько.

Москва 2001 г.

Введение.

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

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

Рассматривая, таким образом, действия сил на объект, используя основные законы механики и соотношения между силами мы можем составить дифференциальное уравнение или систему дифференциальных уравнений, решая которую, можно получить её частное или общее решение (в зависимости от вида системы).

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

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

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

Рассмотрим подводную лодку как материальную точку, которая движется по горизонтали на некоторой глубине, с некоторой постоянной скоростью. Лодка удифферентована, то есть силы, которые действуют на лодку по вертикали, как показано на рис.1, (сила тяжести и выталкивающая сила Архимеда) равны по модулю.
По горизонтали, на лодку действует сила сопротивления, модуль которой примем в виде:

[pic]
Где степень [pic] и коффициент пропорциональности [pic] это некоторые числа, характерные для данной среды, и зависящие от факторов среды, таких как: плотность

Рис. 1 воды, её температура, и величины скорости.
Сила Архимеда, действущая на лодку, зависит от размеров лодки, а именно от её объема, и плотности воды.

[pic]

В этой формуле [pic] – это плотность жидкости, [pic]–объем тела, погруженного в жидкость, [pic] = 9.81 м / c2 – ускорение свободного падения.

Пусть в некоторый момент времени выключены двигатели и сбрасывается балласт. Двигаясь по инерции, а также под действием силы Архимеда, она начнет всплывать по некоторой траектории (рис.2).
Проведем радиус вектор [pic]из начала координат:

[pic]
Вектор скорости также можно разложить на составляющие по осям x и y:

Рис. 2
[pic]
Тогда силу сопротивления мы можем записать так:

, так как вектор скорости всегда направлен по касательной к траектории движения, а сила сопротивления имеет противоположное направление.

По второму закону Ньютона:

[pic], где вектор [pic] - это вектор силы тяжести, действующей на лодку. [pic]- некоторая функция зависящая от времени.

Запишем это векторное уравнение в проекциях на оси.
В проекции на ось [pic]: [pic]
В проекции на ось [pic]: [pic]
В результате получим систему дифференциальных уравнений:

[pic] , где масса [pic]- функция зависящая от времени. Решая эту систему для произвольного значения [pic], и заданных начальных условий, мы получим уравнение траектории движения подводной лодки.
Пусть масса лодки изменяется по линейному закону [pic], где [pic] - масса корпуса, [pic] - это скорость вытеснения воды из цистерн, которую будем считать постоянной, а [pic] - некоторый момент времени, в который вся вода из цистерн вытеснена. Как показано на рис.3, в некоторый момент времени произведение [pic] будет равняться 0, и мы

Рис. 3 получим [pic], то – есть, вся вода из цистерн будет вытеснена.

Решим эту систему для частного случая.
Пусть [pic] = 1. В начальный момент времени лодка находится в начале координат, а вектор её скорости направлен по горизонтали и равен [pic].

Тогда начальные условия будут такими:

[pic]

[pic] .
В рассматриваемом частном случае, система уравнений принимает следующий вид:

[pic] .

Первое уравнение этой системы зависит только от [pic], второе только от
[pic], поэтому их можно разделить. Решим сначала первое уравнение системы.

[pic]

Так как в это уравнение не входит [pic], можно сделать замену [pic]. Решая таким образом полученное уравнение первого порядка с разделяющимися переменными, получим:

[pic]

[pic]

[pic].

[pic].

Решим второе уравнение системы.

[pic]

Делая аналогичную замену, получим линейное неоднородное уравнение, решая которое, получим:

[pic]

[pic]

[pic]

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

[pic]
Траектория движения подводной лодки для заданных начальных условий и
[pic]=1 изображена на рис. 4.
Решим исходную систему для произвольного значения параметра [pic].
На [pic] накладывается ограничение: [pic], так как только при выполнении этого условия, сила сопротивления оказывается прямо

Рис 4. пропорциональна скорости.

Систему

[pic] приведем к нормальной форме Коши, вводя новые переменные.

[pic].

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

Начальные условия для которой имеют вид:

[pic] .
Решения этой системы для нескольких значений параметра [pic] представлены на рис. 5.

Рис. 5 а.

Так как при близких значениях [pic] траектория почти не изменяется и графики сливаются, для большей наглядности изобразим их в более крупном виде.

Рис.5 б.

На рис.5 а,б изображены решения исходной системы для [pic] [pic] [pic]
[pic] [pic][pic]
Найдем значение [pic] для которого время всплытия будет наименьшим и уравнение движения при этом значении параметра. Очевидно, что если [pic] то
[pic], и система принимает следующий вид:
[pic], где [pic] - функция, зависящая от времени.

График решения этой системы представлен на рис.6.
Функция возрастет быстрее, чем в случаях с другим значением [pic]. А это значит, что, при данном значении параметра, она всплывет с определенной глубины за минимальное время.

Рис. 6 При отрицательном значении праметра траектория будет практически совпадать с траекторией
[pic], но, в этом случае, задача теряет физический смысл.

Заключение.

Мы рассмотрели только частные случаи решения задачи. Исходную систему, невозможно решить в общем виде, без использования ЭВМ, или численных методов решения задачи.
Но, уже по частным случаям решений, можно увидеть некоторую закономерность, на основании которых, уже можно делать какие-то выводы.
Сам процесс всплытия подводной лодки – достаточно сложный физический процесс. На всплытие лодки влияют не только несколько сил действующие на неё. Большое значение имеют гидродинамические параметры, которые в построении данной модели не учитывались. Для численных решений системы и построения графиков были взяты реальные размеры и начальная скорость подводной лодки, что позволило как можно больше приблизить рассмотренный процесс к реальному.

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

1. Агафонов С.А., Герман А.Д., Муратова Т.В. Дифференциальные уравнения

М.: Изд-во МГТУ имени Н.Э. Баумана, 2000. - 347 с.

2. Степанов В.В. Курс дифференциальных уравнений

М.: Изд-во технико-теоретической литературы, 1950. - 467 с.

3. Осипенко Л., Жильцов Л., Мормуль Н. Атомная подводная эпопея

М.: “Боргес”, 1994. - 350 c.
-----------------------
[pic]

[pic]

[pic]

[pic]

[pic]

[pic]

[pic]

[pic]

[pic]