标题

Similar documents
Vol. 22 No. 4 JOURNAL OF HARBIN UNIVERSITY OF SCIENCE AND TECHNOLOGY Aug GPS,,, : km, 2. 51, , ; ; ; ; DOI: 10.

标题

T K mm mm Q345B 600 mm 200 mm 50 mm 600 mm 300 mm 50 mm 2 K ~ 0. 3 mm 13 ~ 15 mm Q345B 25

3 : 323 (5cm 200 / m 2 ), ;, 2, [1 ] 1 7 Table 1 The seven typical building thermal walls and their thermal performances [ 1 ]

cm /s c d 1 /40 1 /4 1 / / / /m /Hz /kn / kn m ~

successful and it testified the validity of the designing and construction of the excavation engineering in soft soil. Key words subway tunnel

34 22 f t = f 0 w t + f r t f w θ t = F cos p - ω 0 t - φ 1 2 f r θ t = F cos p - ω 0 t - φ 2 3 p ω 0 F F φ 1 φ 2 t A B s Fig. 1

第 05 期 董房等 : 一种卫星遥测在线状态监测及分析系统的设计 WEB 1 2 总体功能及组成 2.1 总体功能 1 2 3Web 2.2 结构组成 Web WEB WEB 2.3 系统各模块接口关系

2 ( 自 然 科 学 版 ) 第 20 卷 波 ). 这 种 压 缩 波 空 气 必 然 有 一 部 分 要 绕 流 到 车 身 两 端 的 环 状 空 间 中, 形 成 与 列 车 运 行 方 向 相 反 的 空 气 流 动. 在 列 车 尾 部, 会 产 生 低 于 大 气 压 的 空 气 流

[1-3] (Smile) [4] 808 nm (CW) W 1 50% 1 W 1 W Fig.1 Thermal design of semiconductor laser vertical stack ; Ansys 20 bar ; bar 2 25 Fig

Microsoft Word - 专论综述1.doc

Microsoft Word tb 赵宏宇s-高校教改纵横.doc

* CUSUM EWMA PCA TS79 A DOI /j. issn X Incipient Fault Detection in Papermaking Wa

p 3 p 4 p 5 p 6 p 7 p 8 p 9 p 10 p 11 θ 1 θ 2 θ 3 θ 4 θ 5 θ 6 θ 7 θ 8 θ 9 θ d 1 = 0 X c 0 p 1 p 2 X c 0 d pi p j p i p j 0 δ 90

2 3. 1,,,.,., CAD,,,. : 1) :, 1,,. ; 2) :,, ; 3) :,; 4) : Fig. 1 Flowchart of generation and application of 3D2digital2building 2 :.. 3 : 1) :,

f 2 f 2 f q 1 q 1 q 1 q 2 q 1 q n 2 f 2 f 2 f H = q 2 q 1 q 2 q 2 q 2 q n f 2 f 2 f q n q 1 q n q 2 q n q n H R n n n Hessian

θ 1 = φ n -n 2 2 n AR n φ i = 0 1 = a t - θ θ m a t-m 3 3 m MA m 1. 2 ρ k = R k /R 0 5 Akaike ρ k 1 AIC = n ln δ 2

g 100mv /g 0. 5 ~ 5kHz 1 YSV8116 DASP 1 N 2. 2 [ M] { x } + [ C] { x } + [ K]{ x } = { f t } 1 M C K 3 M C K f t x t 1 [ H( ω )] = - ω 2

85% NCEP CFS 10 CFS CFS BP BP BP ~ 15 d CFS BP r - 1 r CFS 2. 1 CFS 10% 50% 3 d CFS Cli

11 : 1345,,. Feuillebois [6]. Richard Mochel [7]. Tabakova [8],.,..,. Hindmarsh [9],,,,,. Wang [10],, (80 µm),.,. Isao [11]. Ismail Salinas [12],. Kaw

31 17 www. watergasheat. com km 2 17 km 15 km hm % mm Fig. 1 Technical route of p

标题

11 25 stable state. These conclusions were basically consistent with the analysis results of the multi - stage landslide in loess area with the Monte

Microsoft Word - 刘 慧 板.doc


增 刊 谢 小 林, 等. 上 海 中 心 裙 房 深 大 基 坑 逆 作 开 挖 设 计 及 实 践 745 类 型, 水 位 埋 深 一 般 为 地 表 下.0~.7 m 场 地 地 表 以 下 27 m 处 分 布 7 层 砂 性 土, 为 第 一 承 压 含 水 层 ; 9 层 砂 性 土

Fig. 1 1 The sketch for forced lead shear damper mm 45 mm 4 mm 200 mm 25 mm 2 mm mm Table 2 The energy dissip

1, : FLUENT 81 CA4113Z,,,,.,,. [ ] : :, 9 9t (<) + div (u<) = div ( < grad<) + S < (1) t, u, <, <, < S <,, 1., u i i, E,. k - RNG, () (3) [ 3-4 ],. 5

Revit Revit Revit BIM BIM 7-9 3D 1 BIM BIM 6 Revit 0 4D 1 2 Revit Revit 2. 1 Revit Revit Revit Revit 2 2 Autodesk Revit Aut

United Nations ~ ~ % 2010

2011_1_核红.indd

~ 4 mm h 8 60 min 1 10 min N min 8. 7% min 2 9 Tab. 1 1 Test result of modified

1. 课 程 负 责 人 情 况 姓 名 蒋 效 宇 性 别 男 出 生 年 月 基 本 信 息 最 终 学 历 研 究 生 职 称 副 教 授 电 话 学 位 博 士 职 务 无 传 真 研 究 方 向 MIS 系 统 整 合 电 子



Mechanical Science and Technology for Aerospace Engineering October Vol No. 10 Web SaaS B /S Web2. 0 Web2. 0 TP315 A

