霍慧 陳國海 王文培 楊迪雄

(大連理工大學(xué)工業(yè)裝備結(jié)構(gòu)分析國家重點實驗室,遼寧大連 116024)

引言

作為工程實際中常見的結(jié)構(gòu)形式,圓柱殼得益于其優(yōu)良的結(jié)構(gòu)特性和力學(xué)性能,被廣泛應(yīng)用于潛艇、火箭、導(dǎo)彈、飛機(jī)、儲液罐等工程結(jié)構(gòu)中[1-5].其中,考慮橫向剪切變形和轉(zhuǎn)動慣量影響的中厚圓柱殼結(jié)構(gòu)不可避免地會承受各類隨機(jī)激勵作用,如隨機(jī)地震動、風(fēng)載、噪聲激勵等[6-7].因此,中厚殼隨機(jī)振動分析和不確定性傳播研究對其設(shè)計和安全服役具有重要的理論意義及應(yīng)用價值.

自Love 建立殼體基本方程以來,已有許多學(xué)者和工程師對中厚圓柱殼結(jié)構(gòu)的自由振動特性進(jìn)行了分析[8-9].Leissa[10]對殼體自由振動理論進(jìn)行了深入廣泛的論述.然而,目前僅有簡支邊圓柱殼能獲得精確解析的自振頻率和振型函數(shù).為此,各種數(shù)值解法被相繼提出,如動力剛度法[11]、傅里葉級數(shù)法[12]、有限元法[13-14]、波函數(shù)法[15]等.

實際工程中,殼體結(jié)構(gòu)經(jīng)常承受環(huán)境動載荷作用.最近幾年,確定性動載荷作用下圓柱殼的動力響應(yīng)研究已經(jīng)得到了廣泛關(guān)注,載荷形式包括黏性流體[16]、移動載荷[17]、溫度或濕度場[18]、磁場作用[19]等.由于服役環(huán)境的復(fù)雜性,動力載荷普遍具有隨機(jī)性,通常采用隨機(jī)過程對其進(jìn)行表征.基于狀態(tài)協(xié)方差分配法,To和Chen[20]給出了非平穩(wěn)隨機(jī)激勵作用下壓電殼體元件的最優(yōu)隨機(jī)振動控制方案.針對邊界隨機(jī)激勵作用下的壓電空心圓柱厚殼,Ying等[21]基于伽遼金法計算了隨機(jī)響應(yīng)的統(tǒng)計矩.Esmailzadeh和Lakis[22]采用有限元方法求解了湍流邊界層激勵作用下的開口薄殼振動響應(yīng).基于FPK (Fokker-Planck-Kolmogorov)方程,Asnafi[23]獲得了隨機(jī)激勵作用下圓柱淺殼的概率密度函數(shù).Li 等[24]利用辛對偶方法研究了層合圓柱殼在軸向壓縮和湍流邊界層激勵下的隨機(jī)振動響應(yīng).

此外,Crandall和Elishkoff[25]以圓柱殼為研究對象,揭示了振型互相關(guān)項對隨機(jī)振動響應(yīng)具有不可忽視的作用.振型互相關(guān)項的重要作用也在球殼[26]、多跨梁[27]中被深入考察.對于離散多自由度系統(tǒng),林家浩等[28]提出了高效的虛擬激勵法,該方法精確計入振型互相關(guān)項.Chen 等[29]將虛擬激勵法推廣到連續(xù)體,獲得了薄殼平穩(wěn)隨機(jī)振動響應(yīng)的解析精確解.相比于薄殼,中厚殼在實際工程中的應(yīng)用更加廣泛,還會受到非平穩(wěn)隨機(jī)激勵作用.因此,本文針對中厚圓柱殼,開展平穩(wěn)、非平穩(wěn)隨機(jī)振動分析,以獲得其隨機(jī)響應(yīng)的精確基準(zhǔn)解.

為了保證結(jié)果的準(zhǔn)確性,本文首先通過自由振動分析,獲得了簡支中厚圓柱殼的精確頻率和解析的振型函數(shù).考慮點、線、面隨機(jī)激勵,構(gòu)造平穩(wěn)和非平穩(wěn)虛擬激勵,結(jié)合振型疊加法,將平穩(wěn)和非平穩(wěn)隨機(jī)振動分別轉(zhuǎn)化為簡諧振動和時程分析,推導(dǎo)了中厚圓柱殼隨機(jī)響應(yīng)功率譜的解析表達(dá)式.將功率譜在頻域上積分即可得到均方響應(yīng)及均方根.此外,在解析解的推導(dǎo)過程中還涉及對振型函數(shù)空間域積分和時間域的Duhamel 積分.基于符號運算的解析法在中厚殼隨機(jī)振動響應(yīng)解析求解中無法實現(xiàn)多空間點同時輸出,且計算難度和效率隨參振振型的增加而顯著增加.為此,進(jìn)一步提出了離散解析法.通過對空間域先解析積分后離散,頻域和時域數(shù)值積分,將解析運算轉(zhuǎn)化為矩陣運算,不僅顯著提高了計算效率,而且能批量獲得殼體隨機(jī)響應(yīng)的分布,便于研究參數(shù)變化對圓柱殼體隨機(jī)響應(yīng)全局影響.此外,基于精確基準(zhǔn)解,深入討論了隨機(jī)激勵類型及結(jié)構(gòu)參數(shù)對中厚殼振動響應(yīng)的影響.

