首页动力电气制冷与低温工程制冷与技术 16 5.1 制冷空调中的计算机仿真.1.ppt

16 5.1 制冷空调中的计算机仿真.1.ppt

页数:145页文件大小:2.85 MB来源: 制冷与技术ppt课件

本课件介绍制冷空调中的计算机仿真,先说明物理仿真与计算机仿真的区别及简单对象建模方法,再展开单级压缩蒸气制冷理论循环的计算机分析与整机模拟,详细给出压缩机、冷凝器、蒸发器、毛细管等部件模型与算法,最后介绍装置优化、计算机辅助设计与常用仿真软件。

第一节制冷空调中的计算机仿真构成说明,对应第5章导览部分的本页内容
仿真技术简介仿真用一个能代表所研概念说明,按照模型性质不同物理仿真计算机仿真
物理仿真用一个与实际系统物理本质相同要点说明,计算机仿真用数学形式表达实际系统
查看本课件的文字版(共 145 页,供检索与复制)

第1页:第一节制冷空调中的计算机仿真

第5章制冷空调中的计算机仿真与控制第一节制冷空调中的计算机仿真第二节制冷空调计算机控制系统的组成

第2页:仿真技术简介仿真用一个能代表所研

5.1.1 仿真技术简介仿真用一个能代表所研究对象的模型去完成的某种实验, 以前常称为模拟。按照模型性质不同物理仿真计算机仿真

第3页:物理仿真用一个与实际系统物理本质相同

物理仿真用一个与实际系统物理本质相同的模型去完成实验。计算机仿真用数学形式表达实际系统的运动规律,数学形式通常是一组微分方程或差分方程,然后用计算机来解这些方程。

第4页:写出实际系统的数学模型

在计算机仿真研究的过程中,一般要经过这样四个步骤

  • (1) 写出实际系统的数学模型。
  • (2) 将它转变成能在计算机上进行运转的数学模型
  • (3) 编出仿真程序
  • (4) 对仿真模型进行修改、校验

第5页:仿真系统有无实物介入实时

仿真系统有无实物介入实时仿真系统非实时仿真系统仿真计算机类型不同用模拟计算机组成的仿真系统用数字计算机组成的数字仿真系统用混合模拟机组成的或用数字-模拟混合计算机组成的混合仿真系统微型机阵列组成的全数字式仿真系统

第6页:简单对象的建模

5.1.2 简单对象的建模

  • 在制冷空调装置仿真中,有些部分在一定假设下,可用一阶微分方程近似描述。下面举例说明。

第7页:例5-1 货物冷却对于货物送入冷藏箱中进

例5-1 货物冷却对于货物送入冷藏箱中进行冷却,如图5-1所示。设冷藏箱中空气温度为;设货物的温度为,质量为M,定容比热为C,与空气传热面积为F,货物与空气的当量传热系统为K。

第8页:货物的蓄热量

货物的蓄热量U为(5-1)传给货物的热量应等于货物蓄热量的变化(5-2)将式(5-1)代入(5-2)并整理得(5-3)上式即是包含对t 求导的一阶微分方程。反映了一定条件下,货物随冷藏室内空气温度的变化规律

第9页:用一阶微分方程描述

用一阶微分方程描述的只能是非常简单与理想化的对象,在制冷空调装置仿真中,如果考虑稍多一些影响参数的话,则必须采用更高阶的方程。例5-2变空气温度下的货物冷却仍然是货物送入冷藏箱中进行冷却的过程计算。与例5-1不同的是,空气温度是变化的,而送入箱内的热量是一定的, 设为Q。设冷藏箱中空气温度为,质量为Ma,定容比热为;设货物的温度为,质量为M,定容比热为C,与空气传热面积为F,货物与空气的当量传热系统为K。货物送入冷藏箱中进行冷却,箱体结构为绝热。

第10页:

第11页:空气的蓄热量

空气的蓄热量U为货物的蓄热量U为传给货物的热量应等于货物蓄热量的变化传给空气的热量与传给货物的热量之和为总热量Q由式(5-6)得(5-4) (5-5) (5-6) (5-7) (5-8)

第12页:将(5-8)代入(5-7)得

将(5-8)代入(5-7)得,(5-9)上面的二阶常微分方程描述了冷藏箱内货物的冷却过程。如果考虑空气与箱体结构的传热,而把箱体结构作为一阶惯性环节,则得到的式子为三阶微分方程。如果对于厚的货物,需要考虑表层与内部温度变化的不一致,则所得到的方程阶数还要高。

第13页:一般地

一般地,描述系统的高阶微分方程可统一用如下形式(5-10) 对于一般的微分方程,难以直接求得分析解,一般采用数值求解方法。对于精度要求较低而速度要求较高的场合,可以采用欧拉法、梯形法;如果精度要求较高,则四阶龙格库塔法是常用的求解方法。

第14页:单级压缩蒸气制冷理论循

