[發明專利]一種特大變形動態網格生成方法

有效

• 申請號: 201510041391.6
• 申請日: 2015-01-27
• 公開（公告）號: CN104572575B
• 公開（公告）日: 2017-10-13
• 發明（設計）人: 郭同慶;陳皓;陸志良
• 申請（專利權）人: 南京航空航天大學
• 主分類號: G06F17/00
• 分類號: G06F17/00
• 代理公司: 南京鐘山專利代理有限公司32252
• 代理人: 戴朝榮

1.一種特大變形動態網格生成方法，其特征在于，包括以下步驟：

步驟一、基于多塊對接結構網格，選擇不同時刻下的變形位置，生成相應的拓撲結構相同的靜態網格；

步驟二、根據步驟一中生成的靜態網格，預估計算出任意時刻下的計算網格；

步驟三、對步驟二中生產的計算網格進行校正；

步驟四、將步驟三中得到的校正網格用于非定常計算，得到隨時間變化的全流場動態網格；

在所述步驟一中，設定選擇n個時刻下的變形位置，將其中第t個時刻下網格點坐標記為：

fx(t,b,i,j,k),fy(t,b,i,j,k),fz(t,b,i,j,k)；

其中，1≤t≤n，b表示網格點所在的網格塊號，i,j,k表示網格點在網格塊三個方向的序號；

所述步驟二中，通過拉格朗日插值算法，預估計算任意時刻下的計算網格，預估計算網格中的網格點坐標表示如下：

fx(n+1,b,i,j,k)＝a₁fx(1,b,i,j,k)+a₂fx(2,b,i,j,k)+…+a_tfx(t,b,i,j,k)+…+a_nfx(n,b,i,j,k)；

fy(n+1,b,i,j,k)＝a₁fy(1,b,i,j,k)+a₂fy(2,b,i,j,k)+…+a_tfy(t,b,i,j,k)+…+a_nfy(n,b,i,j,k)；

fz(n+1,b,i,j,k)＝a₁fz(1,b,i,j,k)+a₂fz(2,b,i,j,k)+…+a_tfz(t,b,i,j,k)+…+a_nfz(n,b,i,j,k)；

其中，a₁,a₂,…a_t,…a_n為拉格朗日插值基函數，且

所述步驟三中，通過基于擾動衰減規律的彈性變形技術對預估計算網格進行校正，校正后的網格點坐標表示如下：

其中，表示預估網格點坐標，表示靜態網格點隨物面邊界作剛性運動的瞬態坐標值，g表示網格點序號的函數；

g=max((i-iwif-iw)2,(j-jwjf-jw)2,(k-kwkf-kw)2);]]>

其中，i_w,j_w,k_w表示物面網格點序號，i_f,j_f,k_f表示相應的遠邊界網格點序號。