1 中厚圓柱殼自由振動精確解

對于殼結(jié)構(gòu),隨著厚徑比(厚度h與曲率半徑R之比)的增大,剪切變形和轉(zhuǎn)動慣量對力學(xué)行為的影響越來越大,Kirchhoff 薄殼理論不再適用.針對中等厚度圓柱殼(如厚徑比≥1/36～ 1/20),分析應(yīng)采用考慮橫向剪切變形和轉(zhuǎn)動慣量影響的中厚殼理論.本文分別討論了如圖1(a)所示的封閉殼(φT=2π)和圖1(b)所示的開口角度φT的開口殼結(jié)構(gòu).R,L和h分別為殼的半徑、長度和厚度,坐標(biāo)系(x,θ,z),x,θ,z分別代表軸向、環(huán)向和徑向.引入5 個廣義獨立位移描述殼體的中面變形:軸向位移u1、環(huán)向位移u2、徑向位移u3、軸向轉(zhuǎn)角u4及環(huán)向轉(zhuǎn)角u5.考慮橫向剪切變形和轉(zhuǎn)動慣量的中厚圓柱殼隨機(jī)偏微分方程可以描述為[10]

圖1 中厚圓柱殼幾何模型及載荷情況Fig.1 Geometric model and load cases for moderately thick cylindrical shell

其中c為等效阻尼系數(shù),ρ為體密度,i=1,2,3 分別代表軸向、環(huán)向和徑向,j=4,5 分別代表軸向轉(zhuǎn)角和環(huán)向轉(zhuǎn)角方向,qi為圖1 所示施加在第i個方向上的隨機(jī)激勵,L1～L5為微分算子

其中,E和υ分別為材料彈性模量和泊松比,κ為剪切修正系數(shù),本文取為5/6[10].薄膜內(nèi)力Nx,Nθ,Nxθ和Nθx,彎矩Mx和Mθ,扭矩Mxθ和Mθx具體形式分別寫為[30]

式中,C=Eh/(1-υ2)表示拉伸剛度,D=Eh3/[12(1-υ2)]為彎曲剛度.

為了得到圓柱殼的自振頻率和振型,需進(jìn)行中厚圓柱殼的自由振動分析.去掉式(1)中的阻尼力項和外力項,即可得到中厚圓柱殼自由振動的微分方程

對于如圖1 (a)所示的簡支封閉圓柱殼,其自由振動精確解可設(shè)為分離變量形式[10]

其中,Us,mn(x,θ)為與第s個廣義位移相應(yīng)的第mn階振型函數(shù),s=1,2,3,4,5 分別代表軸向、環(huán)向、徑向、軸向轉(zhuǎn)角及環(huán)向轉(zhuǎn)角方向,m和n分別為軸向和環(huán)向的半波數(shù),αm=mπ/L,A,B,F,M和N為待求常數(shù).將式(7)代入式(6)可以得到

此外,由式(6)和式(7)可以得到矩陣形式的對稱齊次方程組

式(9)中的各個系數(shù)分別為

式中,Ω=ωR為待求的無量綱頻率,λ=mπR/L,k=h2/(12R2),Iκ=κ(1-υ)/2.為得到式(9)的非平凡解,令其系數(shù)矩陣行列式為0,則可以獲得關(guān)于Ω2的五次方的頻率方程

其中,a1-a5為代數(shù)方程的系數(shù).求解該方程從而獲得殼結(jié)構(gòu)的精確頻率ω.將求得的自振頻率回代到式(9),就可以最終確定式(7)中振型函數(shù)的系數(shù)A,B,F,M和N之間的比值關(guān)系.

類似地,對于如圖1 (b)所示的四邊簡支的開口形式中厚圓柱殼,其自由振動封閉形式解可以表達(dá)為如下形式

其中βn=nπ/φT.采用與封閉圓柱殼類似的計算步驟,將式(12)代入式(6),得到開口中厚圓柱殼自由振動的頻率方程,同樣可以獲得該結(jié)構(gòu)的精確頻率及解析振型函數(shù).

2 中厚圓柱殼隨機(jī)振動響應(yīng)的解析法

由振型疊加法可知,軸向、環(huán)向及徑向上的線彈性中厚圓柱殼位移可以表達(dá)為

其中ηmn(t)為第mn階振型的正則坐標(biāo).將式(13)代入式(1)可以得到

對于i=1,2,3,分別在式(14)兩端乘上第lk階振型Ui,lk,然后將i=1,2,3 的3 個方程在圓柱殼曲面內(nèi)進(jìn)行積分,綜合考慮式(8)及振型正交化,最終得到一系列解耦的單自由度系統(tǒng)

其中 ζmn=c/(2ρhωmn) 為第mn階的中厚殼模態(tài)阻尼比,

