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


 

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

77883. Понятие, виды и исполнение обязательств в гражданском праве 28.5 KB
  Признаки обязательства: относительное гражданское правоотношение его субъектами являются строго определенные лица конкретный должник и конкретный кредитор; объектом обязательства является действие обязанного лица; содержанием являются права и обязанности сторон; опосредуют динамику гражданскоправовых правоотношений процесс перехода имущественных благ от одних лиц к другим; за невыполнение наступает гражданскоправовая ответственность; 6 реализация кредитором своего права возможно только через исполнение обязанности...
77884. Поручительство и банковская гарантия как способы обеспечения исполнения обязательств 30.5 KB
  Поручительство и банковская гарантия как способы обеспечения исполнения обязательств Поручительство способ по которому поручитель обязуется исполнить обязательство за должника. При поручительстве ответственным перед кредитором за неисполнение основного обеспечиваемого обязательства становится наряду с должником еще и другое лицо поручитель...
77885. Залог как способ обеспечения исполнения обязательств 28 KB
  Общие признаки: носят имущественный характер то есть интересы кредитора удовлетворяются за счет имущества должника; направлены на обеспечение удовлетворения интересов кредитора которое состоит в надлежащем исполнении ему обязательства; между должником и кредитором чаще всего существует два обязательства; возникает на основании закона или договора...
77886. Понятие и виды гражданско-правовых договоров 28.5 KB
  Свобода договора: недопущение понуждения к заключению договора; свобода выбора партнера при заключении договора; свобода участников гражданского оборота в выборе вида договора; свобода усмотрения сторон при определении условий договора кроме условий которые предусмотрены законом. Значение договора: договор представляет собой одно из самых уникальных правовых средств в рамках которого интерес каждой стороны в принципе может быть удовлетворен лишь посредством удовлетворения интереса другой стороны; договор...
77887. Формы и способы заключения договоров. Содержание гражданско-правовых договоров 31.5 KB
  Содержание гражданскоправовых договоров Содержание договора условия на которых достигнуто соглашение сторон: существенные условия которые необходимы и достаточны для заключения договора. Во-первых являются условия о предмете договора. Они включаются в текст договора по усмотрению сторон. Однако в отличие от обычных они приобретают юридическую силу лишь в случае включения их в текст договора.
77888. Договор купли-продажи. Обязанности 35 KB
  Продавец обязан передать товар в собственность покупателю надлежащего качества количества комплектности со всеми документами и принадлежностями если отказывается то предоставляется разумный срок не передал покупатель вправе отказаться от договора установленные законном или договором сроки пригодный для целей для которых товар такого рода обычно используется; товар свободный от 3 лиц либо уведомить покупателя об обременении данного товара....
77889. Договор дарения 30.5 KB
  Договор дарения Договором дарения называется договор по которому одна сторона даритель безвозмездно передает или обязуется передать определенное имущество другой стороне одаряемому либо освобождает или обязуется освободить ее от имущественной обязанности перед собой или перед третьим лицом. В отношении договоров дарения с участием юридических лиц ГК также предусматривает ряд специальных ограничений. Другое ограничение касается дарения вещей принадлежащих юридическому лицу на праве хозяйственного ведения или оперативного управления....
77890. Договор ренты 29 KB
  Договор ренты По договору ренты одна сторона получатель ренты передает другой стороне плательщику ренты в собственность имущество а плательщик ренты обязуется в обмен на полученное имущество периодически выплачивать получателю ренту в виде определенной денежной суммы либо предоставления средств на его содержание в иной форме. Реальный односторонний возмездный Возникающие из договора ренты отношения носят длительный стабильный а при пожизненном содержании с иждивением и доверительный характер. Форма подлежит...