import numpy as np
import math
import json
m =(1, 2, 3, 4)#假设为壳程液体流速
shell_resistance = json.dumps({'Density_ave_tem':m})#定义json格式输入
def Shell_resistance(Shell_resistance1):
'''
热交换器壳程流阻
本函数作用于计算管壳式热交换器的壳程流阻
Keyword arguments:
input:
Cor_pre 换热管排列形式对压降的校正系数,当排布方式为站三角形
时,取0.5,正方形排列时取0.3,正方形斜转45°时取0.4
Bundle_center 横过管束中心的管束(根)
Fri_shell 壳程流体摩擦系数
Num_baffles 折流板数目
Density_ave_tem 壳程液体流速
V_tube 折流板间距
Factor_fvc 壳体直径
output:
T_pre_tube 管束阻力损失
Res_way 折流板弓形缺口处的阻力损失
Res_bend 壳程总体阻力损失
calculate process:
利用输入参数分别求出管束阻力损失、折流板弓形缺口处的阻力损失
壳程总体阻力损失为上述阻力损失之和
'''
#定义输入
Cor_pre = float(input('换热管排列形式对压降的校正系数:\t'))
Bundle_center = float(input('横过管束中心的管束:\t'))
Fri_shell = float(input('壳程流体摩擦系数:\t'))
Num_baffles =float(input('折流板数目:\t'))
V_tube = float(input('折流板间距:\t'))
Factor_fvc = float(input('壳体直径:\t'))
#将文件格式还原为python格式
resistance = json.loads(Shell_resistance1)
print(f'壳程流阻:{resistance}')
velocity = np.array(resistance['Density_ave_tem'])
T_pre_tube = (Cor_pre*Fri_shell*Bundle_center*(Num_baffles+1)*
velocity**2)/2#计算管束阻力损失
# print(f'管束阻力损失:{T_pre_tube}\t')
Res_way = (Num_baffles*(3.5-2*V_tube/Factor_fvc)*velocity**2)/2
#print(f'折流板弓形缺口处的阻力损失:{Res_way }\t')
Res_bend = T_pre_tube + Res_way
#print(f'壳程总体阻力损失:{Res_bend}\t')
#将数组转换为元组方便作为字典的值
T_pre_tube = tuple(T_pre_tube)
Res_way = tuple(Res_way)
Res_bend = tuple(Res_bend)
#将输出定义为字典,方便json打包
output_resistance ={'管束阻力损失':T_pre_tube,'折流板弓形缺口处的阻力损失':Res_way,'壳程总体阻力损失':Res_bend}
print(output_resistance)
#json打包
output = json.dumps(output_resistance)
Shell_resistance(shell_resistance)
壳程流阻1.1
最新推荐文章于 2021-03-17 20:08:45 发布