針對隨機(jī)激勵作用下中厚殼的振動分析,基于功率譜的傳統(tǒng)頻域分析方法采用維納-辛欽關(guān)系獲得響應(yīng)的功率譜密度.其中完全二次項組合(CQC)法能夠計及所有參振振型耦合項,但向量乘法運算量過大;而忽略振型耦合項的平方和開平方(SRSS)法計算誤差較大,且僅適用于阻尼比很小、參振頻率稀疏分布的均質(zhì)系統(tǒng).因此,本文利用虛擬激勵法和連續(xù)體結(jié)構(gòu)解析振型函數(shù),在保證結(jié)果準(zhǔn)確性的同時只需要計算一次向量乘法,計算效率得到顯著提高.

基于虛擬激勵法[28],構(gòu)造虛擬激勵

對于平穩(wěn)隨機(jī)激勵qi(x,θ,t)=Γi(x,θ)X(t),其中Γi(x,θ)為作用在第i個方向上的平穩(wěn)隨機(jī)過程X(t)的空間分布函數(shù),式(16)成為

其中SXX(ω) 為X(t)的功率譜密度(PSD)函數(shù).

將虛擬激勵式(17)代入式(15)可以得到

其中

因此,平穩(wěn)隨機(jī)振動分析轉(zhuǎn)化為確定性的簡諧振動分析,從而得到虛擬正則坐標(biāo)

除平穩(wěn)隨機(jī)激勵外,本文還考察統(tǒng)計信息隨時間變化的非平穩(wěn)隨機(jī)激勵.為描述非平穩(wěn)隨機(jī)激勵,Priestley 定義了具有明確物理意義的演變功率譜[31].若平穩(wěn)隨機(jī)過程X(t)表示為:非平穩(wěn)隨機(jī)過程Y(t)可進(jìn)一步描述為

其中,f(ω,t) 為確定性慢變的非均勻調(diào)制函數(shù).Y(t)的方差可表達(dá)為

進(jìn)而得到Y(jié)(t)的演變功率譜密度

然而,通常難以計算式(23)所示的非平穩(wěn)激勵下結(jié)構(gòu)隨機(jī)振動響應(yīng)的演變功率譜密度.當(dāng)激勵功率譜密度函數(shù)中的頻率特性不隨時間變化時,將非均勻調(diào)制函數(shù)f(ω,t)簡化近似為均勻調(diào)制函數(shù)f(t),此時的均勻調(diào)制非平穩(wěn)隨機(jī)過程可以寫為

對于施加在中厚殼第i個方向的均勻調(diào)制非平穩(wěn)激勵qi(x,θ,t)=Γi(x,θ)Y(t)=Γi(x,θ)f(t)X(t),基于虛擬激勵法的非平穩(wěn)虛擬激勵表示為

將式(25)描述的非平穩(wěn)虛擬激勵代入解耦的單自由度系統(tǒng)式(15),可得到

由Duhamel 時域積分,得到虛擬正則坐標(biāo)

其中hmn(t-τ) 為單位脈沖響應(yīng)函數(shù),即

基于平穩(wěn)及非平穩(wěn)隨機(jī)激勵作用下得到的虛擬正則坐標(biāo)如式(20)和式(27)所示,進(jìn)而中厚圓柱殼第i個方向的虛擬位移表示為

進(jìn)一步,可以基于虛擬激勵法得到各個方向位移、速度及加速度響應(yīng)功率譜密度函數(shù)如下

其中,上標(biāo)“*”表示復(fù)共軛.

此外,在虛擬位移響應(yīng)的基礎(chǔ)上,根據(jù)殼結(jié)構(gòu)的幾何關(guān)系及材料本構(gòu)關(guān)系,可以由式(4)和式(5)獲得虛擬薄膜內(nèi)力、彎矩及扭矩

同樣地,基于虛擬激勵法,將虛擬內(nèi)力響應(yīng)的復(fù)共軛乘上自己本身,可以解析獲得中厚圓柱殼薄膜內(nèi)力及彎曲內(nèi)力等待求響應(yīng)的功率譜密度函數(shù).

此節(jié)解析法推導(dǎo)過程沒有引入任何近似,解析地獲得了平穩(wěn)及非平穩(wěn)隨機(jī)激勵作用下中厚圓柱殼的各類響應(yīng)功率譜密度函數(shù)精確解.

3 中厚圓柱殼隨機(jī)振動響應(yīng)的離散解析法

在上一節(jié)的解析推導(dǎo)中,盡管采用了虛擬激勵法,解析解難以滿足求解隨機(jī)振動響應(yīng)多點輸出的需求,且計算難度和效率隨參振振型的增加而增長.為了在保證解析法空間域上精確性的前提下提高計算效率,本節(jié)基于空間域解析積分和微分運算后離散化策略,提出了中厚圓柱殼隨機(jī)振動響應(yīng)分析的離散解析法.

針對解析法推導(dǎo)過程中涉及到的振型函數(shù)在空間域的積分和微分運算(見式(15)、式(19)、式(31)和式(32)),離散解析法同樣采取了解析求解,即

