楊 格， 孫紅碩， 吳 斌， 潘天林， 王 貞

(1.武漢理工大學(xué) 土木工程與建筑學(xué)院，武漢 430070;2.天津大學(xué) 中國(guó)地震局地震工程綜合模擬與城鄉(xiāng)抗震韌性重點(diǎn)實(shí)驗(yàn)室，天津 300350；3.東北電力大學(xué) 建筑工程學(xué)院，吉林 132012)

斜拉索因其結(jié)構(gòu)柔、質(zhì)量輕，在風(fēng)荷載等外力作用下非常容易發(fā)生振動(dòng)。由于拉索振動(dòng)易導(dǎo)致拉索錨固端發(fā)生疲勞破壞等問題，需要采取措施來抑制其大幅度振動(dòng)，工程中廣泛采用的措施是在拉索錨固位置附近安裝阻尼器。為了研究阻尼器對(duì)拉索的減振效果，一般采用數(shù)值模擬[1]或室內(nèi)試驗(yàn)[2]的方式來獲取拉索-阻尼器系統(tǒng)的動(dòng)力響應(yīng)。數(shù)值模擬應(yīng)用方便，但現(xiàn)有的阻尼器數(shù)值模型均存在一定的簡(jiǎn)化，不能真實(shí)反應(yīng)阻尼器在實(shí)際結(jié)構(gòu)中的性能。室內(nèi)試驗(yàn)因試驗(yàn)場(chǎng)地規(guī)模和技術(shù)等問題，很難對(duì)較長(zhǎng)拉索進(jìn)行真實(shí)工況下的模擬。

結(jié)構(gòu)混合試驗(yàn)方法可以將結(jié)構(gòu)中無法準(zhǔn)確模擬的部分作為試驗(yàn)子結(jié)構(gòu)，其他部分采用數(shù)值模型，兩者在線交互完成對(duì)整體結(jié)構(gòu)的模擬。該方法解決了因?qū)嶒?yàn)室場(chǎng)地規(guī)模限制等因素而無法開展大型結(jié)構(gòu)試驗(yàn)的問題，被廣泛應(yīng)用于結(jié)構(gòu)抗震試驗(yàn)研究中[3-6]。其中，實(shí)時(shí)混合試驗(yàn)方法能夠考慮加載速度的影響，可實(shí)現(xiàn)對(duì)速度相關(guān)型試件力學(xué)性能的準(zhǔn)確反映。鑒于實(shí)時(shí)混合試驗(yàn)的這一優(yōu)點(diǎn)，可以將拉索作為數(shù)值子結(jié)構(gòu)，阻尼器作為試驗(yàn)子結(jié)構(gòu)，開展拉索-阻尼器系統(tǒng)的實(shí)時(shí)混合試驗(yàn)。該方法對(duì)研究拉索-阻尼器系統(tǒng)中阻尼器的減振效果有著重要意義。

實(shí)時(shí)混合試驗(yàn)對(duì)數(shù)值模型的計(jì)算效率有著極高的要求。為了保證計(jì)算效率，現(xiàn)有關(guān)于拉索-阻尼器實(shí)時(shí)混合試驗(yàn)的研究大多采用簡(jiǎn)化的數(shù)值模型來模擬拉索。Jung等[7]對(duì)拉索-阻尼器系統(tǒng)在不同風(fēng)速下的響應(yīng)進(jìn)行了實(shí)時(shí)混合試驗(yàn)，其拉索數(shù)值模型采用伽遼金方法得到。伽遼金方法是一種簡(jiǎn)化計(jì)算方法，其在拋物線假設(shè)的基礎(chǔ)上引入一個(gè)參數(shù)來考慮索垂度以減少計(jì)算量；而傳統(tǒng)有限元分析是通過自重作用下索的形狀來得到垂度，因此伽遼金方法計(jì)算精度要比傳統(tǒng)有限元分析計(jì)算精度差。Duan等[8]建立了基于向量式有限元的拉索模型，阻尼器采用線性黏滯阻尼器數(shù)值模型來代替，對(duì)拉索-阻尼器系統(tǒng)的一階模態(tài)振動(dòng)進(jìn)行了實(shí)時(shí)混合試驗(yàn)數(shù)值仿真。相對(duì)于傳統(tǒng)有限元方法，基于向量式有限元的拉索模型具有更高的計(jì)算效率，但其在進(jìn)行質(zhì)點(diǎn)運(yùn)動(dòng)方程求解時(shí)，采用了顯式的中心差分法。中心差分法是一種對(duì)線性體系有條件穩(wěn)定的算法，而拉索-阻尼器系統(tǒng)屬于強(qiáng)非線性系統(tǒng)，一方面拉索振動(dòng)具有較強(qiáng)的幾何非線性;另一方面阻尼器也往往具有很強(qiáng)的非線性，因此無法保證中心差分法在拉索-阻尼器系統(tǒng)動(dòng)力分析中的穩(wěn)定性。當(dāng)拉索存在高頻振動(dòng)時(shí)，通常需要很小的時(shí)間步長(zhǎng)來滿足算法的計(jì)算精度和穩(wěn)定性要求[9]，仍然存在計(jì)算時(shí)間可能超過時(shí)間步長(zhǎng)而導(dǎo)致試驗(yàn)失敗的問題。

