基于MatlabSimulink环境下空调建筑物动态仿真模型的建立.docx

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