由于引入了圓柱殼的封閉振型函數(shù)Ui,mn,如式(7)和式(12)所示,上述積分和微分均可獲得解析解.相比功率譜分析有限元法,采用空間域積分和微分解析運算的離散解析法能保證在空間域的精確性,具有數(shù)值方法無法比擬的優(yōu)勢.

此時,圓柱殼的隨機(jī)偏微分方程便轉(zhuǎn)化為求解式(18)和式(26)所示的常微分方程.在離散化系統(tǒng)中,虛擬激勵法的核心是在響應(yīng)功率譜密度計算中將向量乘法減少到一次,進(jìn)而提高計算效率.而在解析的功率譜密度公式(30)中,只涉及符號運算.

為此,將中厚殼結(jié)構(gòu)空間域進(jìn)行離散,各點位置坐標(biāo)寫為矩陣形式:r=[r1r2···rNG],rα=(xα,θα)T,α=1,2,···,NG,NG為空間離散點總數(shù),則式(13)描述的中厚殼位移公式可寫為向量形式

其中Ui(r)為第i個方向上的振型矩陣,矩陣維度為NG×NF,NF表示參振頻率階數(shù),i=1,2,3,η 代表NF維的正則坐標(biāo)列向量.

作用在中厚殼上的虛擬激勵可以離散化為

因此,對于平穩(wěn)隨機(jī)振動,第i個方向上中厚殼虛擬位移改寫為向量形式

其中,P(ω)為與式(19)對應(yīng)的廣義激勵的振幅矩陣,H(ω)為頻率響應(yīng)函數(shù)矩陣.由于離散解析法中針對振型函數(shù)積分的解析求解(見式(33)中的I1和I2),可知P(ω)在幾何空間域內(nèi)不涉及離散化操作.將多個節(jié)點坐標(biāo)代入式(36)中振型Ui(r)中,使得離散解析法實現(xiàn)了多空間點虛擬響應(yīng)的同時輸出,且虛擬位移響應(yīng)精度與空間離散點多寡無關(guān).

對于非平穩(wěn)隨機(jī)振動,由于Duhamel 積分解析求解通常耗時較長,所以采用高精度的精細(xì)積分法取代Duhamel 積分進(jìn)行時間歷程分析可提高計算效率.將解耦得到的單自由度系統(tǒng)式(26)改寫為狀態(tài)空間形式

式(37)求解過程中需要對時域離散化處理,對任意t∈[tk,tk + 1],采用精細(xì)積分法對式(37)進(jìn)行精確高效計算[32],其中tk=kΔt,k=0,1,2,···,Nt,Δt=T/(Nt),T為激勵持續(xù)時間,Nt為離散時間點數(shù).從而最終得到中厚圓柱殼在非平穩(wěn)隨機(jī)激勵作用下的虛擬位移響應(yīng)向量注意到,針對非平穩(wěn)隨機(jī)振動響應(yīng)分析的離散解析法同樣保證了在空間域的精確性.

基于虛擬激勵法,任一時刻N(yùn)G個空間點第i個方向的位移響應(yīng)功率譜密度矩陣為

需要說明的是,式(39)需要在一系列離散頻率點ωk上進(jìn)行計算,其中ωk=ωL+kΔω,Δω=(ωU-ωL)/(Nω),k=0,1,···,Nω,ωU和ωL為截斷頻率上下限,Nω為離散頻率點數(shù).

響應(yīng)功率譜密度是關(guān)于頻率的單變量光滑函數(shù),通過頻率帶寬離散化操作,對響應(yīng)功率譜密度簡單數(shù)值積分即可獲得響應(yīng)的均方根

同樣地,基于獲得的虛擬位移響應(yīng),亦可容易獲得其他待求響應(yīng)的功率譜密度矩陣及響應(yīng)均方根.需要強(qiáng)調(diào)的是,在基于有限元計算圓柱殼內(nèi)力功率譜密度時,內(nèi)力的計算需要對位移形函數(shù)求導(dǎo),使得連續(xù)性減低,導(dǎo)致內(nèi)力的計算精度比位移低,出現(xiàn)如線性三角形單元中常應(yīng)力情形.而離散解析法采用位移關(guān)于空間坐標(biāo)的偏導(dǎo)數(shù)先解析求導(dǎo),后將導(dǎo)函數(shù)離散化的方式,獲得各類內(nèi)力分布的精確解,避免了空間離散化對結(jié)果精度的影響,所得的內(nèi)力功率譜和均方根與位移是同等精確的.

因此,本節(jié)提出的離散解析法是對解析法的離散化求解,具有極高精度,可作為基準(zhǔn)解驗證中厚圓柱殼隨機(jī)振動分析的其他數(shù)值方法.

4 數(shù)值算例

4.1 中厚圓柱殼自由振動精確解

