主权项 |
一种基于差量采样解的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>Δ</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>‾</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>‾</mo></mover><mo>=</mo><mo>[</mo><mover><mi>y</mi><mo>‾</mo></mover><mo>,</mo><msub><mover><mi>C</mi><mo>‾</mo></mover><mi>p</mi></msub><mo>]</mo><mo>,</mo><mover><mi>y</mi><mo>‾</mo></mover><mo>=</mo><mfrac><mrow><munderover><mi>Σ</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>‾</mo></mover><mi>p</mi></msub><mo>=</mo><mfrac><mrow><munderover><mi>Σ</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。 |