H 2 SO ml ml 1. 0 ml C 4. 0 ml - 30 min 490 nm 0 ~ 100 μg /ml Zhao = VρN 100% 1 m V ml ρ g

图 物理模型

18-陈亚莉.FIT)

12-1b T Q235B ML15 Ca OH Table 1 Chemical composition of specimens % C Si Mn S P Cr Ni Fe

为 止, 以 集 中 式 光 伏 发 电 系 统 为 主, 其 主 要 原 因 是 我 国 政 策 推 动 方 面 以 国 家 主 导 为 主, 这 种 自 上 而 下 的 政 策 和 运 行 方 式, 更 容 易 迅 速 推 动 集 中 式 光 伏 系 统 的 建 设 集 中 式 光 伏 发 电

y 1 = 槡 P 1 1h T 1 1f 1 s 1 + 槡 P 1 2g T 1 2 interference 2f 2 s y 2 = 槡 P 2 2h T 2 2f 2 s 2 + 槡 P 2 1g T 2 1 interference 1f 1 s + n n

m m m ~ mm

Microsoft Word - 专论综述1.doc

untitled

电力信息化2013年第1期.indb

Microsoft Word - A _ doc

Building Technology Experience Center concept air conditioning concept heat pump special energy-saving techniques in hydraulics Concrete core conditio

JOURNAL OF EARTHQUAKE ENGINEERING AND ENGINEERING VIBRATION Vol. 31 No. 6 Dec

m m 5m ~ PTFE Q45B 1 670MPa PMSAP MIDAS SAP2000 MIDAS MIDAS s 0.

Fig. 1 Frame calculation model 1 mm Table 1 Joints displacement mm

国有大型能源企业财务风险内部控制研究

; 3/2, Buck-Boost, 3 Buck-Boost DC-DC ; Y, Fig. 1 1 BBMC The topology of three phase-three phase BBMC 3 BBMC (Study on the control strategy of

~ 10 2 P Y i t = my i t W Y i t 1000 PY i t Y t i W Y i t t i m Y i t t i 15 ~ 49 1 Y Y Y 15 ~ j j t j t = j P i t i = 15 P n i t n Y

<D2BDC1C6BDA1BFB5CDB6C8DAD7CAB8DFB7E5C2DBCCB3B2CEBBE1C3FBB5A52E786C7378>

标题

doc

Microsoft Word D 孙竹森.doc

ph ph ph Langmuir mg /g Al 2 O 3 ph 7. 0 ~ 9. 0 ph HCO - 3 CO 2-3 PO mg /L 5 p

#4 ~ #5 12 m m m 1. 5 m # m mm m Z4 Z5

3-D 王明新1.doc

[1] Liu Hongwei,2013, Study on Comprehensive Evaluation of Iron and Steel Enterprises Production System s Basic Capacities, International Asia Confere

<4D F736F F D20B8BDBCFE3220BDCCD3FDB2BFD6D8B5E3CAB5D1E9CAD2C4EAB6C8BFBCBACBB1A8B8E6A3A8C4A3B0E5A3A92E646F6378>

% 8. 48% 3 80 Alcalase Novozymes Alcalase 2. 4 L Bacillus licheniformis 2. 4 AU /g 1. 2 Hitachi S-4700 JEOL JEM-1200EX Olympus Bu

Microsoft Word - A _ doc

NANO COMMUNICATION 23 No.3 90 CMOS 94/188 GHz CMOS 94/188 GHz A 94/188 GHz Dual-Band VCO with Gm- Boosted Push-Push Pair in 90nm CMOS 90 CMOS 94

1. 前 言 由 於 石 油 價 格 浮 動, 汽 油 價 格 節 節 高 升 及 二 氧 化 碳 等 廢 棄 大 量 排 放 造 成 全 球 環 境 的 改 變, 因 此 世 界 各 國 都 極 力 提 倡 節 能 減 碳 進 而 掀 起 腳 踏 車 城 市 的 風 潮 因 應 目 前 自 行 車

LaDefense Arch Petronas Towers 2009 CCTV MOMA Newmark Hahn Liu 8 Heredia - Zavoni Barranco 9 Heredia - Zavoni Leyva

通 过 厂 变 带 电, 这 种 设 计 减 少 了 机 组 自 带 厂 用 电 负 荷 能 力, 降 低 了 锅 炉 满 足 FCB 时 最 低 稳 燃 工 况, 同 时 造 成 燃 烧 调 整 量 加 大 本 电 厂 在 FCB 试 验 时, 电 泵 不 联 启, 始 终 保 持 汽 泵 运 行

mm ~

风 冷 式 多 联 机 空 调 系 统 的 季 节 经 济 性 作 用 域 Vol. 37ꎬ No. 1 为 名 义 工 况 下 的 经 济 性 作 用 域 )ꎬ 并 指 出 当 多 联 机 的 连 接 管 长 度 小 于 该 经 济 性 作 用 域 时 ꎬ 多 联 机 系 统 的 性 能 优 于

微软用户

热设计网

第 03 期 刘高军等 : 基于 CNONIX 的 XML 与 EXCEL 相互转换技术研究 XML XML CNONIX XML EXCEL EXCEL EXCEL EXCEL CNONIXEXCEL XML EXCEL CNONIX XML EXCEL CNONIX 1 CNONIX 数据元分析

专 技 能 1. 精 通 Matlab/Simulink 平 台 下 的 海 洋 运 载 器 运 动 控 制 系 统 与 仿 真 建 模 设 计 ; 2. 精 通 51 单 片 机 AVR 单 片 机 Arduino 开 源 板 的 开 发 和 设 计 ; 3. 精 通 基 于 Arduino 板

JOURNAL OF EARTHQUAKE ENGINEERING AND ENGINEERING VIBRATION Vol. 31 No. 5 Oct /35 TU3521 P315.

