基于LSTM的振動(dòng)臺(tái)參數(shù)整定方法

張文鵬，紀(jì)金豹，王東岳

(北京工業(yè)大學(xué)，工程抗震與結(jié)構(gòu)診治北京市重點(diǎn)實(shí)驗(yàn)室，北京100124)

0 前言

地震模擬振動(dòng)臺(tái)試驗(yàn)是真正意義上的地震模擬試驗(yàn)，為結(jié)構(gòu)彈塑性地震響應(yīng)研究提供了一種有效的手段[1-4]。地震模擬振動(dòng)臺(tái)一般采用三參量控制算法。三參量控制改善了振動(dòng)臺(tái)系統(tǒng)特性，擴(kuò)展系統(tǒng)加速度響應(yīng)頻寬，增強(qiáng)系統(tǒng)穩(wěn)定性[5]。然而，三參量控制參數(shù)較多，包括前饋增益和反饋增益參數(shù)，在實(shí)際工程應(yīng)用中，三參量控制參數(shù)整定是一個(gè)較困難的環(huán)節(jié)。

目前，在實(shí)際工程中常用手動(dòng)調(diào)參法，通過(guò)逐個(gè)調(diào)整控制增益值，并從中選取控制效果較好的一組增益參數(shù)，以此來(lái)完成三參量的參數(shù)整定。這種方法不僅工作量大，而且參數(shù)很難達(dá)到最優(yōu)。近年來(lái)，國(guó)內(nèi)學(xué)者提出了多種三參量參數(shù)整定方法。欒強(qiáng)利等[6-7]提出一種快速三參量整定法，雖提高了整定速度，但需要復(fù)雜的公式推導(dǎo)；紀(jì)金豹等[8-9]提出基于頻域辨識(shí)和基于專家經(jīng)驗(yàn)的兩種三參量參數(shù)整定算法，但前者需要復(fù)雜的理論計(jì)算，后者需要花時(shí)間獲取和總結(jié)專家經(jīng)驗(yàn)；高春華等[10]利用粒子群算法的尋優(yōu)能力進(jìn)行三參量參數(shù)整定，但其粒子群也只能針對(duì)特定振動(dòng)臺(tái)進(jìn)行調(diào)參，適用性較差。從目前的研究成果可以看出，各種參數(shù)整定算法都能提高參數(shù)整定的效率，但在實(shí)際工程應(yīng)用中還存在各種各樣的局限性。

為了解決目前三參量控制參數(shù)整定的問(wèn)題，提出一種基于LSTM的振動(dòng)臺(tái)三參量控制參數(shù)整定算法。首先，介紹振動(dòng)臺(tái)三參量控制系統(tǒng)的原理；接著詳細(xì)介紹了基于LSTM的振動(dòng)臺(tái)參數(shù)整定算法的詳細(xì)過(guò)程；最后為了驗(yàn)證算法有效性，在Python環(huán)境下進(jìn)行時(shí)頻域的仿真驗(yàn)證。

1 振動(dòng)臺(tái)系統(tǒng)與參數(shù)整定算法原理

1.1 振動(dòng)臺(tái)三參量控制系統(tǒng)原理

振動(dòng)臺(tái)三參量控制系統(tǒng)包括三參量前饋和三參量反饋環(huán)節(jié)，加速度信號(hào)輸入三參量前饋環(huán)節(jié)，生成加速度前饋、速度前饋和位移前饋信號(hào)，后與三參量反饋回來(lái)的加速度反饋、速度反饋和位移反饋信號(hào)一起輸入增益合成器中。在增益合成器中，3個(gè)前饋信號(hào)和3個(gè)反饋信號(hào)各自乘以增益系數(shù)后合并為一個(gè)控制信號(hào)u輸入振動(dòng)臺(tái)系統(tǒng)中，最終得到系統(tǒng)的輸出。在閉環(huán)系統(tǒng)中加入加速度反饋和速度反饋，提高了系統(tǒng)阻尼和系統(tǒng)固有頻率，進(jìn)而提高了系統(tǒng)穩(wěn)定性，前饋環(huán)節(jié)通過(guò)消除系統(tǒng)閉環(huán)傳遞函數(shù)中距離虛軸較近的極點(diǎn)，進(jìn)而拓展系統(tǒng)的頻寬[11]。三參量控制下的系統(tǒng)整體模型如圖1所示。