最常见的制冷装置如家用冰箱、家用空调器等均采用单级蒸气压缩制冷循环5.1.3 单级压缩蒸气制冷理论循环的计算机分析对于单级蒸气压缩制冷理论循环的计算机分析是一种非常简化的制冷循环模拟,可以作为实际制冷装置模拟的基础。

第15页:图5-3 示出了单级蒸气压缩制冷循环

图5-3 示出了单级蒸气压缩制冷循环的lgp–h图。查表可以计算出所要求的各个量,但每次计算都比较复杂。用计算机计算虽然编程需要花时间,但以后每次计算特别快,这对于工况等参数改变时的分析特别能体现出其优势。

第16页:假定输入参数

假定输入参数为4个:蒸发温度Te,冷凝温度Tc,压缩机吸气过热度ΔTe,冷凝器过冷度ΔTc。按理论循环的假设条件,蒸发温度和冷凝温度均为定值,系统的流动阻力忽略不计。压缩过程为等熵过程,节流过程为等焓过程。循环的制冷量(5-11) 单位容积制冷量(5-12) 单位理论热负荷(5-13)制冷系数(5-14)

第17页:图5-4 为计算单级蒸气压缩制冷循环性能

图5-4 为计算单级蒸气压缩制冷循环性能的程序框图。由Te求peT1= Te+ΔTe, p1 =pe由T1, p1求v1, s1, h1由Tc求pcp2 =pc,s2=s1由p2, s2求T2, h2T4= Tc-ΔTc, p4 =pc由T4, p4求h4求q0, qv, qk, w0, ε结束给Te, Tc, ΔTe, ΔTc赋值

第18页:上述程序的用途因为该种计算中只需要知道制

上述程序的用途因为该种计算中只需要知道制冷工质的热力性质,与工质的传输性质以及具体的装置结构均无关所以可以方便地求出当蒸发温度、冷凝温度、压缩机吸气过热度、冷凝器过冷度变化时,理论制冷循环性能的变化现经常被用来比较不同工质的性能

第19页:使用上述方法存

使用上述方法存在的问题因为对于一般的制冷装置来讲,当蒸发温度、冷凝温度变化时,其压缩机吸气过热度、冷凝器过冷度也会变化,定值假定是不符合实际情况的。所以上面分析过程没有牵涉到外界环境对于实际装置的影响方法虽然简单,但同实际装置性能之间是有差距,不能预测外界环境变化时制冷装置的性能变化。

第20页:单级压缩蒸气制冷装置的计算机模拟

5.1.4 单级压缩蒸气制冷装置的计算机模拟5.1.4.1 部件模型在计算机模拟时,并不能任意指定状态,如蒸发温度、冷凝温度、过热度、过冷度,而是应该能把这些参数正确地计算出来。在模型和算法的选取上,应当根据实际需要,在精度、计算稳定性和运算速度之间达到平衡。

第21页:制冷装置部件模型的建立

对于一个简单的单级蒸气压缩制冷装置,设其由往复活塞式压缩机、毛细管、冷凝器与蒸发器这四大件组成。蒸发器与换热器均采用干式换热器,其本身热容可以忽略不计,这两个换热器均采用温度不变的空气冷却。建立各个部件的模型

  • 1.压缩机模型
  • 2.毛细管模型
  • 3. 蒸发器和冷凝器模型
  • 4. 围护结构模型
  • 5. 充注量计算模型要求模拟压缩机开机过程到系统接近稳定的整个过程则主要是要预测制冷剂状态及制冷量随时间的变化

第22页:压缩机模型对于制冷装置来讲

1. 压缩机模型对于制冷装置来讲,活塞在一个运转周期中的流量的变化,是一个频率过高的信号,可以取每个周期的平均值来滤掉该高频信号压缩机进出口状态对于压缩机流量的影响是没有时间迟延的压缩机流量计算的模型可以采用稳态模型,功率则可根据理论功和压缩机的效率确定.所以

第23页:(5-15)(5-16)上面式中

(5-15)(5-16)上面式中,mcom,N分别表示压缩机的制冷剂流量与功率;

  • ,Vh,λ分别为压缩机的输气系数、理论功率;pc,pe,vs,m分别表示冷凝压力、蒸发压力、吸气比容、多变指数。

第24页:毛细管模型毛细管中制冷剂的流速很

2. 毛细管模型毛细管中制冷剂的流速很高,制冷剂流过毛细管所需要的时间也远小于系统的时间常数,因此毛细管进出口状态的影响也可以认为是即时的其模型采用稳态模型因为管内流体流动的高度非线性,各种较为精确的分布参数模型在数值求解时速度较慢且存在计算的稳定性问题所以建立精确,同时又简单、通用的毛细管模型对于实际装置的设计与优化具有重要意义

第25页:对于一维等焓均相流动

对于一维等焓均相流动,有如下控制方程(5-17)式中,p, v, G分别为流体的压力、比容和质流密度,D和L分别为毛细管内径和长度,f为沿程摩阻系数。下面介绍的绝热毛细管的近似积分模型是一种较好的模型。

