顏凝雨，楊自豪*，張潔瓊，林 巍，劉孟源

(1.西北工業(yè)大學(xué) 數(shù)學(xué)與統(tǒng)計(jì)學(xué)院，西安 710072；2.中交懸浮隧道工程技術(shù)聯(lián)合研究組，珠海 519000)

1 引 言

任何非流線型物體在一定的恒定流速下，都會(huì)在物體兩側(cè)交替地產(chǎn)生脫離物理結(jié)構(gòu)表面的周期性漩渦，從而誘發(fā)彈性支撐結(jié)構(gòu)物的渦激振動(dòng)現(xiàn)象VIV[1]。大量實(shí)際工程問題中，許多圓柱結(jié)構(gòu)，如懸浮隧道、深水立管、高空桅桿和新型風(fēng)力發(fā)電等，均可能出現(xiàn)渦激振動(dòng)現(xiàn)象。雖然合理利用渦激振動(dòng)會(huì)為人們帶來便利，但其誘發(fā)的共振也會(huì)對(duì)結(jié)構(gòu)產(chǎn)生較大的疲勞損害。

研究圓柱結(jié)構(gòu)的渦激振動(dòng)主要有試驗(yàn)法[2-4]、數(shù)值方法[5,6]和經(jīng)驗(yàn)?zāi)Ｐ蚚7,8]。試驗(yàn)方法雖然直觀可靠，但成本較高；數(shù)值方法雖然在計(jì)算精度上占有一定優(yōu)勢(shì)，但計(jì)算較為復(fù)雜，尤其是在高雷諾數(shù)的情況下。而經(jīng)驗(yàn)?zāi)Ｐ涂奢^快計(jì)算渦激振動(dòng)的主要特征，還可以在模型的指導(dǎo)下設(shè)計(jì)試驗(yàn)或者研究某個(gè)參數(shù)對(duì)試驗(yàn)結(jié)果的影響，因此經(jīng)驗(yàn)?zāi)Ｐ头ㄕ艿絿?guó)內(nèi)外諸多學(xué)者的關(guān)注，其中尾流振子模型是應(yīng)用較為廣泛的一種模型。

Bishop等[9]認(rèn)為流體的尾流對(duì)圓柱體結(jié)構(gòu)的作用相當(dāng)于一個(gè)非線性振子，首次提出了尾流振子模型。之后，Hartlen等[10]將Vander Pol方程與結(jié)構(gòu)振動(dòng)方程聯(lián)立，使尾流振子模型更加具體化。Skop等[11]進(jìn)一步改進(jìn)了Hartlen的模型，對(duì)模型經(jīng)驗(yàn)參數(shù)進(jìn)行了研究；Landle[12]提出了一個(gè)擁有五階非線性氣動(dòng)阻尼項(xiàng)的無量綱尾流振子模型；Facchinetti等[13]通過分析對(duì)比發(fā)現(xiàn)加速度耦合模型可以更準(zhǔn)確地定量預(yù)測(cè)渦激振動(dòng)的基本特征；Ogink等[14]對(duì)Facchinetti的尾流振子模型建立過程進(jìn)行了具體分析，通過引入頻率耦合變量進(jìn)一步改進(jìn)了模型。但上述對(duì)尾流振子的研究大多針對(duì)橫向單自由度，而實(shí)際的渦激振動(dòng)中橫向振動(dòng)和順流向振動(dòng)往往相互影響，許多學(xué)者也針對(duì)雙自由度的尾流振子模型做了研究。Krenk等[15]根據(jù)流固耦合的能量流平衡原則構(gòu)建了雙振子模型；Srinil等[16]提出了一種改進(jìn)的用于預(yù)測(cè)均勻流中柔性圓柱體的二維耦合尾流雙振子模型。康莊等[17]建立了非線性結(jié)構(gòu)振子和流體振子的耦合模型，可以較準(zhǔn)確預(yù)報(bào)低質(zhì)量比圓柱結(jié)構(gòu)雙自由度渦激振動(dòng)的特征。Yang等[18]提出一種新的渦激振動(dòng)非線性耦合尾流振子模型，其只使用一個(gè)尾流振蕩器來耦合橫向和順流向的運(yùn)動(dòng)。在實(shí)際應(yīng)用中，不是所有模型的參數(shù)都可以通過物理方法測(cè)量得到，尾流振子模型的應(yīng)用往往面臨參數(shù)選取的難題，參數(shù)選取的偏差可能導(dǎo)致計(jì)算結(jié)果的不準(zhǔn)確，從而無法充分發(fā)揮經(jīng)驗(yàn)?zāi)Ｐ驮谘芯繙u激振動(dòng)問題中的優(yōu)勢(shì)。

