Теми рефератів
> Реферати > Курсові роботи > Звіти з практики > Курсові проекти > Питання та відповіді > Ессе > Доклади > Учбові матеріали > Контрольні роботи > Методички > Лекції > Твори > Підручники > Статті Контакти
Реферати, твори, дипломи, практика » Курсовые проекты » Розрахунок стаціонарного теплового поля у двовимірній пластині

Реферат Розрахунок стаціонарного теплового поля у двовимірній пластині





жна визначити за допомогою функцій форм і значень температури в вузлових точках

. (6)

Функціонал (4) можна представити у вигляді суми функціоналів, кожен з яких відображає внесок у функціонал (4) елемента з номером e

. (7)

Мінімум функціонала (4) знаходимо з умови

(8)

Функціонал можна представити у вигляді

(9)

Тут, глобальний вектор температур, - матриця градієнтів, яка для функцій форми (5) прийме вигляд,В  . Локальний вектор температур. Тут матриця геометричних зв'язків має розмірність. Елементи цієї матриці визначаються наступним чином:; всі інші елементи дорівнюють нулю.

Продиференціюємо функціонал (9):

З виразу (8) з урахуванням останнього співвідношення отримуємо, де матриця теплопровідності елемента; вектор навантаження елемента.

У силу особливостей проведеної тріангуляції можна виділити три групи кінцевих елементів. У першу входять трикутники, у яких сторона i - j належить одній з зовнішніх кордонів. У другу - ті, у яких та ж сторона належить одній з внутрішніх кордонів. І, нарешті, третю групу складають елементи, сторони яких лежать всередині розглянутій області.

Залежно від того, до якої групи належить кінцевий елемент з номером e , матриця і вектор визначатимуться кілька різним чином.

Позначимо

.

Поверхневі інтеграли можна порахувати за допомогою відносних координат. Відрізки, що з'єднують будь-яку фіксовану точку P трикутника e c його вершинами, розбивають цей елемент на три трикутні частини площею. Координати визначаються з співвідношень. p> Використовуючи відносні координати, можна отримати такі співвідношення:

В В В 

Якщо кінцевий елемент з номером e належить до першої групи, то. Якщо ко другий, то. Нарешті, якщо елемент належить до третьої групи, то. p> Вектор температур, задовольняє умові (8) мінімуму функціоналу (4), знаходимо рішенням системи лінійних алгебраїчних рівнянь

, (10)

де глобальна матриця теплопровідності K і глобальний вектор навантаження F визначаються за формулами

,. (11)

Для вирішення завдання (10) застосовувався наступний алгоритм:

В· ОбчисленняВ  розкладання матриці ().

В· Оцінка числа обумовленості. Якщо число обумовленості більше (визначається точністю обчислювальної машини), то видається попередження, так як малі відхилення в коефіцієнтах матриці можуть привести до великих відхилень у рішенні.

В·В  . . p> Реалізація описаного вище методу проводилася на мові програмування С + + і FORTRAN в середовищі інтегрованої середовищі розробки Microsoft Visual C + + 6.0. Кінцеві результати даної роботи наведені на рис.4 - 7. table>В 

Рис. 4

В 

Рис.5

В 

Рис.6


Назад | сторінка 2 з 3 | Наступна сторінка





Схожі реферати:

  • Реферат на тему: Анексія Криму, як можна вірішіті Конфлікт України с Россией чі можна его ві ...
  • Реферат на тему: Оцінка логістичної діяльності ТОВ &ВКФ Вектор&
  • Реферат на тему: Опісові композіційно-мовленнєві форми в творах Т. Прохаська &З цього можна ...
  • Реферат на тему: Вектор в просторі. Скалярний твір ненульових векторів
  • Реферат на тему: Основи застосування методу матриці Бостонської консультативної групи на при ...