Учебно-методическое пособие Рекомендовано методической комиссией механико-математического факультета для студентов ннгу, обучающихся по направлению подготовки 010500 «Прикладная математика и информатика» icon

Учебно-методическое пособие Рекомендовано методической комиссией механико-математического факультета для студентов ннгу, обучающихся по направлению подготовки 010500 «Прикладная математика и информатика»


1 чел. помогло.
Смотрите также:
Учебное пособие Рекомендовано методической комиссией механико-математического факультета для...
Учебно-методическое пособие Рекомендовано методической комиссией факультета управления и...
Методические указания по выполнению курсовых работ по курсу...
Учебно-методическое пособие Рекомендовано методической комиссией финансового факультета для...
Учебно-методическое пособие Рекомендовано методической комиссией филологического факультета для...
Практикум Рекомендовано методической комиссией факультета международных отношений для студентов...
Учебно-методическое пособие для студентов экономического факультета специальности 060800...
Учебное пособие Рекомендовано методической комиссией факультета социальных наук для студентов...
Пособие для преподавателей русского языка, ведущих занятия с иностранными студентами...
Учебно-методическое пособие по английскому языку...
Практикум по ценообразованию...
Учебное пособие для студентов ммф томск 2007...



Загрузка...
страницы: 1   ...   5   6   7   8   9   10   11   12   13
вернуться в начало
скачать
^

4. Построение базисных координатных функций в МКЭ

4.1. Основные требования к координатным функциям в МКЭ



Первым и вероятно наиболее ответственным шагом построения конкретных методик решения задач на основе МВН в форме МКЭ является выбор базисных функций используемых конечных элементов. Очевидно, что для получения качественного решения такие функции должны удовлетворять определенным требованиям, обеспечивающим сходимость КЭ решений [2-4]. Часть таких требований была сформулирована в предыдущей главе для базисных функций в МВН, а именно:

- функции должны быть линейно независимы;

-функции должны удовлетворять условиям полноты.

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

Кроме этого к координатным функциям в МКЭ предъявляется еще ряд специфических требований, среди которых важнейшими являются следующие:

- функции должны быть согласованы с формой КЭ, числом и расположением степеней свободы в элементе. Число степеней свободы элемента должно соответствовать числу неопределенных параметров в полиномиальном представлении функции;

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

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

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

^

4.2. Построение базисных функций конечных элементов в обобщенных координатах



Наиболее естественным образом базисные функции элементов могут быть получены на основе решения системы алгебраических уравнений для значений выбранного степенного полинома в узлах элемента. Этот подход является достаточно общим, применим для любых типов элементов и законов распределения функций и носит название «базисные функции в обобщенных координатах» [5]. Основные этапы применения этого подхода удобнее продемонстрировать на примере построения функций формы четырехузлового прямоугольного элемента (см. рис.4.1).

Закон изменения функции u(x,y) в таком элементе может быть записан в виде степенного полинома


,

или

, (4.1)

{aT}=(a0 a1 a2 a3)


Далее составляется система алгебраических уравнений для значений функции в узлах элемента

, (4.2)

,





Рис. 4.1. Схема четырехузлового прямоугольного элемента





Неопределенные параметры выражаются через узловые функции элемента





В результате функции формы элемента могут быть получены в виде


(4.3)

.


Одним из основных недостатков рассмотренного подхода является необходимость обращения матрицы для каждого конечного элемента исследуемой области. В случае использования полиномов выше первого порядка обращение этой матрицы в явном виде представляется весьма сложным и приходится применять численные процедуры. Кроме этого при использовании неполных полиномов выполнение условий межэлементной совместности функций, построенных на основе рассмотренного подхода, не является очевидным. Для выяснения этих условий следует записать закон распределения функций для каждой из сторон (граней) на основе принятого распределения в элементе и установить соответствие между числом неопределенных параметров в получающемся распределении с числом узловых значений функций на рассматриваемой стороне. Например, в представленном на рис.4.2 четырехугольном конечном элементе с билинейным законом распределения функций (4.1) непрерывность функции u(x,y) вдоль сторон 1-2 и 2-3 будет выполняться по разному. Если в элементе ввести локальные координаты и направленные вдоль первой и второй сторон соответственно, то связь этих координат с общими декартовыми координатами x,y можно представить в виде




.





Рис. 4.2. Схема четырехузлового непрямоугольного элемента


Тогда изменения функции и вдоль сторон S1 и S2 могут быть определены подстановкой этих выражений в соотношения (4.1)










.


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

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

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

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





оставить комментарий
страница11/13
Дата02.10.2011
Размер0,69 Mb.
ТипУчебно-методическое пособие, Образовательные материалы
Добавить документ в свой блог или на сайт

страницы: 1   ...   5   6   7   8   9   10   11   12   13
плохо
  2
отлично
  2
Ваша оценка:
Разместите кнопку на своём сайте или блоге:
rudocs.exdat.com

Загрузка...
База данных защищена авторским правом ©exdat 2000-2017
При копировании материала укажите ссылку
обратиться к администрации
Анализ
Справочники
Сценарии
Рефераты
Курсовые работы
Авторефераты
Программы
Методички
Документы
Понятия

опубликовать
Загрузка...
Документы

наверх