因此，在拉索-阻尼器系統(tǒng)的實(shí)時(shí)混合試驗(yàn)中需要對(duì)非線性體系無條件穩(wěn)定的時(shí)間積分方法。然而，目前混合試驗(yàn)中常用的無條件穩(wěn)定時(shí)間積分方法主要是針對(duì)線性體系，如無條件穩(wěn)定的顯式CR方法(Chen and Ricles)[10-12]、顯式Chang方法族[13]、隱式平均加速度法[14]、隱式中點(diǎn)法[15]等。對(duì)于非線性體系，Crisfield等[16]采用平均加速度法求解桁架單元的動(dòng)力方程時(shí)，發(fā)現(xiàn)結(jié)果會(huì)出現(xiàn)發(fā)散。潘天林[17]通過對(duì)具有幾何非線性的桁架體系進(jìn)行分析，發(fā)現(xiàn)隱式平均加速法和隱式中點(diǎn)方法不能保證無條件穩(wěn)定。Kuhl等[18]發(fā)現(xiàn)基于桁架單元的能量耗散積分方法也會(huì)存在計(jì)算結(jié)果不穩(wěn)定的現(xiàn)象。為了實(shí)現(xiàn)對(duì)非線性體系的無條件穩(wěn)定，Wu等[19]提出了能量一致積分方法，并將基于梁柱單元的能量一致積分方法應(yīng)用到足尺鋼框架結(jié)構(gòu)的混合試驗(yàn)中。由于能量一致積分方法為隱式算法需要迭代，目前該積分方法方法尚未應(yīng)用于實(shí)時(shí)混合試驗(yàn)中。考慮到近年來混合試驗(yàn)為應(yīng)對(duì)速度相關(guān)型試件的精細(xì)化模擬需求，逐漸趨于有限元化和實(shí)時(shí)化[20]，因此，將能量一致積分方法應(yīng)用于實(shí)時(shí)混合試驗(yàn)中具有重要意義。

能量一致積分方法是一種隱式方法，將其應(yīng)用于實(shí)時(shí)混合試驗(yàn)時(shí)，會(huì)遇到迭代導(dǎo)致作動(dòng)器加載速度波動(dòng)較大的問題。對(duì)于阻尼器試件，加載速度波動(dòng)較大會(huì)導(dǎo)致測(cè)得的阻尼器出力嚴(yán)重失真，致使試驗(yàn)失敗。為此，本文提出了基于能量一致積分的拉索-阻尼器系統(tǒng)實(shí)時(shí)混合試驗(yàn)方法，一方面通過Jung等[21]提出的固定迭代次數(shù)并對(duì)迭代位移進(jìn)行插值的方式來求解隱式差分方程，實(shí)現(xiàn)平滑加載；另一方面為保證試驗(yàn)過程中拉索-阻尼器系統(tǒng)的能量一致，對(duì)試驗(yàn)測(cè)得的阻尼器出力進(jìn)行恢復(fù)力修正。在不考慮拉索抗彎剛度的情況下，拉索可由若干個(gè)桁架單元模擬，本文將基于桁架單元的能量一致積分方法應(yīng)用于拉索-阻尼器系統(tǒng)的實(shí)時(shí)混合試驗(yàn)中，對(duì)拉索-阻尼器系統(tǒng)進(jìn)行了一階模態(tài)振動(dòng)下的實(shí)時(shí)混合試驗(yàn)數(shù)值仿真，驗(yàn)證了方法的可行性。

