CN110953697A - 基于Modelica的空调风系统初平衡调试方法 - Google Patents
基于Modelica的空调风系统初平衡调试方法 Download PDFInfo
- Publication number
- CN110953697A CN110953697A CN201911125356.7A CN201911125356A CN110953697A CN 110953697 A CN110953697 A CN 110953697A CN 201911125356 A CN201911125356 A CN 201911125356A CN 110953697 A CN110953697 A CN 110953697A
- Authority
- CN
- China
- Prior art keywords
- air
- model
- manual regulating
- flow
- resistance
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F24—HEATING; RANGES; VENTILATING
- F24F—AIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
- F24F11/00—Control or safety arrangements
- F24F11/62—Control or safety arrangements characterised by the type of control or by internal processing, e.g. using fuzzy logic, adaptive control or estimation of values
Abstract
本发明公开了基于Modelica的空调风系统初平衡调试方法,首先使用Modelica语言建立风道系统的数学模型,以手动调节阀的阻力系数作为调节变量输入模型;计算各风口出风量仿真值与设计值的最小绝对误差(LAD)作为优化目标,选取优化算法不断优化,直至误差达到最小值;最后计算优化后仿真结果风量与设计风量的偏差值,并与规范要求进行比较。本发明提出的基于Modelica的空调风系统初平衡调试方法,可以为实际调试工作提供指导,减少实际工程中投入的人力和时间。
Description
技术领域
本发明涉及建筑空调系统调试方法,尤其涉及一种基于Modelica的新建大型公共建筑空调风系统初平衡调试方法。
背景技术
建筑行业中的调试(Commissioning),属于北美建筑行业成熟的管理和技术体系,是一种使得建筑各个系统在方案设计、图纸设计、安装、单机试运转、性能测试、运行和维护的整个过程中确保能够实现设计意图和满足业主使用要求的工作程序和方法。
如图1所示,一套完整的空调风系统,通常由风机、风管、分流三通、弯头连接而成,风机为整个空调风系统的动力源,提供热湿处理过后的空气及克服阻力损失所需的压头,经过热湿处理过后的空气经过风管和分流三通、弯头等管件后(弯头的作用在于改变管道的走向,分流三通的作用在于从主风管中引出一个支管进行流量的分配。在一个风系统中所有的直管段都可以称为风管,通过弯头、三通等管件将风管连接起来),通过风口送入室内,每个风口前都会设置一个手动调节阀,其调试工作的关键在于风平衡的调试,通过调节每个风口前的手动调节阀的阻力系数,保证各风口的出风量满足设计要求。目前,最常用的风平衡调试方法为比例调节法,这是一种迭代方法,需要不断调节手动调节阀以达到调试目标,且调节其中任意一个手动调节阀的阻力系数都会对整个风系统的水力特性产生影响。对于这种方法,调试的质量很大程度上取决于现场操作人员的工程经验,并且需要花费大量的时间。
发明内容
本发明的目的在于克服已有技术的缺点,提供基于Modelica的空调风系统初平衡调试方法,为新建大型公共建筑空调风系统初平衡调试提供指导。
本发明是这样实现的,基于Modelica的空调风系统初平衡调试方法,包括如下步骤:
基于Modelica的空调风系统初平衡调试方法,包括以下步骤:
步骤一、选择位于各个风口处的手动调节阀的阻力系数形成一组手动调节阀阻力系数作为输入空调风系统模型中手动调节阀模型的初始值,使用Modelica语言依照实际工程风系统现场布置情况及尺寸信息建立由风管模型、弯头模型、分流三通模型、手动调节阀模型和风机模型组成一个整体的空调风系统模型;
其中风机模型采用风机的流量-压头曲线表示;
风管模型和弯头模型的计算公式如下:
式中:ΔP1为风管阻力损失,Pa;S为风管阻力系数;Q1为经过风管的空气流量,m3/h,由风机提供,在模型仿真中计算得到;ΔP2为弯头阻力损失,Pa;ζ1为弯头阻力系数;ρ为空气密度,kg/m3,为已知量;v为空气流经弯头的流速,m/s;Q2为经过弯头的空气流量,m3/h,由风机提供,在模型仿真中计算得到;A1为弯头截面积,m2;
分流三通分为直通管段和旁通管段两部分,其数学模型均通过以下公式进行计算:
ζ2=a1R5+a2R4+a3R3+a4R2+a5R+a6
式中:vup为分流三通上游流速,m/s;Qup为分流三通上游流量,m3/h;Aup为分流三通上游截面积,m2;vdown为分流三通下游流速,m/s;Qdown为分流三通上游流量,m3/h;Adown为分流三通下游截面积,m2;R为上下游流速之比;ζ2为分流三通阻力系数;ΔP3为分流三通阻力损失,Pa;ρ为空气密度,kg/m3;[a1…a6]为常系数;
手动调节阀数学模型具体计算公式如下:
式中:ΔP4为手动调节阀阻力损失,Pa;ζ3为手动调节阀阻力系数;ρ为空气密度,kg/m3;v为空气流经手动调节阀的流速,m/s;b1,b2,b3为常系数;Q3为空气流量,m3/h;A3为手动调节阀截面积,m2;
步骤二、将所述的作为初始值的一组手动调节阀阻力系数输入空调风系统模型的手动调节阀数学模型,使用Dymola作为Modelica语言的编译和计算平台对整个空调风系统模型进行仿真,仿真结果输出为每个风口的出风量,并计算仿真结果与风口出风量设计值的最小绝对误差;
步骤三、若最小绝对误差已达到最小值,则进行下一步骤,否则使用粒子群优化算法选择一组新的手动调节阀的阻力系数重复步骤二;
步骤四、计算仿真结果与风口出风量设计值的偏差率,若偏差率满足±10%国家规范要求,则将该组手动调节阀的阻力系数作为最终的输出结果以辅助现场调试,若不满足国家规范要求,使用粒子群优化算法选择一组新的手动调节阀阻力系数,重复步骤二至步骤四直至偏差率满足国家规范要求。
本发明具有以下优点和积极效果:
本发明提出的基于Modelica的空调风系统初平衡调试方法,可以为实际调试工作提供指导,可根据仿真结果得到的手动调节阀阻力系数对现场的阀门进行调整,减少实际工程中投入的人力和时间,并实现良好的管道初平衡效果。
附图说明
图1是简单空调风系统示意图;
图2是本发明基于Modelica的空调风系统初平衡调试方法流程图。
具体实施方式
下面结合实例和附图对本发明做进一步说明。本发明的实例是为了更好地让本领域的技术人员理解本发明,并不对本发明做任何的限制。
如图2所示,本发明的基于Modelica的空调风系统初平衡调试方法,包括以下步骤:
步骤一、选择位于各个风口处的手动调节阀的阻力系数形成一组手动调节阀阻力系数作为输入空调风系统模型中手动调节阀模型的初始值,使用Modelica语言依照实际工程风系统现场布置情况及尺寸信息建立由风管模型、弯头模型、分流三通模型、手动调节阀模型和风机模型组成一个整体的空调风系统模型;
其中风机模型采用风机的流量-压头曲线表示,在实际运行的某一时刻或某一时间段内,风机只能按其曲线上的某一点工作,该工作状态点的压头与空调风系统的阻力损失相等,流量为风机接入空调风系统的流量。风机的流量-压头曲线风机以供应商提供的工频转速下的风机流量-压头曲线作为基准,利用相似原理计算得到风机不同转速下的流量-压头曲线;
风管模型和弯头模型均为具有固定阻力系数的流动阻力模型,风管阻力系数由风管尺寸决定和弯头阻力系数由弯头尺寸决定,风管模型和弯头模型的具体计算公式如下:
式中:ΔP1为风管阻力损失,Pa;S为风管阻力系数,由风管尺寸决定,为已知量;Q1为经过风管的空气流量,m3/h,由风机提供,在模型仿真中计算得到;ΔP2为弯头阻力损失,Pa;ζ1为弯头阻力系数,由弯头尺寸决定,为已知量;ρ为空气密度,kg/m3,为已知量;v为空气流经弯头的流速,m/s;Q2为经过弯头的空气流量,m3/h,由风机提供,在模型仿真中计算得到;A1为弯头截面积,㎡,由弯头尺寸决定,为已知量。
分流三通分为直通管段和旁通管段两部分,其数学模型均通过以下公式进行计算:
ζ2=a1R5+a2R4+a3R3+a4R2+a5R+a6
式中:vup为分流三通上游流速,m/s;Qup为分流三通上游流量,由风机提供,在模型仿真中计算得到,m3/h;Aup为分流三通上游截面积,㎡,由其尺寸决定,为已知量;vdown为分流三通下游流速,m/s;Qdown为分流三通上游流量,由风机提供,在模型仿真中计算得到,m3/h;Adown为分流三通下游截面积,㎡,由其尺寸决定,为已知量;R为上下游流速之比;ζ2为分流三通阻力系数;ΔP3为分流三通阻力损失,Pa;ρ为空气密度,kg/m3,为已知量;[a1…a6]为常系数,由曲线拟合得到,为已知量。
手动调节阀数学模型具体计算公式如下:
式中:ΔP4为手动调节阀阻力损失,Pa;ζ3为手动调节阀阻力系数;ρ为空气密度,kg/m3,为已知量;v为空气流经手动调节阀的流速,m/s;b1,b2,b3为常系数,由曲线拟合得到,为已知量;Q3为空气流量,m3/h,由风机提供,在模型仿真中计算得到;A3为手动调节阀截面积,㎡,由手动调节阀尺寸决定,为已知量。
步骤二、将所述的作为初始值的一组手动调节阀阻力系数输入空调风系统模型的手动调节阀数学模型,使用Dymola作为Modelica语言的编译和计算平台对整个空调风系统模型进行仿真,仿真结果输出为每个风口的出风量,并计算仿真结果与风口出风量设计值的最小绝对误差(LAD);
步骤三、若LAD已达到最小值,则进行下一步骤,否则使用粒子群优化算法选择一组新的手动调节阀的阻力系数重复步骤二。
步骤四、计算仿真结果与风口出风量设计值的偏差率,若偏差率满足国家规范《通风与空调工程施工质量验收规范》中,偏差率±10%的要求,则将该组手动调节阀的阻力系数组合作为最终的输出结果以辅助现场调试,若无法满足偏差率±10%的要求,使用粒子群优化算法选择一组新的阻力系数,重复步骤二至步骤四直至偏差率满足国家规范要求。
具体的计算结果如下表1所示。
阀门编号 | 设计风量m<sup>3</sup>/h | 仿真后风量m<sup>3</sup>/h | 偏差率 | 阀门开度 |
dam | 1000 | 964.3 | -3.57% | 65% |
dam1 | 1000 | 1045.3 | 4.53% | 70% |
dam2 | 1000 | 1032.7 | 3.27% | 70% |
dam3 | 600 | 616.7 | 2.78% | 80% |
dam4 | 1000 | 964 | -3.6% | 75% |
dam5 | 1000 | 1041.6 | 4.16% | 75% |
dam6 | 600 | 587.7 | -2.04% | 80% |
dam7 | 1000 | 1023.7 | 2.37% | 85% |
dam8 | 1000 | 1018.4 | 1.84% | 100% |
需要声明的是,本发明内容及具体实施方式意在证明本发明所提供技术方案的实际应用,不应解释为对本发明保护范围的限定。
本领域技术人员在本发明的精神和原理启发下,可作各种修改、等同替换、或改进。但这些变更或修改均在申请待批的保护范围之内。
Claims (1)
1.基于Modelica的空调风系统初平衡调试方法,其特征在于包括以下步骤:
步骤一、选择位于各个风口处的手动调节阀的阻力系数形成一组手动调节阀阻力系数作为输入空调风系统模型中手动调节阀模型的初始值,使用Modelica语言依照实际工程风系统现场布置情况及尺寸信息建立由风管模型、弯头模型、分流三通模型、手动调节阀模型和风机模型组成一个整体的空调风系统模型;
其中风机模型采用风机的流量-压头曲线表示;
风管模型和弯头模型的计算公式如下:
式中:ΔP1为风管阻力损失,Pa;S为风管阻力系数;Q1为经过风管的空气流量,m3/h,由风机提供,在模型仿真中计算得到;ΔP2为弯头阻力损失,Pa;ζ1为弯头阻力系数;ρ为空气密度,kg/m3,为已知量;v为空气流经弯头的流速,m/s;Q2为经过弯头的空气流量,m3/h,由风机提供,在模型仿真中计算得到;A1为弯头截面积,m2;
分流三通分为直通管段和旁通管段两部分,其数学模型均通过以下公式进行计算:
ζ2=a1R5+a2R4+a3R3+a4R2+a5R+a6
式中:vup为分流三通上游流速,m/s;Qup为分流三通上游流量,m3/h;Aup为分流三通上游截面积,m2;vdown为分流三通下游流速,m/s;Qdown为分流三通上游流量,m3/h;Adown为分流三通下游截面积,m2;R为上下游流速之比;ζ2为分流三通阻力系数;ΔP3为分流三通阻力损失,Pa;ρ为空气密度,kg/m3;[a1…a6]为常系数;
手动调节阀数学模型具体计算公式如下:
式中:ΔP4为手动调节阀阻力损失,Pa;ζ3为手动调节阀阻力系数;ρ为空气密度,kg/m3;v为空气流经手动调节阀的流速,m/s;b1,b2,b3为常系数;Q3为空气流量,m3/h;A3为手动调节阀截面积,m2;
步骤二、将所述的作为初始值的一组手动调节阀阻力系数输入空调风系统模型的手动调节阀数学模型,使用Dymola作为Modelica语言的编译和计算平台对整个空调风系统模型进行仿真,仿真结果输出为每个风口的出风量,并计算仿真结果与风口出风量设计值的最小绝对误差;
步骤三、若最小绝对误差已达到最小值,则进行下一步骤,否则使用粒子群优化算法选择一组新的手动调节阀的阻力系数重复步骤二;
步骤四、计算仿真结果与风口出风量设计值的偏差率,若偏差率满足±10%国家规范要求,则将该组手动调节阀的阻力系数作为最终的输出结果以辅助现场调试,若不满足国家规范要求,使用粒子群优化算法选择一组新的手动调节阀阻力系数,重复步骤二至步骤四直至偏差率满足国家规范要求。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911125356.7A CN110953697A (zh) | 2019-11-18 | 2019-11-18 | 基于Modelica的空调风系统初平衡调试方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911125356.7A CN110953697A (zh) | 2019-11-18 | 2019-11-18 | 基于Modelica的空调风系统初平衡调试方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110953697A true CN110953697A (zh) | 2020-04-03 |
Family
ID=69977676
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911125356.7A Pending CN110953697A (zh) | 2019-11-18 | 2019-11-18 | 基于Modelica的空调风系统初平衡调试方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110953697A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112032860A (zh) * | 2020-08-03 | 2020-12-04 | 四川省建筑科学研究院有限公司 | 一种集中空调系统风力平衡自动调适方法及系统 |
CN114017910A (zh) * | 2021-10-27 | 2022-02-08 | 中国核动力研究设计院 | 一种通风系统风量调试方法、调试系统及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103729336A (zh) * | 2013-12-09 | 2014-04-16 | 中科院成都信息技术股份有限公司 | 一种dae系统的指标约简方法 |
WO2014113751A1 (en) * | 2013-01-21 | 2014-07-24 | Siemens Corporation | Functional debugging of equation-based languages |
CN108151246A (zh) * | 2017-12-13 | 2018-06-12 | 长安大学 | 变风量空调风系统优化节能控制方法及装置 |
CN110107989A (zh) * | 2019-04-30 | 2019-08-09 | 北京工业大学 | 基于冷冻水回水温度最佳设定点的小型定频冷水机组变水温控制方法 |
-
2019
- 2019-11-18 CN CN201911125356.7A patent/CN110953697A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014113751A1 (en) * | 2013-01-21 | 2014-07-24 | Siemens Corporation | Functional debugging of equation-based languages |
CN103729336A (zh) * | 2013-12-09 | 2014-04-16 | 中科院成都信息技术股份有限公司 | 一种dae系统的指标约简方法 |
CN108151246A (zh) * | 2017-12-13 | 2018-06-12 | 长安大学 | 变风量空调风系统优化节能控制方法及装置 |
CN110107989A (zh) * | 2019-04-30 | 2019-08-09 | 北京工业大学 | 基于冷冻水回水温度最佳设定点的小型定频冷水机组变水温控制方法 |
Non-Patent Citations (1)
Title |
---|
孙育英: "Dymola空调风平衡调试虚拟仿真实训平台开发", 《实验室研究与探索》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112032860A (zh) * | 2020-08-03 | 2020-12-04 | 四川省建筑科学研究院有限公司 | 一种集中空调系统风力平衡自动调适方法及系统 |
CN114017910A (zh) * | 2021-10-27 | 2022-02-08 | 中国核动力研究设计院 | 一种通风系统风量调试方法、调试系统及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102519703B (zh) | 空气末端风量平衡检测模拟装置及其模拟检测方法 | |
CN110953697A (zh) | 基于Modelica的空调风系统初平衡调试方法 | |
CN109059217B (zh) | 一种基于运行曲线的变风量空调系统总风量控制方法 | |
CN101655012A (zh) | 一种双洞隧道互补式网络通风的方法 | |
EP3534079B1 (de) | Verfahren zum ermitteln der einzelraum-luftvolumenströme bei zentralen lüftungssystemen und zum pneumatischen abgleich von lüftungssystemen | |
CN112128842A (zh) | 一种供热管网水力平衡快速调节方法 | |
CN204946458U (zh) | 一种全尺寸中央空调与集中供热综合实验平台 | |
CN109719107A (zh) | 一种适用于不同工作模式的有机废气收集方法 | |
CN111397035B (zh) | 一种用于纺织厂空调风系统及其工作方法 | |
CN116467823A (zh) | 一种基于cfd数值模拟的长管道均匀送风设计方法 | |
CN108629115A (zh) | 一种除尘管网阻力平衡方法 | |
CN102841508A (zh) | 分流式气浴风道 | |
CN108506622B (zh) | 一种基于弧形导流片的低阻力三通构件 | |
CN203024323U (zh) | 新型通风均匀限流装置 | |
CN114611428B (zh) | 一种复杂矿井通风网络的模拟分析方法 | |
CN207089622U (zh) | 变风量末端装置和空调系统 | |
CN103982973A (zh) | 一种按风量分布的通风管道系统及设计方法 | |
CN113792443B (zh) | 通风系统管网水力平衡的分析方法、管网设计和调节方法 | |
CN107725945B (zh) | 基于耗散率控制的低阻力三通构件 | |
CN206959224U (zh) | 一种基于bim集成化族块的空调机房 | |
CN107676563B (zh) | 基于仿生学植物分支结构的低阻力三通构件 | |
CN117433125A (zh) | 一种洁净空调系统风量平衡和多房间压差计算调节方法 | |
CN105718631A (zh) | 一种基于数值模拟的通风管道烟气净化系统 | |
CN106446342B (zh) | 一种轴流风机叶片安装角度获得方法 | |
CN113639276B (zh) | 一种防窜烟三通烟道及其控制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200403 |
|
RJ01 | Rejection of invention patent application after publication |