[發(fā)明專(zhuān)利]機(jī)翼相關(guān)變量組和不相關(guān)變量組一階敏感系數(shù)估計(jì)方法

在審

• 申請(qǐng)?zhí)枺?/td>: 202111298682.5
• 申請(qǐng)日: 2021-11-04
• 公開(kāi)（公告）號(hào): CN114065382A
• 公開(kāi)（公告）日: 2022-02-18
• 發(fā)明（設(shè)計(jì)）人: 羅中峰;李洪強(qiáng)
• 申請(qǐng)（專(zhuān)利權(quán)）人: 揚(yáng)州工業(yè)職業(yè)技術(shù)學(xué)院
• 主分類(lèi)號(hào): G06F30/15
• 分類(lèi)號(hào): G06F30/15;G06F30/20;G06F17/16;G06F119/14
• 代理公司: 南京蘇高專(zhuān)利商標(biāo)事務(wù)所(普通合伙) 32204
• 代理人: 柏尚春

1.一種機(jī)翼相關(guān)變量組和不相關(guān)變量組一階敏感系數(shù)估計(jì)方法，其特征在于：包括以下步驟：

(1)建立飛機(jī)機(jī)翼最大應(yīng)力的計(jì)算模型，并以此作為目標(biāo)函數(shù)；

(2)構(gòu)建飛機(jī)機(jī)翼最大應(yīng)力計(jì)算模型中相關(guān)變量的轉(zhuǎn)換矩陣；

(3)根據(jù)轉(zhuǎn)換矩陣，利用服從相關(guān)變量邊緣分布的獨(dú)立隨機(jī)變量，表示飛機(jī)機(jī)翼最大應(yīng)力計(jì)算模型中相關(guān)變量；

(4)根據(jù)獨(dú)立隨機(jī)變量組的一階敏感系數(shù)，獲得飛機(jī)機(jī)翼最大應(yīng)力計(jì)算模型中相關(guān)變量組的一階敏感系數(shù)；

(5)獲得飛機(jī)機(jī)翼最大應(yīng)力計(jì)算模型中不相關(guān)變量組的一階敏感系數(shù)；

(6)根據(jù)獲得的一階敏感系數(shù)，修正飛機(jī)機(jī)翼參數(shù)的分布，降低機(jī)翼最大應(yīng)力的波動(dòng)范圍。

2.根據(jù)權(quán)利要求1所述的機(jī)翼相關(guān)變量組和不相關(guān)變量組一階敏感系數(shù)的估計(jì)方法，其特征在于：步驟(1)中所述目標(biāo)函數(shù)為：

Z＝F(x₁,x₂,…,x_s,x_s+1,…,x_n)

其中：Z為飛機(jī)機(jī)翼的最大應(yīng)力；輸入變量x₁,x₂,…,x_s,x_s+1,…,x_n～N(μ,Σ)；x₁,x₂,…,x_s是飛機(jī)機(jī)翼最大應(yīng)力計(jì)算模型中相關(guān)變量；x_s+1,…,x_n是飛機(jī)機(jī)翼最大應(yīng)力計(jì)算模型中不相關(guān)變量；目標(biāo)函數(shù)擁有有限方差。

3.根據(jù)權(quán)利要求2所述的機(jī)翼相關(guān)變量組和不相關(guān)變量組一階敏感系數(shù)的估計(jì)方法，其特征在于：步驟(2)中包括以下步驟：

(21)獲得輸入變量間的相關(guān)系數(shù)矩陣

其中，ρ_i,j是變量x_i和x_j的相關(guān)系數(shù)，ρ_i,j＝ρ_j,i，ρ_i,i＝1，i＝1,2,…,n；j＝1,2,…,n；

(22)獲得簡(jiǎn)化的相關(guān)系數(shù)矩陣C_ρ-s，所述簡(jiǎn)化的相關(guān)系數(shù)矩陣C_ρ-s只包含所有相關(guān)變量間的相關(guān)系數(shù)

(23)根據(jù)簡(jiǎn)化的相關(guān)系數(shù)矩陣C_ρ-s，利用Cholesky Factorization方法，獲得轉(zhuǎn)換矩陣

其中，d_k,k＝1。

4.根據(jù)權(quán)利要求3所述的機(jī)翼相關(guān)變量組和不相關(guān)變量組一階敏感系數(shù)的估計(jì)方法，其特征在于：步驟(3)中包括以下步驟：

(31)根據(jù)轉(zhuǎn)換矩陣L_ρ-s和相關(guān)系數(shù)的性質(zhì)，利用服從相關(guān)變量邊緣分布的獨(dú)立隨機(jī)變量表示飛機(jī)機(jī)翼最大應(yīng)力計(jì)算模型中相關(guān)變量

其中，變量y_k是服從相關(guān)變量x_k邊緣分布的獨(dú)立隨機(jī)變量，k＝1,2,…,s；μ_k是變量y_k的均值；σ_k是變量y_k的均方差；

(32)將步驟(31)所示公式代入步驟(1)所示目標(biāo)函數(shù)，目標(biāo)函數(shù)轉(zhuǎn)換為

Z＝F_k(y₁,y₂,y₃,...,y_k-1,y_k,y_k+1,...,y_s,x_s+1,...,x_n)。