Динамика идеальной жидкости
Дипломная работа - Физика
Другие дипломы по предмету Физика
Содержание
Определение поля скоростей и вихревого поля2
Нахождение критических точек. Определение обтекаемого контура и линий тока4
Определение распределения давления на обтекаемый контур, направления и величины главного вектора сил давления6
Для заданной функции для плоского деформированного состояния проверить, что - бигармоническая; найти компоненты тензора напряжения и проверить условие нормировки10
Найти гидростатическую и девиаторную составляющие тензора напряжений, интенсивность напряжений13
Построить эпюры напряжений на границах заданной области14
Найти тензор деформации Коши-Грина, полагая материал линейно упругим17
Найти поле перемещений, считая, что точка О(0;0) неподвижна и поворот в плоскости вокруг оси, проходящей через О(0;0) отсутствует18
Определение поля скоростей и вихревого поля
Для определения поля скоростей необходимо в функции отделить действительную и мнимую части:
(1.1.1)
где функции и - функции действительных переменных и связаны между собой условиями Коши-Римана:
(1.1.2)
В данном случае удобнее записать комплексную координату в тригонометрическом виде , а затем отделить мнимую и действительную части:
Тогда и
В данном случае условия Коши-Римана примут вид:
и нетрудно проверить, что они действительно выполняются.
Градиент действительной части и будет полем скоростей, а функция называется функцией тока, которая позволяет определить линии тока и обтекаемый контур.
Найдем поле скоростей. Как было сказано выше, поле скоростей - это градиент от :
(1.1.3)
Оператор в полярных координатах имеет вид:
(1.1.4)
Тогда поле скоростей в полярных координатах будет равно:
(1.1.5а)
Или же в компонентах:
(1.1.5б)
Получив поле скоростей, можем теперь определить наличие вихря. По определению вихревое поле равно:
(1.1.6)
Подставив в (1.1.6) поле скоростей, учитывая (1.1.4), а также правила дифференцирования базисных векторов полярных координат ( ), получим следующее:
(1.1.7)
Учитывая (1.1.5б), получим, что , т.е. при движении не возникает завихрений.
Нахождение критических точек. Определение обтекаемого контура и линий тока
Критической называется точка, в которой скорость течения равно 0. Т.е. для нахождения таких точек необходимо решить систему уравнений:
(1.2.1)
Подставим в (1.2.1) найденные значения из (1.5б) и найдем критические точки:
Поскольку и не могут равняться одновременно нулю (следует из основного тригонометрического тождества), то очевидно, что , а может быть любым действительным числом. В декартовой системе координат такому решению соответствует лишь одна точка: начало координат.
Линии тока и обтекаемый контур, как было сказано выше, будем искать с помощью функции тока .
Линии тока характеризуются тем, что в каждой их точке направление касательной к ним, будут совпадать с направлением вектора скорости. Линии тока имеют вид:
(1.2.2)
Задавая различные константы, будем получать различные линии тока.
Обтекаемый контур будем находить также по формуле (1.2.2), но константу будем определять из следующего условия: линия тока, соответствующая обтекаемому контуру, будет проходить через критические точки. В нашем случае:
И уравнение обтекаемого контура примет вид:
Последнее уравнение характеризует прямые, направленные под углами . Рассмотрим контур, состоящий из прямых, направленных под углами 0 и 15 градусов к горизонтальной оси.
Построим линии тока и обтекаемый контур в математическом пакете MathCAD:
Черные кривые соответствуют линиям тока, а красные - обтекаемому контуру.
Определение распределения давления на обтекаемый контур, направления и величины главного вектора сил давления
обтекаемый контур вихревое поле
Для нахождения закона распределения давления на обтекаемый контур
воспользуемся интегралом Бернулли. При отсутствии массовых сил изменение давления вдоль линии тока связано с изменением скорости следующим выражением:
,(1.3.1)
где - плотность несжимаемой среды.
Применим формулу (3.1) вдоль критической линии тока. Полагая в формуле скоростей (1.1.5б) , получим:
Подставим найденное выражение в (3.1):
(1.3.2)
Определим постоянную C через давление в критической точке r=0:
Распределение нагрузок тогда примет вид:
(1.3.2б)
Построим эпюру давлений в виде для следующих значений:
Поскольку давление зависит только от , на эпюре представлена лишь часть контура (совпадает с горизонтальной осью), так как для оставшейся части контура (прямой под углом 15 градусов) эпюра будет иметь точно такой же вид. Необходимо отметить, что функция давления в данном случае убывает до бесконечности, поэтому мы ограничили контур так, чтобы давление со стороны жидкости было положительно.
Найдем главный вектор сил давления, действующий на горизонтальную составляющую контура. Запишем его проекции на ?/p>