采用文獻(xiàn)[12]中的封閉中厚圓柱殼結(jié)構(gòu),彈性模量E=210 GPa,密度ρ=7800 kg/m3,長度L=0.502 m,半徑R=0.063 5 m,厚度h=0.003 26 m,泊松比υ=0.28.表1 展示了前5 階精確頻率.文獻(xiàn)[12]及有限元(FEM)所得自振頻率也列于表1.有限元結(jié)果基于ABAQUS 軟件,結(jié)構(gòu)模型分別劃分為500,1000,2000,4000,6000 及8000 個S4 R 單元(考慮縮減積分的4 節(jié)點彎曲薄殼或厚殼單元) 或S8 R 單元(考慮縮減積分的8 節(jié)點雙重彎曲厚殼單元).可見有限元計算精度受限于幾何網(wǎng)格疏密程度和單元類型選擇等因素,存在一定誤差,而文獻(xiàn)級數(shù)解[12]與本文所得精確頻率更加接近.

表1 封閉中厚殼的前5 階自振頻率(Hz)Table 1 The first 5 natural frequencies (Hz) of closed moderately thick cylindrical shells

除封閉中厚殼外,考慮一個開口角度為φT=π 的中厚圓柱殼,彈性模量E=27.466 GPa,密度ρ=7850 kg/m3,長度L=14.4 m,半徑R=3.6 m,厚度h=0.36 m,泊松比υ=0.2.圖2 分別給出了開口中厚殼的第1 階、第2 階、第3 階及第5 階自振頻率和振型,有限元結(jié)果由ABAQUS 軟件對結(jié)構(gòu)剖分為4000 個S8 R 單元得到.有限元結(jié)果與本文精確分析獲得的自由振動自振頻率及振型吻合良好.

圖2 開口中厚殼的自振頻率(Hz)及振型(精確解和有限元解)Fig.2 Natural frequencies (Hz) and modal functions of open moderately thick shells (exact solutions and FEM results)

4.2 隨機(jī)點激勵作用下中厚圓柱殼振動響應(yīng)

本節(jié)考慮如下封閉中厚圓柱殼結(jié)構(gòu):楊氏模量E=2.96 GPa、密度ρ=0.733 kg/m3、長度L=14.4 m、半徑R=3.6 m、泊松比υ=0.25、厚度h=0.36 m、阻尼比ζ=0.05.分別在兩個空間點 case 1:(L/2,0);case 2:(L/4,π/4)對中厚殼上施加徑向隨機(jī)點激勵.該隨機(jī)激勵為限帶白噪聲,其頻率范圍為[20,2000] Hz,功率譜密度為1 N2/Hz.

考慮這兩種隨機(jī)點激勵,利用解析法(analytical method,AM)、離散解析法(discrete analytical method,DAM)及ABAQUS 軟件計算,表2 分別給出了(L/2,0)處隨機(jī)響應(yīng)的均方根與各方法相應(yīng)的CPU 運行時間,其中ABAQUS 軟件結(jié)果由基于傳統(tǒng)功率譜方法的隨機(jī)振動分析模塊給出,中厚殼結(jié)構(gòu)被劃分為4000 (40 × 100)個S8 R 單元.為便于比較計算效率,基于離散解析法的中厚殼結(jié)構(gòu)同樣劃分為4000 (40 × 100) 個網(wǎng)格,空間離散點總數(shù)NG為4141 (41 × 101).三類方法均采用激勵帶寬內(nèi)的全部自振模態(tài)參與隨機(jī)振動分析.對本小節(jié)所考慮的中厚殼結(jié)構(gòu),頻帶[20,2000] Hz 內(nèi)存在14 階自振模態(tài),因此式(34)中離散解析法的正則坐標(biāo)η維度NF取為14.

由表2可知,兩類隨機(jī)點激勵作用下,離散解析法所得位移、加速度及彎矩Mx響應(yīng)及均方根與解析法完全吻合,而ABAQUS 軟件結(jié)果與解析法相比存在一定誤差.這驗證了時域、頻域數(shù)值積分的精度是非常高的.由于解析法單次只能計算一個空間點響應(yīng),為了公平比較三類方法的計算效率,將能夠多點運算的離散解析解及ABAQUS 軟件結(jié)果的CPU 運行時間均攤至各節(jié)點.可以看到,離散解析解法在保證計算精度的前提下,計算效率得到了極大提升,展示了該方法的計算優(yōu)越性.值得注意的是,基于空間域離散建模的ABAQUS 軟件結(jié)果計算精度依賴于結(jié)構(gòu)的離散形式,而離散解析法基于精確振型函數(shù)的解析積分或微分后再將節(jié)點坐標(biāo)代入式(39)中,其空間離散點多寡并不影響離散解析法的計算精度.

表2 徑向隨機(jī)點激勵作用下封閉中厚殼(L/2,0)處響應(yīng)均方根Table 2 Response RMSs of (L/2,0) for the closed moderately thick shells under random point excitation

圖3 給出了兩類隨機(jī)點激勵作用下封閉中厚殼徑向位移響應(yīng)功率譜密度曲線.除位移外,中厚圓柱殼的內(nèi)力響應(yīng)同樣值得關(guān)注.隨機(jī)點激勵(case 2)下相應(yīng)的結(jié)構(gòu)薄膜內(nèi)力Nx和彎矩Mx的響應(yīng)功率譜在圖4 中給出.結(jié)果表明,離散解析法和解析法得到的響應(yīng)功率譜密度在[20,2000] Hz 內(nèi)完全吻合,ABAQUS 軟件結(jié)果在低頻范圍內(nèi)與解析解吻合良好,但是隨著激勵頻率增加,ABAQUS 軟件結(jié)果呈現(xiàn)了一定的偏差,且ABAQUS 軟件內(nèi)力計算結(jié)果的精度低于位移計算結(jié)果.由于解析法結(jié)果采用了精確振型,保證了沒有錯根和漏根.因此,這種偏差顯然是由ABAQUS 軟件計算高階自振頻率的誤差引起的.