義 守 大 學 102 學 年 度 第 1 學 期 第 4 次 行 政 會 議 紀 錄 時 間 :103 年 1 月 8 日 ( 星 期 三 ) 下 午 2:00 地 點 : 行 政 大 樓 十 樓 國 際 會 議 廳 主 席 : 蕭 介 夫 校 長 紀 錄 : 楊 育 臻 壹 報 告 事 項 一

Microsoft Word - 11-秦华伟.doc

~ a 3 h NCEP ~ 24 3 ~ ~ 8 9 ~ km m ~ 500 m 500 ~ 800 m 800 ~ m a 200

EL ECTR IC MACH IN ES AND CON TROL Vol113 No11 Jan. 2009,, (, ) :, X 2Y,,,,,,, P ID P ID P ID,, : ; ; ; P ID : TM33 : A : X

/MPa / kg m - 3 /MPa /MPa 2. 1E ~ 56 ANSYS 6 Hz (a) 一阶垂向弯曲 (b) 一阶侧向弯曲 (c) 一阶扭转 (d) 二阶侧向弯曲 (e) 二阶垂向弯曲 (f) 弯扭组合 2 6 Hz

untitled

Microsoft Word doc

SWAN min TITAN Thunder Identification Tracking Analysis SWAN TITAN and Nowcasting 19 TREC Tracking Radar Echo by Correlaction T

5 551 [3-].. [5]. [6]. [7].. API API. 1 [8-9]. [1]. W = W 1) y). x [11-12] D 2 2πR = 2z E + 2R arcsin D δ R z E = πr 1 + πr ) 2 arcsin

Microsoft Word - d27 彭玉柱.doc

Microsoft Word - A doc

山东省招生委员会

Microsoft Word 張嘉玲-_76-83_

京都大学防災研究所年報 第56号(平成24年度)

m 3 m m 84 m m m m m m m

交流活动

