Федеральное государственное бюджетное образовательное учреждение
«Саратовский государственный технический университет имени Гагарина Ю.А.» Кафедра «Промышленная теплотехника"
Контрольная работа
по дисциплине «Математические методы экспериментальных иследований в теплотехнике,теплоэнергетике и теплотехнологии»
Выполнил:
студент группы б1-ТПЭНз41
ИнЭТС
Шифр 181998
Максутов Алексей Владимирович
Проверил:
Серов Дмитрий Юрьевич
Саратов – 2022
Задача № 1
На основе уравнения y = f (x), выбираемого с учётом номера варианта N, составить обратную таблицу зависимости параметра x от y с шагом по y, равным hy в интервале значений y. Метод решения задачи выбирается согласно таблице 1.
Метод решения задачи – метод секущих
тепловой поток, передаваемый через основание шириной 8 мм продольного ребра вогнутого параболического профиля:
Код:
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
Do While f(y, x) > 0.00001
x = x0 - f(y, x0) * (x1 - x0) / (f(y, x1) - f(y, x0))
x0 = x1
Loop
Next
Выполнить балансовые расчёты многоступенчатой выпарной установки (МВУ): изобразить упрощённую схемы выпарной установки, обозначив на ней термодинамические параметры основных материальных потоков; записать уравнения материального баланса; определить концентрации и температуры кипения раствора в ступенях МВУ; записать систему уравнений теплового баланса МВУ, дополнив её уравнением баланса по выпариваемой воде Wi = W; решить полученную систему уравнений, вычислив уточнённое значение производительности каждой ступени МВУ по выпариваемой воде Wi одним из численных методов, выбираемых в соответствии с номером зачётной книжке согласно таблице 3.
Давление, МПа |
Температура, |
Энтальпия, кДж/кг |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Общий вид системы для прямоточной МВУ:
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(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
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, 2) = 1
a(5, 3) = 1
a(5, 4) = 1
a(5, 5) = 1
E = 0.0001
n = 5
M = 0
s = 0
If i <> J Then s = s + a(i, J) * x(J)
Next
v = x(i)
M = Abs(x(i)) - Abs(v)
Next
'Next
' вывод результата
For i = 1 To 5
Range("H" + CStr(i + 3)) = x(i)
Next
End Sub
Решить методом прогонки систему линейных алгебраических уравнений с трёхдиагональной матрицей. Система уравнений выбирается согласно варианту (таблица 5), который определяется по последней цифре номера зачётной книжки.
Решение:
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(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
Next
Range("B8") = x(2)
Range("C8") = x(3)
Range("D8") = x(4)