计算精馏段提馏段塔板数目，终端输出结果和生成cvs文件

软件介绍

在大学的时候，在化工原理里面有怎么一门课，就是计算精馏段提馏段塔板数目，今天我特意把这个制作了可视化界面的代码提供给大家。

源码展示

import pandas as pd

xD = 0.98
xF = 0.4
xW = 0.02
R = 3.5
alpha = 2.46
q = 0.5
FD = (xD - xW) / (xF - xW)
xq = ((q - 1) * xD + (1 + R) * xF) / (q + R)
Yq = (q * xD + R * xF) / (q + R)
Y = xD

result = []

for i in range(1, 200):
    x = Y / (Y + alpha * (1 - Y))
    if x > xq:
        result.append((x, Y))
        print(f"精馏段 第{i}块板 x{i} = {x} Y{i} = {Y}")
        Y = R / (R + 1) * x + (xD / (R + 1))
    else:
        result.append((x, Y))
        print(f"提馏段 第{i}块板 x{i} = {x} Y{i} = {Y}")
        Y = (R + q * FD) / ((R + 1) + (q - 1) * FD) * x - (FD - 1) / ((R + 1) + (q - 1) * FD) * xW
        
    if x < xW:
        break

df = pd.DataFrame(result, columns=['x', 'y'])
df.to_csv('chemical_results.csv', index=False, header=False)

运行展示