CN108120612A - 一种热试验高温峰值的处理方法 - Google Patents
一种热试验高温峰值的处理方法 Download PDFInfo
- Publication number
- CN108120612A CN108120612A CN201611067880.XA CN201611067880A CN108120612A CN 108120612 A CN108120612 A CN 108120612A CN 201611067880 A CN201611067880 A CN 201611067880A CN 108120612 A CN108120612 A CN 108120612A
- Authority
- CN
- China
- Prior art keywords
- temperature
- test
- heat flow
- time
- hot wall
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M99/00—Subject matter not provided for in other groups of this subclass
- G01M99/002—Thermal testing
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明提出一种热试验高温峰值的处理方法，通过确定理论加热温度曲线和理论热壁热流曲线、判断是否进行高温峰值处理、高温峰值处理、判断是否满足总加热量等步骤确立达到试验要求的加热温度曲线和热壁热流曲线用于热试验。本发明将试验条件的温度限制作为约束条件，导入到冷壁热流和恢复焓控制的热试验中，使热环境模拟更为精确，大大提高了试验成功率，提高了工作效率。
Description
技术领域
本发明涉及一种热试验高温峰值的处理方法，属于热试验设计技术领域。
背景技术
伴随着飞行器飞行速度的逐步提高，飞行器所受气动热越来越大，为了模拟飞行器的受热环境，检验热防护效果，需要进行地面加热试验，加热源一般是石英灯或石墨加热器。
在这类地面加热试验中，因为受到加热器材料特质、温度测量以及加热功率等限制，一般都有加热上限。某些试验中若出现模拟的热环境超过了这个加热限制，一般会采用以下几种方式进行处理，一种是在加热上限温度下进行试验，直至试验结束，这样可以会由于热量聚集，导致试验失败；另一种是根据经验确定一个在加热上限温度的时间段，这样设计的温度曲线不能很好的模拟所需的热环境，造成试验结果不准确。
发明内容
本发明的目的在于克服现有技术不足，提供了一种温度模拟精确、提高试验成功率的热试验高温峰值的处理方法。
本发明的技术解决方案：一种热试验高温峰值的处理方法，包括以下步骤：
第一步，确定理论加热温度曲线和理论热壁热流曲线，
A1.1、利用公式(1)确定当前时刻t的热壁焓值Hwt，
其中，qw(t-1)为前一时刻的热壁热流，Hrt为当前时刻的恢复焓，qct为当前时刻的冷壁热流，Hcool为冷壁焓值，Tw(t-1)为前一时刻的温度，Tref为参考温度，ε为热试验材料的灰度、σ为斯忒藩.波尔兹曼常数；试验时间t∈[t0,tJ],t0为试验开始时间，tJ为试验结束时间，试验开始时的取值为试验开始时的冷壁热流试验开始时的温度为试验环境温度；
设计指标给定了的热环境参数冷壁热流qc和恢复焓Hr，均为沿试验时间的曲线，根据试验的时间可以从曲线上插值得到当前时刻的恢复焓Hrt和冷壁热流qct；参考温度Tref也是设计指标给定的已知值，冷壁焓值Hcool为参考温度下得到的焓值。
A1.2、根据步骤A1.1得到的热壁焓值Hwt，利用公式(2)确定当前时刻的热壁热流qwt，
A1.3、判断是否满足|qwt-qw(t-1)|≤Δqw，若满足则转入步骤A1.7,若不满足，则转入步骤A1.4,其中Δqw为预设的热壁热流阈值，热壁热流阈值越小，迭代得到的热壁热流越精确，但迭代次数越多，运算量越大，本领域技术人员根据实际情况进行阈值大小确定，工程中一般取值为冷壁热流qc幅值的0.1％～1％。
A1.4、利用公式(3)确定当前时刻的温度Twt，
其中λ是热试验材料的热导率，x是沿热试验材料法向的位移；
A1.5、根据步骤A1.4确定的当前时刻的温度Twt和步骤A1.1得到的热壁焓值Hwt，利用公式(4)确定当前时刻的热壁热流qwt，
A1.6、判断是否满足|qwt-qw(t-1)|≤Δqw，若满足则转入步骤A1.8,若不满足，则转入步骤A1.4；
A1.7、确定前一时刻的温度Tw(t-1)为当前时刻的温度Twt，转入步骤A1.9；
A1.8、确定此迭代步骤得到的温度为当前时刻的温度Twt，转入步骤A1.9；
A1.9、判断当前时刻是否到试验结束时间，若没到试验结束时间则将下一时刻确定为当前时刻，转入步骤A1.1，若到试验结束时间，得到试验时间内的理论加热温度曲线和理论热壁热流曲线；
第二步，判断是否进行高温峰值处理，若第一步确定的理论加热温度曲线的最高温度K1超过热试验最高加热温度K2则转入第三进行高温峰值处理，若没超过则认为第一步确定的理论加热温度曲线和理论热壁热流曲线符合试验要求，转入第五步；
热试验最高加热温度K2是由热试验所使用的加热设备决定，加热设备能加热的最高温度。
第三步，高温峰值处理，得到高温峰值处理后的加热温度曲线和热壁热流曲线，
A3.1、确定削峰区间初值值T1和T2；
T1为削峰区间开始时间，即第一步得到的理论加热温度曲线中温度刚到K2的对应时间，T2为削峰区间结束时间，即第一步得到的理论加热温度曲线中温度降到K2后的某个时间，T2≥T2′，T2′为理论加热温度曲线中温度刚降低到小于K2的对应时间T2′，本领域技术人员根据实际情况进行选择。
A3.2利用公式(1)确定当前时刻t的热壁焓值Hwt，利用公式(2)确定当前时刻的热壁热流qwt；
A3.3、判断是否满足|qwt-qw(t-1)|≤Δqw，若满足则转入步骤A3.6,若不满足，则转入步骤A3.4；
A3.4、利用公式(3)确定当前时刻的温度Twt，再利用公式(4)确定当前时刻的热壁热流qwt，
A3.5、判断是否满足|qwt-qw(t-1)|≤Δqw，若满足则转入步骤A3.7,若不满足，则转入步骤A3.4；
A3.6、判断当前时刻是否在步骤A3.1确定的削峰区间内，若在则转入步骤A3.8,若不在则确定前一时刻的温度Tw(t-1)为当前时刻的温度Twt，转入步骤A3.12；
A3.7、确定此迭代步骤得到的温度为当前时刻的温度Twt，判断当前时刻是否在步骤A3.1确定的削峰区间内，若在则转入步骤A3.9,若不在则转入步骤A3.12；
A3.8、确定前一时刻的温度Tw(t-1)为当前时刻的温度Twt；
A3.9、判断是否满足Twt-K2≤ΔT，若满足则转入步骤A3.12，若不满足则转入步骤A3.10，ΔT为温度容差；
温度容差ΔT的值越小，确定的削峰区间就越精确，但迭代次数会增加，从而增加运算量，工程中一般取值为加热设备能加热的最高温度K2的0.05％～0.15％，本领域技术人员根据实际情况进行选择。
A3.10、利用公式(5)得到修正后的热壁热流q′w，
q′wt＝qwt-Δqt*ΔT (5)
其中Δqt为当前时刻的热壁热流修正值，且Δqt＝qwt/K2；
A3.11、令qwt＝q′wt，转入步骤A3.5；
A3.12、判断当前时刻是否到试验结束时间，若没到试验结束时间则将下一时刻确定为当前时刻，转入步骤A3.2，若到试验结束时间，得到试验时间内的高温峰值处理后的加热温度曲线和热壁热流曲线；
第四步，判断是否满足|Qw-Q′w|≤ΔQw，若满足，则认为第三步得到的高温峰值处理后的加热温度曲线和热壁热流曲线满足试验要求，转入第五步，若不满足，则调整第三步中的重新设置削峰区间结束时间T2，重复第三步直至满足|Qw-Q′w|≤ΔQw，其中Qw为第一步确定的理论热壁热流曲线得到的在试验时间内的总加热量，Q′w为第三步确定的高温峰值处理后的热壁热流曲线得到的在试验时间内的总加热量，ΔQw为总热量阈值；
对热壁热流曲线进行时间轴上的积分得到总价热量。当不满足|Qw-Q′w|≤ΔQw时，若Qw＞Q′w则将削峰区间结束时间T2往回调即缩短削峰区间；若Qw<Q′w则将削峰区间结束时间T2往后调即延长削峰区间。
总热量阈值ΔQw取值越少，理论和实际的总加热量就越接近，确定的削峰区间就越精确，但迭代次数会增加，从而增加运算量，本领域技术人员可以根据具体情况进行选择，工程中ΔQw的取值一般为理论总加热量的0.1％～1％。
第五步，利用第二步或第四步确认达到试验要求的加热温度曲线和热壁热流曲线进行试验。
本发明与现有技术相比的有益效果：
(1)本发明将试验条件的温度限制作为约束条件，导入到冷壁热流和恢复焓控制的热试验中，使热环境模拟更为精确，大大提高了试验成功率，提高了工作效率；
(2)本发明先进行削峰时间段的确定以及削峰效果预测，减少了试验的测试次数和试验件的数量，节约了成本；
(3)本发明约束控制方法简单可行，为军用、民用领域的高速飞行器的地面模拟加热试验设计提供了一种合理的解决办法。
附图说明
图1为本发明实施例确定的理论加热温度曲线，最高温度K1超过试验上限K2；
图2为本发明实施例削峰后的温度曲线；
图3为本发明流程图。
具体实施方式
下面结合具体实例及附图对本发明进行详细说明。
本发明如图3所示，通过以下步骤实现：
1、确定理论加热温度曲线和理论热壁热流曲线
(1)利用公式(1)确定当前时刻t的热壁焓值Hwt，
其中，qw(t-1)为前一时刻的热壁热流，Hrt为当前时刻的恢复焓，qct为当前时刻的冷壁热流，Hcool为冷壁焓值，Tw(t-1)为前一时刻的温度，Tref为参考温度，ε为热试验材料的灰度、σ为斯忒藩.波尔兹曼常数；试验时间t∈[t0,tJ],t0为试验开始时间，tJ为试验结束时间，试验开始时的取值为试验开始时的冷壁热流试验开始时的温度为试验环境温度；
根据试验的时间从冷壁热流qc和恢复焓Hr曲线上插值得到当前时刻的恢复焓Hrt和冷壁热流qct；本实例中参考温度Tref为280K，冷壁焓值Hcool为280K下得到的焓值。
(2)根据步骤(1)得到的热壁焓值Hwt，利用公式(2)确定当前时刻的热壁热流qwt，
(3)判断是否满足|qwt-qw(t-1)|≤Δqw，若满足则转入步骤(7),若不满足，则转入步骤(4),本实例中热壁热流阈值Δqw取冷壁热流qc幅值的0.1％。
(4)利用公式(3)确定当前时刻的温度Twt，
(5)根据步骤(4)确定的当前时刻的温度Twt和步骤(1)得到的热壁焓值Hwt，利用公式(4)确定当前时刻的热壁热流qwt，
(6)判断是否满足|qwt-qw(t-1)|≤Δqw，若满足则转入步骤(8),若不满足，则转入步骤(4)；
(7)确定前一时刻的温度Tw(t-1)为当前时刻的温度Twt，转入步骤(9)；
(8)确定此迭代步骤得到的温度为当前时刻的温度Twt，转入步骤(9)；
(9)判断当前时刻是否到试验结束时间，若没到试验结束时间则将下一时刻确定为当前时刻，转入步骤(1)，若到试验结束时间，得到试验时间内的理论加热温度曲线和理论热壁热流曲线。
本实例得到如图1所示的理论加热温度曲线。
2、根据图1，理论加热温度曲线的最高温度K1超过热试验最高加热温度K2，对其进行步骤3高温峰值处理。
3、高温峰值处理，得到高温峰值处理后的加热温度曲线和热壁热流曲线
(1)确定削峰区间初值值T1和T2；
T1为削峰区间开始时间，即第一步得到的理论加热温度曲线中温度刚到K2的对应时间，T2为削峰区间结束时间，本实例中取理论加热温度曲线中温度刚降低到小于K2的对应时间T2′。
(2)利用公式(1)确定当前时刻t的热壁焓值Hwt，利用公式(2)确定当前时刻的热壁热流qwt；
(3)判断是否满足|qwt-qw(t-1)|≤Δqw，若满足则转入步骤(6),若不满足，则转入步骤(4)；
(4)利用公式(3)确定当前时刻的温度Twt，再利用公式(4)确定当前时刻的热壁热流qwt；
(5)判断是否满足|qwt-qw(t-1)|≤Δqw，若满足则转入步骤(7),若不满足，则转入步骤(4)；
(6)判断当前时刻是否在步骤(1)确定的削峰区间内，若在则转入步骤(8),若不在则确定前一时刻的温度Tw(t-1)为当前时刻的温度Twt，转入步骤(12)；
(7)确定此迭代步骤得到的温度为当前时刻的温度Twt，判断当前时刻是否在步骤(1)确定的削峰区间内，若在则转入步骤(9),若不在则转入步骤(12)；
(8)确定前一时刻的温度Tw(t-1)为当前时刻的温度Twt；
(9)判断是否满足Twt-K2≤ΔT，若满足则转入步骤(12)，若不满足则转入步骤(10)，本实例中温度容差ΔT为K2的0.1％；
(10)利用公式(5)得到修正后的热壁热流q′w，
q′wt＝qwt-Δqt*ΔT (5)
其中Δqt为当前时刻的热壁热流修正值，且Δqt＝qwt/K2；
(11)令qwt＝q′wt，转入步骤(5)；
(12)判断当前时刻是否到试验结束时间，若没到试验结束时间则将下一时刻确定为当前时刻，转入步骤(2)，若到试验结束时间，得到试验时间内的高温峰值处理后的加热温度曲线(如图2所示)和热壁热流曲线。
4、判断理论热壁热流曲线和高温峰值处理后的热壁热流曲线，是否满足|Qw-Q′w|≤ΔQw，本实例中总热量阈值ΔQw为理论总加热量的1％，满足|Qw-Q′w|≤ΔQw，因此，判定步骤3得到的高温峰值处理后的加热温度曲线和热壁热流曲线满足试验要求。
5、利用步骤4确认的达到试验要求的加热温度曲线和热壁热流曲线进行试验。
本发明未详细说明部分为本领域技术人员公知技术。
Claims (5)
1.一种热试验高温峰值的处理方法，其特征在于，包括以下步骤：
第一步，确定理论加热温度曲线和理论热壁热流曲线，
A1.1、利用公式(1)确定当前时刻t的热壁焓值Hwt，
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>H</mi>
<mrow>
<mi>r</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>w</mi>
<mi>t</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>H</mi>
<mrow>
<mi>r</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>o</mi>
<mi>l</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<msub>
<mi>q</mi>
<mrow>
<mi>c</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&epsiv;</mi>
<mi>&sigma;</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mn>4</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
<mn>4</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，qw(t-1)为前一时刻的热壁热流，Hrt为当前时刻的恢复焓，qct为当前时刻的冷壁热流，Hcool为冷壁焓值，Tw(t-1)为前一时刻的温度，Tref为参考温度，ε为热试验材料的灰度、σ为斯忒藩.波尔兹曼常数；试验时间t∈[t0,tJ],t0为试验开始时间，tJ为试验结束时间，试验开始时的取值为试验开始时的冷壁热流试验开始时的温度为试验环境温度；
A1.2、根据步骤A1.1得到的热壁焓值Hwt，利用公式(2)确定当前时刻的热壁热流qwt，
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>w</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>H</mi>
<mrow>
<mi>r</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>w</mi>
<mi>t</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>H</mi>
<mrow>
<mi>r</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>o</mi>
<mi>l</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<msub>
<mi>q</mi>
<mrow>
<mi>c</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&epsiv;</mi>
<mi>&sigma;</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mn>4</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
<mn>4</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
A1.3、判断是否满足|qwt-qw(t-1)|≤Δqw，若满足则转入步骤A1.7,若不满足，则转入步骤A1.4,其中Δqw为预设的热壁热流阈值；
A1.4、利用公式(3)确定当前时刻的温度Twt，
<mrow>
<msub>
<mi>q</mi>
<mi>w</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mi>&lambda;</mi>
<mfrac>
<mrow>
<msub>
<mi>dT</mi>
<mi>w</mi>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>x</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中λ是热试验材料的热导率，x是沿热试验材料法向的位移；
A1.5、根据步骤A1.4确定的当前时刻的温度Twt和步骤A1.1得到的热壁焓值Hwt，利用公式(4)确定当前时刻的热壁热流qwt，
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>w</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>H</mi>
<mrow>
<mi>r</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>w</mi>
<mi>t</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>H</mi>
<mrow>
<mi>r</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>o</mi>
<mi>l</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<msub>
<mi>q</mi>
<mrow>
<mi>c</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&epsiv;</mi>
<mi>&sigma;</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>w</mi>
<mi>t</mi>
</mrow>
<mn>4</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
<mn>4</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
A1.6、判断是否满足|qwt-qw(t-1)|≤Δqw，若满足则转入步骤A1.8,若不满足，则转入步骤A1.4；
A1.7、确定前一时刻的温度Tw(t-1)为当前时刻的温度Twt，转入步骤A1.9；
A1.8、确定此迭代步骤得到的温度为当前时刻的温度Twt，转入步骤A1.9；
A1.9、判断当前时刻是否到试验结束时间，若没到试验结束时间则将下一时刻确定为当前时刻，转入步骤A1.1，若到试验结束时间，得到试验时间内的理论加热温度曲线和理论热壁热流曲线；
第二步，判断是否进行高温峰值处理，若第一步确定的理论加热温度曲线的最高温度K1超过热试验最高加热温度K2则转入第三进行高温峰值处理，若没超过则认为第一步确定的理论加热温度曲线和理论热壁热流曲线符合试验要求，转入第五步；
第三步，高温峰值处理，得到高温峰值处理后的加热温度曲线和热壁热流曲线，
A3.1、确定削峰区间初值值T1和T2；
T1为削峰区间开始时间，即第一步得到的理论加热温度曲线中温度刚到K2的对应时间，T2为削峰区间结束时间，即第一步得到的理论加热温度曲线中温度降到K2后的某个时间，T2≥T′2，T′2为理论加热温度曲线中温度刚降低到小于K2的对应时间T′2；
A3.2利用公式(1)确定当前时刻t的热壁焓值Hwt，利用公式(2)确定当前时刻的热壁热流qwt；
A3.3、判断是否满足|qwt-qw(t-1)|≤Δqw，若满足则转入步骤A3.6,若不满足，则转入步骤A3.4；
A3.4、利用公式(3)确定当前时刻的温度Twt，再利用公式(4)确定当前时刻的热壁热流qwt；
A3.5、判断是否满足|qwt-qw(t-1)|≤Δqw，若满足则转入步骤A3.7,若不满足，则转入步骤A3.4；
A3.6、判断当前时刻是否在步骤A3.1确定的削峰区间内，若在则转入步骤A3.8,若不在则确定前一时刻的温度Tw(t-1)为当前时刻的温度Twt，转入步骤A3.12；
A3.7、确定此迭代步骤得到的温度为当前时刻的温度Twt，判断当前时刻是否在步骤A3.1确定的削峰区间内，若在则转入步骤A3.9,若不在则转入步骤A3.12；
A3.8、确定前一时刻的温度Tw(t-1)为当前时刻的温度Twt；
A3.9、判断是否满足Twt-K2≤ΔT，若满足则转入步骤A3.12，若不满足则转入步骤A3.10，ΔT为温度容差；
A3.10、利用公式(5)得到修正后的热壁热流q′w，
q′wt＝qwt-Δqt*ΔT (5)
其中Δqt为当前时刻的热壁热流修正值，且Δqt＝qwt/K2；
A3.11、令qwt＝q′wt，转入步骤A3.5；
A3.12、判断当前时刻是否到试验结束时间，若没到试验结束时间则将下一时刻确定为当前时刻，转入步骤A3.2，若到试验结束时间，得到试验时间内的高温峰值处理后的加热温度曲线和热壁热流曲线；
第四步，判断是否满足|Qw-Q′w|≤ΔQw，若满足，则认为第三步得到的高温峰值处理后的加热温度曲线和热壁热流曲线满足试验要求，转入第五步，若不满足，则调整第三步中的重新设置削峰区间结束时间T2，重复第三步直至满足|Qw-Q′w|≤ΔQw，其中Qw为第一步确定的理论热壁热流曲线得到的在试验时间内的总加热量，Q′w为第三步确定的高温峰值处理后的热壁热流曲线得到的在试验时间内的总加热量，ΔQw为总热量阈值；
第五步，利用第二步或第四步确认达到试验要求的加热温度曲线和热壁热流曲线进行试验。
2.根据权利要求1所述的一种热试验高温峰值的处理方法，其特征在于：所述步骤A3.9中温度容差取值为加热设备能加热的最高温度K2的0.05％～0.15％。
3.根据权利要求1所述的一种热试验高温峰值的处理方法，其特征在于：所述第四步中当不满足|Qw-Q′w|≤ΔQw时，若Qw＞Q′w则将削峰区间结束时间T2往回调即缩短削峰区间，若Qw<Q′w则将削峰区间结束时间T2往后调即延长削峰区间。
4.根据权利要求1或3所述的一种热试验高温峰值的处理方法，其特征在于：所述第四步中总热量阈值ΔQw的取值为理论总加热量Qw的0.1％～1％。
5.根据权利要求1所述的一种热试验高温峰值的处理方法，其特征在于：所述步骤A1.3热壁热流阈值Δqw取值为冷壁热流qc幅值的0.1％～1％。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611067880.XA CN108120612B (zh) | 2016-11-29 | 2016-11-29 | 一种热试验高温峰值的处理方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611067880.XA CN108120612B (zh) | 2016-11-29 | 2016-11-29 | 一种热试验高温峰值的处理方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108120612A true CN108120612A (zh) | 2018-06-05 |
| CN108120612B CN108120612B (zh) | 2019-11-15 |
Family
ID=62224269
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611067880.XA Active CN108120612B (zh) | 2016-11-29 | 2016-11-29 | 一种热试验高温峰值的处理方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108120612B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112937913A (zh) * | 2021-02-03 | 2021-06-11 | 中国空气动力研究与发展中心超高速空气动力研究所 | 电弧加热设备上中焓包罩试验状态自动调试方法和装置 |
| CN113252493A (zh) * | 2021-07-13 | 2021-08-13 | 中国飞机强度研究所 | 一种热强度试验系统控制方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8290246B1 (en) * | 2007-12-17 | 2012-10-16 | The United States of America as represented by the Administrator of the National Aeronautics & Space Administration (NASA) | Photogrammetric recession measurements of an ablating surface |
| CN103868750A (zh) * | 2014-03-20 | 2014-06-18 | 航天东方红卫星有限公司 | 适用于星上产品返修后的非对称性热试验方法 |
| CN106508022B (zh) * | 2010-12-31 | 2014-09-10 | 上海机电工程研究所 | 一种烧蚀式防热结构三维温度场分析计算方法 |
| CN104267062A (zh) * | 2014-10-22 | 2015-01-07 | 北京航空航天大学 | 一种气动热模拟试验中冷壁热流转换为热壁热流的方法 |
| CN104461677A (zh) * | 2014-10-30 | 2015-03-25 | 中国运载火箭技术研究院 | 一种基于cfd和fem技术的虚拟热试验方法 |
-
2016
- 2016-11-29 CN CN201611067880.XA patent/CN108120612B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8290246B1 (en) * | 2007-12-17 | 2012-10-16 | The United States of America as represented by the Administrator of the National Aeronautics & Space Administration (NASA) | Photogrammetric recession measurements of an ablating surface |
| CN106508022B (zh) * | 2010-12-31 | 2014-09-10 | 上海机电工程研究所 | 一种烧蚀式防热结构三维温度场分析计算方法 |
| CN103868750A (zh) * | 2014-03-20 | 2014-06-18 | 航天东方红卫星有限公司 | 适用于星上产品返修后的非对称性热试验方法 |
| CN104267062A (zh) * | 2014-10-22 | 2015-01-07 | 北京航空航天大学 | 一种气动热模拟试验中冷壁热流转换为热壁热流的方法 |
| CN104461677A (zh) * | 2014-10-30 | 2015-03-25 | 中国运载火箭技术研究院 | 一种基于cfd和fem技术的虚拟热试验方法 |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112937913A (zh) * | 2021-02-03 | 2021-06-11 | 中国空气动力研究与发展中心超高速空气动力研究所 | 电弧加热设备上中焓包罩试验状态自动调试方法和装置 |
| CN113252493A (zh) * | 2021-07-13 | 2021-08-13 | 中国飞机强度研究所 | 一种热强度试验系统控制方法 |
| CN113252493B (zh) * | 2021-07-13 | 2021-10-01 | 中国飞机强度研究所 | 一种热强度试验系统控制方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108120612B (zh) | 2019-11-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Li et al. | Constitutive equations for high temperature flow stress prediction of Al–14Cu–7Ce alloy | |
| CN103913652B (zh) | 电缆载流量的稳态迭代测定方法 | |
| Akasaka | New fundamental equations of state with a common functional form for 2, 3, 3, 3-tetrafluoropropene (R-1234yf) and trans-1, 3, 3, 3-tetrafluoropropene (R-1234ze (E)) | |
| CN109142083A (zh) | 一种变载荷历程下的蠕变损伤计算方法及模型 | |
| Whittaker et al. | An analysis of modern creep lifing methodologies in the titanium alloy Ti6-4 | |
| CN108120612A (zh) | 一种热试验高温峰值的处理方法 | |
| Liu et al. | A simple analytical crack tip opening displacement approximation under random variable loadings | |
| CN106446390A (zh) | 金属材料稳态蠕变速率拟合方程的计算方法 | |
| CN105468815B (zh) | 一种防雷产品的压敏电阻热致损毁评估方法及系统 | |
| CN105241903B (zh) | 一种预测700℃电站用γ′强化型高温合金持久强度的方法 | |
| CN106771703B (zh) | 霍尔推力器阳极热功率测量方法 | |
| CN109359330B (zh) | 一种描述高温材料蠕变变形三阶段的方法及模型 | |
| KR102341984B1 (ko) | 실리콘 웨이퍼의 휨량의 예측 방법 및 실리콘 웨이퍼의 제조 방법 | |
| CN115034042B (zh) | 一种变物性变压器油对流换热系数的修正方法 | |
| CN101795060B (zh) | 混合集成电路dc/dc变换器可靠性快速评价方法 | |
| JP6430726B2 (ja) | 膨張コンクリートの温度ひび割れ発生リスクの簡易評価方法、および膨張材の使用判定方法 | |
| CN104729909B (zh) | 一种预测金属材料高温松弛应力的方法 | |
| Mei et al. | Piecewise modeling of flow stress of 7075-T6 aluminum alloy in hot deformation | |
| CN103964686B (zh) | 一种用于等离子处理腔室的石英组件及等离子体处理设备 | |
| Sharma et al. | Thermo-mechanical buckling analysis of FGM plate using generalized plate theory | |
| Quan et al. | The flow behavior modeling of as-extruded 3Cr20Ni10W2 austenitic heat-resistant alloy at elevated temperatures considering the effect of strain | |
| CN107002182B (zh) | 热导率得到提高的690合金的有序合金的制备方法及由其制备而成的690合金的有序合金 | |
| KR102690458B1 (ko) | 형상기억합금의 온도측정 방법 및 온도분포의 예측 방법 | |
| Nasriddinov | Investigation of temperature sensors based on Si< P, Ni> | |
| CN105308205A (zh) | 提高了热导率的合金690有序合金的制造方法以及利用该方法来制造的合金690有序合金 |
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 |