第26页:过冷区模型过冷区液体比容和沿程

(1) 过冷区模型过冷区液体比容和沿程摩阻系数可认为不变,对上式积分,得过冷区长度(5-18)式中,ΔpSC表示过冷区压降,下标SC表示过冷区。

第27页:两相区模型用

(2) 两相区模型用p1和v1表示两相区的进口压力和比容,p2和v2表示两相区的出口压力和比容。建立如下经验方程(5-19)因沿程摩阻系数f变化不大,故在积分过程中设为定值,取进出口摩阻系数之算术平均。得二相区长度(5-20)k1是一个仅与边界条件相关的常量(5-21)

第28页:过热区模型对于低压下的过热气体

(3) 过热区模型对于低压下的过热气体,可近似看作理想气体。因此在等焓过程中温度不变(5-22)式中,T和R分别是绝对温度和气体常数。由式(5-22)得(5-23)将式(5-22)和(5-23)代入方程(5-17)并积分,得过热区长度(5-24)式中,下标1和2分别表示过热区的进口和出口参数。

第29页:壅塞流当工质在毛细管出口处的流

在实际计算中,为方便起见,取(4) 壅塞流当工质在毛细管出口处的流速达到当地音速时,毛细管处于壅塞流动。此时毛细管出口压力大于或等于背压背压的降低对毛细管质流率已无影响。此时的质流率GC称为毛细管的壅塞质流率或临界质流率,可按式(5-25)至(5-27)计算

第30页:(5-25

(5-25)(5-26)(5-27)式(5-25)至(5-27)表明毛细管的临界质流量只是当地干度和制冷剂热物性的函数,而与毛细管结构尺寸无关。式(5-26)和(5-27)可以由制冷剂热物性数据拟合成关联式。另外,为了简化计算,若在过冷流动或过热流动中发生壅塞,分别按饱和液体和饱和气体处理。

第31页:其他参数的确定对于毛细管流动

(5) 其他参数的确定对于毛细管流动的沿程摩阻系数f 的计算,采用Churchill关联式:(5-28)上面关联式可覆盖整个Re数区域,且考虑了毛细管内粗糙度的影响,一般毛细管相对粗糙度约为3.27×10μ4。对于两相区的动力粘度μTP按下式计算。(5-29)

第32页:管长计算

(6) 管长计算

  • (7)质流量计算在装置仿真中,毛细管的结构尺寸都是已知的,而需要求得的是流量等参数。其基本实现步骤如下:在进口状态及出口背压已知条件下先要确定进口有无过冷,过冷度有多大一般情况:毛细管进口为过冷,出口为二相管长= 过冷区管长+ 二相区的管长其它情况:先确定存在哪几相总的管长=各相的长度之和

第33页:步骤1:假设毛细管的出口压力等于其背压

步骤1:假设毛细管的出口压力等于其背压,结合进口条件,确定毛细管内是否存在过冷、两相或过热流动区域及存在的各流动区域的进、出口状态,并求出毛细管出口为背压时的壅塞质流率G0。步骤2:假定毛细管的流量为G0,对于存在的各流动区域,计算该区域的长度,并将不同流动区域的计算长度相加后得到毛细管的计算长度。步骤3:将毛细管的计算长度与实际长度比较。若计算长度在误差限之内,则毛细管出口的压力等于背压,质流率等于G0。若计算长度偏长,则说明实际质流率大于G0,毛细管的出口压力高于背压,此时需要重新假定新的出口压力,重复以上的过程。若计算长度偏短,则说明实际质流率小于G0,不出现壅塞,出口压力等于背压,此时只要在小于G0的质流率范围内搜索一个正确的质流率。

第34页:蒸发器和冷凝器模型

3. 蒸发器和冷凝器模型(5-30) (5-31) 建模与求解中忽略蒸发器与冷凝器中制冷剂的阻力损失,制冷剂两相区的温度可近似认为是一致的因此系统不必采用分布参数模型,只要将两器按过冷、二相、过冷分成几个大块即可。对于冷凝器,根据制冷剂的质量和能量守恒方程式,

第35页:其中,M,h,m分别为制冷剂

其中,M,h,m分别为制冷剂的质量、焓和质流率;q为总的热流;下标SH,TP和SC分别表示换热器的过热区、两相区和过冷区。令(5-32) (5-33) 式(5-30) 和(5-31) 在一个短的时间步长内积分得:(5-34) (5-35) 式中,上标1和0分别表示当前时刻和上一时刻的物理量。

第36页:当进出口流量、进口焓值已知时

当进出口流量、进口焓值已知时,冷凝器中其它参数仍然需要通过迭代才能确定。对于上述模型进行求解的一种较为稳定的算法是质量引导法,把质量平衡作为迭代标准。

第37页:冷凝压力估算与能量守恒

估计一个冷凝压力根据能量守恒方程式计算出高压侧制冷剂的状态和质量,从而可得高压侧的制冷剂总质量将该值和由式(5-34) 计算出的质量值进行比较误差小于允许范围yes依次计算出其他状态参数no

第38页:对于蒸发器

对于蒸发器,完全可以采用同样的方法,只是在蒸发器中没有过冷区而已

第39页:围护结构模型制冷装置的性能不仅取

4. 围护结构模型制冷装置的性能不仅取决于制冷系统的特性,而且还跟围护结构的性能密切相关。(1)反应系数法与Z传递系数法计算围护结构特性的原理比如冰箱:

  • 制冷系统在5分钟左右就达到基本稳定,但整个装置基本上没有稳定的时候,主要因素是因为围护结构动态特性的作用。对于一个只有一样材料组成的最简单的围护结构,可以看成如图5-5所示的单层均质平壁热力系统,除导热方程外,还有与热流有关的导热定律:

第40页:(5-36) 平壁两侧表面上有四个时间函

(5-36) 平壁两侧表面上有四个时间函数:内表面温度图5-5 平壁热力系统内表面热流外表面温度外表面热流

第41页:其中两个量给定,另两个量待求

其中两个量给定,另两个量待求。现假定外侧表面上的温度和内表现的热流为已知,内侧温度和外侧为未知,采用过余温度表示,初始状态设为零,数学模型为:(5-37 )对于上述微分方程可通过差分进行数值求解,计算每一时刻的各个参数值,计算量很大。实际的围护结构大多为由多种材料组成,方程更为复杂,求解的量更大。

第42页:围护结构建模方法

适宜于系统仿真的围护结构建模方法

  • 1)谐波法(与正弦传递函数相对应)。
  • 2)反应系数法(与S传递函数相对应)
  • 3)传递系数法(与Z传递函数相对应)这些方法都把扰量和围护结构本身的传递特性分开处理,先求出反映围护结构本身特性的有关参数,最后计算系统的动态响应时,只需要将这些已经计算求得的参数同扰量进行合成。由于对围护结构只计算一次,所以计算量可以大减少。