圖1 振動(dòng)臺(tái)三參量閉環(huán)控制模型

圖1中：E為輸入加速度信號(hào)；G(a)、G(v)、G(d)分別表示三參量發(fā)生器的加速度、速度、位移前饋分量的傳遞函數(shù)；V(s)表示三參量反饋環(huán)節(jié)的速度合成器；Aa、Ad、Av、Afd、Afa、Afv為待整定的控制參數(shù)；af、df、vf為振動(dòng)臺(tái)反饋加速度、位移和速度信號(hào)；u為振動(dòng)臺(tái)輸入信號(hào)；G(s)表示振動(dòng)臺(tái)開(kāi)環(huán)傳遞函數(shù)。各部分表達(dá)式[12]如下：

(1)

(2)

(3)

(4)

(5)

1.2 參數(shù)整定算法原理

1.2.1 網(wǎng)絡(luò)模型

深度學(xué)習(xí)在模式識(shí)別、系統(tǒng)辨識(shí)等領(lǐng)域有很好的應(yīng)用前景[13]。循環(huán)神經(jīng)網(wǎng)絡(luò)(RNN)適用于處理與時(shí)間軸相關(guān)的問(wèn)題，振動(dòng)臺(tái)問(wèn)題就是一種時(shí)間序列處理的問(wèn)題。但普通RNN很難處理時(shí)間序列的長(zhǎng)期依賴問(wèn)題(如距離很遠(yuǎn)的兩步之間的依賴)，因?yàn)殡S著信號(hào)長(zhǎng)度的增加，容易發(fā)生梯度消失或梯度爆炸。在傳統(tǒng)RNN的基礎(chǔ)上，1997年HOCHREITER和SCHMIDHUBER[14]提出LSTM。LSTM的提出解決了傳統(tǒng)RNN的長(zhǎng)期依賴問(wèn)題[15]。后來(lái)隨著深度學(xué)習(xí)的興起，LSTM又經(jīng)過(guò)了許多人的創(chuàng)新發(fā)展[16-17]，已經(jīng)形成比較系統(tǒng)完整的框架，并且在很多領(lǐng)域取得了廣泛的應(yīng)用[18-22]。

LSTM網(wǎng)絡(luò)由輸入層、隱藏層和輸出層構(gòu)成，網(wǎng)絡(luò)結(jié)構(gòu)示意如圖2所示。

圖2 LSTM網(wǎng)絡(luò)結(jié)構(gòu)

如圖2所示，隱藏層是由循環(huán)連接的隱藏單元組成的，神經(jīng)單元中包含3個(gè)門(mén)(輸入門(mén)σi、輸出門(mén)σo和遺忘門(mén)σf)，分別負(fù)責(zé)寫(xiě)入、讀取和選擇的操作。前向傳播公式[15]如下：

(6)

式中：W和b代表網(wǎng)絡(luò)的權(quán)重和偏執(zhí)；xt和ht分別為L(zhǎng)STM的隱藏層輸入和輸出；ct為記憶單元；σ(x)=1/(1-e-x)為sigmoid函數(shù)，它可以將一個(gè)實(shí)數(shù)映射到區(qū)間(0，1)。

輸入門(mén)σi決定了輸入量xt有多少可以保存到記憶中去；遺忘門(mén)σf決定了上一次記憶ct-1有多少可以保留到這次記憶ct中，可以選擇保留或遺忘歷史信息，這是LSTM的關(guān)鍵組成部分，并且可以避免梯度反向傳播時(shí)引發(fā)的梯度消失和爆炸問(wèn)題；輸出門(mén)σo決定了記憶單元ct對(duì)當(dāng)前神經(jīng)元的輸出ht的影響，即判斷當(dāng)前輸出需要調(diào)用多少歷史記憶。