解決復(fù)雜模型參數(shù)反演問題的方法有很多種，其中神經(jīng)網(wǎng)絡(luò)方法以其高速尋優(yōu)和自學(xué)習(xí)等優(yōu)勢(shì)在參數(shù)反演問題中得到廣泛的應(yīng)用，如Spichak等[19]運(yùn)用BP神經(jīng)網(wǎng)絡(luò)求解地電三維逆問題，證實(shí)神經(jīng)網(wǎng)絡(luò)方法可以成功地用于地電三維逆問題參數(shù)反演；Cerdena等[20]在測(cè)定海洋層積云物理性質(zhì)的研究中，運(yùn)用遺傳算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)對(duì)輻射傳輸模型中的參數(shù)進(jìn)行了反演；Zhao等[21]提出一種人工蜂群算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)，并對(duì)巖土工程中的模型參數(shù)進(jìn)行了反演。國(guó)內(nèi)外沒有針對(duì)尾流振子模型參數(shù)反演的研究工作。

本文研究尾流雙振子模型中參數(shù)的反演方法。首先利用降階法將模型變?yōu)橐粋€(gè)常系數(shù)多變量的非線性一階常微分方程組，并利用龍格-庫塔方法進(jìn)行求解。設(shè)計(jì)懸浮隧道物理模型的水槽試驗(yàn)，在給定的多組物理參數(shù)和試驗(yàn)數(shù)據(jù)下，構(gòu)建最優(yōu)化模型，利用遺傳算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)智能算法反演未知參數(shù)。

2 渦激響應(yīng)尾流雙振子模型

柔性圓柱結(jié)構(gòu)雙自由度渦激振動(dòng)模型[16]如 圖1 所示。假定直徑為D的圓柱結(jié)構(gòu)無限長(zhǎng)，水平和豎直方向，即X和Y兩個(gè)方向的彈簧剛度分別為Kx和Ky,水的流速為V,圓柱體受到雙向四個(gè)彈簧的作用力，在X和Y兩個(gè)空間方向上自由振動(dòng)，無量綱形式的渦激響應(yīng)尾流雙振子模型可表示為[16]

(1)

圖1 柔性圓柱結(jié)構(gòu)雙自由度渦激振動(dòng)

λx=2ξxf*+γΩ/μx,λy=2ξy+γΩ/μy

(2)

式中ξx和ξy為結(jié)構(gòu)水平和豎向運(yùn)動(dòng)阻尼比，取值范圍為[0.01,0.08]；f*為X和Y兩方向結(jié)構(gòu)固有頻率之比，取值范圍為[0.5,1.5]，γ為Stall參數(shù)，取值為0.8；Ω=St·Vr為靜水中旋渦脫落頻率與圓柱橫流固有頻率之比，St為Strouhal數(shù)，取值為0.2，Vr=V/(Dfn)為約化速度,fn為Y方向結(jié)構(gòu)固有頻率，μx和μy表達(dá)式為

(3)

此外，(αx,αy,βx，βy)為無量綱幾何非線性系數(shù)，表征位移非線性項(xiàng)對(duì)本身的影響。MD和ML為系統(tǒng)質(zhì)量參數(shù)，

(4)

式中CD 0和CL0分別為靜置圓柱結(jié)構(gòu)的阻力和升力系數(shù)，取值為0.2[16]和0.3；St為斯特魯哈爾數(shù)，取值為0.2；εx和εy為尾流系數(shù)；Λx和Λy為經(jīng)驗(yàn)耦合參數(shù)，均取為12。

對(duì)于質(zhì)量分布均勻的圓柱體，有

(5)

且可取εx=0.3，εy由試驗(yàn)數(shù)據(jù)擬合得到，其擬合表達(dá)式為

εy(m*)=b1eb2 m*

(6)

式中b1=0.00234，b2=0.228。

由式(1)可知尾流雙振子模型是一個(gè)多變量的二階非線性常微分方程組，現(xiàn)運(yùn)用降階法將該模型化為一個(gè)常系數(shù)多變量的一階非線性常微分方程組。引入如下變量

u=(u1,u2,u3,u4,u5,u6,u7,u8)T=

(7)

式中T為矩陣的轉(zhuǎn)置。方程組(1)可轉(zhuǎn)化為一階常微分方程組(8)。