1 基于桁架單元的能量一致積分方法基本原理

對(duì)于桁架單元，連續(xù)的動(dòng)力方程可以表示為

(1)

式中,m、v、r和f分別為質(zhì)量矩陣、速度向量、節(jié)點(diǎn)力向量和外荷載向量。其中節(jié)點(diǎn)力向量r可表示為

r=N[-ee1]T

(2)

式中:N為單元的軸力;e1為桁架單元的軸向單位向量

(3)

式中:x=xQ-xP為節(jié)點(diǎn)坐標(biāo)差，xP與xQ分別為桁架兩個(gè)節(jié)點(diǎn)的坐標(biāo)向量;L為單元的長(zhǎng)度。

桁架單元的節(jié)點(diǎn)力考慮了幾何非線性和材料非線性的影響，其中幾何非線性由向量e1考慮，材料非線性則由N考慮。為了實(shí)現(xiàn)單元的能量一致，需要對(duì)恢復(fù)力進(jìn)行非線性修正，其修正格式為

(4)

ri+β=r[(1-β)ui+βui+1]，

ri+1-β=r[βui+(1-β)ui+1]

(5)

式中：β為修正系數(shù)；ui、ui+1分別為第i、第i+1個(gè)積分點(diǎn)時(shí)刻的位移向量。然后分別對(duì)式(2)中兩項(xiàng)非線性相關(guān)向量進(jìn)行離散，對(duì)e1離散得到

(6)

式中,β1為幾何非線性修正系數(shù)。對(duì)N離散得到

(7)

式中,β2為材料非線性修正系數(shù)。能量一致積分方法采用平均加速度法對(duì)位移和速度的假設(shè)

(8)

(9)

式中：Δt為時(shí)間積分步長(zhǎng)；vi、vi+1分別為第i、i+1個(gè)積分點(diǎn)時(shí)刻速度向量；ai、ai+1分別為第i、第i+1個(gè)積分點(diǎn)時(shí)刻加速度向量。

基于此，式(1)的離散形式可表示為

(10)

式中:fi、fi+1分別為第i、第i+1個(gè)積分點(diǎn)時(shí)刻外荷載向量;εi、εi+1分別為第i、第i+1個(gè)積分點(diǎn)時(shí)刻的應(yīng)變;A0為桁架單元原截面面積;L0為桁架單元原長(zhǎng);σ為工程應(yīng)力;第三式和第四式分別為關(guān)于β1和β2的非線性方程。由于能量一致積分方法為隱式算法，整體節(jié)點(diǎn)力方程的求解通過牛頓迭代實(shí)現(xiàn)。

2 基于能量一致積分的拉索-阻尼器實(shí)時(shí)混合試驗(yàn)方法

2.1 基于能量一致積分的拉索-阻尼器實(shí)時(shí)混合試驗(yàn)方法基本原理

在拉索-阻尼器系統(tǒng)的實(shí)時(shí)混合試驗(yàn)中，以拉索作為數(shù)值子結(jié)構(gòu)，阻尼器作為試驗(yàn)子結(jié)構(gòu)，連續(xù)的動(dòng)力方程可以表示為

(11)

式中:rN為數(shù)值子結(jié)構(gòu)恢復(fù)力向量；rE為試驗(yàn)子結(jié)構(gòu)恢復(fù)力向量，即阻尼器的出力。采用基于桁架單元的能量一致積分方法，式(11)可以離散為

(12)

式中:nN為桁架單元數(shù)目;βN1、βN2分別為幾何非線性修正系數(shù)和材料非線性性系數(shù)。本文假定拉索材料為線彈性，故第四式中所有桁架單元的材料非線性修正系數(shù)βN2=0。由于測(cè)得的試驗(yàn)子結(jié)構(gòu)恢復(fù)力rE不是連續(xù)的，故修正后的等效恢復(fù)力rEβ沒有理論解，需根據(jù)試驗(yàn)數(shù)據(jù)結(jié)合第五式計(jì)算。

