3874

Чисельний метод розв’язання задач стаціонарної теплопровідності

Практическая работа

Энергетика

Чисельний метод розв’язання задач стаціонарної теплопровідності Умова задачі Використовуючи чисельний метод, визначити розподіл температури в поперечному перерізі довгого стального стержня при кроці розрахункової сітки ...

Украинкский

2012-11-09

202.5 KB

8 чел.

Чисельний метод розв’язання задач стаціонарної теплопровідності

Умова задачі

Використовуючи чисельний метод, визначити розподіл температури в поперечному перерізі довгого стального стержня при кроці розрахункової сітки δ=10мм. На основі вирахуваних локальних значень температур накреслити сімейство інтерпольованих ізотерм , для поперечного перерізу стержня провести лінії теплового потоку.

Поперечний переріз стержня має прямокутну форму з розмірами l1=30мм, l2=40мм. Коефіцієнт теплопровідності матеріалу стержня λ=50 Вт/м∙К.

Схема завдання граничних умов в поперечному перерізі стержня показана на рис.1

                                                                  

 

                      

                                                            


     

    

Рис. 1 Поперечний перетин стержня

           

qс10-5, Вт/м2

tж1,0С

tж2,0С

λ, Вт/мК

α1 ,Вт/м2К

α2,Вт/м2К

l1, мм

l2, мм

0,9

50

75

50

1200

2500

30

40

              

2. ПОБУДОВА РОЗРАХУНКОВОЇ СІТКИ

Дана задача двомірна з розмірами 4030мм. Побудуємо розрахункову сітку з постійним кроком розбиття = 10 мм = Dx = Dy. Диференційне рівняння теплопровідності заміняють рівнянням теплопровідності в кінечних різницях. Для того, щоб в розрахунковій сітці отримати вузлові крапки, що лежать безпосередньо на поверхні тіла, розбиття здійснюють з половинним кроком. Тоді в процесі побудови розрахункової сітки досліджуваної  області отримаємо,  що  елементарні об'єми уздовж меж області будуть удвічі менші (в кутках прямокутника - в чотири рази менше) внутрішніх об'ємів . При цьому серед сукупності вузлових точок розрахункової сітки будуть і такі, які належать бічній поверхні і її ребрам. Всім вузловим крапкам з невідомою температурою привласнимо  свій номер числами 1, 2, 3, 20, оскільки температура в кожній точці буде різниця одна від іншої.

 1 2 3 4                      5

   

 

  

 6 7  8 9 10  

 11 12 13 14 15   

    

16 17 18 19                   16

           

Рис. 2 Схема розрахункової сітки

                                                                                                                              

3. СКЛАДАННЯ РІВНЯНЬ ЕНЕРГОБАЛАНСА ДЛЯ ВУЗЛОВИХ

ТОЧОК СІТКИ

Складаємо рівняння  енергобалансу  для  точок  з  невідомою  температурою. Для  цього   для  кожного  із  обємів тіла вибираємо  вузлуву  точку . Після  цього  переходимо  від  диференціального  рівняння  теплопровідності  до  рівняння  теплопровідності  в  кінцевих   різницях  , для  чого  робимо  заміну   диференціалів  на  кінцеві  прирости.

      

 

                       Qж2,2           

1       Q1,2 2 Q3,2       3

                                        Q7,2

 6     7                  8

 

рис. 3 Схема енергобалансу 2 точки

Точка 1

Точка 2

Точка 3

Точка 4

Точка 5

Точка 6

Точка 7

 

Точка 8

Точка 9

Точка 10

Точка 11

Точка 12

Точка 13

Точка  14

Точка 15

Точка 16

Точка 17

Точка 18

Точка 19

Точка 20

З рівнянь енергобалансу отримаємо таку систему лінійних алгебраїчних рівнянь

  1.  
  2.  
  3.  
  4.  
  5.  
  6.  
  7.  
  8.  
  9.  
  10.  
  11.  
  12.  
  13.  
  14.  
  15.  
  16.  
  17.  
  18.  
  19.  
  20.  

4. РОЗВ'ЯЗОК СИСТЕМИ ЛІНІЙНИХ АЛГЕБРАЇЧНИХ РІВНЯНЬ

Для розв’язку даної системи використаємо програмний пакет MS EXEL. Розв’язки рівнянь знаходять за допомогою перемноження матриць.

