[發(fā)明專利]自適應(yīng)步長的順序耦合分析方法

有效

• 申請?zhí)枺?/td>: 201010602510.8
• 申請日: 2010-12-23
• 公開（公告）號: CN102542086A
• 公開（公告）日: 2012-07-04
• 發(fā)明（設(shè)計）人: 譚建榮;張樹有;徐敬華;趙振
• 申請（專利權(quán)）人: 浙江大學(xué)
• 主分類號: G06F17/50
• 分類號: G06F17/50
• 代理公司: 杭州求是專利事務(wù)所有限公司 33200
• 代理人: 陳昱彤

1.一種自適應(yīng)步長的順序耦合分析方法，其特征在于，包括如下步驟：

步驟（1）：建立耦合場分析模型，設(shè)定一個初始迭代步長作為第一個和第二個子步的迭代步長，對瞬態(tài)場耦合系統(tǒng)的各耦合狀態(tài)變量進(jìn)行第一個和第二個子步的順序耦合分析，得到各耦合狀態(tài)變量的第一個和第二個子步的迭代節(jié)點(diǎn)的耦合狀態(tài)變量值；

步驟（2）：提取每個耦合狀態(tài)變量的當(dāng)前迭代節(jié)點(diǎn)及其前兩個迭代節(jié)點(diǎn)的耦合狀態(tài)變量值，利用公式（1）得到各耦合狀態(tài)變量曲線的當(dāng)前迭代節(jié)點(diǎn)的斜率：

???????（1）

式（1）中，i表示迭代節(jié)點(diǎn)的序號且i≥3，j表示耦合狀態(tài)變量的編號，、、分別為各個耦合狀態(tài)變量在當(dāng)前迭代節(jié)點(diǎn)及其前兩個迭代節(jié)點(diǎn)的耦合狀態(tài)變量值，、分別為各個耦合狀態(tài)變量當(dāng)前子步及其前一子步的迭代步長，為各個耦合狀態(tài)變量的狀態(tài)耦合變量曲線在當(dāng)前迭代節(jié)點(diǎn)的斜率；

????步驟（3）：設(shè)定每個耦合狀態(tài)變量的步長調(diào)整區(qū)間，按以下方法得到各個耦合狀態(tài)變量的下一子步的迭代步長：

若步驟（2）所述斜率的絕對值大于對應(yīng)耦合狀態(tài)變量的步長調(diào)整區(qū)間的設(shè)定值上限，則使每個耦合狀態(tài)變量的下一子步的迭代步長與當(dāng)前子步的迭代步長滿足如下關(guān)系式（2）：

??????????????????????????????????????????????（2）

式（2）中，、分別為當(dāng)前子步與下一子步的迭代步長，為步長調(diào)整系數(shù)且；

若所述斜率的絕對值小于對應(yīng)耦合狀態(tài)變量的步長調(diào)整區(qū)間的設(shè)定值下限，則使每個耦合狀態(tài)變量的下一子步的迭代步長與當(dāng)前子步的迭代步長滿足如下關(guān)系式（3）：

?????????????????????????????????????（3）

式（3）中，為步長調(diào)整系數(shù)且；

若所述斜率的絕對值屬于對應(yīng)耦合狀態(tài)變量的步長調(diào)整區(qū)間，則使每個耦合狀態(tài)變量的下一子步的迭代步長與當(dāng)前子步的迭代步長滿足如下關(guān)系式（4）：

???????????????????????????????????????????????（4）

步驟（4）：將步驟（3）所獲得的所有耦合狀態(tài)變量的下一子步的迭代步長的最小值作為瞬態(tài)場耦合系統(tǒng)的下一子步的系統(tǒng)迭代步長，并采用所述系統(tǒng)迭代步長對瞬態(tài)場耦合系統(tǒng)的各耦合狀態(tài)變量進(jìn)行下一子步的順序耦合分析，得到各耦合狀態(tài)變量的下一子步的迭代節(jié)點(diǎn)的耦合狀態(tài)變量值；

步驟（5）：返回執(zhí)行步驟（2），直至達(dá)到所述瞬態(tài)場耦合系統(tǒng)的系統(tǒng)分析終止時間。

?