发明名称 基于差量采样解的GappyPOD翼型反设计方法
摘要 本发明提出一种基于差量采样解的GappyPOD翼型反设计方法,基于特征正交分解方法中基模态的思想,利用GappyPOD降阶模型的数据重构方法,对原始方法中的采样解结构进行重新定义,提出一种差量采样解模式,由已知的翼型表面压力分布对缺失的目标翼型数据进行重构。在建立翼型过程中,采用ΔCST方法对初始翼型进行扰动取样,并用反设计得到翼型替换基础扰动翼型,反复迭代求解。设计结果证明本发明构建的翼型反设计方法的高效性与准确性,可对气动优化设计中的压力分布形态进行精细化的局部修形,通过本发明的反设计方法可得到更精细更符合实际工程要求的飞行器气动外形。
申请公布号 CN104699901A 申请公布日期 2015.06.10
申请号 CN201510106464.5 申请日期 2015.03.11
申请人 西北工业大学 发明人 高正红;夏露;左英桃;李静;赵柯;王超;赵欢
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 陈星
主权项 一种基于差量采样解的GappyPOD翼型反设计方法,其特征在于:包括以下步骤:步骤1:根据给定的翼型目标压力分布<img file="FDA00006803963200000114.GIF" wi="84" he="81" />选取初始翼型,并以初始翼型进入步骤2;步骤2:对于进行本步骤的翼型,采用参数化方法对翼型设计变量进行单独扰动,得到一组翼型;步骤3:对步骤2得到的一组翼型分别进行流场计算,获得每个翼型的采样解矩阵U<sup>i</sup>,i=1,2,…,n,n为步骤2得到的一组翼型的翼型个数,<img file="FDA0000680396320000011.GIF" wi="295" he="80" />y<sup>i</sup>为第i个翼型的坐标向量,<img file="FDA0000680396320000012.GIF" wi="67" he="85" />为流场计算得到的第i个翼型的压力分布;步骤4:根据步骤3得到的n个翼型的采样解矩阵U<sup>i</sup>,进行差量采样解分析过程:取<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>&Delta;</mi><msup><mi>U</mi><mi>i</mi></msup><mo>=</mo><msup><mi>U</mi><mi>i</mi></msup><mo>-</mo><mover><mi>U</mi><mo>&OverBar;</mo></mover><mo>,</mo></mrow>]]></math><img file="FDA0000680396320000013.GIF" wi="310" he="70" /></maths>其中<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><mi>U</mi><mo>&OverBar;</mo></mover><mo>=</mo><mo>[</mo><mover><mi>y</mi><mo>&OverBar;</mo></mover><mo>,</mo><msub><mover><mi>C</mi><mo>&OverBar;</mo></mover><mi>p</mi></msub><mo>]</mo><mo>,</mo><mover><mi>y</mi><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mi>y</mi><mi>i</mi></msup></mrow><mi>n</mi></mfrac><mo>,</mo><msub><mover><mi>C</mi><mo>&OverBar;</mo></mover><mi>p</mi></msub><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>C</mi><mi>p</mi><mi>i</mi></msubsup></mrow><mi>n</mi></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000680396320000014.GIF" wi="823" he="197" /></maths>其中ΔU<sup>i</sup>为第i个翼型的差量采样解矩阵,n个差量采样解矩阵组成一组差量采样解矩阵ΔU;取<img file="FDA0000680396320000015.GIF" wi="445" he="91" />为差量目标压力分布;步骤5:采用特征正交分解方法获取ΔU的一组差量基向量<img file="FDA0000680396320000016.GIF" wi="189" he="84" />其中m为差量基向量个数,第j个差量基向量<img file="FDA0000680396320000017.GIF" wi="438" he="85" />为差量基向量中表示翼型坐标的向量,<img file="FDA0000680396320000018.GIF" wi="82" he="87" />为差量基向量中表示翼型压力分布的向量;步骤6:采用步骤5得到的一组差量基向量<img file="FDA0000680396320000019.GIF" wi="156" he="84" />中的表示翼型压力分布的向量<img file="FDA00006803963200000110.GIF" wi="84" he="80" />拟合差量目标压力分布<img file="FDA00006803963200000111.GIF" wi="128" he="85" />得到拟合系数;根据拟合系数以及<img file="FDA00006803963200000112.GIF" wi="154" he="85" />中的表示翼型坐标的向量得到翼型修正量;步骤7:步骤6得到的翼型修正量加上y得到反设计翼型;对得到的反设计翼型进行流场计算,得到反设计翼型的压力分布;计算反设计翼型的压力分布与目标压力分布<img file="FDA00006803963200000113.GIF" wi="71" he="85" />的差值,若差值满足设计精度要求,则方法结束;否则以得到的反设计翼型进入步骤2。
地址 710072 陕西省西安市友谊西路127号