1.2.2 基于LSTM的參數(shù)整定原理

基于LSTM的振動(dòng)臺(tái)參數(shù)整定主要有以下幾個(gè)部分：選波并獲取輸入輸出數(shù)據(jù)、LSTM網(wǎng)絡(luò)訓(xùn)練、選擇初始參數(shù)、參數(shù)整定、驗(yàn)證。參數(shù)整定結(jié)構(gòu)和流程分別如圖3、4所示。

圖3 基于LSTM的振動(dòng)臺(tái)控制參數(shù)整定結(jié)構(gòu)

圖4 基于LSTM的振動(dòng)臺(tái)控制參數(shù)整定流程

基于LSTM的參數(shù)整定算法的具體過(guò)程如下：

(1)選波并獲取輸入輸出數(shù)據(jù)。在系統(tǒng)辨識(shí)以及參數(shù)整定過(guò)程中都需要選擇合適的地震波加速度時(shí)程信號(hào)。地震波的處理包括：保持所有地震波的采樣頻率一致；取地震波0.1～50 Hz范圍內(nèi)的信號(hào)；將處理好的地震波加速度信號(hào)輸入振動(dòng)臺(tái)系統(tǒng)，并采集得到系統(tǒng)輸入輸出信號(hào)。

(2)LSTM網(wǎng)絡(luò)訓(xùn)練。對(duì)應(yīng)圖3的LSTM系統(tǒng)辨識(shí)模塊，在進(jìn)行系統(tǒng)辨識(shí)時(shí)，要求足夠的振動(dòng)臺(tái)輸入輸出加速度數(shù)據(jù)，并按照6∶1∶3的比例劃分訓(xùn)練集、驗(yàn)證集、測(cè)試集。這里選擇30條地震波進(jìn)行LSTM網(wǎng)絡(luò)訓(xùn)練，5條地震波驗(yàn)證，10條地震波進(jìn)行測(cè)試，當(dāng)系統(tǒng)辨識(shí)的結(jié)果滿足10條測(cè)試波的相關(guān)系數(shù)和峰值誤差分別大于95%、小于10%即認(rèn)為該LSTM網(wǎng)絡(luò)模型已可以正確表達(dá)振動(dòng)臺(tái)開(kāi)環(huán)模型的系統(tǒng)特性。相關(guān)系數(shù)ε和峰值誤差τ的表達(dá)如下：

(7)

(8)

式中：E為加速度控制指令信號(hào)；yi為系統(tǒng)輸出加速度信號(hào)。

(3)選擇初始參數(shù)。完整系統(tǒng)辨識(shí)后，需要選擇一組合適的初始參數(shù)。通過(guò)實(shí)際驗(yàn)證可知，初始參數(shù)的選擇影響訓(xùn)練的效率以及訓(xùn)練結(jié)果的好壞。利用專家經(jīng)驗(yàn)法進(jìn)行初始參數(shù)的選擇，選擇依據(jù)[9]如下：

(9)

(4)參數(shù)整定。進(jìn)行參數(shù)整定時(shí)所需的地震波加速度信號(hào)一般選擇常見(jiàn)的地震波，如EL-centro波、Taft波、CHICHI波等。經(jīng)驗(yàn)證，參數(shù)整定時(shí)地震波的條數(shù)對(duì)于訓(xùn)練結(jié)果的影響不大，甚至3—5條的訓(xùn)練結(jié)果相較8—10條更好。參數(shù)整定采用梯度下降算法，損失函數(shù)選擇MSE，表達(dá)式如下：

(10)

式中：m為信號(hào)長(zhǎng)度；Ei為加速度控制指令信號(hào)；yi為系統(tǒng)輸出加速度信號(hào)；Aa、Ad、Av、Afd、Afa、Afv的整定采用梯度下降法：

(11)

式中：α為學(xué)習(xí)率。采用ADAM優(yōu)化算法，通過(guò)計(jì)算每一輪梯度的一階和二階指數(shù)移動(dòng)加權(quán)平均值，實(shí)時(shí)優(yōu)化調(diào)整學(xué)習(xí)率α，如下式所示：