(13)

(14)

式中:k=0，1，…，n-1,n為每一積分步固定迭代次數(shù),n=Δt/δt，Δt為積分步長(zhǎng)，δt為采樣步長(zhǎng)；系數(shù)m=(k+1)/n；v0、u0分別為初始速度向量和初始位移向量。在計(jì)算第一個(gè)積分步時(shí)，因?yàn)闆]有上一積分點(diǎn)時(shí)刻的位移，位移命令可通過式(14)獲得。

在純時(shí)滯與總時(shí)滯相比非常小的情況下，作動(dòng)器系統(tǒng)的動(dòng)力性能可用不包含純時(shí)滯環(huán)節(jié)的二階傳遞函數(shù)來代替，傳遞函數(shù)可寫為

(15)

式中：s為L(zhǎng)aplace變量；ξ為作動(dòng)器系統(tǒng)的阻尼比；ω為作動(dòng)器系統(tǒng)的頻率。時(shí)滯會(huì)影響實(shí)時(shí)混合試驗(yàn)仿真精度[22]，但在阻尼器作為試驗(yàn)子結(jié)構(gòu)的實(shí)時(shí)混合試驗(yàn)中，時(shí)滯對(duì)仿真精度影響一般較小，故本文不進(jìn)行作動(dòng)器的時(shí)滯補(bǔ)償。

基于能量一致積分的拉索-阻尼器系統(tǒng)實(shí)時(shí)混合試驗(yàn)實(shí)施過程,如圖1所示。第i+1個(gè)積分步的第k次迭代實(shí)施過程可簡(jiǎn)述為：

圖1 拉索-阻尼器系統(tǒng)實(shí)時(shí)混合試驗(yàn)實(shí)施過程示意

重復(fù)上述步驟，直至達(dá)到預(yù)定的固定迭代次數(shù)n，完成第i+1個(gè)積分步的模擬，并依此實(shí)現(xiàn)0～t時(shí)間內(nèi)的仿真。

2.2 恢復(fù)力修正模塊計(jì)算方法

為保證與系統(tǒng)能量一致，需要通過恢復(fù)力修正模塊對(duì)測(cè)得的試驗(yàn)子結(jié)構(gòu)恢復(fù)力進(jìn)行修正，計(jì)算等效恢復(fù)力rEβ。在拉索-阻尼器系統(tǒng)中，第i+1個(gè)積分步中前j次迭代阻尼器出力所做的實(shí)際總功為

(16)

(17)

(18)

通過式(18)可以看出，分子是第i+1個(gè)積分步中前j次迭代的阻尼器出力所做總功的近似值，分母是第i+1個(gè)積分步中前j次迭代的總位移增量。通過式(18)對(duì)測(cè)得的阻尼器出力進(jìn)行修正獲得等效恢復(fù)力，實(shí)現(xiàn)整個(gè)過程系統(tǒng)能量一致。

3 拉索-阻尼器系統(tǒng)實(shí)時(shí)混合模擬

3.1 模型的建立

表1 J26號(hào)拉索結(jié)構(gòu)參數(shù)

圖2 拉索-阻尼器系統(tǒng)模型示意圖

3.2 拉索模型的幾何非線性和正確性驗(yàn)證

對(duì)拉索模型跨中處施加一個(gè)周期的正弦位移，幅值取1 m。通過靜力分析可得拉索跨中節(jié)點(diǎn)恢復(fù)力-位移曲線如圖3所示。可以看出節(jié)點(diǎn)恢復(fù)力與位移并非呈線性關(guān)系，驗(yàn)證了拉索模型的幾何非線性。

圖3 拉索跨中節(jié)點(diǎn)恢復(fù)力-位移曲線

記錄模擬得到的拉索1/4跨處位移時(shí)程，最大振動(dòng)幅值約為2.5 mm。通過圖3可以看出，在小位移下拉索模型的恢復(fù)力-位移曲線基本處于線性階段，此時(shí)可忽略幾何非線性的影響。隨后對(duì)該位移時(shí)程進(jìn)行快速傅里葉變換，可得拉索前三階模態(tài)頻率如圖4所示。與Duan等的研究中基于ANSYS軟件分析得到的拉索前3階模態(tài)頻率對(duì)比，可看出兩者所得前3階模態(tài)頻率基本一致，驗(yàn)證了本文拉索模型的正確性,如表2所示。

