CN110196987B - 基于代理模型的风道结构尺寸优化方法 - Google Patents
基于代理模型的风道结构尺寸优化方法 Download PDFInfo
- Publication number
- CN110196987B CN110196987B CN201910260231.9A CN201910260231A CN110196987B CN 110196987 B CN110196987 B CN 110196987B CN 201910260231 A CN201910260231 A CN 201910260231A CN 110196987 B CN110196987 B CN 110196987B
- Authority
- CN
- China
- Prior art keywords
- optimization
- variable
- air duct
- model
- precision
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Ventilation (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开一种基于代理模型的风道结构尺寸优化方法,包括:(10)原始方案建模仿真:根据风道原始结构尺寸,建立风道CFD模型,计算结构尺寸优化目标变量初始值,(20)数学模型建立:选择结构优化变量,建立优化问题数学模型,(30)结构优化变量采样:在结构优化变量的取值范围内进行结构优化变量采样,(40)样本建模仿真:计算根据结构优化变量样本组合的每个方案的结构尺寸优化目标变量取值;(50)构建代理模型:构建多项式响应面代理模型,拟合每个结构尺寸优化目标变量和结构优化变量之间的函数关系;(60)多目标优化:求解优化问题数学模型,获得风道优化结构尺寸。本发明的风道结构尺寸优化方法,工作量小、周期短。
Description
技术领域
本发明属于通风管道结构设计技术领域,特别是一种工作量小、省时省力、产品开发周期短的基于代理模型的风道结构尺寸优化方法。
背景技术
风道是空调系统用于送风的重要部件,室内温度的变化和出口的风量呈正相关,各 个位置的出风量要尽可能保持一致。实际工作时,由于风道比较狭长,在风道截面形状不作改变的前提下,气流在风道中向后流动的过程中速度衰减较大,靠近进风口和远离 进风口位置的各个出口处气体流速的差异极大。除此以外,风道内气流流动也不可避免 地会产生气动噪声,这又在一定程度上降低了舒适性。噪声会分散注意力,影响人的心 情,降低工作效率。其中,流速是影响噪声的一个重要因素,在风速高的情况下,气动 噪声非常明显。这是目前风道结构设计中存在的两个最主要问题,需要对风道结构进行 优化设计。
经过对现有技术的检索发现,中国专利文献号CN101976276A,公开日为2011年 02月16日,公开了一种空调器风道结构优化设计方法,包括步骤:利用建模软件建立 至少2种不同结构的空调器风道CFD模型;通过CFD软件计算出噪声以及单位时间内 风道模型出口吹出的风量;比对出各种风道模型中风量大且噪声小的模型,并将其作为 预选方案;对所述预选方案做至少一次风道模型修改,再次重复上述步骤选择更优方案; 确定最终风道模型,进行噪声实验和风量实验。但该技术对设计参数采用反复调整择优 的方法,工作量极大且不能准确描述设计参数和性能响应之间的关系,导致风道结构的 最优化设计难以实现。
总之,现有技术存在的问题是:空调风道结构优化设计费时费力、产品开发周期长。
发明内容
本发明的目的在于提供一种基于代理模型的风道结构尺寸优化方法,工作量小、省 时省力、产品开发周期短。
实现本发明目的的技术方案为:
一种基于代理模型的风道结构尺寸优化方法,其特征在于,包括以下步骤:
(10)原始方案建模仿真:根据风道原始结构尺寸,建立风道CFD模型,在最大 风量输出工况下计算得到结构尺寸优化目标变量初始值;
(20)数学模型建立:根据风道原始结构形式,选择入口、出口、弯管内侧倒角、 弯管外侧倒角中至少两处结构尺寸作为结构优化变量,并根据风道整体结构的尺寸限制 确定其变化范围,以结构尺寸优化目标变量最小为优化目标,建立优化问题数学模型;
(30)结构优化变量采样:利用最优拉丁超立方设计方法,在结构优化变量的取值范围内进行结构优化变量采样,结构优化变量样本个数根据代理模型所需的最小值初步拟定;
(40)样本建模仿真:根据结构优化变量样本组合方案分别建立新的CFD模型, 并分别在最大风量输出工况下计算得到每个方案的结构尺寸优化目标变量取值;
(50)构建代理模型:构建多项式响应面代理模型,拟合每个结构尺寸优化目标变量和结构优化变量之间的函数关系,并验证代理模型的精度,精度不满足时,增加多项 式阶数,仍不满足,则增加样本个数;
(60)多目标优化:采用带有精英保留策略的快速非支配多目标优化算法,求解优化问题数学模型,获得使结构尺寸优化目标变量最小的结构优化变量取值,即风道优化 结构尺寸。
与现有技术相比,本发明的显著优点在于:
工作量小、省时省力、产品开发周期短:本发明通过构建多项式响应面代理模型的方式,得到了风道多个结构尺寸优化目标变量与结构优化变量间的数学函数关系,充分 考虑了结构优化变量之间的耦合关系,以及结构尺寸优化目标变量间的关联特性;最优 拉丁超立方设计方法使得样本个数可人为控制且提升了样本点分布的均匀性,有效减少 了仿真分析的次数;采用带有精英保留策略的快速非支配多目标优化算法,实现了风道 结构优化变量的快速寻优。
附图说明
图1为本发明基于代理模型的风道结构尺寸优化方法的主流程图。
图2为实施例风道的三维结构示意图。
图3为实施例风道的二维结构示意图。
图中,1是出口,2是弯管内侧倒角,3是入口。
具体实施方式
如图1所示,本发明基于代理模型的风道结构尺寸优化方法,包括以下步骤:
(10)原始方案建模仿真:根据风道原始结构尺寸,建立风道CFD模型,在最大 风量输出工况下计算得到结构尺寸优化目标变量初始值;
所述结构尺寸优化目标变量包括用以评价出风均匀性的出口风量差和用以评价气 动噪声水平的偶极子噪声声功率级。
可采用CFD软件在最大风量输出工况下计算得到结构尺寸优化目标变量初始值。
(20)数学模型建立:根据风道原始结构形式,选择入口、出口、弯管内侧倒角、 弯管外侧倒角中至少两处结构尺寸作为结构优化变量,并根据风道整体结构的尺寸限制 确定其变化范围,以结构尺寸优化目标变量最小为优化目标,建立优化问题数学模型;
所述(20)数学模型建立步骤中,所述优化问题数学模型的表达式为:
其中,
结构优化变量x={x1,x2,…,xm},m为结构优化变量个数;
V-rnin表示最小化问题,即结构尺寸优化目标变量f(x)=[f1(x),f2(x),…,fn(x)]T中 的所有子目标变量都尽可能达到最小,n为结构尺寸优化目标变量个数;
gi(x)≤0为不等式约束条件,根据结构优化变量或结构尺寸优化目标变量的约束确 定,l为约束条件个数。
(30)结构优化变量采样:利用最优拉丁超立方设计方法,在结构优化变量的取值范围内进行结构优化变量采样,结构优化变量样本个数根据代理模型所需的最小值初步拟定;
可在Isight软件中,利用最优拉丁超立方设计方法,在结构优化变量的取值范围内 进行结构优化变量采样。
所述(30)结构优化变量采样步骤中,结构优化变量样本个数应不少于最高阶构建多项式响应面代理模型所需的最少样本个数,4阶多项式响应面代理模型所需的最少样 本个数计算公式为:
N=(m+1)(m+2)/2+3m
其中,N为最少样本个数,m为结构优化变量个数。
(40)样本建模仿真:根据结构优化变量样本组合方案分别建立新的CFD模型, 并分别在最大风量输出工况下计算得到每个方案的结构尺寸优化目标变量取值;
(50)构建代理模型:构建多项式响应面代理模型,拟合每个结构尺寸优化目标变量和结构优化变量之间的函数关系,并验证代理模型的精度,精度不满足时优先考虑增 加多项式阶数,其次增加样本个数;
所述(50)构建代理模型步骤包括:
(51)代理模型选择:初始可选择公式简单的低阶多项式或精度较高的高阶多项式, 构建多项式响应面代理模型;
(52)代理模型精度验证:根据每个方案的结构尺寸优化目标变量取值,分别拟合每个结构尺寸优化目标变量和结构优化变量之间的函数关系,并用方差分析的方法验证代理模型的精度;
(53)精度判断:将代理模型精度与精度阀值比较,当不满足要求时,优先考虑增加多项式阶数,再用方差分析的方法验证代理模型的精度,如果仍然不满足则增加样本 个数,转至(30)结构优化变量采样步骤;
(60)多目标优化:采用带有精英保留策略的快速非支配多目标优化算法,求解优化问题数学模型,获得使结构尺寸优化目标变量最小的结构优化变量取值,即风道优化 结构尺寸。
下面结合本发明的实施例作详细说明,本实施例在以本发明技术方案为前提下进行 实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
实施例
(10)原始方案建模仿真:根据风道原始结构尺寸,建立风道CFD模型,在最大 风量输出工况下计算得到结构尺寸优化目标变量初始值。建模仿真过程具体包括以下操 作:
1)根据原始风道的CAD图纸,用CATIA软件建立风道的三维模型,保存成.stp 文件导入ANSYS;
2)在ANSYS的Meshing模块对风道进行CFD网格的划分,并检查网格单元质量, 导出.msh文件到Fluent中。
3)在Fluent中设置边界条件,例如流体属性、入口边界、出口边界、壁面边界等。定义监测变量,例如监测点流速、压力等。
4)初始化后进行运算,完成后读取结构尺寸优化目标变量的数值,将计算结果的数据文件保存为.dat文件,供后续研究使用。
如图2所示,风道存在两侧直管段,因此选择以风道左侧直管段的最大风量差x、风道右侧直管段的最大风量差y来分别评价两侧直管段的出风均匀性,选择风道表面偶 极子噪声声功率级最大值z来评价气动噪声水平。最终,结构尺寸优化目标变量确定为 x,y,z三者。
仿真得到的三个结构尺寸优化目标变量的初始值如表1所示。
表1原始结构尺寸方案仿真结果
(20)数学模型建立:如图3所示,风道弯管内侧为垂直拐歪,需对其进行倒角处理,因此选择的三个结构优化变量分别为:出口1的半径a、弯管内侧倒角2的大小b、 入口3的面积c。其中,风道入口选择以面积为优化变量,其长和宽通过等比例缩放来 确定。结构优化变量的取值范围主要根据风道整体结构的尺寸限制来确定,表2所示即 为结构优化变量及其变化范围,以及对应的初始值。
表2结构优化变量的范围和初始值
风道的优化目标是既提高出风均匀性,又降低气动噪声,因此风道左侧直管段的最 大风量差x、风道右侧直管段的最大风量差y和风道表面偶极子噪声声功率级最大值z都要尽可能取到最小值,于是这里的优化问题就变为一个最小化问题。建立的数学模型 表达式为:
其中,s.t.是三个结构优化变量的约束条件。
(30)结构优化变量采样:在Isight软件中利用最优拉丁超立方设计方法在结构优化变量的取值范围内进行采样,软件中计算得到在3个结构优化变量的情况下,构建最 高4阶多项式响应面模型需要的最少样本个数为16个,因此为留有一定余量,初步选 择样本个数为20,20个样本中优化变量a,b,c的组合如表3所示。
(40)样本建模仿真:根据所有20个样本中的结构优化变量组合方案分别建立新的CFD模型,依次通过Fluent在最大风量输出工况下分别计算得到每个方案的结构尺 寸优化目标变量取值,仿真得到的结果如表3所示。
表3所有样本点仿真结果
(50)构建代理模型:提取表3中所有样本点的结构优化变量及其对应的结构尺寸优化目标变量数值,在Isight中将数据进行拟合,此处为保证代理模型的高精度直接建 立4阶多项式响应面模型。于是,得到的三个多项式响应面代理模型分别为:
1)风道左侧直管段最大风量差x
2)风道右侧直管段最大风量差y,
3)风道表面偶极子噪声声功率级最大值z。
代理模型建立好后,采用方差分析的方法来进行模型精度的验证,选取方差决定系 数R2和调整后决定系数Adjusted R2作为代理模型的精度评价指标。这两个指标的数值越接近1,说明模型的精度越高。表4中列出了本实施例的方差分析结果。
表4方差分析结果
由表4数据可知,方差分析的结果均在90%以上,且非常接近1,证明三个响应面模型都达到了很高的精度,完全可以满足后续优化的需要。
(60)多目标优化:采用带有精英保留策略的快速非支配多目标优化算法,对(20)的数学模型进行求解,获得使结构尺寸优化目标变量最小的结构优化变量取值,即风道 优化结构尺寸。其中,优化算法的基本参数配置如表5所示。
表5优化算法基本参数配置
经过软件的迭代计算后,风道优化结构尺寸的取值分别为:出口半径a=34.23mm,弯管内侧倒角b=48.23mm,入口面积c=9.21×104mm2,实际应用时应当根据加工要求 圆整为近似值。
根据风道优化结构尺寸建立风道CFD模型,通过Fluent在最大风量输出工况下计算得到结构尺寸优化目标变量的值,并与(10)中结构尺寸优化目标变量的初始值作对 比,验证优化方法对出风均匀性和气动噪声的改善效果,结构尺寸优化目标变量优化前 后结果对比如表6所示。
表6目标变量优化前后结果对比
由表6数据可知,风道两侧的最大风量差均得到了有效的降低,尤其是左侧风道,最大风量差已经削减到大约原始方案的一半,证明新的结构方案使得风道两侧直管段的出风均匀性得到了显著地提升。偶极子噪声声功率级的最大值也有所降低,证明风道表 面的噪声源强度得到了有效地削弱,气动噪声得到抑制。
Claims (5)
1.一种基于代理模型的风道结构尺寸优化方法,其特征在于,包括以下步骤:
(10)原始方案建模仿真:根据风道原始结构尺寸,建立风道CFD模型,在最大风量输出工况下计算得到结构尺寸优化目标变量初始值;
(20)数学模型建立:根据风道原始结构形式,选择入口、出口、弯管内侧倒角、弯管外侧倒角中至少两处结构尺寸作为结构优化变量,并根据风道整体结构的尺寸限制确定其变化范围,以结构尺寸优化目标变量最小为优化目标,建立优化问题数学模型;
(30)结构优化变量采样:利用最优拉丁超立方设计方法,在结构优化变量的取值范围内进行结构优化变量采样,结构优化变量样本个数根据代理模型所需的最小值初步拟定;
(40)样本建模仿真:根据结构优化变量样本组合方案分别建立新的CFD模型,并分别在最大风量输出工况下计算得到每个方案的结构尺寸优化目标变量取值;
(50)构建代理模型:构建多项式响应面代理模型,拟合每个结构尺寸优化目标变量和结构优化变量之间的函数关系,并验证代理模型的精度,精度不满足时,增加多项式阶数,仍不满足,则增加样本个数;
(60)多目标优化:采用带有精英保留策略的快速非支配多目标优化算法,求解优化问题数学模型,获得使结构尺寸优化目标变量最小的结构优化变量取值,即风道优化结构尺寸。
2.根据权利要求1所述风道结构尺寸优化方法,其特征在于:
所述结构尺寸优化目标变量包括用以评价出风均匀性的出口风量差和用以评价气动噪声水平的偶极子噪声声功率级。
4.根据权利要求1所述风道结构尺寸优化方法,其特征在于,所述(30)结构优化变量采样步骤中,结构优化变量样本个数应不少于构建最高阶多项式响应面代理模型所需的最少样本个数,4阶多项式响应面代理模型所需的最少样本个数计算公式为:
N=(m+1)(m+2)/2+2m
其中,N为最少样本个数,m为结构优化变量个数。
5.根据权利要求1所述风道结构尺寸优化方法,其特征在于,所述(50)构建代理模型步骤包括:
(51)代理模型选择:初始选择公式简单的低阶多项式或精度较高的高阶多项式,构建多项式响应面代理模型;
(52)代理模型精度验证:根据每个方案的结构尺寸优化目标变量取值,分别拟合每个结构尺寸优化目标变量和结构优化变量之间的函数关系,并用方差分析的方法验证代理模型的精度;
(53)精度判断:将代理模型精度与精度阀值比较,当不满足要求时,增加多项式阶数,再用方差分析的方法验证代理模型的精度,如果仍然不满足则增加样本个数,转至(30)结构优化变量采样步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910260231.9A CN110196987B (zh) | 2019-04-02 | 2019-04-02 | 基于代理模型的风道结构尺寸优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910260231.9A CN110196987B (zh) | 2019-04-02 | 2019-04-02 | 基于代理模型的风道结构尺寸优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110196987A CN110196987A (zh) | 2019-09-03 |
CN110196987B true CN110196987B (zh) | 2022-09-06 |
Family
ID=67751879
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910260231.9A Active CN110196987B (zh) | 2019-04-02 | 2019-04-02 | 基于代理模型的风道结构尺寸优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110196987B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110765569A (zh) * | 2019-09-05 | 2020-02-07 | 珠海格力电器股份有限公司 | 一种空调管路结构减振的多目标优化方法、计算机可读存储介质及终端 |
CN112800496A (zh) * | 2020-11-06 | 2021-05-14 | 大唐环境产业集团股份有限公司 | 一种烟风道结构设计智能计算模块和方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103699720B (zh) * | 2013-12-09 | 2016-05-04 | 浙江大学 | 基于区间约束违反度的高速压力机滑块机构尺寸优化方法 |
CN107357962A (zh) * | 2017-06-19 | 2017-11-17 | 西安电子科技大学 | 一种基于自适应代理模型的天线罩肋截面尺寸优化方法 |
-
2019
- 2019-04-02 CN CN201910260231.9A patent/CN110196987B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110196987A (zh) | 2019-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108153998B (zh) | 离心鼓风机叶轮的全三维气动优化设计方法 | |
CN111709096B (zh) | 一种强化自然对流换热的异型翅片结构设计方法 | |
CN105677964A (zh) | 一种基于cfd仿真和网格自适应的阀门流量系数计算方法 | |
CN110196987B (zh) | 基于代理模型的风道结构尺寸优化方法 | |
CN109766589B (zh) | 一种管翅式换热器非均匀迎面风速下性能评价方法 | |
CN110750855B (zh) | 一种外形定尺寸限制下的蜗壳型线设计方法 | |
CN111767663A (zh) | 一种基于cfd仿真的阀门流量系数便捷计算方法 | |
CN112417773B (zh) | 多级轴流膨胀机的多学科优化设计方法、装置及设备 | |
CN107450471A (zh) | 一种基于三次ph曲线插值实现刀具轨迹参数弧长化的方法 | |
CN114117877B (zh) | 一种基于等几何粒子描述的拓扑优化方法 | |
CN108446452B (zh) | 一种混流泵叶轮鲁棒优化设计方法 | |
CN111241752B (zh) | 一种基于数字孪生和强化学习的离心式叶轮综合优化方法 | |
CN113536640B (zh) | 一种基于正交试验的布风器内部流道结构的优化设计方法 | |
CN112651191B (zh) | 换热器的表面的风速分布确定方法、装置及空调设备 | |
Malanichev et al. | Pressure loss reduction in ventilation ducts by shape optimization of the removable profiled components | |
Michaeli et al. | Approach of an automatic extrusion die optimization | |
Jing et al. | A variable gradient descent shape optimization method for transition tee resistance reduction | |
CN113792443B (zh) | 通风系统管网水力平衡的分析方法、管网设计和调节方法 | |
CN104679953A (zh) | 一种用于户内变电站主变室流场温度场数值模拟的快速计算方法及系统 | |
CN113722856B (zh) | 一种管道泵进口导流叶片的自动建模及优化设计方法 | |
CN113435133B (zh) | 一种螺旋槽管优化方法 | |
CN112163288B (zh) | 基于大涡模拟的流体机械叶片优化设计方法 | |
CN117034584A (zh) | 一种转子叶片几何偏差影响压气机稳定性的分析方法 | |
CN109241585B (zh) | 一种高低压涡轮过渡流道型面反问题设计方法 | |
CN117436209B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |