摘要：為實(shí)現艦船縱搖和升沉運動(dòng)的解耦，基于保結構同譜流算法，提出一種解耦變換的尋找方法，將尋找解耦變換的非線(xiàn)性問(wèn)題轉化為Sylvester方程的求解問(wèn)題，并利用矩陣Kronecker積的相關(guān)知識快速找到解耦變換.基于水池實(shí)驗獲得的縱向運動(dòng)數據進(jìn)行的數值實(shí)驗仿真結果表明該方法確實(shí)可行.

關(guān)鍵詞：艦船縱向運動(dòng)；保結構同譜流；水動(dòng)力系數；二階微分系統

中圖分類(lèi)號：U661.32 文獻標志碼：A

Numerical decoupling of ship vertical motion system

WANG Shu-juana，SHEN Ji-honga，LI Ji-deb

（a.College of Science;b.College of Shipbuilding Engineering，Harbin Engineering University，Harbin 150001，China)

Abstract:To realize the ship vertical motion decoupling，a method to find out decoupling transforms was proposed based on structure preserving isospectral flows (SPIF)，which converting nonlinear problem in the process of finding out decoupling transforms into the solution of Sylvester equation，and matrix Kronecker product knowledge was used for finding decoupling transforms out quickly.Numerical experiments based on pool experiment data show the feasibility of the proposed method.

Key words:ship vertical motion;structure preserving isospectral flows(SPIF);hydrodynamic parameters;quadratic system

0 引 言

船舶在海上的運動(dòng)往往是幾種簡(jiǎn)單運動(dòng)的疊加，可以概括為六個(gè)自由度的搖蕩運動(dòng).船舶各自由度的運動(dòng)是相互耦合的，因此，在研究船舶運動(dòng)特性時(shí)，通常假設六個(gè)自由度運動(dòng)是相互獨立的.但實(shí)際上，船舶運動(dòng)方程可以分解為兩組耦合方程，即縱向運動(dòng)一一升沉、縱搖和縱蕩，橫向運動(dòng)一一橫蕩、橫搖搖和艏搖.縱向運動(dòng)中，縱蕩對升沉和縱搖的耦合作用較小，通常忽略.縱搖和升沉通常是在惡劣氣候里限制航速的主要因素，在大浪中對船體結構有重大影響，因此研究艦船縱向運動(dòng)具有重要意義[1-2].但研究艦船縱向運動(dòng)需要研究其解耦問(wèn)題，以去除縱搖和升沉運動(dòng)的相互耦合影響.

基于船舶水動(dòng)力理論建立的縱向運動(dòng)方程為典型的二階微分系統，為此，本文引入二階系統解耦理論研究該問(wèn)題.數值代數領(lǐng)域通過(guò)保持Lancaster結構來(lái)研究二階系統的解耦問(wèn)題[3-6]，通過(guò)尋找等價(jià)變換來(lái)實(shí)現Lancaster結構中塊陣的對角化，但該變換的數值求解涉及非線(xiàn)性方程組求解問(wèn)題，難以實(shí)現.在文獻[3-4]的基礎上，文獻[5]從理論上證明幾乎對所有的二階系統均存在等價(jià)變換將系統解耦，但并未給出等價(jià)變換的數值求解方法.文獻[6]提出應用保結構同譜流方法研究二階系統的解耦問(wèn)題，通過(guò)一系列的保結構、保譜變換實(shí)現二階系統的解耦.但該方法只能給出解耦后系統的形式，不能給出相應的解耦變換，這使得系統還原及系統分析無(wú)法進(jìn)行，而且該方法所定義的保譜流的保譜性質(zhì)還有待進(jìn)一步完善.

本文提出一種解耦變換的尋找方法，將尋找解耦變換的非線(xiàn)性問(wèn)題轉化為Sylvester方程求解問(wèn)題，并利用矩陣Kronecker積的相關(guān)知識快速而方便地給出二階系統的解耦變換.對水池試驗獲得的艦船縱向運動(dòng)數據所建立的運動(dòng)方程進(jìn)行解耦，數值試驗結果表明該方法確實(shí)可行.

1 二階系統解耦方法

1.1 艦船縱向運動(dòng)方程

根據船舶水動(dòng)力理論，在波浪中航行的船舶在水平舵作用下縱向運動(dòng)方程可表示為：

其中：z表示垂蕩；θ表示縱搖.將式(1)表示成矩陣形式：

其中：M、C、K和F分別為質(zhì)量、阻尼、剛度和外力矩陣，且

則艦船縱向運動(dòng)方程為典型的二階微分系統[1-2].

1.2 基于保結構變換的二階系統解耦

設式(2)的齊次解x(t)有如下形式：

x(t)=eλtu （3）

則數值λ和U向量為二階特征值問(wèn)題的非平凡解.

Q(λ)u=(λ2M+λC+K)u=0 （4）

文獻[3-4]中提出了實(shí)現三矩陣的同時(shí)對角化的一種方法.很容易證明式(4)所描述的二階特征值問(wèn)題等價(jià)與廣義特征值問(wèn)題：

其中，L(λ)為Lancaster結構，

很明顯，M為非奇異時(shí)，有

如果存在非奇異的2n×2n矩陣Πt和Πr表示的等價(jià)變換保持式(6)的Lancaster結構，即

使得MD、CD和KD均為對角矩陣，則式(3)表述的二階特征值問(wèn)題等價(jià)于完全解耦的系統，即

（λ2MD+λCD+KD）z=0 （9）

在MD，M均為非奇異情況下，特征向量u和z具有如下關(guān)系

通過(guò)保持Lancaster結構，該方法將多自由度的系統直接與單自由度的系統鏈接起來(lái).根據Garvey等的思路，文獻[5-6]中給出將原始的n自由度系統解耦為n個(gè)單自由度系統的集合的實(shí)值變換幾乎對所有的二階系統均存在，并通過(guò)保結構同譜流的數值方法來(lái)研究二階系統的解耦.

1.3 基于同譜流的二階系統解耦

若令

定義兩個(gè)隨時(shí)間變化的保結構變換TL(t)，TR(t)∈R2n×2n,且

其中，t∈R，且TL(0)=TR(0)=I2n.若TL(t)和TR(t)非奇異，則（A(t)，B(t)）與(A0，B0)同譜.此時(shí)一類(lèi)具有特殊形式的TL(t)和TR(t)可定義為

其中，Lij(t)，Rij(t)(i，j=1，2)為n×n階矩陣.則有

若保結構由式(11)有

式(13)構成了一個(gè)具有5n2個(gè)方程、8n2個(gè)未知數的線(xiàn)性系統，系統的解具有3n2個(gè)自由度，即三個(gè)n×n階自由參數矩陣.在此按照文獻[6]的方式引入參數矩陣D、NL、NR，使得

則此時(shí)定義了系統參數矩陣M、C、K隨時(shí)間的發(fā)展方向集合，可以利用數值積分的方法來(lái)求解該微分系統的解.但該方法只能求出解耦后的系統參數MD、CD和KD，卻無(wú)法求得解耦變換TL和TR.

2 基于保結構同譜流的系統對角化

2.1 保結構同譜流的實(shí)現

系統參數矩陣M，C，K的同時(shí)對角化可用一個(gè)目標函數來(lái)表示，即

其中，‖·‖F表示矩陣的Frobenius范數；offdiag(M)為矩陣M非對角線(xiàn)上的部分.式(15)前半部分為三個(gè)系數矩陣非對角線(xiàn)上元素的平方和(全局最優(yōu)值為0)，對應于三矩陣的對角形式；后半部分為系數矩陣對角線(xiàn)上元素的平方和.

式(15)描述目標數下降最快的方向為其負梯度方向：

系統參數矩陣M，C，K沿著(zhù)最接近目標函數負梯度的方向發(fā)展變化，則在一定的迭代步數后便可實(shí)現系統的對角化.按照矩陣的Kronecker積相關(guān)知識，式(14)可表示為

其中，vec(X)表示矩陣X的按列向量化.為使M，C，K沿著(zhù)接近目標函數負梯度的方向發(fā)展變化，給出自由參數矩陣D，NL，NR的最小二乘估計：

其中，X+表示矩陣X的Moore-Penrose廣義逆.

2.2 解耦變換的求解

若令

則尋找一對非奇異的Πl，和Πr，滿(mǎn)足

顯然，當M與MD非奇異時(shí)，B與BD非奇異，且

將式(16)變形，得

式(17)中方程2可化為Sylvester方程的一般形式：

其中：為已知矩陣；X=Πr為待求矩陣.式(18)可轉化為齊次線(xiàn)性方程組求解問(wèn)題，其方程式為

因為解耦前后系統具有相同的譜信，式(19)必有非零解，且不難找到非奇異的Πr=unvec(X).其中，vec表示矩陣的向量化函數；unvec表示向量的矩陣化函數[7].

3 船舶運動(dòng)系統解耦實(shí)例

對船模水池實(shí)驗獲得的船舶縱向運動(dòng)方程數據進(jìn)行解調算法仿真(其中，船舶重量為425 t，船長(cháng)為60 m).實(shí)驗中船模航速為18 kn.表1為2個(gè)頻率下的運動(dòng)水動(dòng)力參數.

表1 艦船縱向運動(dòng)水動(dòng)力參數

利用Matlab編譯代碼，并調用其ode函數實(shí)現保結構同普流算法及相應的解耦變換方法.當遭遇頻率為1.308時(shí)，解耦變換及對應的解耦后系統分別為：

圖1為三個(gè)參數矩陣對角線(xiàn)上元素的平方和與的目標函數的變化曲線(xiàn).由圖1可知，曲線(xiàn)最終趨于平穩.圖2為非對角線(xiàn)元素平方和的變化曲線(xiàn).由圖2可知，曲線(xiàn)最終趨于零，實(shí)現了三個(gè)參數矩陣的同時(shí)對角化.

當遭遇頻率為1.068時(shí)，解耦變換及對應的解耦后系統分別為：

圖3為三個(gè)參數矩陣對角線(xiàn)上元素平方和與目標函數變化曲線(xiàn)；圖4為非對角線(xiàn)元素平方和的變化曲線(xiàn).

從實(shí)驗結果可以看出，該方法在極小誤差下實(shí)現了艦船縱向運動(dòng)系統的數值解耦，保證了系統的同譜性質(zhì).給出的相應的解耦變換進(jìn)一步完善了系統的解耦理論.

4 結 論

本文將保結構同譜流方法引入艦船縱向運動(dòng)的解耦研究中，通過(guò)選定合理的系統參數來(lái)實(shí)現縱搖和升沉運動(dòng)的解耦.對保結構同譜方法進(jìn)行改進(jìn)，將求解尋找解耦變換的非線(xiàn)性問(wèn)題轉化為Sylvester方程求解，并利用矩陣的卡式積理論給出解耦變換.根據水池實(shí)驗獲得的縱向運動(dòng)數據進(jìn)行算法仿真，給出浪向角為0°，航速為18 kn時(shí)兩個(gè)頻率下的艦船縱向運動(dòng)系統解耦結果.數值試驗結果表明，該方法可將原始艦船縱向運動(dòng)系統解耦，并給出相應的解耦變換.

參考文獻(References)：

[1]李積德.船舶耐波性[M].哈爾濱：哈爾濱工程大學(xué)出版社，2003.

[2]李積德，王淑娟，李焱，等.基于灰色動(dòng)態(tài)MGM(1，n)模型的艦船縱搖-升沉運動(dòng)預報[J].船舶力學(xué)，2008，12(1):31-36.

[3]GARVEY S D，FRISWELL M I，PRELLS U.Co-ordinate transforms for second order systems，I:General transforms[J].Journal of Sound and Vibration，2002，258(5): 885-909.

[4]GARVEY S D，FRISWELL M I，PRELLS U.Co-ordinate transforms for second order systems，II:Elementary structure-preserving transforms [J].Journal of Sound and Vibration，2002，258(5):911-930.

[5]CHU M T，BUONO N D.Total decoupling of a general quadratic pencil，Part I:Theory [J].Journal of Sound and Vibration，2008，309(1-2):96-111.

[6]CHU M T，BUONO N D.Total decoupling of a general quadratic pencil，Part II:Structure preserving isospectral flows[J]. Journal of Sound and Vibration，2008，309(l­2):112-128.

[7]張賢達.矩陣分析與應用[M].北京:清華大學(xué)出版社，2004.

作者：王淑娟，沈繼紅，李積德  來(lái)源：大連海事大學(xué)學(xué)報