第43页:在反应系数法中

在反应系数法中,假定室外温度变化引起室内温度和室外热流变化的反应系数分别为、,而室内热流变化引起室内温度与室外热流变化的反应系数分别为、。计算的时间步长为则在第时刻的室内温度与室外热流分别为的输出值为(5-38)(5-39)

第44页:Z传递函数的定义

Z传递函数的定义(5-40) 为保证分子、分母的系数唯一,取定分母多项式的首项恒为。输出函数的Z变换/输入函数的Z变换对于平壁热力系统,其Z传递函数记作

第45页:(5-41

(5-41)(5-42) 反应系数法项数得取较多而Z传递函数所取系数少得多。如果只考虑室内热量引起温度变化的关系, 只要先求出对应此两个参数输入输出关系的Z传递函数, 确定了此函数的各个系数, 则有:

第46页:状态空间法求反应系数在状态空间

(2)状态空间法求反应系数在状态空间法中使用标准形式的状态方程和输出方程,如下所示。状态方程(5-43) 输出方程(5-44)

第47页:对于平壁围护结构

对于平壁围护结构,为了建立状态空间,将平壁适当分层,作为一个n层的集中热容系统处理(见图5-6),从而可建立起一个m维(m= n+ 1)的状态空间。图中,阴影部分分别表示内外边界及内部的控制体。

第48页:(5-45)由能量守衡知控制体

(5-45)由能量守衡知控制体的内能变化等于进出控制体的热流量的代数和,由此可以列出一组常微分方程,即状态方程。

第49页:所要求的内表面热流为

所要求的内表面热流为(5-46) 上面式中的符号为热容,热阻,:厚度,:材料导热系数,:密度,:比热,:墙体外侧和内侧的空气换热系数。

第50页:当我们进行吸热反应计算时

当我们进行吸热反应计算时,内表面的过余温度,将输入输出关系整理成标准的状态方程和输出方程,如式(5-43)、(5-44)。主要的参数为各状态点温度,各状态点温度变化率,内表面热流,室外温度变化,

第51页:

第52页:根据前面的状态方程和输出方程

根据前面的状态方程和输出方程,可以求解在一定扰量作用下的系统的参数输出。在计算单个反应系数时,系统的输入是单个量,输出也是单个量。只要能构造与前面定义的反应系数相一致的输入,则所得到的输出即为相应的反应系数对于三角波反应,可以由斜坡反应构成。设U为一个斜坡扰量,如能根据前面的状态方程和输出方程,求出时间间隔为时的响应系数列,根据线性迭加原理,即可求出三角波反应系数(5-47)

第53页:对于状态方程

对于状态方程(5-45),其解的一般形式为:(5-48)式中,称为矩阵指数,与其相关的积分在本书中统称为矩阵指数的积分。若设为离散化时间步长,并在上式中分别令则可以求得(5-49)

第54页:上面的解中既有自由项,又有强制项

上面的解中既有自由项,又有强制项,计算复杂。如能把控制量增广到状态量中去使状态方程变成齐次的,求解就简便多了。对于为斜坡函数的情况,增广是能够实现的。令(5-50)(5-51) 从而构成齐次的增广状态方程(5-52)可简记为(5-53)

第55页:对于齐次状态方程

对于齐次状态方程(5-54) 其解的形式为(5-55) 取时间步长为, 可得如下的递推方程:(5-56) 矩阵指数及其积分的计算方法有很多,下式所示的直接级数展开程序实现比较简单。(5-57) 式中I为单矩矩阵。

第56页:状态空间法求

(3) 状态空间法求Z传递系数对于方程(5-43),其解的离散形式为式(5-49),但是除了一些特殊输入函数外,该式无法直接用于计算。因此,有必要在保证一定精度的条件下,采取一些近似方法。常用的近似方法有两类保持器法中零阶保持和一阶保持比较简单,高阶保持比较复杂,而且对于不同的扰量输入,精度并不与保持器的阶数成正比,因此常用零阶保持和一阶保持。对于零阶保持,数学上表述为(5-58)保持器法数值积分法

第57页:代入式(5-49)可得(5-59

代入式(5-49)可得(5-59) 式中对于一阶保持,数学上可表示为(5-60)

第58页:代入式(5-59

代入式(5-59),可得(5-61) 式中,数值积分法就是对式(5-59)右边第二项直接进行数值积分。这里给出梯形公式的结果:(5-62) 式中,式(5-61)与式(5-62)形式相同

第59页:要完成上面的计算,需要计算下面三个矩阵指

要完成上面的计算,需要计算下面三个矩阵指数。(5-63)(5-64)(5-65)

第60页:完成状态方程的求解后

完成状态方程的求解后,结合输出方程的离散化形式,得平壁的离散状态空间模型:(5-66)(5-67)(5-68) 式(5-67)和(5-68)可以通过线性变换转化为式(5-66)的形式。

第61页:对于一个n+1阶的系统

对于一个n+1阶的系统,Z传递函数的形式如下:(5-69) 问题归结为如何确定分子与分母中的系数与。以标准离散状态空间模型(5-66)为例,取Z变换,(5-70)

第62页:整理得

整理得,(5-71) 式中为矩阵行列式为伴随矩阵为阶常数阵,及由Leverrier-Faddeeva 算法确定, 这样也就确定了

第63页:(5-72

(5-72)(5-73) 如果已知反应系数序列Y(i),则亦可由下述关系简捷地求得:(5-74)具体计算按式(5-72)、(5-73):

第64页:充注量计算模型制冷剂充注量与制冷

5. 充注量计算模型制冷剂充注量与制冷装置的工作特性是紧密相关的,对于制冷装置,适宜的制冷剂充注量是非常重要的。

第65页:对于一个典型的小型制冷装置,制冷剂量可一

对于一个典型的小型制冷装置,制冷剂量可一般地表示成如下形式:(5-75) 上式中等式右边各项分别对应蒸发器二相区,蒸发器过热区(包括回气管),冷凝器二相区,冷凝器过热区,冷凝器过冷区,压缩机空腔,干燥过滤器和润滑油。

第66页:空泡系数与质量计算

为什么研究空泡系数: 质量计算需要对于单相区的制冷剂密度容易确定,但对要计算二相区的制冷剂密度,则必须计算空泡系数。空泡系数跟干度有关。式中,A是流道内截面积,LTP是两相区长度。

第67页:空泡系数定义又称为截面含气率或真实含气

空泡系数定义

  • 又称为截面含气率或真实含气率, 指两相混合物在任一流动截面内气相所占的总面积份额A、Ag分别表示流道面积与气体流通面积

第68页:干度定义也叫质量含气率

干度定义

  • 也叫质量含气率, 是指单位时间内流过流道截面的两相流总质量中, 气相质量所占的份额, 其定义式为
  • 式中, M、Mg、Ml 分别表示总的两相流质量流率以及气相、液相的质量流率。

第69页:在传热计算中

在传热计算中,制冷剂质量的计算不能直接利用干度来进行,而需要由空泡系数来确定二相区制冷剂的密度可用下式来表示

第70页:干度和空泡系数关系

干度和空泡系数关系

  • S为滑动比。

第71页:空泡系数模型分为均相模型滑动比修正

空泡系数模型分为均相模型滑动比修正Xtt修正考虑质流率的模型

第72页:均相模型两相均匀混合,滑动比

均相模型两相均匀混合,滑动比为1

第73页:Zivi滑动比模型

Zivi滑动比模型

  • 导出条件: 无流体夹带的环状流, 在管壁摩擦为零, 熵增为零

第74页:Smith滑动比修正模型式中K为夹带系

Smith滑动比修正模型

  • 式中K为夹带系数,推荐值为0.64。
  • 导出条件: 基于均匀混合物核心与环状液相具有相等的速度头的假设

第75页:Xtt 修正模型

Xtt 修正模型(5-85)(5-86) 式中,(5-87)(5-88)

第76页:考虑质流率的修正模型-

考虑质流率的修正模型-Tandon模型式中

第77页:考虑质流率的修正模型-

考虑质流率的修正模型-Premoli模型

  • 这是一个经验修正模型, 它是能过滑动比的计算进行的, 滑动比的计算过程如下

第78页:考虑质流率的修正模型-

考虑质流率的修正模型-Hughmark模型

  • 式中KH=f(Z),其具体关系式见教材
  • 适用范围: 很广, 许多文献中均采用这些关系式

第79页:.2 稳态仿真制冷空调装置的系统仿真

5.1.4.2 稳态仿真制冷空调装置的系统仿真,是将部件模型组合一个有机的整体,以表现实际装置的特性。根据不同的对象和不同的研究目的,可以对部件模型进行不同的组合。

第80页:图5-8

对于稳态仿真,以空调器设计企业设计需要为目的仿真包括两种算法:第一种算法是已知蒸发器过热度、冷凝器过冷度(毛细管的内径和并联数给定,其他结构参数与环境参数也已知),求整机的充注量和毛细管长度,以及制冷量、压缩机功率等。算法如图5-7所示。第二种算法是已知系统充注量和毛细管长度(毛细管的内径和并联数给定,其他结构参数与环境参数已知),求系统性能(制冷量、压缩机功率、蒸发器过热度、冷凝器过冷度等),算法如图5-8

第81页:图5-7 系统稳态仿真算法流程图一

图5-7 系统稳态仿真算法流程图一(过热度、过冷度为输入,充注量、毛细管长度为输出)

第82页:图5-8 系统稳态仿真算法流程图二

图5-8 系统稳态仿真算法流程图二(充注量、毛细管长度为输入,过热度和过冷度为输出)

第83页:.3 动态仿真下面结合电冰箱

5.1.4.3 动态仿真下面结合电冰箱动态仿真进行介绍。电冰箱中,各个参数间的相互影响关系可分成两类通过制冷剂质量流动发生的各部件间的参数联系通过热量的传递发生的各部件参数的联系

第84页:图5-9制冷系统进出口参数耦合图毛细管蒸

图5-9制冷系统进出口参数耦合图毛细管蒸发器冷凝器压缩机mcaphcaphcoheopcpehcommcom

第85页:箱内参数联系图从压缩机进口断开进行分析

图5-10 箱内参数联系图从压缩机进口断开进行分析

第86页:图5-11简化示意图

图5-11简化示意图图5-12函数的交点可以归结为求某一输入x,使得y=x,即求出函数y=y(x)与y=x的交点。

第87页:当已知两点

当已知两点(x1, y1)、(x2, y2)时, 根据线性插值, 可得图5-12中的x3值。(5-100) 当已知三点后,可以通过拉格朗日插值公式,确定一条二次曲线,用它和y= x的交点作为新的估计值.经过推导,得到该点为(5-101) 式中,

第88页:实际使用时

实际使用时,x为估计的箱内空气温度值,y为在此估计值下经过蒸发器、试验包、箱体这几部分计算后所得的箱内空气温度值。当用程序求得y=x的点时,则找到了正确的箱内空气温度值。

第89页:制冷装置优化与计算机辅助设计简介

5.1.5制冷装置优化与计算机辅助设计简介5.1.5.1 优化的含义制冷空调装置的优化首先要使装置设计最佳,其次要保证系统能够工作在最优的工作状态下。制冷空调装置的优化包括最优设计最优控制

第90页:制冷装置的优化原则

制冷装置的优化原则

  • l首先要确定优化的原则:1)优化目标,2)优化参数,3)优化计算的约束条件,
  • l然后才是优化的方法的确定。

第91页:优化目标的确定对于不同的装置

1. 优化目标的确定对于不同的装置,不同的人员,所选择的优化目标都会有所不同。一般来讲,优化的目标应该包括:装置能够正常工作,达到其功能要求效率与经济性最高

第92页:优化参数的选择

优化参数的选择

优化参数
优化计算中的可变量。改变这些参数,寻找其最佳组合,即是优化计算过程。
  • 连续取值的优化参数:毛细管的管长,管板式换热器的散热面积等;
  • 不连续取值的优化参数:只能在有限个类型中进行选择,如压缩机的容量大小,冷凝器与蒸发器的管径与外表面的面积,膨胀阀的容量等。
  • 如果选择太多的参数作为优化参数必然使得计算十分复杂,在参数的选择上,要兼顾各种因素。

第93页:约束条件的选取适当选择约束有二个作用

约束条件的选取

  • 适当选择约束有二个作用
  • 1)实际装置各参数值的优化都必须在一定范围内进行,超过这个范围得到的优化值是毫无意义的。
  • 2)当参数可变化范围增大时,可能出现多个极值,寻优过程在不为最值的某一极值处停止。数学模型的准确性有一定范围,如超出适用范围,模型的精确度就要降低,因此在优化计算时,有时还需要人为地定一些约束条件,以使优化计算有效地搜索。
  • 对于第1)类约束条件,它的存在会使得计算时间变大、迭代次序增加。而第2)类约束条件有利的。