(12)

式中：β1、β2(取值區(qū)間0～1)為衰減率；gi為各分量的梯度；mi為一階指數(shù)移動(dòng)加權(quán)平均值；ui為二階指數(shù)移動(dòng)加權(quán)平均值。偏差修正公式為

(13)

(14)

(5)驗(yàn)證。參數(shù)整定完畢后需要在振動(dòng)臺(tái)三參量控制系統(tǒng)模型中進(jìn)行驗(yàn)證。首先，將參數(shù)整定的結(jié)果輸入測(cè)試系統(tǒng)的增益合成器中，采用10條與訓(xùn)練波不同的地震波加速度信號(hào)進(jìn)行測(cè)試，記錄測(cè)試波的相關(guān)系數(shù)ε和峰值誤差τ，若ε>95%且τ<10%，則認(rèn)為該組參數(shù)滿足要求，整定完畢，否則，執(zhí)行步驟(4)。

2 仿真對(duì)比與分析

2.1 模型參數(shù)選取

以某大學(xué)3 m×3 m水平雙向振動(dòng)臺(tái)的x向?yàn)檠芯繉?duì)象。圖5(a)為該振動(dòng)臺(tái)的控制系統(tǒng)，圖5(b)為該振動(dòng)臺(tái)的臺(tái)面系統(tǒng)。振動(dòng)臺(tái)的具體性能指標(biāo)如表1所示。

表1 3 m×3 m振動(dòng)臺(tái)性能指標(biāo)

圖5 3 m×3 m地震模擬振動(dòng)臺(tái)系統(tǒng)

首先根據(jù)專家調(diào)參經(jīng)驗(yàn)進(jìn)行手動(dòng)調(diào)參，然后輸入測(cè)試地震波進(jìn)行測(cè)試；接下來(lái)用文中所提出的基于LSTM的振動(dòng)臺(tái)控制參數(shù)整定算法進(jìn)行調(diào)參。LSTM選擇4層網(wǎng)絡(luò)：輸入層、隱藏層1、隱藏層2、輸出層，每個(gè)隱藏層取60個(gè)隱藏單元，初始學(xué)習(xí)率取0.01，訓(xùn)練次數(shù)設(shè)置為1 500次。

手動(dòng)調(diào)參采用白噪聲加速度信號(hào)，依據(jù)輸入輸出的頻譜特性對(duì)參數(shù)逐個(gè)調(diào)整；基于LSTM的振動(dòng)臺(tái)參數(shù)自整定采用5條天然地震波加速度信號(hào)進(jìn)行參數(shù)整定。兩種方法的參數(shù)整定結(jié)果如表2所示。

表2 振動(dòng)臺(tái)三參量控制參數(shù)整定結(jié)果

由表2可以看出：手動(dòng)參數(shù)整定值與自整定結(jié)果差異較大，其控制效果通過(guò)測(cè)試整定后的時(shí)頻域信號(hào)以及相關(guān)系數(shù)、峰值誤差等控制指標(biāo)對(duì)比后可知。

2.2 對(duì)比驗(yàn)證

此次仿真不考慮試件對(duì)振動(dòng)臺(tái)的影響，測(cè)試采用5條地震波加速度信號(hào)和1條白噪聲信號(hào)，采樣周期為0.005 s，運(yùn)行時(shí)間30 s，輸入信號(hào)的頻率為0.1～50 Hz。仿真得到的振動(dòng)臺(tái)系統(tǒng)輸入、輸出如圖6—9所示(這里僅展示2條測(cè)試波)。

圖6 白噪聲時(shí)程響應(yīng)曲線

測(cè)試波形的具體對(duì)比數(shù)據(jù)如表3所示，采用相關(guān)系數(shù)和峰值誤差兩個(gè)指標(biāo)來(lái)對(duì)比。

表3 LSTM整定值與手動(dòng)整定值性能對(duì)比

結(jié)合參數(shù)整定過(guò)程和圖6—9及表3可得以下結(jié)論：