(8)

進(jìn)一步，式(7，8)可以改寫為

(9)

式中M和g的定義如下

(10)

g(t,u,αx,αy,βx,βy)=(g1,g2,g3,g4,g5,g6,g7,g8)T

式中

g1=u5，g2=u6，g3=u7，g4=u8

MDΩ2u2-2πMLΩ2(u4u7/Vr)

MLΩ2u4+2πMDΩ2(u2u7/Vr)

(11)

給定初值u0,該一階非線性常微分方程系統(tǒng)可使用經(jīng)典的龍格-庫塔方法進(jìn)行求解。

3 圓柱結(jié)構(gòu)渦激振動(dòng)響應(yīng)水槽試驗(yàn)

3.1 試驗(yàn)?zāi)康?/h3>

尾流雙振子模型中(αx,αy,βx，βy)四個(gè)參數(shù)表征位移非線性項(xiàng)的影響，無法通過試驗(yàn)直接測(cè)得，其他參數(shù)均可通過試驗(yàn)或經(jīng)驗(yàn)獲得。若想得到上述四個(gè)參數(shù)的值，由模型(9)可知，一種可行的方案是通過圓柱體渦激振動(dòng)時(shí)的位移數(shù)據(jù)來反演這四個(gè)未知參數(shù)。因此，本節(jié)設(shè)計(jì)圓柱結(jié)構(gòu)的水槽試驗(yàn)，測(cè)得圓柱結(jié)構(gòu)在不同水流力(流速)作用下水平和豎直方向的位移時(shí)間序列及運(yùn)動(dòng)軌跡。

3.2 試驗(yàn)設(shè)計(jì)

現(xiàn)有的圓柱結(jié)構(gòu)渦激振動(dòng)試驗(yàn)不考慮懸浮隧道等實(shí)際工程結(jié)構(gòu)試驗(yàn)節(jié)段的真實(shí)剛度，試驗(yàn)結(jié)果無法反映真實(shí)結(jié)構(gòu)的剛度和自振頻率，試驗(yàn)結(jié)果與實(shí)際值誤差較大。本節(jié)給出新型的圓柱結(jié)構(gòu)水槽試驗(yàn)設(shè)計(jì)方案。

首先運(yùn)用等效固有振頻法計(jì)算懸浮隧道等實(shí)際工程結(jié)構(gòu)某一節(jié)斷的真實(shí)剛度值k0,計(jì)算公式為

k0=(4π2m)/T2

(12)

式中m為某一節(jié)段的質(zhì)量，T為該節(jié)段在水中的固有頻率。其次根據(jù)相似準(zhǔn)則，將真實(shí)剛度k0按照比尺換算成試驗(yàn)?zāi)Ｐ偷膭偠萲1,并據(jù)此配置水平和豎向彈簧的剛度。最后將試驗(yàn)?zāi)Ｐ头湃胨刂羞M(jìn)行試驗(yàn)，浮筒式試驗(yàn)裝置如圖2所示。

圖2 試驗(yàn)裝置設(shè)計(jì)

試驗(yàn)采用直徑為8 cm和長(zhǎng)度為70 cm的透明樹脂材料制成的空心圓柱體模擬深海懸浮隧道等實(shí)際工程結(jié)構(gòu)的節(jié)段(等剛度模型)，模型重 3.3 kg，浮重比1.07，試驗(yàn)比尺為 1∶157.5。圓柱體使用4根彈簧與透明固定板鏈接，其中橫向和豎向彈簧剛度分別取為30.6 N/m和21.6 N/m，彈簧原始長(zhǎng)度均為100 mm，圓柱體、彈簧和透明固定板之間均采用L型鋼扣栓接，如圖3(a)所示。靜止?fàn)顟B(tài)下，彈簧長(zhǎng)度為192 mm，如圖3(b)所示。

選取天津水運(yùn)工程科學(xué)研究院水槽廳的波流水槽為試驗(yàn)場(chǎng)地，水槽寬70 cm，長(zhǎng)30 m，水深 57 cm，造流范圍為0.02 m/s～0.5 m/s，水槽尾端采用斜坡式消波框防止波浪反射。將試驗(yàn)?zāi)Ｐ头湃?.5倍水深處進(jìn)行試驗(yàn)，并在模型截面貼上紅色方塊標(biāo)志點(diǎn)，如圖3(b)所示。不斷改變水的流速，利用相機(jī)對(duì)標(biāo)志點(diǎn)進(jìn)行動(dòng)態(tài)捕捉，然后利用 Matlab 編寫的視頻識(shí)別程序，抓取每一幀中標(biāo)志點(diǎn)中心的位置，并將前后兩幀位置數(shù)據(jù)進(jìn)行比較，從而得到圓柱體在不同物理流速下水平和垂向位移隨時(shí)間變化的數(shù)據(jù)。試驗(yàn)共進(jìn)行9組，每組試驗(yàn)所對(duì)應(yīng)的水的流速和約化速度列入表1。

圖3 試驗(yàn)裝置

表1 9組試驗(yàn)中水的流速和約化速度Tab.1 Water velocity and reduced velocity in 9 sets of tests

3.3 試驗(yàn)結(jié)果

不同流速下圓柱結(jié)構(gòu)的運(yùn)動(dòng)軌跡經(jīng)過無量綱處理后如圖4所示?？梢钥闯?，運(yùn)動(dòng)軌跡呈8字形，且隨著流速的增大，模型在豎直方向的運(yùn)動(dòng)幅度增大，在第五組達(dá)到最大，考慮可能達(dá)到了鎖定流速區(qū)間，即當(dāng)流速在此區(qū)間內(nèi)時(shí)，圓柱結(jié)構(gòu)振動(dòng)頻率近乎不變，振幅維持較高水平；且隨著流速增加，模型所受橫向拖拽力增加，無量綱橫向位移一直增大，這與以往研究保持一致[8]。

4 參數(shù)反演的BP神經(jīng)網(wǎng)絡(luò)模型

4.1 參數(shù)反演的最優(yōu)化模型

根據(jù)第3節(jié)中每組試驗(yàn)，可以得到對(duì)應(yīng)的位移試驗(yàn)數(shù)據(jù)

(i=1,…,n;ti∈[0,tmax](13)

而每給定一組參數(shù)(αx,αy,βx，βy),求解微分方程系統(tǒng)(9)，可得到對(duì)應(yīng)的位移計(jì)算數(shù)據(jù)，記作

x(t,αx,αy,βx,βy),y(t,αx,αy,βx,βy)

(t∈[0,tmax])

(14)

下面構(gòu)建參數(shù)反演的最優(yōu)化模型，即求(αx,αy,βx,βy),使

(15)

由于本文關(guān)注的核心為振子的幅值和均方根，因此可引入向量函數(shù)

opt(x,y)=[rms(x),rms(y),max(x),max(y),

min(x),min(y)]T

式中max(·)和min(·)分別為向量的最大和最小分量，rms(·)為均方根。

N為向量x的長(zhǎng)度。目標(biāo)函數(shù)可以取為

(16)

求解最優(yōu)化問題(15，16)即可得到未知參數(shù)的一組最優(yōu)值。

圖4 不同流速下圓柱的運(yùn)動(dòng)軌跡

4.2 BP神經(jīng)網(wǎng)絡(luò)的構(gòu)建

人工神經(jīng)網(wǎng)絡(luò)是一種由大量處理單元互聯(lián)組成的非線性、自適應(yīng)信息處理系統(tǒng)，具有快速求解最優(yōu)化解的能力，其最小最重要的單元叫神經(jīng)元，每個(gè)神經(jīng)元都有輸入連接和輸出連接，這些連接模擬了大腦中突觸的行為，每一個(gè)連接都有對(duì)應(yīng)權(quán)重。BP神經(jīng)網(wǎng)絡(luò)由輸入層、隱藏層和輸出層構(gòu)成，是一種按誤差反向傳播訓(xùn)練的多層前饋網(wǎng)絡(luò)。BP神經(jīng)網(wǎng)絡(luò)求解優(yōu)化問題的基本思想是利用BP算法，通過信息的正向傳播和誤差的反向傳播，不斷調(diào)整網(wǎng)絡(luò)中的權(quán)值和閾值，使輸出和樣本的實(shí)際標(biāo)簽之間的誤差達(dá)到最小。渦激響應(yīng)尾流雙振子模型的參數(shù)反演思路為，生成多組未知參數(shù)數(shù)據(jù)，通過求解尾流雙振子模型(9)得到每組未知參數(shù)下的位移時(shí)間序列數(shù)據(jù)；分別計(jì)算每組豎直和水平位移數(shù)據(jù)的最大值、最小值和均方根，并以這6個(gè)數(shù)據(jù)作為神經(jīng)網(wǎng)絡(luò)的輸入，以未知的4個(gè)參數(shù)為輸出，優(yōu)化目標(biāo)為式(15)，即可訓(xùn)練出對(duì)應(yīng)的神經(jīng)網(wǎng)絡(luò)模型；將試驗(yàn)得到的一組位移數(shù)據(jù)代入已訓(xùn)練好的神經(jīng)網(wǎng)絡(luò)，即可輸出一組未知參數(shù)的值。

4.3 基于遺傳算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)

BP神經(jīng)網(wǎng)絡(luò)存在局部極值，誤差較大。現(xiàn)采用遺傳算法對(duì)BP神經(jīng)網(wǎng)絡(luò)進(jìn)行優(yōu)化，即使用遺傳算法對(duì)BP神經(jīng)網(wǎng)絡(luò)的初始權(quán)值和閾值進(jìn)行優(yōu)化，進(jìn)一步提高參數(shù)反演的精度。遺傳算法以優(yōu)勝劣汰的自然選擇和基因階段的遺傳變異原理作為基礎(chǔ)，是一種具有很強(qiáng)隨機(jī)搜索能力的智能算法，可以取得很好的搜索效果。遺傳算法對(duì)BP神經(jīng)網(wǎng)絡(luò)的優(yōu)化，實(shí)際上是對(duì)BP神經(jīng)網(wǎng)絡(luò)的初始權(quán)值和閾值進(jìn)行優(yōu)化。首先確立BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)并得到一組初始權(quán)值和閾值，遺傳算法將對(duì)其編碼，并再次通過神經(jīng)網(wǎng)絡(luò)訓(xùn)練，得到的誤差作為適應(yīng)度。然后不斷經(jīng)過選擇、交叉和變異操作使適應(yīng)度逐漸增加，最終滿足精度條件，得到最優(yōu)權(quán)值和閾值。最后再代回BP神經(jīng)網(wǎng)絡(luò)進(jìn)行求解。具體流程如圖5所示。

4.4 參數(shù)設(shè)置及誤差定義

尾流雙振子模型的參數(shù)取值列入表2。

反演數(shù)據(jù)與試驗(yàn)數(shù)據(jù)的誤差計(jì)算公式為

(17)

表2 參數(shù)取值

圖5 求解參數(shù)反演問題分析流程

5 參數(shù)反演結(jié)果分析

利用遺傳算法、BP神經(jīng)網(wǎng)絡(luò)算法和基于遺傳算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)算法反演參數(shù)。遺傳算法中初始種群數(shù)目取為20，進(jìn)化代數(shù)取為80，交叉概率取為0.95，變異概率取為0.001。BP神經(jīng)網(wǎng)絡(luò)由一個(gè)輸入層、一個(gè)隱藏層和一個(gè)輸出層構(gòu)成，網(wǎng)絡(luò)模型參數(shù)選取如下，輸入層神經(jīng)元個(gè)數(shù)取為6，隱藏層神經(jīng)元個(gè)數(shù)取為8，輸出層神經(jīng)元個(gè)數(shù)取為2，數(shù)據(jù)庫樣本個(gè)數(shù)取為8000，訓(xùn)練集個(gè)數(shù)取為7200，驗(yàn)證集個(gè)數(shù)取為800，訓(xùn)練誤差目標(biāo)取為10-4，學(xué)習(xí)效率取為0.4，傳遞函數(shù)取為{tansig,purelin}，未知參數(shù)αx和αy的取值范圍為(0,1)，未知參數(shù)βx和βy的取值范圍為(0,10)。此外，每組試驗(yàn)數(shù)據(jù)分別訓(xùn)練50次，取平均值作為參數(shù)的最終預(yù)測(cè)結(jié)果。需要指出的是，由于試驗(yàn)誤差，試驗(yàn)得到的位移數(shù)據(jù)沒有呈現(xiàn)穩(wěn)定后的周期震蕩模型，但大致穩(wěn)定在某個(gè)固定區(qū)間之內(nèi)。為了減小參數(shù)反演誤差，現(xiàn)對(duì)神經(jīng)網(wǎng)絡(luò)的輸入?yún)?shù)，即豎直和橫向位移的最大和最小實(shí)驗(yàn)數(shù)據(jù)做以下處理，以第3組試驗(yàn)的水平方向位移最大值為例，利用Matlab程序找到數(shù)據(jù)的峰值，再求出峰值的平均值作為神經(jīng)網(wǎng)絡(luò)位移最大值的輸入，如圖6所示?；谶z傳算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)算法流程如圖5所示，其中遺傳算法和BP神經(jīng)網(wǎng)絡(luò)的初始參數(shù)如前所述。

圖6 試驗(yàn)數(shù)據(jù)輸入最大值的處理

利用遺傳算法、BP神經(jīng)網(wǎng)絡(luò)算法和基于遺傳算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)算法反演參數(shù)的誤差Ret列入表3，誤差隨水的流速變化如圖7所示。由表3可知，三種方法所得結(jié)果的平均誤差分別為10.37%、7.53%和5.50%，基于遺傳算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)算法的平均誤差相比遺傳算法和BP神經(jīng)網(wǎng)絡(luò)算法分別降低了46.96%和26.96%。此外，由圖7可知，遺傳算法的穩(wěn)定性最差，BP神經(jīng)網(wǎng)絡(luò)次之，遺傳算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)的穩(wěn)定性最好。遺傳算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)參數(shù)參演結(jié)果列入表4。

圖7 三種模型方法的誤差比較

圖8為利用遺傳算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)的反演結(jié)果與第3、4和第6組試驗(yàn)無量綱位移數(shù)據(jù)的對(duì)比，其中圖8左圖分別表示第3、4和第6組試驗(yàn)方法與反演方法水平和豎直方向無量綱位移數(shù)據(jù)的對(duì)比，圖8右圖分別表示第3、4和第6組試驗(yàn)數(shù)據(jù)與反演數(shù)據(jù)水平和豎直方向運(yùn)動(dòng)軌跡的對(duì)比?？梢钥闯觯谶z傳算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)算法反演結(jié)果與試驗(yàn)結(jié)果吻合較好。因此，基于遺傳算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)智能算法能夠精確穩(wěn)定實(shí)現(xiàn)尾流雙振子模型的參數(shù)確定。

表4 遺傳算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)參數(shù)反演結(jié)果Tab.4 Parameter inversion results of BP neural network optimized by genetic algorithm

最后以第3組參數(shù)反演問題為例進(jìn)行靈敏度分析，研究反演參數(shù)的微小偏差對(duì)反演精度的影響。將兩個(gè)參數(shù)值分別在參數(shù)反演結(jié)果附近微小擾動(dòng)，圖9給出了誤差Ret的值，其中灰色點(diǎn)代表遺傳算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)參數(shù)反演結(jié)果。

可以看出，參數(shù)反演結(jié)果處于在局部最優(yōu)解附近，進(jìn)一步驗(yàn)證了反演方法的正確性。此外，通過計(jì)算誤差沿兩個(gè)參數(shù)的變化率可發(fā)現(xiàn)，誤差Ret對(duì)參數(shù)βx和βy更加敏感。在反演結(jié)果(0.8214±0.2，10.0674±0.5)附近的參數(shù)取值(灰色實(shí)線區(qū)域)對(duì)應(yīng)的誤差Ret約為6%，在反演結(jié)果(0.8214±0.5，10.0674±1)附近的參數(shù)取值(灰色虛線區(qū)域)對(duì)應(yīng)的誤差Ret約為10%。該靈敏度分析結(jié)果可為類似問題的參數(shù)取值提供可靠的理論分析基礎(chǔ)。

圖8 第3、4和第6組試驗(yàn)位移數(shù)據(jù)與反演數(shù)據(jù)對(duì)比

圖9 第3組參數(shù)反演問題中的參數(shù)取值域

6 結(jié) 論

本文主要研究了尾流雙振子模型中未知參數(shù)的反演方法，給出了基于物理試驗(yàn)和BP神經(jīng)網(wǎng)絡(luò)模型的參數(shù)反演方法。利用降階法將尾流雙振子模型轉(zhuǎn)換為一階常微分方程組，并使用龍格-庫塔方法實(shí)現(xiàn)模型的求解。給出了一種新型的圓柱結(jié)構(gòu)水槽試驗(yàn)設(shè)計(jì)方案，獲取了圓柱結(jié)構(gòu)在水平和豎直方向的位移試驗(yàn)數(shù)據(jù)?；谠囼?yàn)數(shù)據(jù)，構(gòu)建了最優(yōu)化模型，并利用基于遺傳算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)算法對(duì)尾流雙振子模型中四個(gè)未知參數(shù)進(jìn)行了反演，誤差處于可接受范圍之內(nèi)。本文的工作將為懸浮隧道等實(shí)際工程結(jié)構(gòu)的渦激振動(dòng)響應(yīng)分析提供理論基礎(chǔ)。