Изучение безвихревого течения идеальной жидкости

Страницы работы

Содержание работы

2.1.1.Кубичный четырёхугольный элемент.

  До сих пор применение метода конечных элементов были связаны с использованием одномерных линейных элементов , двумерных треугольных элементов и трёхмерного тетраэдра . Теперь рассмотрим новую группу элементов: двумерный четырёхугольник.

  Четырехугольный элемент представляет собой мультиплекс-элемент . Границы такого элемента должны быть параллельны координатным линиям для сохранения непрерывности при переходе от одного элемента к другому. Прямоугольный элемент является специальным случаем четырёхугольника . Свойства прямоугольного элемента служат основой для применения криволинейной системы координат, необходимой при использовании четырёхугольного элемента .

  Рассмотрим четырехугольный элемент, содержащий 12 узлов. Такой элемент называются соответственно кубичным элементом, так как его интерполяционный полином является кубичной функцией вдоль линий x=const и h=const.

  Интерполяционный полином соответственно для кубичного элемента (рис.2.1) записывается в виде

(2.1).

                                     

Рис.2.1. Кубичный четырехугольный элемент.

Функции формы для этого элемента представляют собой полиномы , идентичные по форме (2.1). Функции формы для двумерных элементов равны нулю во всех узлах , за исключением узла , номер которого совпадает с номером соответствующей функции формы; кроме того , они принимают узловые значения вдоль всех границ элемента , которые не содержат указанного узла . Например , функция формы N1 для элемента (рис.2.1 ) обращается в ноль во всех узлах , за исключением первого узла .

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

 (2.2)

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

f1=(1+h),

                             f2=(1-x),                      (2.3)

f3=(1-h),

f4=(1+x).

Каждая из них обращается в нуль на одной из границ элемента.

  Введем еще множество функций Fi , i=1, 2, 3, 4:

              fk , если узел b не принадлежит стороне k ,

 Fk=              k=1,2,3,4,                                                             (2.4)

              1 ,если узел b принадлежит стороне k .

Функция формы для кубичного элемента дается формулой

     (2.5)

Степень многочлена в (2.5) определяется числом имеющихся узловых условий. Его коэффициенты определяются приравниванием Nb единице в узле b и нулю во всех других узловых точках, которые не входят в произведение .

2.1.2.Определение функций формы .

  Рассмотрим определение функции формы N2 для элемента, показанного на рис.2.1.

  Узел 2 принадлежит первой стороне элемента, поэтому F1=1. Остальные три функции следующие: F2=f2=(1-x), F3=f3=(1-h) и F4=f4=(1+x).Произведение в (2.5) равно

 .

Общее выражение для N2 имеет вид

N2=(1-x2)(1-h)(a1+a2 x+a3h+a4x2+a5h2).

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

N2=1 при x= -1/3 , h= -1

 и

N2=0 при x=1/3 , h= -1.

Коэффициенты a4 и a5 должны быть вычеркнуты, поскольку члены вида x4 и x2h2 не входят в формулу (2.1) . Сохранение члена a3h приводит к системе уравнений с нулевым определителем, поэтому этот член тоже должен быть зачёркнут. Таким образом, для N2 имеем

N2=(1-x2)(1-h)(a1-a2x).

Используя условия в узлах, получаем систему

(1-1/9)2(a1-a2/3)=1,

(1-1/9)2(a1+a2/3)=0,

откуда находим a1=9/32 и a2= - 27/32. Окончательное выражение для функции формы имеет вид

N2=9/32(1-x2)(1-h)(1-3x).

  Приведём     расчёт значений функций формы для кубичного элемента (рис.2.1):

         

Простым сложением можно убедиться, что эти функции формы удовлетворяют критерию   для элемента.

2.1.3.Вычисление производных функций формы.

 Матрица Якоби определяется соотношением

  (2.6)

которое можно обратить, чтобы получить частные производные, по x и по у.

  Матрица Якоби является функцией x и h даже для простейших четырехугольных элементов. Эта зависимость легко обнаруживается при рассмотрении преобразования координат:

x=R1C1+R2C2+R3C3+R4C4 ,(2.6)

где C1 ,C2 ,C3 ,C4¾x-координаты четырех вершин. Rb используется здесь для обозначения функций, определяющих форму элемента. Дифференцируя x по x, получаем

где

(2.7)

После подстановки и перемножения имеем

Отсюда видно, что коэффициенты матрицы [J], являются функциями x и h.

   Определим частные производные  и  .

Похожие материалы

Информация о работе