ГДЗ №2 Вариант 3 (1051185)

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

Текст из файла

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ

ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ

ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ

АЭРОКОСМИЧЕСКИЙ ФАКУЛЬТЕТ

КАФЕДРА

«Вычислительная математика и математическая физика»

ДОМАШНЕЕ ЗАДАНИЕ №1

по дисциплине

«Методы конечных элементов »

Вариант № 3

Студент: ________________ Волков М.Н.

Группа АК3-71

Руководитель: ___________________ Шпакова Ю. В.

МОСКВА 2015

Содержание

1. Постановка задачи стр. 2

2. Математическая модель стр. 3

3. Решение задачи теплопроводности стр. 5

4. Численная реализация задачи теплопроводности методом конечных элементов стр.7

5. Алгоритм решения задачи с использованием ЭВМ стр. 11

6. Метод Холецкого для решения СЛАУ. стр. 15

7. Реализация метода на ЭВМ. стр. 16

8. Результаты решения задачи теплопроводности стр. 18

9. Выводы стр.20

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

1. Постановка задачи.

Дана исследуемая область (Рис. 1) с граничными условиям (Tсреды или Tср равносильно теплообмену со средой). Геометрические параметры области A, B, L, a, b, l [см] задаются самостоятельно. Воздействие теплового потока принять равным , коэффициент теплоотдачи от стенки к среде ; T – заданная температура стенки, 150 ; - температура окружающей среды, 22 . Кроме того, задан коэффициент теплопроводности . Требуется найти распределение температуры в исследуемой области: свести задачу к решению СЛАУ и решить СЛАУ, используя метод Холецкого для разреженных матриц.

Рис.1 – Исследуемая область

1. Математическая модель

Рассмотрим двумерную краевую задачу:

, (1)

С граничными условиями вида:

-подвод тепла к телу

-теплообмен со средой (граничные условия 3 рода)

где область, изображенная на Рис. 2.

Рис.2

Для решения задачи составим невязку для уравнения (1):

, (2)

Проверим, что выполняется условие:

, (3)

Подставим (2)в (3):

, (4)

Запишем формулу Гаусса-Остроградского:

(5)

Проинтегрируем по частям:

(6)

Подставим (6) в (5):

(7)

Подставим (7) в (4), с учётом :

(8)

Введем оператор (9)

Кроме того, учтём ,что на некоторой части границы Г заданы краевые условия:

(10)

Перепишем (8) с учётом (9)и (10):

(11)

Далее разобьем область Ω на треугольники и перенесем задачу на область одного такого элемента:

(12)

1. Решение задачи теплопроводности

Рассмотрим частный случай двумерной задачи (1) – задачу теплопроводности:

, (13)

Краевые условия принимают вид:

,

Где T-функция температуры,

h-коэффициент теплообмена,

Тср-температура среды,

-коэффициенты теплопроводности.

Запишем (8) для задачи теплопроводности, учитывая, что и что внутренние источники тепла в нашем случае отсутствуют.

Границу рассматриваемой области, к которой подводится тепловой поток, обозначим , а границу, на которой происходит теплообмен со средой- :

(14)

Сведем интегральную задачу к СЛАУ вида

, (15)

где

(16)

-заданные условия

- локальная матрица характеристик элемента,

-матрица формы,

.

Так как матрицы B и D не зависят от координат элемента:

, (17)

где - площадь одного конечного элемента

1. Численная реализация задачи теплопроводности методом конечных элементов

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

Рассмотрим отдельный конечный элемент (Рис.3).

Р
ис.3

Аппроксимацию основной переменной на конечном элементе зададим в виде полного полинома первой степени

, (18)

где коэффициенты полинома.

Представим линейную интерполяцию (6) так, чтобы в качестве неизвестных выступали не , а значения функции в точках 1, 2, 3, т.е. , , . Для этого подставим координаты узлов в (6), в результате чего получим систему

Используя метод Крамера, получим:

, (19)

где

-удвоенное значение площади конечного элемента.

Перегруппируем в (19) слагаемые относительно .

(20)

Основное свойство матрицы формы N:

Подставим матрицу в (16):

(21)

Тогда условие подведения тепла на границе примут вид:

(22)

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

(23)

Условие конвективного теплообмена :

(24)

(25)

(26)

= (27)

= (28)

= (29)

Условие, если на границах задан тепловой поток : = (30-32)

Схема составления глобальной матрицы жёсткости.

Локальная матрица жесткости i-го элемента имеет вид:

, (33)

Тогда

, (34)

где -номера узлов i-го элемента,m-число элементов.

Схема составление глобальной правой части.

Локальная правая часть i-го элемента имеет вид:

(35)

Тогда

, (36)

где -номера узлов i-го элемента,m-число элементов.

1. Алгоритм решения задачи с использованием ЭВМ

1. Разбиение области на конечные элементы (Рис. 4):

Рис. 4 – Дискретизация области

1. Нумеруем узлы элементов.

1. Составление локальной матрицы жесткости и локальной правой части.

Приведем пример для элемента №5 и №40 (Реализация на С++)

Характеристики

Тип файла документ

Документы такого типа открываются такими программами, как Microsoft Office Word на компьютерах Windows, Apple Pages на компьютерах Mac, Open Office - бесплатная альтернатива на различных платформах, в том числе Linux. Наиболее простым и современным решением будут Google документы, так как открываются онлайн без скачивания прямо в браузере на любой платформе. Существуют российские качественные аналоги, например от Яндекса.

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

Файлы такого типа обычно разбиты на страницы, а текст может быть форматированным (жирный, курсив, выбор шрифта, таблицы и т.п.), а также в него можно добавлять изображения. Формат идеально подходит для рефератов, докладов и РПЗ курсовых проектов, которые необходимо распечатать. Кстати перед печатью также сохраняйте файл в PDF, так как принтер может начудить со шрифтами.

Список файлов домашнего задания