ЛБ7

Инженерная школа природных ресурсов

Направление подготовки Химическая технология

Отделение химической инженерии

ОСНОВЫ ПРОГРАММИРОВАНИЯ НА PYTHON

Отчет по лабораторной работе № 7

Функции и функциональные возможности языка Python

Выполнил студент гр. ХИМ54 А.В. Чижова

(Подпись)

01.01.2026 г.

Отчет принят:

Преподаватель

доцент ОХИ ИШПР, к.т.н. В.А. Чузлов

(Подпись)

_____ _____________ 2026 г.

Томск 2026 г.

Задание 1

Бензиновая фракция ( ) нагревается в теплообменном аппарате от 30 до 50 ⁰С.

1. Определить относительную плотность этой фракции при температуре из указанного диапазона с шагом 5 ⁰С.

2. Определить относительную плотность для следующих фракций:

№	Относительная плотность	Температурная поправка
1	0.7123	8.84e-4
2	0.7834	7.92e-4
3	0.8099	7.65e-4
4	0.8254	7.38e-4
5	0.999	5.15e-4

Программная реализация:

def Mendeleev_equation(

t: float,

rho_204: float,

alpha: float

) -> float:

return rho_204 - alpha * (t - 20)

rho_204 = 0.7486

alpha = 0.0008444

result = {}

for t in range(30, 56, 5):

result[t] = Mendeleev_equation(t, rho_204, alpha)

print(result)

rho_204 = [0.7123, 0.7834, 0.8099, 0.8254, 0.999]

alpha = [0.000884, 0.000792, 0.000765, 0.000738, 0.000515]

for i in range(0, 5):

print(Mendeleev_equation(50, rho_204[i], alpha[i]))

Ответ:

{30: 0.740156, 35: 0.7359340000000001, 40: 0.731712, 45: 0.7274900000000001, 50: 0.723268, 55: 0.7190460000000001}

0.6857800000000001

0.75964

0.7869499999999999

0.80326

0.98355

Задание 2

Составьте функцию для расчета средней молярной массы нефтепродукта по формуле Б.М. Воинова и рассчитайте молярные массы прямогонных бензиновых фракций, если их средние молярные температуры кипения t_{ср.м .} равны 115 и 132 °С, соответственно.

Программная реализация:

def Voinov_equation(

t: float,

) -> float:

return 60 + 0.3*t + 0.001*t*t

t1 = 115

t2 = 132

M1 = Voinov_equation(t1)

M2 = Voinov_equation(t2)

print(f'M1 = {M1}')

print(f'M2 = {M2}')

Ответ:

M1 = 107.725

M2 = 117.024

Задание 3

Компонент дизельного топлива имеет среднюю молярную температуру кипения 274 °С, его характеризующий фактор 10.8. Рассчитайте молярную массу компонента, реализовав соответствующую функцию.

Программная реализация:

def Voinov_with_K_factor_equation(

t: float,

K: float,

) -> float:

return (7*K - 21.5) + (0.76-0.04)*t + (0.0003*K-0.00245)*t*t

t = 274

K = 1038

M = Voinov_with_K_factor_equation(t, K)

print(f'M = {M}')

Ответ:

M = 30636.510199999993

Задание 4

Реализуйте функцию для расчета плотности чистого компонента для состояния насыщенной жидкости и определите с ее помощью плотность жидкого пропана в состоянии насыщения для t ∈ [30,50] °F с шагом 5 °F . Критические параметры пропана: T_c = 206.02 °F ; p_c = 616.12 psia

Программная реализация:

def fahrenheit_to_renkin(f):

return f + 459.67

def rho_sat(

t: float,

tc:float,

pc:float,

ZRA: float

):

return 1/((10.732*tc/pc)*ZRA**(1.0+(1.0-t/tc)**(2/7)))

Tc = 206.02

pc = 616.12

ZRA = 0.2766

t_range = range(30, 51, 5)

for t in t_range:

result = rho_sat(fahrenheit_to_renkin(t),fahrenheit_to_renkin(Tc), pc, ZRA )

print(f'Для T={t} плотность {result}')

Ответ:

Для T=30 плотность 0.750812499763869

Для T=35 плотность 0.7454215113988283

Для T=40 плотность 0.7399566770062249

Для T=45 плотность 0.7344143627682101

Для T=50 плотность 0.7287906361002597