(1)手動(dòng)參數(shù)整定和參數(shù)自整定均可使振動(dòng)臺(tái)系統(tǒng)復(fù)現(xiàn)地震波的大概走勢(shì)，但手動(dòng)整定需要技術(shù)人員反復(fù)調(diào)整參數(shù)，費(fèi)時(shí)費(fèi)力；參數(shù)自整定所需操作少，只需設(shè)置好訓(xùn)練次數(shù)和初始學(xué)習(xí)率即可，實(shí)測(cè)訓(xùn)練1 000次需要60 min左右，減少了整定時(shí)間，不需要技術(shù)人員反復(fù)調(diào)試，并且整定精度遠(yuǎn)好于手動(dòng)參數(shù)整定。

(2)從圖6和圖7的白噪聲對(duì)比曲線可看出：白噪聲覆蓋頻率為0.1～50 Hz，手動(dòng)整定的控制效果在頻率曲線圖的10～50 Hz有一些偏低，時(shí)程圖的峰值上明顯有一些跟不上；基于LSTM 的整定算法在時(shí)程的峰值上雖也略低，但是誤差很小，頻率響應(yīng)曲線也基本可以吻合。

圖7 白噪聲頻率響應(yīng)曲線

(3)從圖8可以看出，手動(dòng)參數(shù)整定的控制效果超調(diào)量大，響應(yīng)滯后；參數(shù)自整定的響應(yīng)曲線峰值和響應(yīng)都達(dá)到很高的精度。圖8的后半段波形都沒(méi)有完全吻合，對(duì)比手動(dòng)整定和自整定可明顯看出自整定的擬合效果更好。從圖9可以看出：自整定值在0.1～50 Hz全頻段都可以很好地?cái)M合目標(biāo)信號(hào)，不足之處在0.1～1 Hz左右有一點(diǎn)偏差；手動(dòng)整定在15～30 Hz頻段的擬合效果好，但在0.1～15 Hz和20～30 Hz頻段的擬合效果較差。

圖8 VAH-1波時(shí)程響應(yīng)曲線

圖9 VAH-1波頻率響應(yīng)曲線

(4)從表3可以明顯看出：手動(dòng)整定的相關(guān)系數(shù)可達(dá)90%左右，峰值誤差在10%～20%；參數(shù)自整定的相關(guān)系數(shù)除VAH-1波外都達(dá)到99.5%以上，峰值誤差都在7%以下，完全符合工程試驗(yàn)要求。

(5)相較于手動(dòng)參數(shù)整定，參數(shù)自整定無(wú)論是在時(shí)域響應(yīng)和頻譜上的表現(xiàn)，還是相關(guān)系數(shù)與峰值誤差等控制指標(biāo)的值都更優(yōu)，這說(shuō)明了基于LSTM參數(shù)自整定算法的優(yōu)越性。

3 結(jié)論

針對(duì)振動(dòng)臺(tái)控制參數(shù)整定問(wèn)題，提出一種基于LSTM的振動(dòng)臺(tái)控制參數(shù)整定算法，并基于Python進(jìn)行對(duì)比驗(yàn)證。根據(jù)仿真結(jié)果，得到如下結(jié)論：

(1)該參數(shù)整定算法需要的時(shí)間相比目前手動(dòng)調(diào)參少很多，并且不需要技術(shù)人員在旁反復(fù)調(diào)參，省去了大量人力物力。

(2)相比手動(dòng)調(diào)參的結(jié)果，該算法得到的控制參數(shù)超調(diào)量小，相關(guān)系數(shù)由平均90%提高到了平均99%，峰值誤差由平均16%降低到了4%，控制精度有很大的提高。

未來(lái)的研究應(yīng)對(duì)該算法進(jìn)行實(shí)機(jī)驗(yàn)證，且仍需要解決的問(wèn)題包括：考慮試件影響的振動(dòng)臺(tái)參數(shù)整定；多自由度振動(dòng)臺(tái)不同方向之間的耦合；振動(dòng)臺(tái)開(kāi)環(huán)模型的系統(tǒng)辨識(shí)等。以上問(wèn)題是下一步研究的重點(diǎn)。