VLBI2010 [2] 1 mm EOP VLBI VLBI [3 5] VLBI h [6 11] VLBI VLBI VLBI VLBI VLBI GPS GPS ( ) [12] VLBI 10 m VLBI 65 m [13,14] (referen

ANSYS在高校《材料力学》教学中的应用

u d = R s i d - ωl q i q u q = R s i q + ωl d i d + ωψ 1 u d u q d-q i d i q d q L d L q d q ψ f R s ω i 1 i 5th i th 5 θ 1 θ θ 3 5 5

mm 5 1 Tab 1 Chemical composition of PSB830 finishing rolled rebars % C Si Mn P S V 0 38 ~ 1 50 ~ 0 80 ~ ~


untitled

132 包 装 工 程 2016 年 5 月 网 产 品 生 命 周 期 是 否 有 与 传 统 产 品 生 命 周 期 曲 线 相 关 的 类 似 趋 势 旨 在 抛 砖 引 玉, 引 起 大 家 对 相 关 问 题 的 重 视, 并 为 进 一 步 研 究 处 于 不 同 阶 段 的 互 联 网

240 生 异 性 相 吸 的 异 性 效 应 [6] 虽 然, 心 理 学 基 础 研 [7-8] 究 已 经 证 实 存 在 异 性 相 吸 异 性 相 吸 是 否 存 在 于 名 字 认 知 识 别 尚 无 报 道 本 实 验 选 取 不 同 性 别 的 名 字 作 为 刺 激 材 料, 通

Transcription:

文章编号 :0253-339(2018) 06-0070 - 07 doi:10 3969 / j issn 0253-339 2018 06 070 1 杨怀毅 1 任滔 1 丁国良 2 陈绍楷 2 刘忠民 (1 上海交通大学制冷与低温工程研究所上海 20020ꎻ 2 海信科龙电器股份有限公司佛山 528303) 摘要智能控制技术与互联网技术在房间空调系统应用的基础是对其性能的在线实时监测 ꎮ 本文提出一种房间空调系统性能在线监测的方法 ꎮ 首先将空调器布置温度测点采集到的换热器铜管表面温度和电流电压等电路参数结合 ꎬ 转化为制冷循环中制冷剂的状态参数 ꎻ 再通过压缩机流量半经验计算模型计算制冷剂的质量流量 ꎬ 进而得出室内机中制冷剂侧的换热量 ꎻ 最后结合电路参数求得的整机功率 ꎬ 实现房间空调系统制冷量 制热量 功率 能效的实时监测 ꎮ 结果表明 : 采用此方法计算房间空调器性能参数的最大偏差小于 15% ꎬ 主要误差来源为蒸发压力的计算误差 ꎮ 关键词制冷系统 ꎻ 性能监测 ꎻ 压缩式制冷循环 ꎻ 制冷量中图分类号 :TB61 + 5ꎻ TU831. 3 文献标识码 : A Performance Monitoring Method of Air conditioning System based on Limited Measuring Points Yang Huaiyi 1 Ren Tao 1 Ding Guoliang 1 Chen Shaokai 2 Liu Zhongmin 2 (1. Institute of Refrigeration and Cryogenicsꎬ Shanghai Jiao Tong Universityꎬ Shanghaiꎬ 20020ꎬ Chinaꎻ 2. Hisense Kelon Electrical Holdings Co. ꎬ Ltd. ꎬ Foshanꎬ 528303ꎬ China) Abstract The alication of intelligent control technologies and internet technologies to room air conditioning systems is based on the erformance monitoring of air conditioning systems. In the studyꎬ the method of real time monitoring of the erformance of room air condi tioning system was roosed. Firstlyꎬ temerature measuring oints were arranged in the air conditionerꎬ and the temerature and circuit data were transformed into refrigerant state arameters. Subsequentlyꎬ the mass flow rate of refrigerant was calculated via a semi emirical comressor modelꎬ and the heat exchange of the indoor unit was calculated via the enthaly difference. Finallyꎬ the ower was calculated by using circuit data. Thusꎬ the mass flow rateꎬ caacityꎬ ower and energy efficiency were monitored. The test results indicate that the average deviation of the calculated erformance obtained by the method in the study is less than 15%. The deviation is mainly attributed to the calculation of evaorating ressure that significantly influences the calculation of the mass flow rate. Keywords refrigerating systemꎻ erformance monitoringꎻ comression refrigeration cycleꎻ refrigerating caacity 随着人们对房间空气舒适性能要求的日益提高 ꎬ 传统空调因控制方法陈旧 舒适度差 温度调节速度 慢等原因已不能满足用户需求 [1] ꎬ 而智能空调通过 对智能传感 互联网等技术的开发 ꎬ 使室内空气的智 能调节成为可能 [2 - ] ꎮ 智能空调的传感技术与互联 网技术等应用的基础是对空调的运行数据进行采集 与传输 ꎬ 而数据的采集是基于对制冷系统运行参数的 实时监测 [5-6] ꎮ 根据标准 GB / T 7725200 [7] ꎬ 房间空调器制冷 量和热泵制热量性能的测试方法主要为房间型量热 [8] [9 - 计法和空气焓值法 10] ꎮ 但在用户实际使用过程 中 ꎬ 没有条件采集风量 温湿度数据进行产品性能判 定 ꎮ 只能通过测量制冷剂侧的参数 ꎬ 求出进 出口制 冷剂的焓值和流量 ꎬ 计算系统制冷剂侧的换热量来获取系统的制冷量 制热量 ꎮ 采用制冷剂侧的参数计算性能 ꎬ 需要获取制冷剂的质量流量 ꎮ 目前有研究通过测量载冷剂的流量进而计算出系统制冷量等参数 ꎬ 实现对中央空调系统性能的在线监测 [11] ꎬ 但该方法不适用于风冷式空调系统的监测 ꎮ 对于风冷式空调系统的监测 ꎬ 额外增加测量装置测量制冷剂流量会影响系统性能 ꎮ 目前已有很多学者研究了半经验压缩机流量模 [12 - 型 1] ꎬ 可以通过测量压缩机的吸 排气压力 ꎬ 结合 [15 - 实验数据计算得到制冷剂流量 16] ꎮ 由于压力测点相较于温度测点成本较高 ꎬ 因此可以通过测量换热器铜管温度和进风温度 ꎬ 结合实验数据估算对应的饱和 收稿日期 :2017 年 9 月 6 日 70

压力 ꎮ 本文开发了基于有限测点的空调系统在线监测方法 ꎬ 包括空调系统的测点布置方法及根据测点的数据计算系统制冷剂质量流量 制冷量 压缩机功率的算法 ꎮ 1 技术路线 空调系统的主要性能参数包括制冷量 制热量 整机功率 能效比等 ꎮ 由于室内换热存在析湿现象 ꎬ 测量空气侧的热量变化对传感器的要求较高 ꎬ 而测量制冷剂侧的参数 ꎬ 计算制冷剂质量流量及室内机制冷剂进出口焓差 ꎬ 可以较为准确地得到换热量 : Q R = m r (h rꎬout - h rꎬin ) (1) 式中 :Q R 为室内机换热量 ꎬWꎻm r 为制冷剂质量流量 ꎬkg / sꎻ h rꎬout 为室内机制冷剂出口焓值 ꎬ kj / kgꎻ h rꎬin 为室内机制冷剂入口焓值 ꎬkJ / kgꎮ 综合考虑性能参数的计算精度和测量成本 ꎬ 对家用空调系统性能参数进行检测的思路为 : 在空调系统中布置温度 电流 电压测点 ꎬ 并建立测点数据和实验测得的温度 压力 质量流量 功率因数 室外侧电流等参数转化的模型 ꎬ 计算时依据模型逐步求解得到系统的性能参数 ꎮ 压力 焓的转换 : 根据实验数据和物理特性对测点数据进行拟合 ꎬ 建立测点数据与系统中实际制冷循环参数的对应关系 ꎻ2) 质量流量 制冷量 功率等性能参数的计算 : 建立压缩机的质量流量模型 ꎬ 结合第一步得到的压力 焓 功率因数等数据计算空调系统的制冷量 制热量及能效比 ꎮ 本文认为蒸发压力和冷凝压力分别与吸 排气压力近似相等 ꎮ 2 测点布置 图 2 所示为空调部件与对应测点示意图 ꎬ 系统部件由压缩机 室内机 膨胀阀 室外机组成 ꎬ 实际产品中共布置 5 个温度测点 ꎻ 另外 ꎬ 需要测量室外侧的电流 电压及压缩机的运行频率 ꎮ 图 2 空调部件与对应测点示意图 Fig. 2 Schematic diagram of air conditiming comonents and corresonding measuring oints 实际产品中测点的布置方案主要考虑以下因素 : 1) 测点能反应整个系统的运行情况 ꎬ 可用于求解制 冷剂质量流量 室内机进出口焓差 功率等 ꎻ2) 测点 不影响空调系统的性能及可靠性 ꎬ 减少不必要的成 本 ꎻ3) 换热器温度测点布置在流路中部 ꎬ 使其能测得 两相区的换热器铜管温度 ꎮ 各个测点的测量精度如表 1 所示 ꎮ 表 1 测量精度 Tab. 1 Uncertainties of direct measurements 测点 测点说明 测量精度 T 1 室外机液管温度 ± 1 T 2 室外机进风温度 ± 1 T 3 压缩机排气温度 ± 1 T 室内机液管温度 ± 1 T 5 室内机进风温度 ± 1 U 输入电压 ± 1 V 图 1 系统性能在线监测技术路线 Fig. 1 Road ma for air conditioner erformance monitoring 系统性能参数求解可分为两步 :1) 测点数据向 I 室外侧输入电流 ± 0 2 A F 压缩机运行频率 71