Вихідна таблиця з коефіцієнтами при відповідних невідомих температурах та вільними членами рівняння.

Табл .1

Коефіцієнти в рівняннях для відповідних вузлових точок при температурах

Вільний член рівняння

t1

t2

t3

t4

t5

t6

t7

t8

t9

t10

t11

t12

t13

t14

t15

t16

t17

t18

t19

t20

1

-5

2

0

0

0

2

0

0

0

0

0

0

0

0

0

0

0

0

0

0

-111

2

1

-5

1

0

0

0

2

0

0

0

0

0

0

0

0

0

0

0

0

0

-75

3

0

1

-5

1

0

0

0

2

0

0

0

0

0

0

0

0

0

0

0

0

-75

4

0

0

1

-5

1

0

0

0

2

0

0

0

0

0

0

0

0

0

0

0

-75

5

0

0

0

25

-68,5

0

0

0

0

25

0

0

0

0

0

0

0

0

0

0

-1237,5

6

1

0

0

0

0

-4

2

0

0

0

1

0

0

0

0

0

0

0

0

0

-36

7

0

1

0

0

0

1

-4

1

0

0

0

1

0

0

0

0

0

0

0

0

0

8

0

0

1

0

0

0

1

-4

1

0

0

0

1

0

0

0

0

0

0

0

0

9

0

0

0

1

0

0

0

1

-4

1

0

0

1

1

0

0

0

0

0

0

0

10

0

0

0

0

25

0

0

0

50

-112

0

0

0

0

25

0

0

0

0

0

-600

11

0

0

0

0

0

1

0

0

0

0

-4

2

0

0

0

1

0

0

0

0

0

12

0

0

0

0

0

0

1

0

0

0

1

-4

1

0

0

0

1

0

0

0

0

13

0

0

0

0

0

0

0

1

0

0

0

1

-4

1

0

0

0

1

0

0

0

14

0

0

0

0

0

0

0

0

1

0

0

0

1

-4

1

0

0

0

1

0

0

15

0

0

0

0

0

0

0

0

0

25

0

0

0

50

-112

0

0

0

0

25

-600

16

0

0

0

0

0

0

0

0

0

0

1

0

0

0

0

-2

1

0

0

0

-18

17

0

0

0

0

0

0

0

0

0

0

0

2

0

0

0

1

-4

1

0

0

-36

18

0

0

0

0

0

0

0

0

0

0

0

0

2

0

0

0

1

-4

1

0

-36

19

0

0

0

0

0

0

0

0

0

0

0

0

0

2

0

0

0

1

-4

1

0

20

0

0

0

0

0

0

0

0

0

0

0

0

0

0

25

0

0

0

25

-56

-300

При розв’язуванні даної системи отримаємо такі результати:

      табл. 2

t1

115,13

t2

103,66

t3

96,01

t4

89,23

t5

82,47

t6

128,66

t7

116,09

t8

106,09

t9

96,33

t10

87,25

t11

131,35

t12

125,93

t13

115,91

t14

102,77

t15

91,72

t16

144,87

t17

140,39

t18

128,83

t19

107,13

t20

94,13

5. АНАЛІЗ ТЕМПЕРАТУРНОГО ПОЛЯ

На основі розрахованих локальних значень температур, зобразимо температурне поле у поперечному перерізі стержня (рис.3). Для цього накреслимо у масштабі поперечний переріз стержня, відмітимо положення вузлових точок і покажемо значення температур в них, а потім побудуємо з певним вибраним інтервалом сім’ю інтерпольованих ізотерм.

На графіку температурного поля (рис.3) нанесемо лінії теплового потоку ортогонально до ізотерм.

З рисунка видно, що температурне поле спадає справа наліво і зверху  вниз, тому що на правій частині нижньої грані та на нижній половині лівої грані немає теплообміну з навколишнім середовищем, а зверху і справа відбувається теплообмін з рідиною, температура якої менша температури стінки. Також ізотерми густішають в лівому нижньому куті, тому що на малому інтервалі розміщується велика різниця температур. Лінії теплового потоку рідкі у правому верхньому куті, так як ізотерми там менш густіші.

qc=const

15

15

qc=const

qc=const

qc=const


 

А также другие работы, которые могут Вас заинтересовать

