CN106909747B - 一种热对流扩散系统中模糊参数隶属度函数识别方法 - Google Patents
一种热对流扩散系统中模糊参数隶属度函数识别方法 Download PDFInfo
- Publication number
- CN106909747B CN106909747B CN201710128743.0A CN201710128743A CN106909747B CN 106909747 B CN106909747 B CN 106909747B CN 201710128743 A CN201710128743 A CN 201710128743A CN 106909747 B CN106909747 B CN 106909747B
- Authority
- CN
- China
- Prior art keywords
- interval
- fuzzy
- temperature response
- variable
- variables
- 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/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/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种热对流扩散系统中模糊参数隶属度函数识别方法,步骤如下:引入模糊变量表征系统不确定参数;统计测量数据,得到模糊温度响应的测量值;选取截集水平,将待识别模糊变量转化为区间变量,将模糊温度响应测量值转化为区间温度响应测量值;利用顶点组合方法,快速计算各截集水平下区间温度响应的预测值;利用区间温度响应的测量值和预测值,建立误差函数;基于优化理论建立区间变量的识别模型;选用优化算法对区间变量识别模型进行编程计算,得到最优解;利用模糊分解定理重组各截集水平下区间变量的识别结果,得到模糊参数的隶属度函数。本发明在保证系统参数识别精度的前提下,可有效提高计算效率,这是一般商用软件所不能实现的。
Description
技术领域
本发明属于机械工程领域,具体涉及一种热对流扩散系统中模糊参数隶属度函数识别方法。
背景技术
在各种生产技术领域中,由温度差异引起的热能传递是一种极其普遍的物理现象。小到电子元件的开发,大到飞行器结构系统的设计,如何更为有效的实现热量传递是工程师所面临的重要问题。而在实际工程中,由于环境的复杂,某些系统参数无法或很难直接测量得到。以现场测量得到的响应信息为基础,通过反演技术来估计这些系统参数,是解决这类问题的有效方法。现在,反问题的研究已经遍及定向设计、无损探伤、扫描成像等现代化生产、生活的各个领域。与正问题相比,传热反问题的研究起步较晚,发展还远不成熟。
现有关于热分析的许多研究都是针对确定性模型而进行的,没有考虑模型输入参数的不确定性。由于制造工艺的限制、测量误差以及认知的局限,结构的材料属性、外部载荷和边界条件等物理参数不可避免的受到多种不确定因素的影响,使得传热系统的温度响应也表现出一定的波动。由于反问题与不确定性问题的双重复杂性,使得不确定传热反问题的研究更具挑战性。用随机变量对不确定性因素进行定量化描述的研究至今已经取得了许多成果,但概率模型的建立需要大量样本信息来事先确定其概率密度函数。而获得足够的样本数据往往花费较大或代价过高,这就限制了概率模型和分析方法的进一步推广。而在模糊不确定性分析中,尽管某些事物的概念或参数的数值是难以确定的,但可以根据实验数据或主观经验确定一个大致的范围。如此一来,模糊模型在不确定性建模方面表现出了很强的方便性和经济性。因此,在温度测量信息具有模糊不确定性的情况下,如何建立准确高效的计算方法,通过反演技术识别出系统模糊参数的隶属度函数,是目前学术领域的一个研究热点,对于弥补现有传热分析方法和反演技术的不足,具有重要的理论和工程应用价值。
发明内容
本发明所要解决的技术问题为:克服现有技术在传热系统参数识别问题中存在的不足,充分考虑温度测量信息具有的模糊不确定性,基于反演技术提出了一种更加先进的模糊参数识别方法。
本发明采用的技术方案为:一种热对流扩散系统中模糊参数隶属度函数识别方法,该方法包括:
步骤一:针对热对流扩散物理模型,引入模糊变量表征系统待识别的不确定参数;
步骤二:统计测量数据,得到模糊温度响应的测量值;
步骤三:选取截集水平,利用截集运算将步骤一中待识别模糊变量转化为区间变量,将步骤二中的模糊温度响应测量值转化为区间温度响应测量值;
步骤四:利用顶点组合方法,快速计算步骤三各截集水平下区间温度响应的预测值;
步骤五:利用步骤三和步骤四中区间温度响应的测量值和预测值,建立误差函数;
步骤六:以步骤三中区间变量上下界为设计变量,利用步骤五中的误差函数,基于优化理论建立区间变量的识别模型;
步骤七:选用优化算法对步骤六中的区间变量识别模型进行编程计算,得到最优解,作为区间变量的识别结果;
步骤八:利用模糊分解定理重组步骤七中各截集水平下的区间变量识别结果,最终得到模糊参数的隶属度函数。
其中,所述步骤三中,截集水平的选取并不是固定不变的;根据模糊温度响应测量值隶属度函数的分布类型来确定所需截集水平的数量规模和数值大小。
其中,所述步骤七中采用的优化算法并不是固定不变的;根据优化模型的类型和问题的规模来确定合适的优化算法。
该方法具体包括以下步骤:
步骤二:统计测量数据,得到模糊温度响应的测量值Ti F i=1,2,...,n,其中n为观测点的数量。
步骤三:在0到1范围内选取截集水平λ,利用截集运算可以将步骤一中待识别模糊变量转化为区间变量记为向量的形式其中α i,λ,为转化后区间变量的下界和上界,α λ,为转化后区间向量的下界和上界。同样的,利用截集运算可以将步骤二中的模糊温度响应测量值Ti F转化为区间温度响应测量值其中T i,λ,为的下界和上界。
其次,对于m个区间变量而言,可以得到M=2×2×…×2=2m个顶点组合,记为:
然后,利用有限元方法对上述M个顶点组合处的温度响应进行求解,通过筛选可以快速得到温度响应的下界和上界:
步骤五:利用步骤三和步骤四中区间温度响应的测量值和预测值,建立误差函数:
s.t.α λ,l≤α λ≤α λ,u
步骤七:选用优化算法对步骤六中的区间变量识别模型进行编程计算,定义最大循环次数Itermax和收敛因子ε,当如下三个条件中的任一个得到满足时,计算终止:
(1)循环迭代次数q>Itermax;
(2)在连续两次迭代过程中,目标函数相对变化量满足:
其中|| ||2表示向量的2范数。
对选定的所有截集水平重复上述操作,进而可得到各截集水平对应的区间变量识别结果。
本发明与现有技术相比的优点在于:
(1)与传统的参数识别方法相比,所建立的参数识别模型充分考虑到温度测量信息的模糊不确定性,计算结果对传热问题的反演分析具有更重要的指导意义。
(2)利用顶点组合方法对区间温度响应进行预测,计算精度高,实施方便,可有效降低计算耗费。
(3)本发明提出的模糊参数识别方法抗噪性能强,当温度测量信息受到不同程度噪声影响时,保证了系统参数的识别精度。
附图说明
图1为本发明的热对流扩散系统模糊参数隶属度函数识别流程;
图2为本发明的三维空气冷却系统模型示意图;
图3为模糊参数cF的隶属度函数曲线图;
图4为模糊参数uF的隶属度函数曲线图;
图5为模糊参数QF的隶属度函数曲线图。
具体实施方式
下面结合附图和实施例对本发明做进一步说明。
本发明适用于含有模糊不确定性的热对流扩散系统模糊参数隶属度函数识别问题。本发明实施方式以某三维空气冷却系统为例,具体说明所述的热对流扩散系统模糊参数隶属度函数识别方法。另外,此三维空气冷却系统的模糊参数隶属度函数识别方法可以推广到其他含有模糊不确定性的传热问题反演分析中。
此三维空气冷却系统模糊参数隶属度函数识别过程如图1所示,引入模糊变量表征系统不确定参数,统计模糊温度响应的测量值,选取截集水平,将模糊变量转化为区间变量,将模糊温度响应测量值转化为区间温度响应测量值,利用顶点组合方法计算区间温度响应的预测值,基于区间温度响应测量值和预测值间的误差函数,建立区间变量的识别模型,并选用优化算法对其进行编程计算,利用模糊分解定理重组各截集水平下区间变量的识别结果,得到模糊参数的隶属度函数。可分为如下几个步骤进行:
步骤一:如图2所示的空心圆柱体式三维空气冷却系统,长100cm,内外直径分别为10cm和20cm,实体结构7用640个六面体单元来离散,内部管道8用960个六面体单元来离散。实体结构有容积热产生,而入口处9有冷却空气以一定速度流经此管道。在管道中心线和实体结构外沿上分别选取节点1、2、3和4、5、6作为此热对流扩散系统的温度响应观测点。由于环境的复杂,空气的比热容c,流动速度u和结构容积热密度Q均无法直接测量,且均含有一定的模糊不确定性。引入模糊变量表征系统待识别的3个不确定模糊参数,记为向量的形式
步骤三:在0到1范围内选取11个截集水平λj=(j-1)×0.1j=1,...,11,为了表示方便,将其统一记为λ。利用截集运算可以将步骤一中待识别模糊变量转化为区间变量记为向量的形式其中α i,λ,为转化后区间变量的下界和上界,α λ,为转化后区间向量的下界和上界。同样的,利用截集运算可以将步骤二中的模糊温度响应测量值Ti F转化为区间温度响应测量值其中T i,λ,为的下界和上界。
其次,对于3个区间变量而言,可以得到M=23=8个顶点组合,记为:
然后,利用有限元方法对上述8个顶点组合处的温度响应进行求解,通过筛选可以快速得到温度响应的下界和上界:
步骤五:利用步骤三和步骤四中区间温度响应的测量值和预测值,建立误差函数:
步骤七:选用模拟退火算法对步骤六中的区间变量识别模型进行编程计算,定义最大循环次数Itermax=2000和收敛因子ε=10-4,当如下三个条件中的任一个得到满足时,计算终止:
(1)循环迭代次数q>Itermax;
(2)在连续两次迭代过程中,目标函数相对变化量满足:
其中|| ||2表示向量的2范数。
本实施算例中,当截集水平λ=0.4时,经过206次迭代计算,达到了上述第2条所示的终止条件,设计变量的最优解为c λ=922.4J/(kg·℃), u λ=4.38m/s, Q λ=4379.3W/m3,即λ=0.4截集水平下区间变量的识别结果。
对选定的所有截集水平重复上述操作,进而可得到各截集水平对应的区间变量识别结果。
本实施算例中,得到的3个模糊参数cF,uF,QF的隶属度函数曲线分别如图3-5所示。
以上所述的仅为本发明的较佳实施例而已,本发明不仅仅局限于上述实施例,凡在本发明的精神和原则之内所作的局部改动、等同替换、改进等均应包含在本发明的保护范围之内。
Claims (2)
1.一种热对流扩散系统中模糊参数隶属度函数识别方法,其特征在于包括以下步骤:
步骤一:针对热对流扩散物理模型,引入模糊变量表征系统待识别的不确定参数;
步骤二:统计测量数据,得到模糊温度响应的测量值;
步骤三:选取截集水平,利用截集运算将步骤一中待识别模糊变量转化为区间变量,将步骤二中的模糊温度响应测量值转化为区间温度响应测量值;
步骤四:利用顶点组合方法,快速计算步骤三各截集水平下区间温度响应的预测值;
步骤五:利用步骤三和步骤四中区间温度响应的测量值和预测值,建立误差函数;
步骤六:以步骤三中区间变量上下界为设计变量,利用步骤五中的误差函数,基于优化理论建立区间变量的识别模型;
步骤七:选用优化算法对步骤六中的区间变量识别模型进行编程计算,得到最优解,作为区间变量的识别结果;
步骤八:利用模糊分解定理重组步骤七中各截集水平下的区间变量识别结果,最终得到模糊参数的隶属度函数;
该方法具体包括以下步骤:
步骤二:统计测量数据,得到模糊温度响应的测量值Ti F i=1,2,...,n,其中n为观测点的数量;
步骤三:在0到1范围内选取截集水平λ,利用截集运算可以将步骤一中待识别模糊变量转化为区间变量记为向量的形式其中α i,λ,为转化后区间变量的下界和上界,α λ,为转化后区间向量的下界和上界,同样的,利用截集运算可以将步骤二中的模糊温度响应测量值Ti F转化为区间温度响应测量值其中T i,λ,为的下界和上界;
其次,对于m个区间变量而言,可以得到M=2×2×…×2=2m个顶点组合,记为:
然后,利用有限元方法对上述M个顶点组合处的温度响应进行求解,通过筛选可以快速得到温度响应的下界和上界:
步骤五:利用步骤三和步骤四中区间温度响应的测量值和预测值,建立误差函数:
s.t.α λ,l≤α λ≤α λ,u
步骤七:选用优化算法对步骤六中的区间变量识别模型进行编程计算,定义最大循环次数Itermax和收敛因子ε,当如下三个条件中的任一个得到满足时,计算终止:
(1)循环迭代次数q>Itermax;
(2)在连续两次迭代过程中,目标函数相对变化量满足:
其中|| ||2表示向量的2范数;
对选定的所有截集水平重复上述操作,进而可得到各截集水平对应的区间变量识别结果;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710128743.0A CN106909747B (zh) | 2017-03-06 | 2017-03-06 | 一种热对流扩散系统中模糊参数隶属度函数识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710128743.0A CN106909747B (zh) | 2017-03-06 | 2017-03-06 | 一种热对流扩散系统中模糊参数隶属度函数识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106909747A CN106909747A (zh) | 2017-06-30 |
CN106909747B true CN106909747B (zh) | 2020-06-30 |
Family
ID=59186232
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710128743.0A Active CN106909747B (zh) | 2017-03-06 | 2017-03-06 | 一种热对流扩散系统中模糊参数隶属度函数识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106909747B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109408926B (zh) * | 2018-10-12 | 2023-04-07 | 大连理工大学 | 求解复杂结构多维瞬态非线性热传导反问题的方法 |
CN114018972B (zh) * | 2021-11-04 | 2024-05-28 | 株洲国创轨道科技有限公司 | 基于分散模糊推理机制的固液相变材料表面热流测量方法及系统 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH04356611A (ja) * | 1991-06-03 | 1992-12-10 | Matsushita Electric Ind Co Ltd | 温風暖房機の制御装置 |
US5410890A (en) * | 1994-01-27 | 1995-05-02 | Fujitsu General Limited | Control apparatus of air-conditioner |
FR2836205A1 (fr) * | 2002-02-15 | 2003-08-22 | Renault | Procede de determination de la temperature instantanee d'un gaz sous pression, notamment d'un gaz pour une pile a combustible |
CN101916394A (zh) * | 2010-08-25 | 2010-12-15 | 清华大学 | 一种基于知识融合的在线软测量方法 |
CN102175460A (zh) * | 2011-01-28 | 2011-09-07 | 中国人民解放军海军工程大学 | 基于模糊Petri网的船用燃气轮机气路故障诊断方法 |
CN102425863A (zh) * | 2011-09-19 | 2012-04-25 | 河海大学 | 一种dsg槽式太阳能集热器出口蒸汽温度控制方法 |
CN104950674A (zh) * | 2015-06-12 | 2015-09-30 | 河海大学常州校区 | 基于反演模糊滑模控制的悬臂梁振动控制方法 |
CN105677995A (zh) * | 2016-01-12 | 2016-06-15 | 北京航空航天大学 | 一种基于全网格配点理论的模糊稳态热传导问题数值求解方法 |
CN105760586A (zh) * | 2016-02-03 | 2016-07-13 | 北京航空航天大学 | 一种基于配点理论的模糊温度响应隶属度函数求解方法 |
-
2017
- 2017-03-06 CN CN201710128743.0A patent/CN106909747B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH04356611A (ja) * | 1991-06-03 | 1992-12-10 | Matsushita Electric Ind Co Ltd | 温風暖房機の制御装置 |
US5410890A (en) * | 1994-01-27 | 1995-05-02 | Fujitsu General Limited | Control apparatus of air-conditioner |
FR2836205A1 (fr) * | 2002-02-15 | 2003-08-22 | Renault | Procede de determination de la temperature instantanee d'un gaz sous pression, notamment d'un gaz pour une pile a combustible |
CN101916394A (zh) * | 2010-08-25 | 2010-12-15 | 清华大学 | 一种基于知识融合的在线软测量方法 |
CN102175460A (zh) * | 2011-01-28 | 2011-09-07 | 中国人民解放军海军工程大学 | 基于模糊Petri网的船用燃气轮机气路故障诊断方法 |
CN102425863A (zh) * | 2011-09-19 | 2012-04-25 | 河海大学 | 一种dsg槽式太阳能集热器出口蒸汽温度控制方法 |
CN104950674A (zh) * | 2015-06-12 | 2015-09-30 | 河海大学常州校区 | 基于反演模糊滑模控制的悬臂梁振动控制方法 |
CN105677995A (zh) * | 2016-01-12 | 2016-06-15 | 北京航空航天大学 | 一种基于全网格配点理论的模糊稳态热传导问题数值求解方法 |
CN105760586A (zh) * | 2016-02-03 | 2016-07-13 | 北京航空航天大学 | 一种基于配点理论的模糊温度响应隶属度函数求解方法 |
Non-Patent Citations (3)
Title |
---|
"Adaptive neuro-fuzzy based inferential sensor model for estimating the average air temperature in space heating systems";S.Jassar等;《Building and Environment》;20081021;第44卷(第8期);第1609-1616页 * |
"Hybrid uncertain analysis for temperature field prediction with random, fuzzy and interval parameters";Chong Wang等;《International Journal of Thermal Sciences》;20150801;第98卷;第124-134页 * |
"基于边界元和分散式模糊推理的圆筒内壁形状识别算法";姜曙等;《世界科技研究与发展》;20150430;第37卷(第2期);正文第122-126页摘要,第1-4节 * |
Also Published As
Publication number | Publication date |
---|---|
CN106909747A (zh) | 2017-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106897520B (zh) | 一种含有模糊参数的传热系统可靠性分析方法 | |
CN113051831B (zh) | 机床热误差自学习预测模型建模方法及热误差控制方法 | |
CN108038325A (zh) | 一种3d打印技术制造的多孔框架结构宏观弹性性能可靠性预测方法 | |
CN112580855A (zh) | 基于自适应变异pso-bp神经网络的电缆群稳态温升预测方法 | |
CN106909747B (zh) | 一种热对流扩散系统中模糊参数隶属度函数识别方法 | |
Peng et al. | Physics-informed graph convolutional neural network for modeling fluid flow and heat convection | |
Yan et al. | Bi-directional prediction of structural characteristics and effective thermal conductivities of composite fuels through learning from finite element simulation results | |
Lee et al. | Machine learning-based inverse design methods considering data characteristics and design space size in materials design and manufacturing: a review | |
CN117521518A (zh) | 一种基于机器学习的镁合金热处理工艺优化方法 | |
CN108717505A (zh) | 一种基于k-rvfl的固化热过程时空建模方法 | |
Wen et al. | Cost reduction for data acquisition based on data fusion: Reconstructing the surface temperature of a turbine blade | |
Lu et al. | Three-dimensional temperature field inversion calculation based on an artificial intelligence algorithm | |
Decke et al. | DADO–low-cost query strategies for deep active design optimization | |
CN117195743B (zh) | 一种热障涂层裂纹结构的喷涂参数优化方法 | |
JP5457737B2 (ja) | プラント制御情報生成装置及び方法、並びにそのためのコンピュータプログラム | |
Deka et al. | Geometric size optimization of annular step fin using multi-objective genetic algorithm | |
Wang et al. | A Review of Intelligent Airfoil Aerodynamic Optimization Methods Based on Data-Driven Advanced Models | |
Song et al. | Topology optimization for cold plate using neural networks as proxy models | |
Ahrari et al. | Multi-objective optimization of cellular scanning strategy in selective laser melting | |
CN112084723B (zh) | 光纤预制棒一次拉伸工艺仿真方法及装置 | |
Xiao et al. | Isomorphism identification of kinematic chains using novel evolutionary approaches | |
CN115169265A (zh) | 基于数值分析的搅混系数分析方法、系统、设备和介质 | |
Gonzales et al. | Evaluation of a differential evolution/constructal design algorithm for geometrical optimization of double T-shaped cavity intruded into a heat generating wall | |
CN112182739A (zh) | 一种飞行器结构非概率可信可靠性拓扑优化设计方法 | |
Ding et al. | Online Dynamic Modelling for Digital Twin Enabled Sintering Systems: An Iterative Update Data‐Driven Method |
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 |