第94页:建立在动态仿真基础上的制冷装置优

5.1.5.3 制冷装置优化方法1.建立在动态仿真基础上的制冷装置优化对优化方法的要求一般说来,利用函数梯度信息的优化方法的寻优速度较快。但在实际应用中,此类方法往往受到一定的限制。2.多维寻优方法的选择在直接法优化方法中,坐标轮换法最简易。但是坐标轮换法的效能,很大程度上取决于目标函数的性质。

第95页:另一种较为简单的方法是模式搜索法

另一种较为简单的方法是模式搜索法。模式搜索法的应用范围很广,对变量的极值问题分析是较有效的,程序也较方便,算法收敛速度同步长选择有较大的关系。步长加速法在寻优开始阶段应用,可获得较快的逼近速度,但在后期搜索中的收敛速度不是最理想。Powell方法则是目前多变量寻优直接法中较好的一种方法。

第96页:一维优化方法的选择二次插值法优点

3.一维优化方法的选择二次插值法优点:比较简单,在最优点附近收敛速度很快缺点:要求初始知道高─低─高三点成功失败法优点:在最优点所在区间的寻找上是有效的缺点:最后的收敛速度不是太高相结合先用成功失败法寻找高─低─高三点,然后用二次插值法找出最优解,可使一维寻优快速可靠。

