发明名称 一种考虑加工误差的高超声速机翼鲁棒优化设计方法
摘要 本发明公开了一种考虑加工误差的高超声速机翼鲁棒优化设计方法,属于优化设计技术领域。充分考虑高超声速机翼设计中存在的几何加工误差,在加工误差系数概率密度函数未知的情况下,利用区间向量实现加工误差系数的定量化表征。通过参数化方法建立高超声速机翼气动外形,进行非结构表面网格划分。将区间参数顶点法与气动力/热工程算法相结合,计算机翼驻点总加热量及升阻比的区间上下界。在此基础上,建立多目标区间鲁棒优化模型,应用遗传算法对机翼外形进行优化设计。数值结果表明,本发明方法在保持机翼升阻比约束的前提条件下,降低了所设计机翼的驻点总加热量,同时减小了驻点总加热量的波动范围,为高超声速机翼外形设计提供了新思路。
申请公布号 CN106126860A 申请公布日期 2016.11.16
申请号 CN201610554691.9 申请日期 2016.07.14
申请人 北京航空航天大学 发明人 邱志平;郑宇宁;王晓军;王磊;王睿星;许孟辉;李云龙;仇翯辰;姜南
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 北京科迪生专利代理有限责任公司 11251 代理人 杨学明;顾炜
主权项 一种考虑加工误差的高超声速机翼鲁棒优化设计方法,其特征在于实现步骤如下:步骤(1)、首先,确定高超声速机翼的平面轮廓参数,包括翼根弦长C<sub>r</sub>、副翼宽度C<sub>a</sub>、机翼展弦比AR及稍根比TR;步骤(2)、利用类函数/形函数转换方法建立高超声速机翼翼型的参数化表达式,翼型的几何曲线可用下列函数表示:<maths num="0001"><math><![CDATA[<mrow><mfrac><mi>y</mi><mi>c</mi></mfrac><mrow><mo>(</mo><mfrac><mi>x</mi><mi>c</mi></mfrac><mo>)</mo></mrow><mo>=</mo><mi>C</mi><mrow><mo>(</mo><mfrac><mi>x</mi><mi>c</mi></mfrac><mo>)</mo></mrow><mi>S</mi><mrow><mo>(</mo><mfrac><mi>x</mi><mi>c</mi></mfrac><mo>)</mo></mrow><mo>+</mo><mfrac><mi>x</mi><mi>c</mi></mfrac><mfrac><msub><mi>z</mi><mrow><mi>t</mi><mi>e</mi></mrow></msub><mi>c</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001048911450000011.GIF" wi="1349" he="135" /></maths>式中,x/c为翼型弦向的无量纲坐标值,y/c为翼型法向的无量纲坐标值,C(x/c)和S(x/c)分别为类型函数和形状函数,z<sub>te</sub>/c为翼型后缘点的无量纲坐标值,C(x/c)可表示为:<maths num="0002"><math><![CDATA[<mrow><mi>C</mi><mrow><mo>(</mo><mfrac><mi>x</mi><mi>c</mi></mfrac><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><mfrac><mi>x</mi><mi>c</mi></mfrac><mo>)</mo></mrow><msub><mi>N</mi><mn>1</mn></msub></msup><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><mi>x</mi><mi>c</mi></mfrac><mo>)</mo></mrow><msub><mi>N</mi><mn>2</mn></msub></msup><mo>,</mo><mn>0</mn><mo>&le;</mo><mfrac><mi>x</mi><mi>c</mi></mfrac><mo>&le;</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001048911450000012.GIF" wi="1415" he="144" /></maths>对于Clark‑Ys翼型,取指数N<sub>1</sub>=0.5,N<sub>2</sub>=1,S(x/c)可表示为:<maths num="0003"><math><![CDATA[<mrow><mi>S</mi><mrow><mo>(</mo><mfrac><mi>x</mi><mi>c</mi></mfrac><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mi>n</mi></munderover><mrow><mo>&lsqb;</mo><mrow><msub><mi>&lambda;</mi><mi>i</mi></msub><mo>&CenterDot;</mo><mfrac><mrow><mi>n</mi><mo>!</mo></mrow><mrow><mi>i</mi><mo>!</mo><mrow><mo>(</mo><mrow><mi>n</mi><mo>-</mo><mi>i</mi></mrow><mo>)</mo></mrow><mo>!</mo></mrow></mfrac><mo>&CenterDot;</mo><msup><mrow><mo>(</mo><mfrac><mi>x</mi><mi>c</mi></mfrac><mo>)</mo></mrow><mi>i</mi></msup><mo>&CenterDot;</mo><msup><mrow><mo>(</mo><mrow><mn>1</mn><mo>-</mo><mfrac><mi>x</mi><mi>c</mi></mfrac></mrow><mo>)</mo></mrow><mrow><mi>n</mi><mo>-</mo><mi>i</mi></mrow></msup></mrow><mo>&rsqb;</mo></mrow><mo>,</mo><mn>0</mn><mo>&le;</mo><mfrac><mi>x</mi><mi>c</mi></mfrac><mo>&le;</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001048911450000013.GIF" wi="1606" he="165" /></maths>S(x/c)与翼型前缘半径R<sub>le</sub>/c及后缘倾角β满足以下关系:<maths num="0004"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>S</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo>=</mo><msub><mi>&lambda;</mi><mn>0</mn></msub><mo>=</mo><msqrt><mfrac><mrow><mn>2</mn><msub><mi>R</mi><mrow><mi>l</mi><mi>e</mi></mrow></msub></mrow><mi>c</mi></mfrac></msqrt><mo>,</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>S</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msub><mi>&lambda;</mi><mi>n</mi></msub><mo>=</mo><mi>t</mi><mi>a</mi><mi>n</mi><mi>&beta;</mi><mo>+</mo><mfrac><msub><mi>z</mi><mrow><mi>t</mi><mi>e</mi></mrow></msub><mi>c</mi></mfrac></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001048911450000014.GIF" wi="1300" he="279" /></maths>在形状函数S(x/c)中取n=3,这样共包含<img file="FDA0001048911450000015.GIF" wi="583" he="63" />等8个几何设计参数,其中,<img file="FDA0001048911450000016.GIF" wi="234" he="55" />表示翼型前缘半径,β<sub>1</sub>和β<sub>2</sub>表示上、下翼面的后缘倾角,<img file="FDA0001048911450000017.GIF" wi="219" he="63" />为翼型后缘点的无量纲坐标值,λ<sub>1</sub>、λ<sub>2</sub>、λ<sub>1</sub>'、λ<sub>2</sub>'为上、下翼面形状函数多项式的加权系数;步骤(3)、将机翼展弦比AR、稍根比TR及翼型前缘半径<img file="FDA0001048911450000018.GIF" wi="58" he="55" />作为优化变量,记为:<img file="FDA0001048911450000019.GIF" wi="325" he="55" />其余设计参数均视为常量,给出设计变量的初始值及取值范围;步骤(4)、考虑机翼加工制造过程中存在的表面加工误差,引入表面加工误差系数<img file="FDA00010489114500000110.GIF" wi="59" he="47" />并利用区间向量对<img file="FDA00010489114500000111.GIF" wi="33" he="46" />进行定量化,记为<img file="FDA00010489114500000112.GIF" wi="257" he="63" />步骤(5)、分别在<img file="FDA00010489114500000113.GIF" wi="106" he="63" />和<img file="FDA00010489114500000114.GIF" wi="114" he="60" />处得到翼型的参数化表达式,表示如下:<img file="FDA0001048911450000021.GIF" wi="1462" he="127" /><img file="FDA0001048911450000022.GIF" wi="1470" he="119" />步骤(6)、根据式(5)和(6)中翼型的参数化表达式及步骤(3)中机翼平面轮廓参数,建立高超声速机翼气动外形的几何模型;步骤(7)、对于步骤(6)中生成的气动外形,对机翼表面进行非结构网格的自由划分;步骤(8)、根据步骤(7)中生成的机翼网格文件,在给定飞行工况下,利用高超声速工程算法计算机翼驻点热流密度Q及机翼升阻比K;步骤(9)、根据不同时刻得到的机翼驻点热流密度值,利用驻点总加热量工程计算方法得到所给机翼的总加热量;步骤(10)、利用区间参数顶点法得到机翼驻点总加热量及升阻比的区间上下界,然后基于区间数学理论计算机翼驻点总加热量及升阻比的区间中心值和区间半径;步骤(11)、在保持机翼升阻比约束条件下,以机翼驻点总加热量的区间中心值和半径最小化为优化目标,建立多目标区间鲁棒优化模型;步骤(12)、利用遗传算法,在设计变量取值范围内对高超声速机翼进行鲁棒优化设计;步骤(13)、判断优化目标是否满足收敛条件,若不满足,转到步骤(3),更新设计变量,重复步骤(4)~(12);步骤(14)、直至设计目标相邻两次迭代值的变化小于设定容许偏差时,完成高超声速机翼的鲁棒优化设计;步骤(15)、将通过鲁棒优化得到的高超声速机翼与原始机翼进行对比,比较两种机翼的驻点总加热量。
地址 100191 北京市海淀区学院路37号