一种低阻力非烧蚀高超声速折叠舵计算方法及结构

专利主权项内容

1.一种低阻力非烧蚀高超声速折叠舵计算方法，其特征在于：包括如下步骤：Step1：采用简化的后掠圆柱驻点热流计算舵面前缘热流；根据后掠圆柱驻点热流计算公式得到舵面前缘热流，对前缘曲率半径和前缘后掠角参数有：其中：
为等半径的球形端头驻点热流；q为舵面前缘热流，单位为kW/m；s2r为前缘曲率半径，单位为m；cμ为自由来流处粘性系数，单位为kg/m.s，由Sutherland公式计算；∞
其中T为来流温度，根据弹道热环境直接获得；∞λ为前缘后掠角；
为湍流放大因子和层流加热放大因子之比；T为舵表面温度，T默认值为300K；wwh为滞止焓，h为壁面气体焓值，h为室温300K下的焓值，因此h＝h，此时sW300W300以空气密度ρ、飞行速度v、前缘曲率半径r、前缘后掠角λ为输入参数，得到等半径的球形端头驻点热流q；按照简化的后掠圆柱驻点热流计算公式得到前缘热流q；∞∞cds前缘曲率半径的一组取值为：前缘后掠角的一组取值为：ω，ω+5°、ω+10°，ω+15°；根据3个前缘曲率半径和4个前缘后掠角的输入组合，得到12组热流条件；Step2：根据Step1获得的12组热流条件，选取折叠舵前缘展向中截面作为传热计算简化模型，将前缘三维热响应温度计算简化为二维瞬态传热计算：边界条件为前缘外表面的热量交换，热量交换包括气动加热、热辐射以及向前缘结构内部的热传导：除与气流存在热量交换的表面外均按理想绝热壁计算：式中，ρ为材料密度，C为材料定压比热容，K为材料热传导系数，T为前缘驻点表温，x为空气舵展向坐标，y为空气舵厚度方向坐标，q气动热环境计算提供的随飞行时间变化的冷壁热流密度，h为壁面气体焓值，h为气流恢复焓，ε为涂层材料热辐射系数，σ＝5.6704×10W/mK为Stephan-Boltzmann常数，T为舵表面温度；计算获得导弹飞行全程折叠舵前缘驻点表温；pwWr-824wStep3：将Step2计算得到的多组前缘驻点表温与前缘材料的许用温度进行对比分析，若计算得到的多组前缘驻点表温均小于前缘材料的许用温度，则选择最小的前缘曲率半径作为最终设计结果，本次迭代结束；若计算得到的多组前缘驻点表温部分结果小于前缘材料的许用温度，计算得到的多组前缘驻点表温部分结果大于前缘材料的许用温度，则在前缘驻点表温部分结果小于前缘材料的许用温度的结果中选择最小的前缘曲率半径作为最终设计结果，本次迭代结束；若计算得到的前缘驻点表温均高于或等于前缘材料许用温度，则改变前缘曲率半径前缘曲率半径/>增加1mm，跳转至Step1，将改变后的前缘曲率半径/>值按照Step1、Step2重新计算分析，直到前缘驻点表温满足前缘材料的许用温度。