第97页:约束条件的处理对不同的约束类型可

4. 约束条件的处理对不同的约束类型可以用不同的处理方法,通常对不等式约束用内点法构造惩罚项,而对等式约束用外点法构造惩罚项。对于一般同时有等式与不等式约束的优化问题,可以用混合罚函数法,其惩罚函数具体形式为∶(5-102) 式中,gi(x)为不等式约束,hj(x)为等式约束,r为罚因子,是一个递减的无穷正数数列。尽管混合罚函数法是一种比较成熟的方法,但在实际使用中仍有一些需要注意的地方。在式(5-102)中,必须保证为正,否则,不等式惩罚项所起的作用正好远离最优点,因此在每一维的寻优中都必须检验不等式约束是否满足要求。

第98页:优化目标对于冰箱

5.1.5.4优化设计实例下面以冰箱为例,对优化过程加以进一步的说明。1.优化目标对于冰箱,在性能可靠的前提下,要求制造成本低,使用费用即耗电量低。在设计时主要是尽可能降低耗电量。冰箱工作过程可分为初始打冷工况和常规开停工况,装置的绝大多数时间工作于开停工况(图5-13),选择此工况的耗电量最小为优化目标比较合理。

第99页:制冷装置工作过程

图5-13 制冷装置工作过程