圖3 隨機(jī)點激勵下封閉中厚殼(L/2,0)處徑向位移響應(yīng)功率譜密度曲線Fig.3 Deflection response PSD curves of (L/2,0) for the closed moderately thick shells under random point excitation

圖4 隨機(jī)點激勵(case 2)下封閉中厚殼(L/2,0)處薄膜內(nèi)力Nx和彎矩Mx 功率譜密度曲線Fig.4 Response PSD curves of membrane force Nx and bending moment Mx of (L/2,0) for the closed moderately thick shells under random point excitation (case 2)

4.3 隨機(jī)環(huán)形線激勵作用下中厚圓柱殼振動響應(yīng)

本節(jié)對于上節(jié)中討論的半徑R=3.6 m 的封閉中厚圓柱殼,分別取厚度h為0.1 m,0.18 m,0.36 m,也即將厚徑比(h/R)分別取為1/36,1/20,1/10.考慮在封閉中厚殼的環(huán)線 (L/2,0～2π)上施加徑向隨機(jī)線激勵,該激勵同樣是限帶寬高斯白噪聲,頻帶為[20,2000] Hz,功率譜密度為1 N2/Hz,針對圓柱殼的三種厚徑比情形,采用激勵帶寬內(nèi)的全部振型參與隨機(jī)振動分析(NF分別取為46,27,14).針對上述三種厚徑比,表3 分別給出了中厚殼在(L/2,0)處的徑向位移、速度和加速度均方根.結(jié)果表明,在徑向隨機(jī)線激勵下,隨著中厚殼厚徑比的增加,結(jié)構(gòu)的抗彎剛度增大,從而導(dǎo)致徑向位移、速度及加速度響應(yīng)均方根減小.

表3 隨機(jī)環(huán)形線激勵作用下封閉中厚殼(L/2,0)處徑向響應(yīng)均方根Table 3 Radial response RMSs of (L/2,0) for the closed moderately thick shells under random annular line excitation

而且,也討論了徑向隨機(jī)線激勵下中厚殼結(jié)構(gòu)的厚徑比對彎矩均方根的影響.圖5 分別給出了厚徑比為1/36,1/20,1/10 下的中厚殼的彎矩Mx和Mθ的功率譜密度函數(shù).由式(5)可知,彎曲剛度D的表達(dá)式中包含厚度的三次方項.在虛擬激勵法中,彎矩響應(yīng)Mx和Mθ的功率譜密度可由相應(yīng)的虛擬彎矩與其自共軛相乘,即與D2呈正比.因此,隨著中厚殼厚徑比的增加,彎矩響應(yīng)Mx和Mθ的功率譜密度隨之增大.

圖5 隨機(jī)環(huán)形線激勵下中厚殼(L/2,0)處彎矩響應(yīng)功率譜密度曲線Fig.5 Bending moment response PSD curves of (L/2,0) for the closed moderately thick shells under random annular line excitation

4.4 分布隨機(jī)激勵下中厚圓柱殼振動響應(yīng)

針對4.1 小節(jié)討論的開口中厚圓柱殼結(jié)構(gòu),在其表面施加徑向分布的隨機(jī)激勵.考慮激勵具有如下兩類空間分布形式case 1:p1(x,θ)=x;case 2:p2(x,θ)=5θ.隨機(jī)激勵為[0,100] Hz 內(nèi)的限帶寬白噪聲,其功率譜密度Spp=1 Pa2/Hz,采用激勵帶寬內(nèi)的全部振型參與隨機(jī)振動分析(NF=32).阻尼比ζ=0.05.圖6 給出了這兩類分布激勵在開口中厚圓柱殼上的具體空間分布形式.

圖6 兩類隨機(jī)分布載荷的空間分布形式Fig.6 Spatial distribution forms of two random distributed loads

