在MatlabSimulink环境下建立空调建筑物的动态仿真模型.doc

基于Matlab/Simulink环境下空调建筑物动态仿真模型的建立 哈尔滨工业大学 陈建成 姜永成 王森 唐勇辉 摘要 分析了空调建筑物的特性，介绍了运用状态空间法建立空调建筑物数字模型的过程，以某大学空调系统试验台的空调房间作为建模对象，利用Matlab软件中的Simulink工具箱建立了空调建筑物的动态仿真模型，并对模型进行了动态特性仿真，通过改变不同干扰通道的扰量，分析得出了房间的不同干扰通道的过程特性参数：放大系数K和时间常数T等，最后总结了该模型的应用范围。 关键词 Matlab/Simulink 空调建筑物 动态仿真模型 0 引言 空调建筑物是一个复杂的热力模型，受室内外等因素的影响，具有多干扰性。因此，建立精确的数学模型十分困难，近年来，许多学者对此做了研究，文献[1]在建立VAV空调房间数学模型之前进行了过多的简化(如不考虑室外温度的变化、室内负荷、太阳辐射等的影响)，所建的数学模型与实际的房间特性相差太大；文献[2]、[3]、[4]均用一个惯性环节和纯滞后环节串联来表示房间模型，这种方法不能充分反映空调建筑物的动态特性。随着我国经济的高速发展和节能意识的提高，建立建筑物精确的数学模型和动态仿真模型，对于研究建筑物围护结构的特性及其暖通空调系统的特性，确定合理的控制系统具有非常重要的意义。 1 空调建筑物数字模型的建立 空调房间的示意图如图1所示，房间模型建立的基本方法：状态空间法。整个空调房间是一个复杂的热力系统，想用完全精确的数学模型来描述几乎不可能，因此，在建立数学模型前先作如下少量简化：1)墙体、房间及各独立环节各看成一个单容对象；2）忽略房间内部各物体的蓄热量；3)忽略房间内部气体的流动；4)设房间内部温度为分布参数。  假设组成建筑物的围护结构的各层材料各向均质同性，则墙体的导热可用偏微分方程描述[5]： 图1 空调房间的示意图（1） 其中， 为材料的密度，；为材料的比热， ；为材料的导热系数，；t为壁体内的温度，℃。 该方程加上墙体两侧的边界条件和整个墙体的初始条件就构成一个定解问题，采用空间上离散而时间上保持连续的半离散化技术，处理该方程得到一个常微分方程组，由于建筑物有一系列的围护结构，可以得到一系列常微分方程组。各围护结构表面之间的长波辐射以及各房间之间的空气对流渗透，使上述方程组相互关联成一个矩阵微分方程，描述了整个建筑物的热平衡。 单一房间由围护结构、室内家具和室内空气构成，围护结构包括外墙、屋顶、内墙、楼板、楼地、门和窗，其中窗户由于具有热容较小、可能内含空气层等特点，因而与其他围护结构略有不同，因此在围护结构热平衡部分，窗户和其他围护结构的温度节点离散方法和热平衡方程不同。 图2 多层材料组成的单一围护结构离散温度节点示意图 图3 窗户的温度节点示意图 1.1 单一外墙、内墙、楼板、楼地的热平衡方程 对于由多层材料组成的单一围护结构，如图2所示，内外表面分别是1节点和n+1节点。由于室外温度采用室外空气综合温度，所以相当于n+1节点的温度已知，则内表面温度节点、内部温度节点和n温度节点的热平衡方程如下[5]： （2） （3） （4） 以上三式中：、、、为第i个差分层的比热容，KJ/(kg·K)、密度kg/m3、导热系数，W/(m·K)、厚度，m；为房间内的空气温度，℃；为室外空气综合温度，℃；为围护结构内表面的换热系数，W/(m2·K)；为第j面墙内表面温度，℃；为与第j个表面的辐射换热系数，W/(m2·K)；为围护结构内表面获得的辐射热量，包括太阳辐射热量以及照明、人体和设备的辐射得热，W。 1.2 窗户的热平衡方程 建立双层玻璃窗户的内、外层玻璃温度节点的热平衡方程如下[5]： （5） （6）以上两式中，为窗户内层和外层玻璃的热容，KJ/(m2·K)；，分别是窗户内层和外层玻璃的温度；与分别为各层玻璃对太阳辐射的吸收得热，W；为内层温度节点获得室内照明、人体和设备的辐射得热，W。 1.3 房间空气的热平衡方程 这里所谓的房间空气热平衡，包括室内家具、设备，因为他们同样影响着房间整体的热特性。房间空气热平衡方程式用文字表示为： 与各壁面的对流换热量+各种对流得热量+空气渗透得热量+空调系统显热除热量=单位时间内房间空气中显热量的增值[6]，其数学式表示为： （7）式中：为照明、人体显热和设备显热等的对流散热量，W， ；为空气渗透量，m3/h；为