小型制冷系统用压缩机稳态性能仿真模型.docxVIP

小型制冷系统用压缩机稳态性能仿真模型 目前，对冷压缩缩缩缩缩缩缩缩缩缩缩缩缩缩缩缩缩差的研究可分为两类：一类是通用模型，主要用于空气压缩和系统仿真。另一类是对其研究感兴趣的部分建立模型,多用于压缩机的设计和改进。通用模型的建模方法一般有两种:一种是以压缩机性能曲线为基础,在压缩机性能测试工况下对压缩机性能参数如功率、制冷量等进行曲面回归,然后推广到任意工况。这种方法免去了物理模型参数的迭代求解,比较简单,但由于完全是在性能曲线上回归得出,模型没有物理意义,适用面较窄,在偏离测试标准工况时误差较大。另一种是对压缩机的控制容积建立能量方程、动量方程和质量方程,再辅以阀片的流量方程和进出气管及壳体的热量交换方程,从联立方程组求解压缩机各性能参数。这种方法物理意义明确,普适性强,但计算复杂,动态模型中对时间步长要求严格,否则难以收敛,因此难以应用到整个系统仿真中去。 本文将对小型制冷系统用全封闭往复式压缩机建立一个稳态和动态模型,并通过实验进行验证。 1 数学模型 1.1 热容计算模型 全封闭往复式压缩机内部热力过程如图1所示。为简化模型,对模型作如下假定: ① 压缩机转速为定值; ② 流经吸气管和排气管的压降忽略不计,故有p1=p2,p3=p4; ③ 忽略润滑油的传热; ④ 压缩机壳内制冷剂状态均匀。 通过上述假定,可以按如下方法计算压缩机的流量、功率、出口温度以及各部件的表面温度。 压缩机的质量流量按下式计算: 式中输气系数λ是压缩机流量计算是否准确的关键。 容积系数λV和压力系数λP按如下计算: 对全封闭小型往复式制冷压缩机,式中各常数根据实验取值为:C=1.43%,m=1.01,a=0.1,b=0.05。由于模型中质量流量是以吸气阀前比容v2计算,因此,温度系数λT和泄漏系数λl都取为常数:λT=0.97,λl=0.99。由制冷剂物性方程计算吸气阀前2点的熵s2,压缩机的等熵压缩功率Pis为: 其中,h3=f(p3,s3),s3=s2。因此,压缩机的实际功率为: 式中ηis,ηs,ηe分别是指示效率、轴效率和电机效率,取值分别为0.66,0.80,0.95。前述的排气阀后状态3点的焓h3是根据等熵压缩由物性方程计算,实际压缩过程是非等熵的多变过程,因此该点焓值应重新计算: 对稳态问题,压缩机进出口质量流量不变,并忽略各部件和壳内制冷剂的热容,压缩机吸气阀前状态2的比焓h2可由热平衡得出: 式中Qmr,Qcr,Qdr分别为电机、气缸和压缩机内部排气管传给壳内制冷剂的热量。 通常假定压缩机排气管在压缩机机壳内的传热与压缩机的功耗有关。 式中ζ为常数,其值与吸气温度有关。Qrs为压缩机壳内制冷剂传给壳体的热量,可通过壳体的热平衡决定。 Qrs=Qsa(12) 其中,Qrs=αiFi(Tr-Ts) Qsα=αoFo(Ts-Tα)+σεFo(Ts4-Tα4) 式中α为压缩机壳体对流传热系数,与壳体表面温度Ts有关,因此需迭代求解壳体温度Ts和制冷剂对壳体的传热量Qrs。同理,已知电机和气缸对制冷剂的传热Qmr和Qcr,可以迭代求解电机和气缸表面温度Tm和Tc。电机和气缸以及压缩机壳体均当作圆柱体处理,其对流传热系数按圆柱体对流传热计算。 由Qdr最终可计算压缩机排气管出口4的焓值h4: 从而根据制冷剂状态方程计算其出口温度T4。 由于在计算压缩机质量流量时需知道吸气阀前状态2的比容,而2点状态又要根据压缩机热平衡才能决定,显然需要迭代求解,其步骤如下:首先假定2点比焓h2,由p2根据状态方程可计算其比容v2和比熵s2;从而计算压缩机质量流量和耗功,进一步计算电机、气缸、排气管以及壳体与壳内制冷剂的传热及其表面温度;根据压缩机内部热平衡式(8)可计算出一个新的焓值h′2,与前述假定值比较,如果两者之差大于收敛标准,则调整先前假定值,重新计算,直至两者之差小于收敛标准。这样在已知压缩机进、出口压力和进口温度以及必须的几何尺寸后,可以按上述模型计算出压缩机稳定状态下的质量流量、消耗功率、出口温度以及电机、气缸和壳体的表面温度。模型中所有状态方程和物性方程采用NIST发布的商业计算程序REFPROP(V5.0ac)计算。 1.2 热容系数和排气热容量的确定 考虑到系统仿真中只关心压缩机的热力性能,并且由于压缩机的转速远远大于内部热交换的速率,丁国良等就指出,即使对压缩机的开机特性也可以忽略其流量的间歇性和压差的建立过程,只考虑其内部热平衡的建立过程。因此,笔者建立的动态模型只考虑压缩机内部热平衡的动态交换过程,着重压缩机出口温度动态特性的计算,并改善一般动态模型时间步长不能过大这一问题,以满足系统动态仿真的需要。 对于动态模型,由于前述的假定,压缩机质量流量和耗功的计算与稳态模型相同,但在计算壳内制冷剂与周围部件的换热时必须考虑各部件及制冷剂