針對這兩類隨機(jī)分布載荷,圖7～圖9 分別給出了開口中厚殼徑向位移、速度及加速度響應(yīng)的均方根空間分布.對于分布形式case 1 的隨機(jī)載荷,其關(guān)于θ=π/2 的軸對稱,激勵峰值為x=L=14.4 m.由圖7 (a)、圖8 (a)及圖9 (a)可以看到,位移、速度及加速度響應(yīng)均方根峰值逐漸向激勵峰值方向移動,這些均方根關(guān)于θ=π/2 的軸是對稱的,并且位移、速度響應(yīng)均方根關(guān)于x=L/2 的軸呈現(xiàn)一定的對稱性,而加速度響應(yīng)均方根關(guān)于x=L/2 的軸沒有對稱性.對于隨機(jī)載荷case 2,其關(guān)于x=L/2 的軸對稱,激勵峰值為θ=φT=π.由圖7 (b)、圖8 (b)及圖9 (b)可知,位移、速度及加速度響應(yīng)均方根峰值逐漸向θ最大值位置移動,這些均方根關(guān)于x=L/2 的軸也是對稱,并且位移、速度響應(yīng)均方根關(guān)于θ=π/2 的軸呈現(xiàn)一定的對稱性,而加速度響應(yīng)均方根關(guān)于θ=π/2 的軸沒有對稱性.其原因在于,與其他空間位置相比,激勵峰值作用位置往往能夠激起結(jié)構(gòu)更高階振型.由式(30)可知,高階頻率對加速度響應(yīng)影響最大,其次是速度響應(yīng).由此可知,由于高階參振振型的影響,激勵峰值作用位置附近的加速度響應(yīng)均方根出現(xiàn)較大值,導(dǎo)致圖9 (a)、圖9 (b)所示的加速度響應(yīng)均方根分別關(guān)于x=L/2,θ=π/2 的軸呈現(xiàn)非對稱分布.分布隨機(jī)激勵的空間分布形式對中厚圓柱殼的加速度響應(yīng)均方根分布影響較大.由于隨著響應(yīng)均方根的增大,結(jié)構(gòu)響應(yīng)的變異性隨之增大.因此在實際工程中,為了保證中厚殼結(jié)構(gòu)在隨機(jī)環(huán)境激勵作用下的安全性,除徑向位移及速度響應(yīng)外,應(yīng)尤其關(guān)注結(jié)構(gòu)的加速度響應(yīng)分布.

圖7 兩類分布載荷作用下開口圓柱殼的位移響應(yīng)均方根分布Fig.7 RMS distributions of deflection for open shell under two distributed random loads

圖8 兩類分布載荷作用下開口殼的速度響應(yīng)均方根分布Fig.8 RMS distributions of velocity for open shell under two distributed random loads

圖9 兩類分布載荷作用下開口殼的加速度響應(yīng)均方根分布Fig.9 RMS distributions of acceleration for open shell under two distributed random loads

4.5 均勻調(diào)制非平穩(wěn)激勵下中厚圓柱殼振動響應(yīng)

除平穩(wěn)隨機(jī)激勵外,均勻調(diào)制非平穩(wěn)激勵下中厚圓柱殼的隨機(jī)振動響應(yīng)同樣值得重視.對于如下開口中厚圓柱殼(φT=π,E=27.466 GPa,ρ=7850 kg/m3,L=14.4 m,R=3.6 m,υ=0.2),考慮在開口殼上施加均勻調(diào)制非平穩(wěn)的隨機(jī)徑向均布載荷q3(x,θ,t)=f(t)X(t).隨機(jī)激勵X(t)假設(shè)為[0,100] Hz 內(nèi)的限帶寬高斯白噪聲激勵,激勵PSDSxx(ω)=1 Pa2/Hz,采用激勵帶寬內(nèi)的全部振型參與隨機(jī)振動分析.

采用如圖10 所示的時間調(diào)制函數(shù)f(t)[33]

圖10 時間調(diào)制函數(shù)f(t)Fig.10 Time-modulated function f(t)

其中,A為時間調(diào)制函數(shù)幅值,?(t)為Heaviside 函數(shù).本文中,取A=105,激勵持續(xù)時間T=9.6 s.可知,t=4.8 s 時,時間調(diào)制函數(shù)f(t)取最大值.

表4 給出了t=4.8 s 時開口殼中心點(L/2,π/2)處的徑向位移、速度及加速度響應(yīng)均方根,其中結(jié)構(gòu)厚徑比分別采用1/20和1/10,對應(yīng)的參振頻率階數(shù)NF分別為65和32.為了驗證結(jié)果,基于10 000次并行蒙特卡洛模擬(MCS)結(jié)果也被給出.利用譜表達(dá)方法,隨機(jī)激勵描述為

表4 非平穩(wěn)隨機(jī)均布激勵下中厚殼中心點響應(yīng)均方根 (t=4.8 s)Table 4 Stochastic response RMSs of the center of moderately thick shell under nonstationary random distributed excitations(t=4.8 s)

其中,Sqq(x,θ,t,ωk)為非平穩(wěn)隨機(jī)載荷在第k個頻率點的功率譜密度函數(shù),φk是[0,2π]內(nèi)獨立均勻分布的隨機(jī)相位角.對于頻帶[0,100] Hz 的離散頻率點數(shù)Nω,離散解析法和蒙特卡洛模擬均將其取為500 以保證公平比較計算效率.此外,若假設(shè)0～9.8 s 內(nèi)f(t)均等于A,則該均勻調(diào)制非平穩(wěn)激勵轉(zhuǎn)變?yōu)槠椒€(wěn)隨機(jī)激勵.表4 同時給出了中厚殼在該平穩(wěn)均布激勵下的解析解.可以看到,非平穩(wěn)隨機(jī)振動的離散解析解法在采用精細(xì)積分法進(jìn)行時程分析時,取了較大的時間步長(Δt=0.2 s),得到的計算結(jié)果與平穩(wěn)隨機(jī)響應(yīng)的解析解吻合良好,且計算效率很高.而蒙特卡洛模擬即使時間步長已經(jīng)足夠小,其計算誤差仍大于離散解析解,這主要歸因于蒙特卡洛的隨機(jī)收斂性.此外,從耗費的CPU 計算時間,也可看出,即使蒙特卡洛模擬采用了并行算法,計算效率仍遠(yuǎn)低于離散解析法.