43264. Себестоимость продукции, ее учет и пути снижения 236.5 KB
  Затраты живого и овеществленного труда в процессе производства составляют издержки производства. В условиях товарно-денежных отношений и хозяйственной обособленности предприятия неизбежно сохраняются различия между общественными издержками производства и издержками предприятия. Общественные издержки производства - это совокупность живого и овеществленного труда, находящая выражение в стоимости продукции. Издержки предприятия состоят из всей суммы расходов предприятия на производство продукции и ее реализацию.
43265. Процесс сравнения двух четырёх разрядных чисел 132.5 KB
  Подсказанная Максимом схема. Схема сравнения или компаратор от англ compre сравнивать операционный узел ЦВМ предназначенный для сравнения двух чисел А и В. Схема.
43266. Проектирование и реализация информационно-поисковой системы Виртуальный Деканат с помощью программы MySQL-Front 458 KB
  База данных представляет собой структурированную совокупность данных. Эти данные могут быть любыми - от простого списка предстоящих покупок до перечня экспонатов картинной галереи или огромного количества информации в корпоративной сети. Для записи, выборки и обработки данных, хранящихся в компьютерной базе данных, необходима система управления базой данных, каковой и является ПО MySQL. Поскольку компьютеры замечательно справляются с обработкой больших объемов данных, управление базами данных играет центральную роль в вычислениях. Реализовано такое управление может быть по-разному - как в виде отдельных утилит, так и в виде кода, входящего в состав других приложений.
43267. СУЩНОСТЬ, ФОРМЫ И ВИДЫ ПРЕДПРИНИМАТЕЛЬСКОЙ ДЕЯТЕЛЬНОСТИ 176 KB
  Неотвратимо осуществляются экономические реформы. Результатом реформ является становление и развитие новых экономических, финансовых, социальных и других отношений, базирующихся на формировании рыночной экономики, в которой ведущим субъектом хозяйствования являются предприниматели (коллективные и индивидуальные).
43268. Методы и методики перевода 226 KB
  Современные технологии политика наука и конечно литература – без перевода все это было бы недоступно для многих наций. Когда в мире появилось как минимум два языка наверное тогда человек впервые столкнулся с проблемой адекватности перевода. Когда знания людей о каком-либо виде деятельности стали обобщаться превращаясь в науку перевод не стал исключением.
43269. Сестринский процесс при ревматоидном артрите 982.5 KB
  Ревматоидный артрит – аутоиммунное заболевание неизвестной этиологии, для которого характерны симметричный эрозивный артрит (синовиит) и широкий спектр внесуставных (системных) проявлений. Ревматоидный артрит является одним из весьма распространенных заболеваний, причем заболеваемость им за последние 50—60 лет увеличилась. Через 5 лет от начала заболевания, несмотря на лечение „базисными“ препаратами, 16% пациентов теряют трудоспособность, а через 20 лет – 90%, треть всех пациентов становятся полными инвалидами.
43270. Історія, види та характеристики веб-браузерів 579 KB
  Але доступ до нього неможливий без багатьох факторів одним з яких є Webбраузери. Webбраузери це програмне забезпечення для комп’ютера під’єднаного до Інтернет що дає можливість працювати з текстом малюнками або іншою інформацією на вебсторінці. Поняття веббраузера Браузер це програма що представляє в зручному для сприйняття вигляді інформацію що отримується з Інтернету.
43271. Решение задач оптимизации 730.5 KB
  Цель задания: Найти минимум функции методами прямого поиска и градиентными методами. Цель работы отработка навыков решения задач безусловной оптимизации функции нескольких переменных методами прямого поиска и отработка навыков решения задач безусловной оптимизации градиентными методами. Решена задача безусловной оптимизации функции нескольких переменных методами прямого поиска и градиентными методами. Нахождение стационарной точки Целевая функция: ; Частные производные f по x1 и x2: ∂f ∂x1=2x1 x2 –6; ∂f ∂x2=x1 2x2 –4; Приравниваем...
43272. Синтез линейной дискретной антенны с равномерным амплитудным распределением 892.5 KB
  Акустическая антенна обычно состоит из электроакустических преобразователей элементов антенны акустических экранов несущей конструкции акустических развязок амортизаторов и линий электрокоммуникаций. По способу создания пространственной избирательности антенны можно подразделять на интерференционные фокусирующие рупорные и параметрические. Интерференционные антенны можно подразделять на непрерывные и дискретны. И непрерывные и дискретные антенны подразделяют по конфигурации геометрического образования объединяющего активные...