圖4 位移功率譜密度函數(shù)

表2 拉索模態(tài)頻率對(duì)比

3.3 實(shí)時(shí)混合試驗(yàn)的數(shù)值仿真

本文假定作動(dòng)器系統(tǒng)頻率和阻尼比分別為ω=3.14 rad/s、ξ=0.8[24],以式(15)所示傳遞函數(shù)來模擬實(shí)際作動(dòng)器系統(tǒng)。在積分步長(zhǎng)Δt=0.05 s下，取固定迭代次數(shù)n=15，對(duì)拉索-阻尼器系統(tǒng)進(jìn)行實(shí)時(shí)混合試驗(yàn)數(shù)值仿真。為了盡可能真實(shí)地再現(xiàn)實(shí)際試驗(yàn)的情況，阻尼器出力采用迭代點(diǎn)時(shí)刻的實(shí)際速度來計(jì)算，在第i+1個(gè)積分步第k個(gè)迭代點(diǎn)時(shí)刻速度為

(19)

(20)

將頻率為拉索第1階模態(tài)頻率的正弦荷載300×sin(2π×0.5t)N施加于拉索的23個(gè)結(jié)點(diǎn)上，荷載方向均垂直于拉索，激勵(lì)100 s后釋放荷載。在阻尼器作用下，振動(dòng)自由衰減，可獲得前150 s的拉索跨中位移時(shí)程，其仿真結(jié)果如圖5和圖6所示。從圖5、圖6可以看出，本文所提實(shí)時(shí)混合試驗(yàn)方法的仿真結(jié)果與整體模擬的數(shù)值解吻合較好。計(jì)算得到拉索跨中位移的均方根誤差為1.96%，滿足精度要求。

圖5 拉索跨中位移時(shí)程曲線對(duì)比(1階模態(tài))

通過圖6可以看出，本文方法計(jì)算幅值比數(shù)值解偏大。因此，本文方法的計(jì)算結(jié)果偏于保守，有利于保證結(jié)構(gòu)的安全性。

在計(jì)算效率方面，該實(shí)時(shí)混合試驗(yàn)數(shù)值仿真時(shí)長(zhǎng)為150 s，計(jì)算總耗時(shí)為9.45 s，每一積分步平均計(jì)算耗時(shí)0.003 15 s?？梢钥闯霰疚乃岱椒ㄓ?jì)算效率較高，滿足實(shí)時(shí)要求。其中，MATLAB軟件版本為MATLAB2017b，計(jì)算機(jī)配置如表3所示。

表3 計(jì)算機(jī)配置

為了觀察基于能量一致積分的拉索-阻尼器系統(tǒng)實(shí)時(shí)混合試驗(yàn)方法的迭代效果，取第二個(gè)積分步發(fā)給作動(dòng)器模型的位移命令，如圖7所示。可以看出發(fā)給作動(dòng)器的位移命令時(shí)程曲線較光滑，因此本文所提方法可保證對(duì)阻尼器試件的平滑加載。

該仿真中阻尼器出力與速度關(guān)系曲線如圖8所示。阻尼器出力與速度關(guān)系曲線滿足本文阻尼器力學(xué)模型。

4 結(jié) 論

本文提出了基于能量一致積分的實(shí)時(shí)混合試驗(yàn)方法，并應(yīng)用于拉索-阻尼器系統(tǒng)的實(shí)時(shí)混合試驗(yàn)仿真中，主要結(jié)論如下：

(1) 基于桁架單元的能量一致積分方法可考慮拉索的幾何非線性，實(shí)現(xiàn)對(duì)拉索-阻尼器系統(tǒng)的振動(dòng)過程數(shù)值模擬。

(2) 基于能量一致積分的實(shí)時(shí)混合試驗(yàn)方法可以實(shí)現(xiàn)速度相關(guān)型試件的光滑加載，可應(yīng)用于拉索-阻尼器系統(tǒng)的實(shí)時(shí)混合試驗(yàn)中。