壳程流阻1.1

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)