第100页:从理论上讲

从理论上讲,当环境条件不变、系统工作完全稳定时,每一个周期的工作过程都应该相等。实际状况有些偏差,数值仿真是以一定的步长进行的,每个周期都有些差异,因此不宜仅以一个周期的平均功耗最小作为最后的优化目标,而适当多取几个周期。写成数学表达式为(5-103) 一般来说n取3或4就够了。

第101页:优化参数对家用冰箱进行优化计算

2. 优化参数对家用冰箱进行优化计算,可选择以下四个可连续变化参数作为优化参数1) 系统充注量2) 冷凝管的长度3) 毛细管的管长4)冷藏室蒸发器的传热面积,或当肋化系数一定时的流道长度。

第102页:约束条件在冰箱优化计算中选择的几

3.约束条件在冰箱优化计算中选择的几个主要约束条件为:

  • ①毛细管的长度应大于最小布置长度。
  • ②冷藏室蒸发器应该小于最大可布置的面积。
  • ③冷凝器的传热面积应小于最大可能布置面积。
  • ④冷冻室空气温度应该达到国标要求。

第103页:优化方法这是一个约束优化问题

4.优化方法这是一个约束优化问题。需要把上面这些约束条件分别处理。把约束条件④这类非结构参数的约束条件通过修改仿真部分的程序,使其作用在仿真程序中体现出来。这样在优化部分的约束中,都是清一色的结构参数,可以用相近的方法处理,带来许多方便之处。

第104页:由于上面的几个约束条件均为不等式约束

由于上面的几个约束条件均为不等式约束,所以可以取消优化程序中罚函数循环收敛这一层次,借用无约束优化的计算方法来解决此类有约束的优化问题,只要在一维寻优过程中检验不等式约束条件是否满足,这样可使计算时间可大大减少。多维无约束优化采用POWELL方法。一维优化采用成功失败法寻找高低高三点,再用二次插值法找出最优解。

第105页:优化设计步骤

图5-14 优化设计步骤

第106页:计算机辅助设计的基本概念计算机辅

5.1.5.5制冷装置计算机辅助设计入门1.计算机辅助设计的基本概念计算机辅助设计(Computer Aided Design)–CAD技术是近年来得到迅速发展的科技新领域。一个CAD系统一般应该包括有专业计算、分析、优化程序,数据库系统,以及自动化绘图系统。