3 循环状态参数的转换 由于温度测点无法直接测量制冷剂的温度 ꎬ 需要对测点的温度值进行转换以得到相应的冷凝温度 蒸发温度等参数 ꎻ 测点的电流值为室外侧的总电流 ꎬ 由于变频空调中存在电感 ꎬ 需要进行转换才能得到整机的有效电流 ꎻ 而节流前后没有布置测点 ꎬ 需要根据实验数据及饱和点的焓值对该点的焓值进行估算 ꎮ 3 1 电流的转换测点的电流值为室外侧电流 ꎬ 加上室内侧的风机电流即为流入空调系统的总电流 ꎬ 可用线性公式表示测点电流与整机电流的关系 ꎮ 由于变频器的存在 ꎬ 总电流不能直接用于计算 ꎬ 应结合功率因数转换为有效电流 ꎮ 当运行频率较高时 ꎬ 功率因数矫正装置 ( PFC) 开启 ꎬ 功率因数 PF (ower factor) 维持在 0 9 以上 ꎻ 当运行频率较低时 ꎬ PFC 关闭 ꎬ 功率因数较低 ꎮ 根据空调系统中 PFC 的开关状态 ꎬ 可分别用关于电流的二次多项式对 PF 进行拟合 ꎮ 总的有效电流可以拟合为以下形式 : I all = I PF = (C 1 I r + C 2 )(C 3 I r 2 + C I r + C 5 ) (2) 式中 :I all 为有效电流 ꎬAꎻI 为室外侧测点的电流 ꎬ AꎻI r 为测点电流 ꎬAꎻPF 为对应的功率因数 ꎮ 式中 C 1 ~ C 18 均为拟合系数 ꎮ 3 2 蒸发 冷凝压力的转换冷凝压力和蒸发压力的转换精度直接影响质量流量的计算精度 ꎬ 制冷工况下的蒸发压力 制热工况下的冷凝压力还影响焓差的计算 ꎮ 出于成本考虑 ꎬ 本文通过测量进气温度与铜管表面温度拟合冷凝温度与蒸发温度 ꎬ 再根据制冷剂物性得出冷凝压力和蒸发压力 ꎮ 布置在换热器流路中部的温度测点 ꎬ 可反映两相段的壁温 ꎬ 通过拟合 ꎬ 可减少计算蒸发和冷凝温度的系统性偏差 ꎮ 如图 3 所示 ꎬ 测点为铜管外表面的温度 ꎬ 铜管与空气 管内制冷剂分别进行换热 ꎮ 图 3 换热器铜管温度测点示意图 Fig. 3 Coil temerature measurement oints diagram 当温度稳定时 ꎬ 忽略铜内部的导热 ꎬ 温度可以表示为 : α 1 A 1 (T 空气 - T 测点 ) + α 2 A 2 (T 测点 - T 制冷剂 ) = 0 (3) 72 即 : T 制冷剂 = α 1A 1 α 2 A 2 T 空气 + (α 2A 2 - α 1 A 1 ) α 2 A 2 T 测点 () 式中 :α 为表面传热系数 ꎬW / (m 2 )ꎻA 为传热 面积 ꎬm 2 ꎮ 表面传热系数的计算较为复杂 ꎬ 这里可采 用换热器铜管表面测点温度及进风测点温度来拟合 制冷剂的冷凝温度 ꎬ 再计算得到该温度对应的饱和 压力 : T C = C 6 T L + C 7 T A + C 8 (5) 式中 :T C 为冷凝温度 ꎬ ꎻT L 为对应的换热器铜管 温度 ( 制冷时为室外机液管温度 T 1 ꎬ 制热时为室内机 液管温度 T )ꎬ ꎻT A 为对应的进风温度 ( 制冷时为室 外机进风温度 T 2 ꎬ 制热时为室内机进风温度 T 5 )ꎬ ꎮ 实践表明 ꎬ 采用换热器铜管温度及进风温度对冷 凝温度进行二维线性拟合 ꎬ 可以得到较好的拟合精度 ( 相关系数 R 2 > 0 99)ꎻ 但蒸发温度的拟合效果不理 想 ꎬR 2 < 0 8ꎮ 考虑到蒸发压力与质量流量的变化关 系很大 ꎬ 因此增加运行频率进行三维线性拟合 ꎬR 2 > 0 99ꎬ 拟合效果较好 : T E = C 9 T L + C 10 T A + C 11 F + C 12 (6) 式中 :T E 为蒸发温度 ꎬ ꎻF 为压缩机的运行频 率 ꎬHzꎻT L 为对应的换热器铜管温度 ( 制冷时为室内 机液管温度 T ꎬ 制热时为室外机液管温度 T 1 )ꎬ ꎻT A 为对应的进风温度 ( 制冷时为室内机进风温度 T 5 ꎬ 制 热时为室内机进风温度 T 2 )ꎬ ꎮ 3 3 吸气温度 节流焓的转换 系统中没有压缩机吸气状态和节流前后状态相 应位置的测点 ꎬ 因此需要采用其他参数进行拟合 ꎮ 压缩机的吸气温度用于计算压缩机的吸气密度 及制冷工况下的室内机出口焓 ꎮ 由于没有吸气温度 测点 ꎬ 本文采用排气过热度拟合吸气过热度 ꎬ 结合蒸 发温度得到吸气温度 ꎮ 当排气过热度小于 10 时 ꎬ 吸气过热度通常小于 2 ꎬ 而排气过热度较大时 ꎬ 吸 气过热度才相应增加 ꎮ 本文采用排气过热度对吸气 过热度进行二次拟合 ꎬ R 2 > 0 95ꎬ 拟合结果如图 所示 ꎮ T S = C 13 (T 3 - T C ) 2 + C 1 (T 3 - T C ) 2 + C 15 + T E (7) 制冷剂经过节流装置时的焓 ꎬ 在制冷工况下为室 内机的进口焓 ꎬ 在制热工况下为室内机的出口焓 ꎮ 由 空调系统的实际运行状态可知 ꎬ 节流前制冷剂应有一 定的过冷度 ꎬ 节流后变为两相态 ꎬ 即节流前后的焓值 应小于冷凝压力对应的汽化焓 ꎬ 应大于蒸发压力对应 的汽化焓 ꎮ 结合实验数据 ꎬ 当运行频率较高时 ꎬ 过冷 度较大 ꎻ 运行频率较低时 ꎬ 过冷度较小 ꎮ

