[發(fā)明專利]一種兩自由度機翼的顫振邊界預測方法

在審

• 申請?zhí)枺?/td>: 201510320937.1
• 申請日: 2015-06-11
• 公開（公告）號: CN104881585A
• 公開（公告）日: 2015-09-02
• 發(fā)明（設計）人: 周麗;楊秉才;李揚
• 申請（專利權）人: 南京航空航天大學
• 主分類號: G06F19/00
• 分類號: G06F19/00
• 代理公司: 南京瑞弘專利商標事務所(普通合伙) 32249
• 代理人: 楊曉玲

1.一種兩自由度機翼的顫振邊界預測方法，其特征在于，該方法包括以下步驟：

1)設定待測風速為v_i,測得風速v_i下機翼的加速度響應信號，用帶通濾波器濾掉低頻噪聲和高頻模態(tài)，保留信號中的顫振耦合兩個模態(tài)；

2)用遞推最小二乘法計算出步驟1)所得到的每個風速下機翼響應信號ARMA模型的AR項系數(shù)，將AR項系數(shù)設為{α₀,α₁,α₂,α₃,α₄}；

3)用步驟2)給出的AR項系數(shù)構造穩(wěn)定性參數(shù)，根據(jù)穩(wěn)定性準則進行判定：

G(1)＝1+α₁+α₂+α₃+α₄＞0

G(-1)＝1-α₁+α₂-α₃+α₄＞0

F⁺(1)＝1+α₄＞0,F^-(1)＝1+α₄＞0

F⁺(3)＝det(X+Y)＞0

F^-(3)＝det(X-Y)＞0

其中，

X=α0α1α20α0α100α0,Y=α2α3α4α3α40α400;]]>

若在風速v_i下能滿足以上所有方程式，則判定系統(tǒng)是穩(wěn)定的，進行步驟4)；如不穩(wěn)定，則說明發(fā)生顫振；

4)用步驟2)給出的AR項系數(shù)構造顫振邊界判據(jù)，定義顫振邊界判據(jù)F，計算公式如下：

F=det(X-Y)(1-α4)2]]>

5)計算不同風速v_i下的顫振判據(jù)F，用每個風速下得到的顫振判據(jù)F與風速v_i進行二次函數(shù)擬合，通過曲線外推，當顫振判據(jù)F下降為零時的風速為顫振邊界。