Контрольная работа по дисциплине «Математические методы экспериментальных иследований в теплотехнике,теплоэнергетике и теплотехнологии» - файл

скачать (46.4 kb.)


Министерство образования и науки Российской Федерации

Федеральное государственное бюджетное образовательное учреждение

«Саратовский государственный технический университет имени Гагарина Ю.А.» Кафедра «Промышленная теплотехника"

Контрольная работа

по дисциплине «Математические методы экспериментальных иследований в теплотехнике,теплоэнергетике и теплотехнологии»

Выполнил:

студент группы б1-ТПЭНз41

ИнЭТС


форма обучения заочная

Шифр 181998

Максутов Алексей Владимирович

Проверил:

Серов Дмитрий Юрьевич

Саратов – 2022

Задача № 1

На основе уравнения y = f (x), выбираемого с учётом номера варианта N, составить обратную таблицу зависимости параметра x от y с шагом по y, равным hy в интервале значений y. Метод решения задачи выбирается согласно таблице 1.

Метод решения задачи – метод секущих

тепловой поток, передаваемый через основание шириной 8 мм продольного ребра вогнутого параболического профиля:















Код:


Function f(y, x As Double) As Double

f = 153.6 * (1 + 40 * x ^ 2) ^ 0.5 - y * x - 153.6

End Function
Sub mSec()

Dim y As Double

Dim x As Double

Dim x0 As Double

Dim x1 As Double

Dim iAs Integer

For i = 1 To 9

y = Range("A" + CStr(i + 1)).Value

x0 = 1

x1 = 2


x = x0 - f(y, x0) * (x1 - x0) / (f(y, x1) - f(y, x0))

Do While f(y, x) > 0.00001

x = x0 - f(y, x0) * (x1 - x0) / (f(y, x1) - f(y, x0))

x0 = x1


x1 = x

Loop


Range("B" + CStr(i + 1)).Value = x

Next


EndSub




Задача № 2

Выполнить балансовые расчёты многоступенчатой выпарной установки (МВУ): изобразить упрощённую схемы выпарной установки, обозначив на ней термодинамические параметры основных материальных потоков; записать уравнения материального баланса; определить концентрации и температуры кипения раствора в ступенях МВУ; записать систему уравнений теплового баланса МВУ, дополнив её уравнением баланса по выпариваемой воде Wi = W; решить полученную систему уравнений, вычислив уточнённое значение производительности каждой ступени МВУ по выпариваемой воде Wi одним из численных методов, выбираемых в соответствии с номером зачётной книжке согласно таблице 3.



Метод Зейделя



























































По давлению паров в секциях найдем температуры секций и энтальпии:

Давление, МПа

Температура,

Энтальпия, кДж/кг





































Общий вид системы для прямоточной МВУ:






Dim a(1 To 5, 1 To 5) As Double

Dim b(1 To 5) As Double

Dim x(1 To 5) As Double

Dim alf, bet, W, G As Double

Dim i As Integer

Sub Zeidel()

alf = 0.974

bet = 0.062

G = 5

W = 4.6


x(1) = 0.8364

x(2) = 0.8782

x(3) = 0.92

x(4) = 0.9618

x(5) = 1.0036

For i = 1 To 4

b(i) = 0.01 * alf - G * bet

Next


b(5) = W

a(1, 1) = alf - bet

a(1, 2) = -1

a(1, 3) = 0

a(1, 4) = 0

a(1, 5) = 0

a(2, 1) = -bet

a(2, 2) = alf - bet

a(2, 3) = -1

a(2, 4) = 0

a(2, 5) = 0

a(3, 1) = -bet

a(3, 2) = -bet

a(3, 3) = alf - bet

a(3, 4) = -1

a(3, 5) = 0


a(4, 1) = -bet

a(4, 2) = -bet

a(4, 3) = -bet

a(4, 4) = alf - bet

a(4, 5) = -1


a(5, 1) = 1

a(5, 2) = 1

a(5, 3) = 1

a(5, 4) = 1

a(5, 5) = 1
E = 0.0001

n = 5


Do

M = 0


For i = 1 To n

s = 0


For J = 1 To n

If i <> J Then s = s + a(i, J) * x(J)

Next

v = x(i)


x(i) = (b(i) - s) / a(i, i)

M = Abs(x(i)) - Abs(v)

Next

'Next


Loop Until M < E

' вывод результата

For i = 1 To 5

Range("H" + CStr(i + 3)) = x(i)

Next

End Sub





Задача № 3

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



Решение:


Dim a(1 To 4) As Double

Dim b(1 To 4) As Double

Dim c(1 To 4) As Double

Dim d(1 To 4) As Double

Dim e(1 To 4) As Double

Dim Ab(0 To 4) As Double

Dim Bb(0 To 4) As Double

Dim x(1 To 5) As Double

Dim i As Integer

Sub Progon()

Ab(0) = 0

Bb(0) = 0

x(5) = 0

a(1) = 0


a(2) = 0.0408

a(3) = 0.0477

a(4) = 0.0546

b(1) = 13.4

b(2) = 12.5

b(3) = 11.6

b(4) = 10.7

c(1) = 0.0581

c(2) = 0.065

c(3) = 0.0718

c(4) = 0

d(1) = 17.7828

d(2) = 19.0599

d(3) = 19.9744

d(4) = 20.5261

For i = 1 To 4

Ab(i) = -c(i) / (b(i) + a(i) * Ab(i - 1))

Bb(i) = (d(i) - a(i) * Bb(i - 1)) / (b(i) + a(i) * Ab(i - 1))

Next

For i = 4 To 1 Step -1



x(i) = -c(i) * x(i + 1) / (b(i) + a(i) * Ab(i - 1)) + (d(i) - a(i) * Bb(i - 1)) / (b(i) + a(i) * Ab(i - 1))

Next


Range("A8") = x(1)

Range("B8") = x(2)

Range("C8") = x(3)

Range("D8") = x(4)



End Sub





Учебный материал
© nashaucheba.ru
При копировании укажите ссылку.
обратиться к администрации