组成 ꎮ 对于定频压缩机 ꎬ 质量流量可用容积效率计 算 ꎬ 而容积效率可由实验数据拟合得到 ꎮ 即 : m = λ V V th ν suc = V th λ V λ ν suc λ Vꎬref Vꎬref 图 吸气过热度拟合结果 Fig. Fitting results of suction gas suerheat 图 5 节流焓示意图 Fig. 5 Refrigerant enthaly diagram at the throttle 节流焓示意图如图 5 所示 ꎬ 节流时的过冷度与冷凝压力 蒸发压力有关 ꎬ 当冷凝压力较高时 ꎬ 过冷度较大 ꎻ 蒸发压力较低时 ( 此时运行频率较高 ꎬ 质量流量较大 )ꎬ 过冷度较大 ꎮ 另外 ꎬ 随着风机的档位不同 ꎬ 略有变化 ꎮ 本文采用运行频率 蒸发压力下干度为 0 的焓 冷凝压力下干度为 0 的焓对过冷焓进行拟合 : h sc = C 16 h min + C 17 h max + C 18 F (8) 式中 :h min 为蒸发压力下干度为 0 的焓值 ꎬkJ / kgꎻ h max 为冷凝压力下干度为 0 的焓值 ꎬkJ / kgꎮ 流量模型与性能参数计算 1 在线监测压缩机流量模型 已有很多学者对半经验压缩机质量流量模型进行了研究 ꎬ 本文在前人研究的基础上 ꎬ 针对空调在线监测进行了改动 ꎮ 在线监测中 ꎬ 蒸发压力与冷凝压力不是直接测量得到 ꎬ 存在一定误差 ꎬ 因而对模型拟合实验数据的精度要求更高 ꎬ 以减少输入量变化对总误差的影响 ꎮ 对于变频压缩机 ꎬ 当工作在特定频率时 ꎬ 其性能与工作在该频率下的定频压缩机相似 ꎬ 因此可将变频压缩机视为无数个工作在不同频率下的定频压缩机 = V th ν suc f(f)f( C ꎬ E ) (9) 式中 :λ V 为变频压缩机的输气系数 ꎻλ Vꎬref 为压缩 机在定速下的参考输气系数 ꎻV th 为压缩机的理论容 积输气量 ꎬm 3 / sꎻν suc 为开式压缩机环节吸气口的制 冷剂气体比容 ꎬm 3 / kgꎻ C 为冷凝压力 ꎬkPaꎻ E 为蒸 发压力 ꎬkPaꎮ 质量流量可由吸气温度 蒸发压力计 算得到 ꎮ 根据容积效率比与定频部分质量流量的 拟合形式中是否包含压力项一次 二次项和二次频 率项 ꎬ 对以下 6 种流量拟合形式进行对比 ꎬ 得到相 应拟合系数 ꎬ 并对比公式的拟合偏差 ꎬ 结果如表 2 所示 ꎮ f(f)f( C ꎬ E ) = ì F [c 1 F 2 æ C ö ç 3 ] (10 1) F [c 1 F 2 æ C ö ç 3 C 5 E 6 ] (10 2) 2 F [c 1 F 2 C ö 3 æ C ö ç è E 5 C + c 6 2 C 7 E 8 2 E 9 ] (10 3) í [c 1( F F ref ) 2 F 2 F 3 C ö 5 ] (10 ) [c 1( F F ref ) 2 F 2 F 3 C ö + c 5 C 6 E 7 ] (10 5) [c 1( F F ref ) 2 2 F 2 F 3 C ö C ö 5 + è E îc 6 C 7 2 C 8 E 9 2 E 10 ] (10 6) 式中 :c 1 ~ c 10 为拟合系数 ꎮ 表 2 不同流量模型的拟合偏差对比 Tab. 2 Fitting errors for each mass flow rate equation 公式平均偏差 / % 最大偏差 / % 式 (10 1) 8 9 29 5 式 (10 2) 6 6 23 3 式 (10 3) 3 8 12 0 式 (10 ) 6 3 17 7 式 (10 5) 2 5 5 2 式 (10 6) 1 3 0 73

