发明名称 一种探空火箭气动外形优化目标函数的确定方法
摘要 为克服探空火箭气动外形设计优化过程现有目标函数确定方法在优化效率与优化结果性能上的不足,本发明提供一种探空火箭气动外形优化目标函数的确定方法。本发明以气动阻力沿弹道所致动量损失最小为出发点,提出了包含以下步骤的探空火箭气动外形优化目标函数确定方法:(1)给出火箭上升段初始弹道数据;(2)初始弹道数据预处理,确定弹道特征点;(3)计算各特征点加权系数;(4)以弹道特征点加权阻力系数最小作为探空火箭气动外形优化目标函数。本发明显著提高了探空火箭气动外形优化的效率、改善了优化结果的性能,适用于各类探空火箭、各个设计阶段的外形设计,且可推广应用于其它飞行器气动外形优化的目标函数确定问题。
申请公布号 CN105160108A 申请公布日期 2015.12.16
申请号 CN201510574018.7 申请日期 2015.09.10
申请人 中国人民解放军国防科学技术大学 发明人 江振宇;彭科;张士峰;胡凡;向敏
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 北京中济纬天专利代理有限公司 11429 代理人 胡伟华;张小雪
主权项 一种探空火箭气动外形优化目标函数的确定方法,其特征在于,包括以下步骤:第一步:给出火箭上升段初始弹道数据,火箭上升段初始弹道数据包括时间、高度、马赫数、攻角、动压参数;第二步:火箭上升段初始弹道数据预处理,剔除弹道飞行后期动压较小的弹道数据,确定弹道特征点;探空火箭的弹道动压q在火箭上升段[0,t<sub>s</sub>]后期弹道动压数值相对较小,动压数值较小的飞行段火箭动量损失也较小,在确定探空火箭气动外形优化目标函数过程中可不予考虑;设时刻t<sub>c</sub>满足下式:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mo>&Integral;</mo><mn>0</mn><msub><mi>t</mi><mi>c</mi></msub></msubsup><mrow><mi>q</mi><mi>d</mi><mi>t</mi></mrow><mo>/</mo><msubsup><mo>&Integral;</mo><mn>0</mn><msub><mi>t</mi><mi>s</mi></msub></msubsup><mi>q</mi><mi>d</mi><mi>t</mi><mo>=</mo><mi>C</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000799624350000011.GIF" wi="1065" he="115" /></maths>式中,C取0.99~0.999间的常数;时刻t<sub>c</sub>后的弹道动压数值相对较小,对应火箭动量损失也较小,因此在确定探空火箭气动外形优化目标函数时仅考虑[0,t<sub>c</sub>]时段弹道数据;将该时段等分为份,得分段时刻t<sub>0</sub>,t<sub>1</sub>,t<sub>2</sub>…,t<sub>n‑2</sub>,t<sub>n‑1</sub>,t<sub>n</sub>,时刻(t<sub>i‑1</sub>+t<sub>i</sub>)/2对应的弹道[0,t<sub>c</sub>]点即为弹道特n征点,其中<sub>i</sub>为1,2,3…,n;第三步:由步骤二得到的[0,t<sub>c</sub>]时段弹道时间、动压参数计算各弹道特征点加权系数w<sub>i</sub>;各弹道特征点加权系数根据下式计算:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>w</mi><mi>i</mi></msub><mo>=</mo><msubsup><mo>&Integral;</mo><msub><mi>t</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><msub><mi>t</mi><mi>i</mi></msub></msubsup><mi>q</mi><mi>d</mi><mi>t</mi><mo>/</mo><msubsup><mo>&Integral;</mo><mn>0</mn><msub><mi>t</mi><mi>c</mi></msub></msubsup><mi>q</mi><mi>d</mi><mi>t</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000799624350000012.GIF" wi="1069" he="127" /></maths>式中,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msub><mi>w</mi><mi>i</mi></msub><mo>=</mo><mn>1</mn><mo>;</mo></mrow>]]></math><img file="FDA0000799624350000013.GIF" wi="254" he="93" /></maths>步骤四:取C<sub>xi</sub>为步骤三中各弹道特征点的火箭阻力系数,得目标函数为<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>min</mi><mi> </mi><mi>f</mi><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msub><mi>w</mi><mi>i</mi></msub><msub><mi>C</mi><mrow><mi>x</mi><mi>i</mi></mrow></msub><mo>.</mo></mrow>]]></math><img file="FDA0000799624350000014.GIF" wi="499" he="99" /></maths>
地址 410073 湖南省长沙市开福区德雅路109号