非平穩(wěn)均布激勵下,圖11 分別繪制了厚徑比為1/20和1/10 的兩種情況下位移響應(yīng)的均方根曲線.可以看到,兩種方法的計算結(jié)果吻合良好,但蒙特卡洛模擬曲線存在一定數(shù)值振蕩.對于均勻調(diào)制非平穩(wěn)激勵下的中厚殼隨機(jī)振動響應(yīng),其徑向位移、速度及加速度響應(yīng)均方根隨時間的變化趨勢與圖10所示的時間調(diào)制函數(shù)趨勢一致.

圖11 非平穩(wěn)隨機(jī)激勵下中厚殼中心點位移響應(yīng)均方根Fig.11 Stochastic deflection responses RMSs of the center of moderately thick shell under nonstationary random excitations

圖11 非平穩(wěn)隨機(jī)激勵下中厚殼中心點位移響應(yīng)均方根(續(xù))Fig.11 Stochastic deflection responses RMSs of the center of moderately thick shell under nonstationary random excitations(continued)

為了討論開口中厚圓柱殼的厚徑比對位移響應(yīng)功率譜密度的影響,圖12 分別繪制了厚徑比h/R=1/36,1/20,1/15 下不同時刻(t=4.8 s,t=7 s)的中厚殼中心點處的位移響應(yīng)功率譜密度曲線.由圖12可知,徑向位移響應(yīng)的功率譜密度曲線的第一個峰值出現(xiàn)位置從約7.4 Hz 后移到約14.2 Hz.這是由于中厚圓柱殼的基頻隨著中厚殼厚徑比的增加而增大.這表明,在線性隨機(jī)振動分析中,結(jié)構(gòu)的自振頻率精確與否將直接影響響應(yīng)功率譜密度峰值出現(xiàn)的位置.此外,在徑向激勵作用下,中厚殼徑向位移響應(yīng)功率譜密度的幅值整體減小.因此,隨著結(jié)構(gòu)厚徑比的增加,徑向位移的均方根逐漸減小,也即響應(yīng)的變異性減小,圓柱殼的安全性得到提高.

圖12 非平穩(wěn)隨機(jī)激勵下中厚殼中心點位移響應(yīng)功率譜密度Fig.12 Stochastic deflection responses PSD curves of the center of moderately thick shell under nonstationary random excitations

圖12 非平穩(wěn)隨機(jī)激勵下中厚殼中心點位移響應(yīng)功率譜密度(續(xù))Fig.12 Stochastic deflection responses PSD curves of the center of moderately thick shell under nonstationary random excitations(continued)

5 結(jié)論

本文考察了考慮轉(zhuǎn)動慣量和剪切變形的中厚圓柱殼的隨機(jī)振動,獲得了平穩(wěn)、非平穩(wěn)隨機(jī)振動響應(yīng)的基準(zhǔn)解.通過引入中厚圓柱殼的精確自振頻率和封閉振型函數(shù),基于虛擬激勵法,解析推導(dǎo)了中厚圓柱殼在平穩(wěn)、非平穩(wěn)隨機(jī)激勵作用下的響應(yīng)功率譜密度函數(shù).為了充分發(fā)揮虛擬激勵法在矩陣運算中的高效性優(yōu)勢,將空間域積分解析求解,頻域和時域積分?jǐn)?shù)值求解,提出了離散解析法高效計算了中厚殼在各類隨機(jī)激勵下響應(yīng)均方根的基準(zhǔn)解.通過與解析解、ABAQUS 軟件結(jié)果及蒙特卡洛模擬結(jié)果的比較,展示了離散解析法的準(zhǔn)確性和高效性.通過算例分析,得到以下主要結(jié)論:

(1) 徑向隨機(jī)線激勵下,隨著圓柱殼厚度的增加,結(jié)構(gòu)徑向抗壓強(qiáng)度增大,位移、速度及加速度響應(yīng)均方根等比例降低;結(jié)構(gòu)彎曲剛度增大,導(dǎo)致彎矩響應(yīng)均方根增加.

(2) 分布隨機(jī)激勵的空間分布形式對中厚圓柱殼的加速度響應(yīng)均方根分布影響較大.為保證結(jié)構(gòu)安全性,中厚殼結(jié)構(gòu)的加速度響應(yīng)也應(yīng)受到關(guān)注.

(3) 非平穩(wěn)隨機(jī)激勵下,中厚圓柱殼位移響應(yīng)均方根隨時間變化與時間調(diào)制函數(shù)相一致.此外,隨著圓柱殼厚度的增加,位移響應(yīng)功率譜的幅值整體呈減小趨勢,相應(yīng)的均方根(即變異性)減小,因此結(jié)構(gòu)安全性得到提高.