由表 2 可知 ꎬ 含拟合系数最多的式 (10 6) 平均偏差和最大偏差均最小 ꎬ 分别为 1 % 和 3 0% ꎬ 拟合效果最好 ꎮ 考虑到在线监测中 ꎬ 蒸发压力和冷凝压力一般存 在一定偏差 ꎬ 需要质量流量模型在蒸发压力和冷凝压力略有变动时尽量保持稳定 ꎬ 表 3 对比了输入的蒸发压力和冷凝压力在 100 kpa 范围内变化时 ꎬ 式 (10) 对总流量的平均拟合偏差 ꎮ 表 3 压力偏移时不同流量模型的平均拟合偏差对比 Tab. 3 Average fitting errors for each mass flow rate equation with ressure offsets 蒸发压力 偏移 / kpa 冷凝压力 偏移 / kpa 平均拟合偏差 / % 式 (10 1) 式 (10 2) 式 (10 3) 式 (10 ) 式 (10 5) 式 (10 6) 0 0 8 9 6 6 3 8 6 3 2 5 1 50 0 6 2 5 3 2 7 2 3 2 3-50 0 6 1 3 3 0 5 0 2 6 2 1 100 0 6 3 5 0 3 7 9 3 2 5 1-100 0 6 1 5 3 8 5 5 3 8 3 8 0 100 6 1 3 9 2 8 5 1 2 6 2 3 100 100 6 1 5 5 3 9 6 0 5 5 3 由表 3 可知 ꎬ 式 (10 5) 和式 (10 6) 蒸发压力在 ± 50 kpa 范围内变动时 ꎬ 流量计算结果变化较小 ꎻ 式 (10 3) 由于不含有频率的二次项 ꎬ 蒸发压力在 ± 100 kpa 范围内变动时 ꎬ 流量计算结果变化较小 ꎻ 冷凝压力的变化对流量的影响较小 ꎮ 实际应用中 ꎬ 可根据蒸发压力与冷凝压力的转换精度 ꎬ 采用不同的拟合公式 : 当蒸发压力转换精度较高时 ꎬ 选用含有频率 压力一次 二次项的式 (10 5) 或式 (10 6)ꎻ 当蒸发压力及冷凝压力的转换精度较低时 ꎬ 选用不含频率二次项的式 (10 3)ꎮ 2 室内换热量及能效的计算依据测点数据转换得到的室内机进出口焓值及压缩机流量公式计算得到的质量流量 ꎬ 可由式 (1) 求得室内侧的换热量 ꎮ 功率可由测点数据转换得到的有效电流及电压计算得到 ꎮ 实验中发现测点的电压与外接电压计测得的电压存在 0 ~ 6 V 的压降 ꎬ 计算时可以加上一定的压降修正提高功率的计算准确性 ꎮ 图 6 换热量计算值与实验值对比 Fig. 6 Comarison between calculation and exeriment results of indoor unit heat exchange 5 实验验证及误差分析 5 1 实验验证根据海信科龙空调公司提供的空调在 16 个制冷工况及 11 个制热工况下的测点数据 ꎬ 按本文的性能参数算法 ꎬ 得到的性能参数与实验数据的偏差如图 6 和图 7 所示 ꎮ 由图 6 和图 7 可知 ꎬ 本文开发的各个系统性能参数的算法与实验实测值的最大偏差小于 15% ꎬ 满足在线监测的精度要求 ꎮ 7 图 7 功率计算值与实验值对比 Fig. 7 Comarison between calculation and exeriment results of ower consumtion 5 2 误差分析为降低成本 ꎬ 空调系统中的温度测点的精度有限

