1.一种内嵌管式围护结构供冷房间热过程仿真模拟方法，其特征在于，该方法包括以下步骤：
步骤一、内嵌管式围护结构供冷房间模型的建立及设定：
步骤101、设定供冷房间由四面墙体和顶部墙体组成，其中，供冷房间朝向南的墙体为南墙体，南墙体为外墙，且南墙体上设置一个外窗户，供冷房间的其余三个墙体均为内墙，内嵌管式供冷管埋设在顶部墙体内，且顶部墙体称作辐射冷板；
步骤102、设定三个内墙的壁面温度和表面发射率相同，则三个内墙等效为一个内墙；
步骤二、构建墙体的RC有限差分模型:
步骤201、构建外墙的RC网络有限差分模型：
步骤2011、采用数据处理器将房间中外墙的总厚度记作 将外墙的导热系数记作λw，
3
单位为W/(m·K)，将外墙的密度记作ρw，单位为kg/m ，将外墙的比热容记作cw，单位为J/(kg·K)，采用数据处理器根据公式 得到外墙的热阻Rw，单位为m2·K/W,采用数据处理器根据公式 得到外墙的热容Ci，单位为J/m2·K；
步骤2012、采用数据处理器将外墙沿厚度方向分割成N层外墙；其中，N层外墙的厚度均相同，N层外墙按照从内到外依次标记为第1层外墙，第2层外墙，...，第j层外墙，...，第N层外墙；其中，j和N均为正整数，且1≤j≤N，N≥3；
步骤2013、采用数据处理器根据公式 得到外墙中第j层外墙的热阻Rw,j；采用数据处理器根据公式 得到外墙中第j层外墙的热容Cw,j；
步骤2014、采用数据处理器对第j层外墙构建T型中心热容RC网络模型结构，得到第j层外墙构的T型中心热容RC网络模型结构，第j层外墙的T型中心热容RC网络模型结构包括第j层外墙的第一电阻Rw,j,1、第j层外墙的第二电阻Rw,j,2和第j层外墙的热容Cw,j，且第j层外墙的第一电阻Rw,j,1的一端、第j层外墙的第二电阻Rw,j,2的一端和第j层外墙的热容Cw,j的一端连接；其中，
步骤2015、按照步骤2014所述的方法，得到第j-1层外墙构的T型中心热容RC网络模型结构和第j+1层外墙构的T型中心热容RC网络模型结构；其中，第j层外墙的第一电阻Rw,j,1的另一端与第j-1层外墙的第二电阻Rw,j-1,2的另一端连接，第j层外墙的第二电阻Rw,j,2的另一端与第j+1层外墙的第一电阻Rw,j+1,1的另一端连接；
步骤2016、按照步骤2011至步骤2015所述的方法，得到内墙的RC有限差分模型；其中，第i层内墙的T型中心热容RC网络模型结构包括第i层内墙的第一电阻Rn,i,1、第i层内墙的第二电阻Rn,i,2和第i层内墙的热容Cn,i，且第i层内墙的第一电阻Rn,i,1的一端、第i层内墙的第二电阻Rn,i,2的一端和第i层内墙的热容Cn,i的一端连接，内墙划分为Nn层，Nn和i均为正整数，1≤i≤Nn，Nn≥3； 第1层内墙的第一电阻Rn,1,1的另一端连接墙体内壁温度，Rn,i表示外墙中第i层内墙的热阻；
步骤三、构建供冷构件的简化RC网络模型：
步骤301、设定内嵌管式供冷管内水温相同和水温分布均匀；其中，内嵌管式供冷管为冷盘管；将内嵌管式供冷管上部的顶部墙体记作上顶部墙体，将内嵌管式供冷管下部的顶部墙体记作下顶部墙体，并按照步骤2012的方法将上顶部墙体划分为Ns层，下顶部墙体划分为Nx层，Ns和Nx均为正整数；
步骤302、采用数据处理器根据公式 得到内嵌管式供冷管的核心层热阻RL；其中，L表示内嵌管式供冷管的埋管间距，单位为m，λup表示内嵌管式供冷管所在的顶部墙体的导热系数，单位为W/(m·K)；
步骤303、采用数据处理器根据公式 得到内嵌管式供冷管的核心
层热容CL；其中，ρup表示内嵌管式供冷管所在的顶部墙体的密度，单位为kg/m3，cup表示内嵌管式供冷管所在的顶部墙体的比热容，单位为J/(kg·K)，中间变量
x1表示上顶部墙体的厚度，单位为m，x2表示下顶部墙体的厚度，单位为m，比例系数d2表示内嵌管式供冷管的外径，单位为m；
步骤304、采用数据处理器根据公式 得到第j′层上顶部墙体的热阻
Rup,j′；采用数据处理器根据公式 得到第j′层上顶部墙体的热容Cup,j′；
采用数据处理器根据公式 得到第j″层下顶部墙体的热阻Rx,j″；采用数据处理器根据公式 得到第j″层下顶部墙体的热容Cx,j″；其中，j′和j″均为正整数，且1≤j′≤Ns，1≤j″≤Nx；
步骤305、按照步骤2014至步骤2015所述的方法，得到上顶部墙体的T型中心热容RC网络模型和下顶部墙体的T型中心热容RC网络模型；其中，上顶部墙体中第j′层上顶部墙体的T型中心热容RC网络模型结构包括第j′层上顶部墙体的第一上电阻Rup,j′,1、第j′层上顶部墙体的第二上电阻Rup,j′,2和第j′层上顶部墙体的热容Cup,j′，且第j′层上顶部墙体的第一上电阻Rup,j′,1的一端、第j′层上顶部墙体的第二上电阻Rup,j′,2的一端和第j′层上顶部墙体的热容Cup,j′的一端连接，
下顶部墙体中第j″层下顶部墙体的T型中心热容RC网络模型结构包括第j″层下顶部墙体的第一下电阻Rx,j″,1、第j″层下顶部墙体的第二下电阻Rx,j″,2和第j″层下顶部墙体的热容Cx,j″，且第j″层下顶部墙体的第一下电阻Rx,j″,1的一端、第j″层下顶部墙体的第二下电阻Rx,j″,2的一端和第j″层下顶部墙体的热容Cx,j″的一端连接，
步骤306、当内嵌管式供冷管的管内流体的状态为层流时，采用数据处理器根据公式得到内嵌管式供冷管的管内流体与管壁面间的对流换热系数h，单位为W/
(m2·K)，其中，d表示内嵌管式供冷管的内径，单位为m，λsf表示内嵌管式供冷管的管内流体的导热系数；
当内嵌管式供冷管的管内流体的状态为过渡流时，采用数据处理器根据公式
得到内嵌管式供冷管的管内流体与管
壁面间的对流换热系数h；其中，Re表示雷诺数，Pr表示内嵌管式供冷管的管内流体温度下的普朗特数，Prb表示内嵌管式供冷管管壁温度下的普朗特数；
当内嵌管式供冷管的管内流体的状态为紊流时，采用数据处理器根据公式
得到内嵌管式供冷管的管内流体与管壁面间的对流换热系数h；
步骤307、采用数据处理器根据公式 得到内嵌管式供冷管的管内流体与管壁面间的对流换热热阻Rw；
步骤四、构建房间热力系统RC网络模型：
步骤401、设定内墙的辐射换热热阻为Rr，顶部墙体与外墙之间的辐射换热热阻为Rr,1，顶部墙体与内墙之间的辐射换热热阻为Rr,2，外墙与内墙之间的辐射换热热阻为Rr,3，顶部墙体的辐射换热热阻为Rr,4，且Rr＝Rr,1＝Rr,2＝Rr,3＝Rr,4，设定内嵌管式供冷管表面对流换热热阻Rc,c-a、外墙内表面与室内空气间对流换热热阻Rc,ew-a、内墙内表面与室内空气之间的对流换热热阻Rc,iw-a，1、顶部墙体表面与室内空气间对流换热热阻Rc,iw-a，2、外墙外表面与室外空气间的对流换热热阻Rc,ew-oa和设定内墙另一个内表面与室内空气间对流换热热阻Rc,iw-a；
步骤402、采用数据处理器根据公式Q2＝3.6×Cxf×ρxf×G×(T1a-Ts)，得到新风承担的显热负荷Q2，单位为W；其中，Cxf表示空气定压比热容，单位为kJ/(kg·K)，ρxf表示新风的密度，单位为kg/m3，G表示新风送风量，单位为m3/h，T1a表示室内空气温度，单位为K，Ts表示新风送风温度，单位为K；新风送风量G不小于36.3m3/h；
采用数据处理器根据公式Q3＝Q3，1+Q3，2+Q3，3，得到供冷房间内部得热量Q3，单位为W；其中，Q3，1表示供冷房间内部人员作用下得热量，单位为W；Q3，2表示供冷房间内部灯光作用下得热量，单位为W；Q3，3表示供冷房间内部设备作用下得热量，单位为W；
设定通过外窗户进入房间的日射得热按照面积比例平均分配到房间外墙、内墙和顶部墙体，设定房间日射得热量分配在顶部墙体内表面的单位面积日射得热量为Qdq、分配在内墙内表面的单位面积日射得热量为Qnq、分配在外墙内表面的单位面积日射得热量为Qwq，且Qdq＝Qnq＝Qwq，单位为W/m2，外墙外表面吸收的单位面积太阳辐射得热量为Q4，单位为W/m2；
窗户传热得热量为Q5，单位为W；
步骤五、构建供冷房间综合RC简化模型：
步骤501、设定内嵌管式供冷管的核心层第一热阻为RL，1，内嵌管式供冷管的核心层第二热阻为RL，2，且 内嵌管式供冷管的核心层第一热容为CL，1，内嵌
管式供冷管的核心层第二热容为CL，2，且
步骤502、采用数据处理器将内嵌管式供冷管的管内流体与管壁面间的对流换热热阻Rw的一端、内嵌管式供冷管的核心层第一热阻RL，1的一端和内嵌管式供冷管的核心层第一热容CL，1的一端连接，形成节点1；将内嵌管式供冷管的核心层第热阻RL，1的另一端、内嵌管式供冷管的核心层第二热阻RL，2的一端和内嵌管式供冷管的核心层第二热容CL，2的一端连接，形成节点2；将内嵌管式供冷管的核心层第二热阻RL，2的另一端、第1层上顶部墙体的第一上电阻Rup,1,1的另一端和第1层下顶部墙体的第一下电阻Rx,1,1的另一端连接，形成节点3；
步骤503、采用数据处理器将第1层上顶部墙体的第一上电阻Rup,1,1的一端、第1层上顶部墙体的第二上电阻Rup,1,2的一端和第1层上顶部墙体的热容Cup,1的一端连接，形成节点4；
步骤504、按照步骤503所述的方法，采用数据处理器将第j′层上顶部墙体的第一上电阻Rup,j′,1的一端、第j′层上顶部墙体的第二上电阻Rup,j′,2的一端和第j′层上顶部墙体的热容Cup,j′的一端连接，形成节点3+j′；
步骤505、多次重复步骤504，直至采用数据处理器将第Ns层上顶部墙体的第一上电阻的一端、第Ns层上顶部墙体的第二上电阻 的一端和第Ns层上顶部墙体的热容的一端连接，形成节点3+Ns；
步骤506、采用数据处理器将第Ns层上顶部墙体的第二上电阻 的另一端、顶部墙体表面与室内空气间对流换热热阻Rc,iw-a，2的一端和顶部墙体的辐射换热热阻Rr,4的一端连接，形成节点4+Ns；
步骤507、采用数据处理器将第1层下顶部墙体的第一下电阻Rx,1,1的一端、第1层上顶部墙体的第二上电阻Rx,1,2的一端和第一层下顶部墙体的热容Cx,1的一端连接，形成节点5+Ns；
步骤508、按照步骤507所述的方法，采用数据处理器将第j″层下顶部墙体的第一下电阻Rx,j″,1的一端、第j″层下顶部墙体的第二下电阻Rx,j″,2的一端和第j″层下顶部墙体的热容Cx,j″的一端连接，形成节点4+Ns+j″；
步骤509、多次重复步骤508，直至采用数据处理器将第Nx层下顶部墙体的第一下电阻的一端、第Nx层下顶部墙体的第二下电阻 的一端和第Nx层下顶部墙体的热容的一端连接，形成节点4+Ns+Nx；
步骤5010、采用数据处理器将第Nx层下顶部墙体的第二下电阻 的另一端、顶部墙体与外墙之间的辐射换热热阻Rr,1的一端、顶部墙体与内墙之间的辐射换热热阻Rr,2的一端和内嵌管式供冷管表面对流换热热阻Rc,c-a的一端连接，形成节点5+Ns+Nx；其中，节点5+Ns+Nx考虑顶部墙体内表面的单位面积日射得热量Qdq；
步骤5012、采用数据处理器将内嵌管式供冷管表面对流换热热阻Rc,c-a的另一端、内墙内表面与室内空气之间的对流换热热阻Rc,iw-a，1的一端和外墙内表面与室内空气间对流换热热阻Rc,ew-a的一端连接，形成节点6+Ns+Nx；其中，节点6+Ns+Nx中考虑新风承担的显热负荷Q2、供冷房间内部得热量Q3和窗户传热得热量Q5；
步骤5013、采用数据处理器将顶部墙体与外墙之间的辐射换热热阻Rr,1的另一端、外墙内表面与室内空气间对流换热热阻Rc,ew-a的另一端、外墙与内墙之间的辐射换热热阻为Rr,3的一端和第1层外墙的第一电阻Rw,1,1的另一端连接，形成节点7+Ns+Nx；其中，节点7+Ns+Nx考虑外墙内表面的单位面积日射得热量Qwq；
步骤5014、按照步骤503所述的方法，采用数据处理器将第1层外墙的第一上电阻Rw,1,1的一端、第1层外墙的第二上电阻Rw,1,2的一端和第1层外墙的热容Cw,1的一端连接，形成节点
8+Ns+Nx；
步骤5015、按照步骤5014所述的方法，采用数据处理器将第j层外墙的第一电阻Rw,j,1、第j层外墙的第二电阻Rw,j,2和第j层外墙的热容Cw,j，且第j层外墙的第一电阻Rw,j,1的一端、第j层外墙的第二电阻Rw,j,2的一端和第j层外墙的热容Cw,j的一端连接，形成节点7+Ns+Nx+j；
步骤5016、多次重复步骤5015，直至采用数据处理器将第N层外墙的第一电阻Rw,N,1的一端、第N层外墙的第二电阻Rw,N,2的一端和第N层外墙的热容Cw,N的一端连接，形成节点7+Ns+Nx+N；
步骤5017、采用数据处理器将第N层外墙的第二电阻Rw,N,2的另一端和外墙外表面与室外空气间的对流换热热阻Rc,ew-oa的一端连接，形成节点8+Ns+Nx+N；其中，节点8+Ns+Nx+N考虑外墙外表面吸收的单位面积太阳辐射得热量Q4，外墙外表面与室外空气间的对流换热热阻Rc,ew-oa的另一端连接室外空气温度To；
步骤5018、采用数据处理器将顶部墙体与内墙之间的辐射换热热阻Rr,2的另一端、外墙与内墙之间的辐射换热热阻为Rr,3的另一端、内墙内表面与室内空气之间的对流换热热阻Rc,iw-a，1的另一端和第1层内墙的第一电阻Rn,1,1的另一端连接，形成节点9+Ns+Nx+N；其中，节点9+Ns+Nx+N考虑内墙内表面的单位面积日射得热量Qnq；
步骤5015、采用数据处理器将第1层内墙的第一电阻Rn,1,1的一端、第1层内墙的第二电阻Rn,1,2的一端和第1层内墙的热容Cn,1的一端连接形成节点10+Ns+Nx+N；
步骤5016、按照步骤5015所述的方法，采用数据处理器将第i层内墙的第一电阻Rn,i,1的一端、第i层内墙的第二电阻Rn,i,2的一端和第i层内墙的热容Cn,i的一端连接，形成节点9+Ns+Nx+N+i；
步骤5017、多次重复步骤5016，直至采用数据处理器将第Nn层内墙的第一电阻 的一端、第Nn层内墙的第二电阻 的一端和第Nn层内墙的热容 的一端连接，形成节点9+Ns+Nx+N+Nn；
步骤5018、采用数据处理器将第Nn层内墙的第二电阻 的另一端、内墙内表面与室内空气之间的对流换热热阻Rc,iw-a的一端和内墙的辐射换热热阻Rr的一端连接，形成节点
10+Ns+Nx+N+Nn；其中，内墙内表面与室内空气之间的对流换热热阻Rc,iw-a的另一端连接室内空气温度T1a，内墙的辐射换热热阻Rr的另一端连接室内平均辐射温度T1m。
步骤六、依据kirchhoff电流定律对供冷房间综合RC网络简化模型建立数学模型：
步骤601、采用数据处理器建立各个节点的数学模型，具体过程如下：
节点1的数学模型：
其中，T1表示节点1的温度，T2表示节点2的温度，Tw表示内嵌
管式供冷管内流体的温度，τ表示时间；
节点2的数学模型：
其中，T3表示节点3的温度；
节点3的数学模型：
其中，T4表示节点4的温度， 表示节点5+Ns的温度；
节点4～节点3+Ns中节点3+j′的数学模型：
其中，T2+j′表示节点2+j′的温度，T3+j′表示节点3+
j′的温度，T4+j′表示节点4+j′的温度；
节点4+Ns的数学模型：
其中， 表示节点3+Ns的温度， 表示节
点4+Ns的温度；
节点5+Ns～节点4+Ns+Nx中节点4+Ns+j″的数学模型：
当j″＝1时， 其中， 表示节点5+Ns的温度，
表示节点6+Ns的温度；
当j″大于1时， 其中， 表示节
点3+Ns+j″的温度， 表示节点4+Ns+j″的温度， 表示节点5+Ns+j″的温度；
节点5+Ns+Nx的数学模型：
；其
中， 表示节点4+Ns+Nx的温度， 表示节点5+Ns+Nx的温度， 表示节点6+Ns+
Nx的温度， 表示节点7+Ns+Nx的温度， 表示节点9+Ns+Nx+N的温度；
节点6+Ns+Nx的数学模型：
其中，
表示节点7+Ns+Nx的温度，Q2为新风承担的显热负荷，Q3为供冷房间内部得热量，Q5为窗户传热得热量；
节点7+Ns+Nx的数学模型：
其中， 表示节点8+Ns
+Nx的温度，Qwq为外墙单位面积日射的热量；
节点8+Ns+Nx～节点7+Ns+Nx+N中节点7+Ns+Nx+j的数学模型：
其中， 表示节点
6+Ns+Nx+j的温度， 表示节点7+Ns+Nx+j的温度， 表示节点8+Ns+Nx+j的温
度；
节点8+Ns+Nx+N中的数学模型：
其中， 表示节点8+Ns+Nx+N的温
度，Q4为外墙外表面吸收的单位面积太阳辐射得热量；
节点9+Ns+Nx+N中的数学模型：
其中， 表示节点9+Ns+Nx+N的温度， 表示节点10+Ns+Nx+N的温度；
节点10+Ns+Nx+N～节点9+Ns+Nx+N+Nn中节点9+Ns+Nx+N+i的数学模型：
其中，
表示节点8+Ns+Nx+N+i的温度， 表示节点9+Ns+Nx+N+i的温度，
表示节点10+Ns+Nx+N+i的温度；
节点10+Ns+Nx+N+Nn的数学模型：
其中，
表示节点10+Ns+Nx+N+i的温度；
步骤602、采用数据处理器联立各个节点的数学模型求解，得到各个节点的温度、以及室内空气温度T1a和室内平均辐射温度T1m。
2.按照权利要求1所述的一种内嵌管式围护结构供冷房间热过程仿真模拟方法，其特征在于：步骤2012中，当外墙的导热系数λw小于0.5W/(m·K)时，第j层外墙的厚度小于
0.02m；当外墙的导热系数λw大于等于0.5W/(m·K)时，第j层外墙的厚度小于0.03m；
步骤303中内嵌管式供冷管的埋设几何结构满足 比例系数0.06＜k2＜
0.2。
3.按照权利要求1所述的一种内嵌管式围护结构供冷房间热过程仿真模拟方法，其特征在于：步骤306中判断内嵌管式供冷管的管内流体的状态，具体过程如下：
步骤3061、采用数据处理器根据公式 得到雷诺数Re；其中，um表示内嵌管式供冷管内供冷流体的断面平均流速，单位为m/s，ν表示内嵌管式供冷管内流体的运动粘度,单位为m2/s；
步骤3062、采用数据处理器判断，当Re＜2300时，内嵌管式供冷管的管内流体状态为层流；
4
当2300≤Re≤10时，内嵌管式供冷管的管内流体状态为过渡流；
当Re＞104时，内嵌管式供冷管的管内流体状态为紊流。
4.按照权利要求1所述的一种内嵌管式围护结构供冷房间热过程仿真模拟方法，其特征在于：步骤302中供冷房间内部人员作用下得热量Q3，1、供冷房间内部灯光作用下得热量Q3，2和供冷房间内部设备作用下得热量Q3，3的获取具体过程如下：
步骤3021、采用数据处理器根据公式Q3，1＝(qx+qe)×Se×Sz×ηe，得到供冷房间内部人员作用下得热量Q3，1；其中，Se表示供冷房间内单位面积的人员数量，Sz表示供冷房间的面积，单位为m2，qx表示供冷房间内部单个人员的显热，且qx＝61W，qe表示供冷房间内部单个人员的潜热，且qe＝73W，ηe表示供冷房间内部人员的在室率，且ηe的取值范围为10％～
95％；
步骤3022、采用数据处理器根据公式Q3，2＝qd×Sz×ηd，得到供冷房间内部灯光作用下得热量Q3，2；其中，qd表示灯光功率密度，且qd＝11W/m2，ηd表示供冷房间内部照明灯的使用率，且ηd的取值范围为10％～95％；
步骤3023、采用数据处理器根据公式Q3，3＝qs×Sz×ηs，得到供冷房间内部设备作用下得热量Q3，3；其中，qs表示设备功率密度，且qs＝20W/m2，ηz表示供冷房间内部设备的使用率，且ηz的取值范围为10％～95％。