第107页:计算机辅助设计系统的组成及基本功

2.计算机辅助设计系统的组成及基本功能一个完整的计算机辅助设计系统是由一系列硬件系统和软件系统组成的。作为一个计算机辅助设计系统应包括以下几个功能;

  • (1)计算功能。
  • (2)存储功能。
  • (3)输入功能。
  • (4)输出功能。

第108页:制冷装置计算机辅助设计的内容

3.制冷装置计算机辅助设计的内容一个完整的制冷装置计算机辅助设计系统应该包括从初步规划到最后图纸输出的这样一个功能强大的系统, 大致可以分为

  • (1)结构规划。
  • (2)系统初步分析计算。
  • (3)仿真与优化。
  • (4)自动图纸绘制。

第109页:部分仿真软件介绍

5.1.6部分仿真软件介绍1)空调器仿真软件2)电冰箱仿真软件3)冷水机组仿真软件

第110页:房间空调器仿真软件实现以下的功能模拟

房间空调器仿真软件实现以下的功能

  • 模拟房间空调器(包括窗式空调器和壁挂式空调器)在制冷和制热运行模式下的整机的变工况性能,工况范围覆盖常见的制冷和空调工况。
  • 预测空调器所需的合理的制冷剂充注量。
  • 模拟房间空调器所用的部件特性。包括单独适用于压缩机、冷凝器、蒸发器和毛细管四大主要部件的仿真子系统。可以从部件库中推荐合适的部件用于系统的匹配。

第111页:仿真软件的主界面

仿真软件的主界面仿真软件的主界面空调器仿真软件的主要参数输入界面

第112页:仿真结果一:列表形式

仿真结果一:列表形式

第113页:仿真软件的主界面仿真结果二:压焓图形式

仿真软件的主界面仿真结果二:压焓图形式

第114页:仿真结果三:流程图形式

仿真结果三:流程图形式

第115页:敏感性分析

敏感性分析

第116页:蒸发器部件的仿真界面

蒸发器部件的仿真界面

第117页:毛细管部件的仿真界面

毛细管部件的仿真界面

第118页:电冰箱仿真软件实现以下

电冰箱仿真软件实现以下的功能

  • 模拟电冰箱在国标规定的六种实验工况下的动态过程和性能指标。
  • 模拟电冰箱在自定义实验工况下的动态过程和性能指标。
  • 预测电冰箱所需的合理的制冷剂充注量。
  • 模拟电冰箱压缩机、毛细管、冷凝器、蒸发器、箱体等部件的动态特性。

第119页:电冰箱仿真软件的启动画面

电冰箱仿真软件的启动画面

第120页:仿真软件的主输入界面

仿真软件的主输入界面

第121页:冰箱结构参数显示

冰箱结构参数显示

第122页:箱体结构参数的输入界面

箱体结构参数的输入界面

第123页:空气温度变化曲线

空气温度变化曲线

第124页:冷量与功率变化曲线

冷量与功率变化曲线

第125页:蒸发与冷凝温度变化曲线

蒸发与冷凝温度变化曲线

第126页:制冷剂流量变化曲线

制冷剂流量变化曲线

第127页:冷却速度的仿真演示画面

冷却速度的仿真演示画面

第128页:负载温度回升的仿真演示画面

负载温度回升的仿真演示画面

第129页:冷冻能力的仿真演示画面

冷冻能力的仿真演示画面

第130页:制冰能力的仿真演示画面

制冰能力的仿真演示画面

第131页:耗电量实验的仿真演示画面

耗电量实验的仿真演示画面

第132页:冷水机组仿真软件实现以下的功能模拟

冷水机组仿真软件实现以下的功能

  • 模拟冷水机组的变工况稳态性能指标(包括制冷量、输入功率、冷媒水和冷却水的出口温度、蒸发温度、冷凝温度等)。
  • 模拟冷水机组的变结构稳态性能指标。
  • 模拟冷水机组的开机动态特性。
  • 模拟冷水机组的停机动态特性。
  • 模拟冷水机组的变负荷动态特性。

第133页:冷水机组仿真软件的启动画面

冷水机组仿真软件的启动画面

第134页:仿真的类型选择界面

仿真的类型选择界面

第135页:仿真软件的主界面

仿真软件的主界面

第136页:参数输入界面

参数输入界面

第137页:稳态性能敏感性分析

稳态性能敏感性分析(图形)

第138页:稳态性能敏感性分析

稳态性能敏感性分析(表格)

第139页:制冷剂温度变化曲线

制冷剂温度变化曲线

第140页:冷却水和冷媒水温度变化曲线

冷却水和冷媒水温度变化曲线

第141页:制冷量和输入功率变化曲线

制冷量和输入功率变化曲线

第142页:动态COP曲线

动态COP曲线

第143页:制冷剂质量分布变化曲线

制冷剂质量分布变化曲线

第144页:制冷剂流量变化曲线

制冷剂流量变化曲线

第145页:仿真曲线的局部放大分析

仿真曲线的局部放大分析

顶部目录