( ± 1 )ꎬ 在转换压力 焓等参数时也存在一定偏差 ꎮ 根据企业提供的实验数据 ꎬ 各中间变量与测量值的偏差如表 所示 ꎮ 表 循环状态参数的计算偏差 Tab. Deviations of the calculated refrigerant state arameters 制冷工况制热工况参数平均偏差最大偏差平均偏差最大偏差蒸发压力 28 kpa - 87 kpa 9 kpa 18. 5 kpa 冷凝压力 20 kpa - 5 kpa 6 kpa 188 kpa 体积流量 1. 6% - 5. 2% 1. 5% 7. 0 % 吸气密度 2. 7% 8. 3% 1. 0% 2. 5% 节流焓 0. 9% 1. 8% 1. 3% - 2. 7% 室内出口焓 0. 8% - 1. 7% 室内入口焓 0. 7% - 1. 9% 注 : 由于蒸发 冷凝压力计算偏差与制冷量不存在直接的对应关系 ꎬ 因此采用差值表示 ꎮ 由表 可知 ꎬ 换热量的误差主要来源于质量流量计算 ꎬ 而流量的误差主要有两方面 :1) 温度测点的误差 ꎬ 产生蒸发 冷凝压力的计算误差 ꎬ 从而在计算体积流量时产生较大误差 ꎻ2) 蒸发压力误差与拟合吸气过热度的误差 ꎬ 导致计算吸气密度时产生误差 ꎮ 因此蒸发压力是本文换热量监测方法的主要误差来源 ꎬ 蒸发压力对体积流量 吸气密度和制冷工况下的室内机出口焓均有影响 ꎮ 增加蒸发压力测点可以有效提高换热量预测精度 ꎬ 但同时也会显著增加在线监测的成本 ꎮ 功率的计算误差主要来自电流测量的误差与功率因数的拟合误差 ꎬ 改善电流测量的精度 ꎬ 可有效减少这部分误差 ꎮ 6 结论 1) 对于房间空调 ꎬ 实际使用过程中可以通过采集温度值及电路参数 ꎬ 结合实验数据 ꎬ 拟合制冷剂状态和其他电路参数的方法 ꎬ 求得制冷剂侧的室内换热量及功率 ꎬ 对其性能进行监测 ꎮ 2) 本文开发的性能参数算法与实验测量值的最大偏差均小于 15% ꎬ 可较为准确的监测系统性能 ꎮ 3) 换热量计算的误差 ꎬ 主要来源于蒸发压力的计算误差导致的体积流量和吸气密度的计算误差 ꎻ 功率计算的误差 ꎬ 主要来源于有效电流的计算误差 ꎮ 参考文献 [1] 曹恩国 ꎬ 张宇希 ꎬ 张宇红 ꎬ 等. 基于立体送风模式的智能空调创新设计研究 [ J]. 机械设计 ꎬ 2015ꎬ 32 (5): 105-108. ( CAO Enguoꎬ ZHANG Yuxiꎬ ZHANG Yu hongꎬ et al. Innovative design of intelligent air conditioner based on stereoscoic air suly technology[ J]. Journal of Machine Designꎬ 2015ꎬ 32(5): 105-108. ) [2] 郭旭红 ꎬ 芮延年 ꎬ 李军涛. 基于遗传算法的模糊智能变频空调仿真系统的研究 [ J]. 系统仿真学报 ꎬ 2005ꎬ 17 (5):1237-120. (GUO Xuhongꎬ RUI Yannianꎬ LI Jun tao. Research on the simulation system of fuzzy intelligence inverter air conditioner based on genetic algorithm [ J ]. Journal of System Simulationꎬ 2005ꎬ 17 ( 5 ): 1237-120. ) [3] HAVES Pꎬ SALSBURY T Iꎬ WRIGHT J A. Condition mo nitoring in HVAC subsystems using first rinciles models [ J]. ASHRAE Transactionsꎬ 1996ꎬ 102(1):519-527. [] BUSWELL R A. Uncertainty in the first rincile model based condition monitoring of HVAC systems[ J]. Lough borough:loughborough Universityꎬ 2001. [5] PI HAI LINꎬ BROBERG H L. Internet based monitoring and controls for HVAC alications[ J]. IEEE Industry A lications Magazineꎬ 2002ꎬ 8(1):9-5. [6] TAYLOR Kꎬ WARD Jꎬ GERASIMOV Vꎬ et al. Sensor / actuator networks suorting agents for distributed energy management[ C] / / IEEE International Conference on Local Comuter Networks. Tama: IEEE Comuter Societyꎬ 200:63-70. [7] 房间空气调节器 :GB / T 7725200[ S]. 北京 : 中国标准出版社 ꎬ200. ( Room air conditioners: GB / T 7725 200[ S]. Beijing: Standards Press of Chinaꎬ200. ) [8] 柳胜耀 ꎬ 李瑛 ꎬ 黄彩凤 ꎬ 等. 平衡环境型房间量热计凝结水量测量设计 [ J]. ꎬ 201ꎬ 35 (2):106-110. ( LIU Shengyaoꎬ LI Yingꎬ HUANG Caifengꎬ et al. The design of condensate measurement for balanced ambient room tye calorimeter[ J]. ꎬ 201ꎬ 35(2):106-110. ) [9] 徐志亮. 空气焓差法的试验分析 [ J]. 发电与空调 ꎬ 2005ꎬ 26(6):6-67. ( XU Zhiliangꎬ Exerimental study on air enthaly difference method [ J]. Refrigeration Air Conditioning & Electric Power Machineryꎬ 2005ꎬ 26 (6): 6-67. ) [10] 曹小林 ꎬ 李雄林 ꎬ 喻首贤 ꎬ 等. 空调器性能测试平台测量不确定度的研究 [ J]. ꎬ 2009ꎬ 30(5):58-62. ( CAO Xiaolinꎬ LI Xionglinꎬ YU Shouxianꎬ et al. In vestigation on measurement uncertainty in air conditioner erformance testing bench [ J ]. ꎬ 2009ꎬ 30(5):58-62. ) 75

[11] 罗林聪 ꎬ 王刚 ꎬ 李明 ꎬ 等. 浅谈能源计量系统用于中央空调能效在线检测 [ J]. 能源与节能 ꎬ 2016 (7):7-75. ( LUO Lincongꎬ WANG Gangꎬ LI Mingꎬ et al. Intro duction to the energy measurement system alication in en ergy efficiency insection on line to central air conditioning [ J]. Energy and Conservationꎬ 2016(7):7-75. ) [12] SHAO Shuangquanꎬ SHI Wenxingꎬ LI Xiantingꎬ et al. Performance reresentation of variable seed comressor for inverter air conditioners based on exerimental data [ J]. International ꎬ 200ꎬ 27(8):805-815. [13] LI Wenhua. Simlified steady state modeling for variable seed comressor [ J ]. Alied Thermal Engineeringꎬ 2013ꎬ 50(1):318-326. [1] 孙浩然 ꎬ 任滔 ꎬ 丁国良 ꎬ 等. 一种产品数据交互式的变频压缩机理论模型 [J]. ꎬ 2015ꎬ 36(3):73-78. (SUN Haoranꎬ REN Taoꎬ DING Guoliangꎬ et al. A roduct data in teractive theoretical model of variable seed comressor[ J]. ꎬ 2015ꎬ 36(3):73-78. ) [15] 丁国良 ꎬ 张春路. 制冷空调装置仿真与优化 [ M]. 北京 : 科学出版社 ꎬ 2001. ( DING Guoliangꎬ ZHANG Chunlu. Simulation and otimization of refrigeration and air condi tioning units[ M]. Beijing:Science Pressꎬ2001. ) [16] 张春路. 制冷空调系统仿真原理与技术 [ M]. 北京 : 化学工业出版社 ꎬ2013. ( ZHANG Chunlu. Fundamentals of vaor comression refrigeration and air conditioning system modeling[ M]. Beijing: Chemical Industry Pressꎬ2013. ) 通信作者简介丁国良 ꎬ 男 ꎬ 教授 ꎬ 博士生导师 ꎬ 上海交通大学机械与动力工程学院 ꎬ(021) 3206378ꎬE mail: glding @ sjtu. edu. cnꎮ 研究方向 : 制冷空调装置的仿真 优化与新工质应用 ꎮ About the corresonding author Ding Guoliangꎬ maleꎬ rofessorꎬ Ph. D. suervisorꎬ School of Mechanical Engineeringꎬ Shanghai Jiao Tong Universityꎬ + 86 21 3206378ꎬ E mail: glding @ sjtu. edu. cn. Research fields: simulation and otimization research for room air conditioner and utilization of new refrigerant. 76