CN109883660A - 一种热模拟试验控制方法 - Google Patents
一种热模拟试验控制方法 Download PDFInfo
- Publication number
- CN109883660A CN109883660A CN201711258050.XA CN201711258050A CN109883660A CN 109883660 A CN109883660 A CN 109883660A CN 201711258050 A CN201711258050 A CN 201711258050A CN 109883660 A CN109883660 A CN 109883660A
- Authority
- CN
- China
- Prior art keywords
- flow density
- heat flow
- heat
- testpieces
- test
- 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.)
- Granted
Links
Landscapes
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
一种热模拟试验控制方法,本发明属于航空航天飞行器环境模拟技术领域。针对热模拟试验中温度或热流密度载荷谱与真实热载荷存在误差的问题，本发明提出一种热模拟试验控制方法：全方程热流密度控制，即通过测量试验件的表面温度实时核算出来热流密度命令值，将其与测得的试验件表面真实热流进行比较，计算机将相应的误差信号输出到功率控制器，进行热试验控制。本发明将气动加热与结构热响应的耦合效应和高温热力学参数随温度的变化对热流密度命令值的影响计入到了计算当中，更加真实的模拟了试验件的瞬态受热情况。
Description
技术领域
本发明涉及航空航天飞行器环境模拟技术领域，特别是用于高速飞行器瞬态地面热模拟试验控制方法。
背景技术
随着高速飞行器速度的不断提高，在超声速飞行时出现的“热障”问题变得越来越突出，严重的气动加热所产生的高温会降低材料的强度极限和飞行器结构的承载能力，使结构产生热变形，破坏部件的气动外形并影响飞行器的安全飞行。为保证高速飞行器的安全，确认飞行器的材料和结构是否能经得起高速飞行时所产生的热冲击及高温热应力破坏，须建立高速飞行器瞬态地面热模拟试验系统，模拟飞行器材料和结构在高速飞行时的真实受热状况，对高速飞行器热强度进行验证试验。
对飞行器材料和结构在高速飞行时的受热状况的模拟是否真实，主要取决于试验中所控制的物理量及控制方法。热模拟试验常用的控制方法有两种：温度控制和热流密度控制，这两种方法的共同点是需要预先给定温度或热流密度的控制曲线，而温度或热流密度控制曲线一般是根据气动加热参数和结构件简单的一维热分析模型计算得到的，无法考虑气动加热与结构热响应的耦合效应，此外，由于缺乏某些材料的高温热力学参数，使计算得到的温度或热流密度载荷谱与结构表面所承受的真实热载荷有比较大的差别。
本发明为一种新的热模拟试验控制方法，即全方程热流密度控制。所谓全方程热流密度控制即试验过程中热流密度命令值不是在试验开始前提前输入的，而是通过测量试验件的表面温度实时核算出来的，它是针对飞行器飞行状态实时模拟的一种控制方法。该控制方法可以弥补前两种控制方法的不足，可以模拟气动加热与结构热响应的耦合效应，以及材料的高温热力学参数随温度的变化。
发明内容
本发明的目的是提供一种热模拟试验控制方法，精确实现高速飞行器瞬态地面热模拟，使试验模拟更加准确，以达到真实模拟试验件(飞行器)受热情况的目的。
本发明的技术方案是：将试验委托方提供的飞行器轨道参数和相关的空气动力加热数据预先输入计算机，试验时用热电偶(温度传感器)测出试验件表面温度，并将此数据传给计算机，计算机就对任何瞬间按热流密度控制方程进行运算，得到该瞬间试验件表面应该达到的加热热流密度，即热流密度命令值。将此热流密度命令值与热流密度计所测得的试验件表面真实热流进行比较，如有差异，则计算机将相应的电信号(误差信号)输出到功率控制器中，进而改变辐射式加热器当中的石英灯的电压，从而进行热试验控制。
与现有技术相比，本发明在地面模拟试验中，将气动加热与结构热响应的耦合效应对热流密度命令值的影响、高温热力学参数随温度的变化对热流密度命令值的影响计入到了计算当中，全方程热流密度控制方法针对飞行器飞行状态实时模拟的一种控制方法，因此更加真实的模拟了试验件的瞬态受热情况。
附图说明
图1为本发明全方程热流密度控制方法框图；
具体实施方式
下面通过实施例对本发明做进一步的说明：
一种热模拟试验控制方法针对某飞行器结构件进行全方程热流密度控制，其具体步骤为：
步骤1：确定飞行器(试验件)气动加热热流密度控制方程，由试验委托方给定。由气动加热理论可知，在热流密度控制方程中除黑度系数ε和波尔兹曼常数σ这两个常数外，其它各参数(包括附面层底部气体与飞行器壁面间的对流换热系数h、附面层底部气体在阻尼温度下的定压比热Cp、附面层底部气体的恢复焓Je、壁面温度下的热焓Jw)均为飞行器表面温度的函数，所以输入飞行器的热流密度最终可表示为飞行器表面温度的函数，因此只要测得飞行器 (试验件)表面温度，就可确定热流密度值；
步骤2：标定热流密度损失项qloss。地面热模拟试验采取石英灯辐射加热模拟飞行器的热流密度输入，由于存在空气对流和试件本身的热辐射损失，试验运行过程中必须对热流密度理论计算值进行修正，补偿由于空气对流和试件本身的热辐射产生的热流密度损失，而热流密度损失与试验件与空气自然对流换热系数及试验件表面温度等多种因素有关，通过计算直接得到总的热损失非常困难，所以热流密度损失项通过单独热流密度损失标定试验方法来测定；
步骤3：将热流密度损失项及飞行器气动加热参数引入控制系统。将热流密度损失项qloss、黑度系数ε、波尔兹曼常数σ、附面层底部气体与飞行器壁面间的对流换热系数h、附面层底部气体在阻尼温度下的定压比热Cp、附面层底部气体的恢复焓Je、壁面温度下的热焓Jw等参数引入控制系统。这个过程为载荷谱的编制过程，对于载荷谱的编制分两步进行：1.为每一个需要引入的参数各自添加一个通道，然后排列对应通道的控制载荷谱，即输入同一时间点内各通道对应的载荷，若采用美国MTS控制系统，控制载荷谱在Load Table中编制；2.排列飞行谱，即排列两级加载段之间的运行时间，飞行谱在Profile中编制。对于以温度为自变量的物理量可如热电偶分度表插值一样分段进行线性插值，可通过“基于信号的命令”功能(SBC，Signal Based Command)进行输入，然后将SBC表格在各自通道的控制载荷谱选定。
步骤4：用贴在飞行器表面的温度传感器测得飞行器表面温度TS。
步骤5：将飞行器气动加热参数、飞行器(试验件)表面温度及热流密度损失项代入热流密度控制方程中计算出气动加热热流密度值。对应试验程序中，先添加一个通道，在此通道中输入热流密度控制方程，再将其他通道中所计算得到的、所测量得到的物理量及引入的物理量引用到该通道的控制方程中，进而计算出热流密度值qsta。
步骤6：用贴在飞行器表面的热流密度计测得当前时刻热流密度反馈值。
步骤7：将计算得到的热流密度值qsta作为下一时刻热流密度命令值，与当前时刻测得的热流密度反馈值进行比较。
步骤8：根据二者之差调整功率控制器电压，实现对结构热试验的全方程热流密度控制。热流密度命令值大于热流密度反馈值，控制系统给功率控制器正电压信号，功率控制器会加大功率，进而石英灯加热器也会加大功率；热流密度命令值小于热流密度反馈值，控制系统给功率控制器负电压信号，功率控制器会减小功率，进而石英灯加热器也会减小功率，以此来实现地面热模拟试验的全方程热流密度控制。
Claims (4)
1.一种热模拟试验控制方法，其特征在于：包括以下步骤：
步骤1：确定试验件气动加热热流密度控制方程，根据气动加热理论得到试件气动加热热流密度控制方程，在热流密度方程中除黑度系数ε和波尔兹曼常数σ这两个常数外，其它各参数均为试验件表面温度的函数，即输入试验件的热流最终可表示为表面温度的函数，测得试验件表面温度，就可确定热流密度载荷值；
步骤2：标定热损失项，通过单独热损失标定试验方法来测定热损失项，在试验过程中补偿由于对流和热辐射产生的热损失；
步骤3：将试验件轨道参数及相关的空气动力加热数据引入计算机，即控制载荷谱的编制过程，首先为每一个需要引入的参数各自添加一个通道，排列对应通道的控制载荷谱，即输入同一时间点内各通道对应的载荷；其次排列飞行谱，即排列两级加载段之间的运行时间；
步骤4：用贴在试验件表面的温度传感器测得试验件表面温度；
步骤5：将试验件轨道参数、相关空气动力加热数据、试验件表面温度及热损失项代入热流密度控制方程计算出气动加热热流密度值，对应试验程序中，先添加一个通道中，在此通道中输入热流密度控制方程，再将其他通道中所计算得到的、所测量得到的物理量及输入的物理量引入到该通道中，进而计算出热流密度值；
步骤6：用贴在试验件表面的热流密度计测得当前时刻热流密度反馈值；
步骤7：将计算得到的热流密度值作为下一时刻热流密度命令值，与当前时刻热流密度反馈值进行比较；
步骤8：根据二者比较结果调整功率控制器电压，实现对结构热试验的全方程热流密度控制，热流密度命令值大于热流密度反馈值，控制系统给功率控制器正电压信号，功率控制器会加大功率，进而试验加热器也会加大功率；热流密度命令值小于热流密度反馈值，控制系统给功率控制器负电压信号，功率控制器会减小功率，进而试验加热器也会减小功率。
2.根据权利要求1所述的一种热模拟试验控制方法，其特征在于：步骤2标定热损失项时，选用标定件为与试验件为同一种材料的小板。
3.根据权利要求1所述的一种热模拟试验控制方法，其特征在于：步骤3中控制系统选用美国MTS控制系统，弹道参数的引入选用“基于信号的命令”功能实现。
4.根据权利要求1所述的一种热模拟试验控制方法，其特征在于：步骤5利用热流密度方程计算气动加热热流密度值是在虚拟控制通道中通过编程实现。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711258050.XA CN109883660B (zh) | 2017-12-01 | 2017-12-01 | 一种热模拟试验控制方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711258050.XA CN109883660B (zh) | 2017-12-01 | 2017-12-01 | 一种热模拟试验控制方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109883660A true CN109883660A (zh) | 2019-06-14 |
| CN109883660B CN109883660B (zh) | 2021-09-14 |
Family
ID=66923055
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711258050.XA Active CN109883660B (zh) | 2017-12-01 | 2017-12-01 | 一种热模拟试验控制方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109883660B (zh) |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111581720A (zh) * | 2020-04-30 | 2020-08-25 | 中国飞机强度研究所 | 一种评定飞行器全方程控热试验温度数据不确定度的方法 |
| CN112937913A (zh) * | 2021-02-03 | 2021-06-11 | 中国空气动力研究与发展中心超高速空气动力研究所 | 电弧加热设备上中焓包罩试验状态自动调试方法和装置 |
| CN112986038A (zh) * | 2021-05-11 | 2021-06-18 | 中国飞机强度研究所 | 一种热强度试验系统控制方法 |
| CN113155885A (zh) * | 2021-03-30 | 2021-07-23 | 中国飞机强度研究所 | 一种石英灯辐射加热试验热损失标定方法及标定装置 |
| CN113252493A (zh) * | 2021-07-13 | 2021-08-13 | 中国飞机强度研究所 | 一种热强度试验系统控制方法 |
| CN113504064A (zh) * | 2021-07-09 | 2021-10-15 | 哈尔滨工业大学 | 一种在线仿真驱动的飞行器结构热力联合试验系统及方法 |
| CN114674546A (zh) * | 2022-05-30 | 2022-06-28 | 中国飞机强度研究所 | 空天飞机测试用复杂热场下曲面结构高温热强度实验方法 |
| CN114706295A (zh) * | 2022-06-07 | 2022-07-05 | 中国飞机强度研究所 | 用于空天飞机强度测试的热试验中快时变热载荷控制方法 |
| CN114722543A (zh) * | 2022-06-09 | 2022-07-08 | 中国飞机强度研究所 | 一种高超声速飞行器结构热强度试验中热反射屏设计方法 |
| CN114721450A (zh) * | 2022-06-10 | 2022-07-08 | 中国飞机强度研究所 | 空天飞机结构强度测试中温度控制用热试验迭代控温方法 |
| CN114818407A (zh) * | 2022-03-11 | 2022-07-29 | 西北工业大学 | 一种用于结构强度分析的虚拟热试验方法 |
| CN114815931A (zh) * | 2022-06-23 | 2022-07-29 | 中国飞机强度研究所 | 一种极端高温环境下飞机构件热试验温度控制方法 |
| CN114940266A (zh) * | 2021-12-31 | 2022-08-26 | 成都流体动力创新中心 | 一种复杂防冰腔所能维持的蒙皮表面温度预测方法及系统 |
| CN115356372A (zh) * | 2022-10-24 | 2022-11-18 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种新型材料在飞行试验中的时变热响应测试方法及系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104015942A (zh) * | 2014-06-16 | 2014-09-03 | 北京卫星环境工程研究所 | 航天器真空热试验超高温度热流模拟系统 |
| CN104535438A (zh) * | 2014-12-31 | 2015-04-22 | 北京航空航天大学 | 一种试验件高温高低周复合疲劳裂纹扩展试验系统及测量方法 |
| CN105548250A (zh) * | 2016-01-13 | 2016-05-04 | 北京机电工程研究所 | 飞行器气动热试验的热流控制方法、装置及系统 |
| CN106557597A (zh) * | 2015-09-29 | 2017-04-05 | 中国飞机强度研究所 | 一种热虚拟试验方法 |
-
2017
- 2017-12-01 CN CN201711258050.XA patent/CN109883660B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104015942A (zh) * | 2014-06-16 | 2014-09-03 | 北京卫星环境工程研究所 | 航天器真空热试验超高温度热流模拟系统 |
| CN104535438A (zh) * | 2014-12-31 | 2015-04-22 | 北京航空航天大学 | 一种试验件高温高低周复合疲劳裂纹扩展试验系统及测量方法 |
| CN106557597A (zh) * | 2015-09-29 | 2017-04-05 | 中国飞机强度研究所 | 一种热虚拟试验方法 |
| CN105548250A (zh) * | 2016-01-13 | 2016-05-04 | 北京机电工程研究所 | 飞行器气动热试验的热流控制方法、装置及系统 |
Non-Patent Citations (1)
| Title |
|---|
| 秦强等: "临近空间高超声速飞行器地面热防护试验技术", 《飞机设计》 * |
Cited By (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111581720A (zh) * | 2020-04-30 | 2020-08-25 | 中国飞机强度研究所 | 一种评定飞行器全方程控热试验温度数据不确定度的方法 |
| CN112937913A (zh) * | 2021-02-03 | 2021-06-11 | 中国空气动力研究与发展中心超高速空气动力研究所 | 电弧加热设备上中焓包罩试验状态自动调试方法和装置 |
| CN112937913B (zh) * | 2021-02-03 | 2022-07-19 | 中国空气动力研究与发展中心超高速空气动力研究所 | 电弧加热设备上中焓包罩试验状态自动调试方法和装置 |
| CN113155885A (zh) * | 2021-03-30 | 2021-07-23 | 中国飞机强度研究所 | 一种石英灯辐射加热试验热损失标定方法及标定装置 |
| CN112986038A (zh) * | 2021-05-11 | 2021-06-18 | 中国飞机强度研究所 | 一种热强度试验系统控制方法 |
| CN113504064A (zh) * | 2021-07-09 | 2021-10-15 | 哈尔滨工业大学 | 一种在线仿真驱动的飞行器结构热力联合试验系统及方法 |
| CN113252493A (zh) * | 2021-07-13 | 2021-08-13 | 中国飞机强度研究所 | 一种热强度试验系统控制方法 |
| CN113252493B (zh) * | 2021-07-13 | 2021-10-01 | 中国飞机强度研究所 | 一种热强度试验系统控制方法 |
| CN114940266A (zh) * | 2021-12-31 | 2022-08-26 | 成都流体动力创新中心 | 一种复杂防冰腔所能维持的蒙皮表面温度预测方法及系统 |
| CN114818407A (zh) * | 2022-03-11 | 2022-07-29 | 西北工业大学 | 一种用于结构强度分析的虚拟热试验方法 |
| CN114818407B (zh) * | 2022-03-11 | 2024-03-12 | 西北工业大学 | 一种用于结构强度分析的虚拟热试验方法 |
| CN114674546A (zh) * | 2022-05-30 | 2022-06-28 | 中国飞机强度研究所 | 空天飞机测试用复杂热场下曲面结构高温热强度实验方法 |
| CN114706295A (zh) * | 2022-06-07 | 2022-07-05 | 中国飞机强度研究所 | 用于空天飞机强度测试的热试验中快时变热载荷控制方法 |
| CN114706295B (zh) * | 2022-06-07 | 2022-08-26 | 中国飞机强度研究所 | 用于空天飞机强度测试的热试验中快时变热载荷控制方法 |
| CN114722543A (zh) * | 2022-06-09 | 2022-07-08 | 中国飞机强度研究所 | 一种高超声速飞行器结构热强度试验中热反射屏设计方法 |
| CN114722543B (zh) * | 2022-06-09 | 2022-08-12 | 中国飞机强度研究所 | 一种高超声速飞行器结构热强度试验中热反射屏设计方法 |
| CN114721450A (zh) * | 2022-06-10 | 2022-07-08 | 中国飞机强度研究所 | 空天飞机结构强度测试中温度控制用热试验迭代控温方法 |
| CN114721450B (zh) * | 2022-06-10 | 2022-08-26 | 中国飞机强度研究所 | 空天飞机结构强度测试中温度控制用热试验迭代控温方法 |
| CN114815931B (zh) * | 2022-06-23 | 2022-09-13 | 中国飞机强度研究所 | 一种极端高温环境下飞机构件热试验温度控制方法 |
| CN114815931A (zh) * | 2022-06-23 | 2022-07-29 | 中国飞机强度研究所 | 一种极端高温环境下飞机构件热试验温度控制方法 |
| CN115356372A (zh) * | 2022-10-24 | 2022-11-18 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种新型材料在飞行试验中的时变热响应测试方法及系统 |
| CN115356372B (zh) * | 2022-10-24 | 2023-03-10 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种新型材料在飞行试验中的时变热响应测试方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109883660B (zh) | 2021-09-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109883660A (zh) | 一种热模拟试验控制方法 | |
| CN112987566B (zh) | 一种气动热的超螺旋非线性分数阶滑模无模型控制方法 | |
| CN105173128B (zh) | 模拟飞行器飞行过程中热环境的试验方法 | |
| CN112965364B (zh) | 一种高超声速飞行器气动热地面模拟系统建模及iPID方法 | |
| Patil et al. | Inverse determination of temperature distribution in partially cooled heat generating cylinder | |
| CN112859617B (zh) | 一种iPI无模型自适应全局非奇异快速终端滑模控制方法 | |
| CN113625564B (zh) | 一种结构热试验基于系统模型的滑模控制方法 | |
| CN112987569B (zh) | 一种扰动上界自适应分数阶全局终端滑模无模型控制方法 | |
| CN114637277B (zh) | 空天飞机测试实验的全方程热流密度测控系统及测控方法 | |
| CN114355779B (zh) | 一种结构热试验非线性扩展状态观测器全局滑模控制方法 | |
| CN112965384B (zh) | 一种自适应神经网络全局非奇异快速终端滑模无模型方法 | |
| CN112965382A (zh) | 一种基于神经网络观测器非线性全局滑模无模型控制方法 | |
| CN113900373B (zh) | 一种结构热试验的全局非奇异快速终端滑模控制方法 | |
| CN115236975B (zh) | 一种导弹气动热地面有限元分析及ipd非线性滑模控制方法 | |
| CN113495486A (zh) | 一种结构热试验基于扩展状态观测器的模型预测控制方法 | |
| CN114139420A (zh) | 一种石英灯辐射加热虚拟试验方法 | |
| CN114706295B (zh) | 用于空天飞机强度测试的热试验中快时变热载荷控制方法 | |
| CN117922841A (zh) | 一种提升高超声速飞行器结构热试验模拟精度的控制方法 | |
| Nguyen et al. | Sensitivity of predicted temperature in a fillet weld T-joint to parameters used in welding simulation with prescribed temperature approach | |
| Sulamanidze et al. | Verification and validation of multi-physics numerical analysis of thermomechanical fatigue test conditions under induction heating and forced convection | |
| CN114721450A (zh) | 空天飞机结构强度测试中温度控制用热试验迭代控温方法 | |
| Lalande et al. | Aerothermal Computations for Laminar–Turbulent Transition Onset Measurement Using Infrared Imaging Technique | |
| US10801925B2 (en) | Test bypass for a cooling apparatus, having a liquid vessel with a variable pressure level | |
| Chu et al. | Dynamic Compensation Method for Temperature Measurement Error of High Temperature and High-speed Airflow | |
| Lua et al. | A Multiphysics Coupling for Evaluation of Effects of Local Boundary Conditions on Autoclave-Cured Composite |
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 |