謝皓宇，潘飛，朱翊洲，鄭萬山，仉文崗

(1. 橋梁工程結(jié)構(gòu)動力學國家重點實驗室，重慶 400067；2. 招商局重慶交通科研設(shè)計院有限公司，重慶 400067；3.重慶大學 土木工程學院，重慶 400045；4. 上海核工程研究設(shè)計院有限公司，上海 200233)

地震災(zāi)害對核電站造成的破壞會導(dǎo)致非常嚴重的后果，放射性物質(zhì)泄漏對于人類和自然環(huán)境都是無法承受的。鑒于此，世界各國對于核電廠及核電廠樓面設(shè)備的抗震可靠度都給予了高度重視?；谛阅艿目拐鹪O(shè)計理念要求對目標結(jié)構(gòu)的有限元模型進行時程動力分析或符合相似原則的模型振動臺試驗，因此，設(shè)計地震動的選擇被認為至關(guān)重要[1-2]。對于核電廠設(shè)備的設(shè)計地震動，中國的《核設(shè)備抗震鑒定試驗指南》(HAF·J0053)[3]、《核電廠抗震設(shè)計標準》(GB 50267—2019)[4]以及美國的NUREG-0800號《標準審查大綱(SRP)》[5-6]、IEEE 344-2013[7]均有相關(guān)規(guī)定，即輸入地震波的頻域信息在匹配典型反應(yīng)譜(Representative Response Spectrum)的同時，還需要在指定的頻域范圍內(nèi)(0.3～24 Hz)包絡(luò)標準功率譜密度(Power Spectral Density)曲線的70%或80%。然而，無論是HAF·J0053或IEEE 344-2013，對于核設(shè)備標準PSD的規(guī)定都停留在定性的層面而非定量的描述。因此，針對核電廠設(shè)備的抗震設(shè)計，目標PSD包絡(luò)的要求通常被列為低于RRS匹配的次要要求，在2014版SRP頒布之前都并未得到足夠重視。孫渝剛等[8]討論了基于2014版SRP推薦的核電廠廠址抗震設(shè)計的標準PSD擬合程序，但并未針對核設(shè)備的RRS做進一步的優(yōu)化，趙鳳新等[9]、張郁山等[10]則對包絡(luò)標準PSD的人工波擬合方法做出了基于時域的研究，但他們提出的算法中需要根據(jù)具體的譜型調(diào)整的參數(shù)較多，且也未針對核設(shè)備的RRS做出優(yōu)化。此外，對于同時滿足RRS匹配及PSD包絡(luò)的輸入地震動和只有RRS匹配的輸入地震動，時程動力分析中核設(shè)備被激發(fā)的實際響應(yīng)特性會有何區(qū)別，目前也缺乏相關(guān)試驗研究與效果驗證。

筆者基于2014版SRP建議的生成標準PSD的方法，做出了提升迭代效率的修改，并在包絡(luò)PSD的情況下使用傳統(tǒng)的頻域法擬合了匹配RRS的人工地震動，再以典型的核設(shè)備壓力容器的模型作為研究對象進行振動臺試驗，對安全停堆地震(SSE)的RRS匹配及PSD包絡(luò)的人工波性能進行驗證，以期為核電廠設(shè)備動力分析輸入地震動的選擇提供更可靠的依據(jù)。

1 RRS匹配和PSD包絡(luò)

1.1 標準PSD的生成

SRP 3.7.1附錄B.1中建議使用一種對RRS頻域控制點逐點迭代的方法生成標準PSD。這個方法的步驟可以簡要地概括：1)對于作為目標的RRS用RSrep表示，使用NUREG/CR-6728時程數(shù)據(jù)庫中的平均PSD作為迭代的初始PSD；2)對于第M次迭代，以M-1步迭代生成的PSD為基礎(chǔ)通過快速傅里葉逆變換及隨機相位譜生成10×M條人工地震波；3)將生成的人工地震波時程乘以一個包含上升段、平臺段和下降段的時程包絡(luò)函數(shù)；4)計算10×M條人工地震波所對應(yīng)5%阻尼比的反應(yīng)譜，并求出他們的算數(shù)平均RSavg；5)將M-1步迭代生成的PSD逐個頻率控制點都乘以(RSrep/RSavg)2，再使用新生成的PSD作為下一次迭代的基礎(chǔ)，直到滿足規(guī)定的收斂條件。孫渝剛等[11]在研究中已經(jīng)實現(xiàn)了這個方法。本質(zhì)上這個方法是求出了匹配RRS的人工波的一個“均值”，并要求輸入地震動的PSD高于這個均值的PSD。

SRP建議的迭代方法中涉及到頻域法擬合人工地震波的內(nèi)容。為了優(yōu)化該方法、提升迭代效率，參考胡聿賢等[12]、謝皓宇等[13]的研究工作，提出考慮迭代相關(guān)及隨機相位譜對迭代收斂效率的影響，修正上一段中SRP方法第5)步的迭代公式，具體內(nèi)容包括：1)對于頻率ω的線性振子，相同頻率的傅里葉分量激勵下的響應(yīng)加速度aω(ω)與人工波激勵下的響應(yīng)加速度ag(ω)在ag(ω)峰值出現(xiàn)的時刻t反向，那么將功率譜密度迭代的指數(shù)修改為負數(shù)；2)當人工波的平均反應(yīng)譜大于RRS并且t時刻aω(ω)與ag(ω)同向，或人工波平均反應(yīng)譜小于RRS并且t時刻aω(ω)與ag(ω)反向，那么將頻率ω對應(yīng)的隨機相位加上π。可以用公式來描述以上兩項修改。

(1)

(2)

1.2 人工波擬合

核設(shè)備設(shè)計地震動的擬合應(yīng)以RRS匹配和包絡(luò)標準PSD的80%為目標。前述方法生成的標準PSD是基于RRS的，即標準PSD與RRS所對應(yīng)的人工時程在頻域上的特性及包含的信息相似。因此，初始時程加速度xinitial直接由標準PSD所計算出的傅里葉幅值譜和隨機的相位譜通過傅里葉逆變換之后乘以一個時程包絡(luò)函數(shù)I(t)得到[14]。

(3)

(4)

式中：Sx(ω)為標準PSD；A(ω)為標準PSD對應(yīng)的傅里葉幅值譜；Δω為頻域采樣間隔；φ(ω)為隨機的傅里葉相位譜，隨機相位譜由馬特賽特回旋算法(Mersenne Twister)[15]生成。

由于PSD包絡(luò)的目標，需要修正傳統(tǒng)頻域法的迭代過程，令迭代過程中的傅里葉相位譜始終高于0.8×A(ω)，即

An(ω)>0.8·A(ω)→An(ω)=An(ω)

(5)

An(ω)<0.8·A(ω)→An(ω)=0.8·A(ω)

(6)

式中：An(ω)為第n次迭代后的傅里葉幅值譜。由于標準PSD和RRS之間的高相性，通常3～5次迭代就能夠使人工地震動的反應(yīng)譜收斂。

1.3 壓力容器算例

以壓力容器SSE的RRS作為算例，先計算出相應(yīng)的標準PSD，再擬合出匹配RRS的同時還包絡(luò)80%標準PSD的人工地震動。RRS的控制頻率區(qū)間從1 Hz到100 Hz，其中，1.0～5.0 Hz為反應(yīng)譜幅值的上升段；5.0～8.6 Hz為峰值平臺段，幅值為6.23g；8.6～42.5 Hz為下降段；42.5～100 Hz為下降段之后的平臺段，加速度幅值為0.92g左右。按照前述計算方法，生成了10條人工地震波，持時為40 s，綜合振動臺硬件的優(yōu)化需求以及人工波擬合過程中的效率，決定使用1 024 Hz作為采樣頻率。10條人工波反應(yīng)譜與RRS的匹配情況見圖1，規(guī)范要求0.3～24 Hz的頻域里人工地震動包絡(luò)80%標準PSD，而算例中RRS的控制頻率并不包含0.3～1 Hz的控制點，因此，僅考慮1～24 Hz頻域內(nèi)PSD的包絡(luò)情況，結(jié)果見圖2，其中1條人工地震動的加速度波形示意見圖3。

圖1 10條人工地震動反應(yīng)譜匹配Fig.1 Response spectra compatibility of 10 artificial ground motions

圖2 10條人工地震動標準PSD包絡(luò)Fig.2 Standard PSD envelopment of 10 artificial ground motions

圖3 人工地震動的加速度時程Fig.3 Acceleration time-history of artificial ground motion

從圖1～圖3中可知，人工地震動時程的峰值加速度(PGA)為14.2 m/s2，PSD完全包絡(luò)了80%的標準PSD，反應(yīng)譜的匹配在高頻部分無法實現(xiàn)的原因是設(shè)計譜中高頻部分的譜加速度與峰值平臺段對應(yīng)的加速度幅值相比有6～7倍的差距，高頻振子在人工波的作用下主要模態(tài)為剛體運動，因此，受人工地震波峰值加速度影響較大，在中間頻率平臺段擬合的條件下高頻部分難以完全收斂。

2 振動臺試驗

2.1 試驗設(shè)備及模型設(shè)計

試驗使用的雙點振動臺試驗機可實現(xiàn)3向6自由度的運動，單點最大試樣重量可達35 t，最大傾覆力矩為70 t·m，最大扭轉(zhuǎn)力矩為35 t·m，工作頻率范圍是0.1～100 Hz，最大工作位移、速度及加速度分別為±150 mm，±800 mm/s以及±2.5g。

壓力容器作為比較典型的核電廠樓面設(shè)備被選作試驗的目標系統(tǒng)。作為機理性試驗，選用模型與實際的反應(yīng)堆壓力容器有一定區(qū)別，試件以較為簡單的低壓壓力容器為原型進行以1∶1的比例制作，試樣主要部件包括桶體、蓋板、內(nèi)件、氣閥、底座等，材質(zhì)為Q235。試樣總高度為1 464 mm，桶內(nèi)徑為550 mm，壁厚3 mm，設(shè)計承受壓力0.8 MPa。底板與桶體采用厚20 mm的45°三角形肋板加固。壓力容器上有桶蓋，能夠拆卸便于安裝內(nèi)件。試樣與底座之間使用螺栓連接，底座與振動臺采用連接螺栓進行錨固。利用錘擊法得到壓力容器外件1階模態(tài)頻率為64.9 Hz，壓力容器內(nèi)件1階模態(tài)頻率為10.3 Hz。壓力容器模型裝配構(gòu)造見圖4，振動臺模型實際布置見圖5。

圖4 壓力容器模型構(gòu)造圖Fig.4 Detail structure drawing of pressure vessel

圖5 壓力容器內(nèi)件及外件模型Fig.5 Model of pressure vessel trim partand exterior part

2.2 測點布置

為了測量輸入地震載荷作用下試樣內(nèi)件與外件的動態(tài)響應(yīng)及應(yīng)力變化規(guī)律，模型試驗采用了加速度傳感器、位移傳感器及應(yīng)變傳感器。在模型上的測點布置見圖6。

圖6 模型測點布置圖Fig.6 Layout of measuring points on the model

由于內(nèi)件和外件相對獨立，因此，傳感器主要依據(jù)分別測量內(nèi)件及外件沿豎向分布的地震動響應(yīng)而布置。另外，通過經(jīng)驗判斷，外件最大變力的發(fā)生位置大概在桶底，因此，在相應(yīng)位置布置了應(yīng)變傳感器。如圖6所示，加速度傳感器沿試樣外側(cè)面從低至高布設(shè)，依次為從隔板第2、第4、第6層分別對應(yīng)桶壁外側(cè)的高度位置，桶體上的傳感器測點編號為A01～A03，并在外桶頂端布置加速度傳感器，測點編號A04，內(nèi)件的第2層隔板、第4層隔板以及最上層隔板上布設(shè)加速度傳感器，測點編號依次為A05～A07。位移傳感器在壓力容器模型外件的頂端布設(shè)，測試試樣頂端的位移，測點編號為D01。另外，地震臺臺面的加速度信號從振動臺獲取，測點編號為A00，臺面的位移從振動臺獲取，測點編號為D00。桶底距下法蘭20 mm處截面布置了2個應(yīng)變測點，測點編號為S01、S02，每個測點布設(shè)1個應(yīng)變花。

2.3 輸入地震動

以前述10條隨機人工地震波作為試驗輸入波的實驗組。對照組僅針對RRS匹配進行傳統(tǒng)的頻域法擬合人工地震波，即不使用標準PSD作為初始傅里葉幅值譜的依據(jù)，并且計算過程中沒有式5、式6所實現(xiàn)的包絡(luò)標準PSD的迭代修正條件。對照組總共也擬合了10條隨機人工地震動，比照實驗組的信號參數(shù)，持時取40 s，采樣頻率取1 024 Hz。對照組10條人工波反應(yīng)譜與RRS的匹配情況見圖7，PSD見圖8。

圖7 對照組人工地震動反應(yīng)譜匹配Fig.7 Response spectra compatibility of artificial ground motions in control group

圖8 對照組人工地震動PSDFig.8 Standard PSDs of artificial ground motions in control group

從圖8可見，80%的標準PSD處于10條對照組人工波PSD的中間位置。而對比圖7和圖1可以發(fā)現(xiàn)，由于對照組沒有包絡(luò)PSD的限制，對照組RRS的擬合精度略優(yōu)于實驗組。

2.4 模型動力響應(yīng)

著重關(guān)注桶頂最大響應(yīng)加速度a1(A04測點)、桶頂最大響應(yīng)位移d(D01測點)、內(nèi)件頂層隔板最大響應(yīng)加速度a2(A07測點)以及桶底最大應(yīng)變s(S01測點)4項動力響應(yīng)參數(shù)，這4項參數(shù)能夠相對全面地反映試件模型在地震過程中的響應(yīng)特點，數(shù)據(jù)見表1。

從表1可以看到，外件最大響應(yīng)加速度與輸入地震動的PGA接近，說明壓力容器外件在輸入地震動作用下做近似剛體運動，而內(nèi)件頂層的最大響應(yīng)加速度是外件頂部最大響應(yīng)加速度的4倍左右，這是因為其模態(tài)頻率比外件更低，所以內(nèi)件沿高度放大了人工波的PGA。實驗組與對照組內(nèi)外件最頂層最大加速度、外件底部應(yīng)變以及外件頂層最大位移4項數(shù)據(jù)平均值歸一化之后的對比見圖9，圖10為外件及內(nèi)件沿高度方向上各測點平均最大響應(yīng)加速度的包絡(luò)圖，其中，外件4個測點，內(nèi)件3個測點。

表1 模型地震動力響應(yīng)統(tǒng)計表Table 1 Statistical table of model’s seismic responses

圖9 實驗組與對照組地震響應(yīng)對比Fig.9 Comparison of seismic responses between experimental group and control group

圖10 實驗組與對照組響應(yīng)加速度對比Fig.10 Comparison of response accelerations between experimental group and control group

由圖9可知，對于外件頂部最大加速度、頂層最大位移、內(nèi)件頂層最大加速度以及桶底最大應(yīng)變4項地震動響應(yīng)參數(shù)，對照組的平均值均小于實驗組的平均值，分別小3.3%、2.8%、1.2%及3.4%；由圖10可以發(fā)現(xiàn)，在模型沿豎直高度方向上，內(nèi)件與外件的最大響應(yīng)加速度都是實驗組包絡(luò)了對照組，平均高2.0%。除了對于地震動響應(yīng)強度峰值的比較，圖11還對比了實驗組與對照組總共20個工況下，內(nèi)件頂層隔板響應(yīng)的加速度時程的0～50 Hz的功率譜密度曲線，對于他們在頻域中所包絡(luò)的平均能量(曲線所包絡(luò)的面積)，實驗組為1.601，對照組為1.551。可以看到，實驗組激勵的加速度響應(yīng)在0～50 Hz的能量高于對照組的3.3%。

以上結(jié)果均說明，實驗組的人工地震動激勵下模型的實際響應(yīng)高于對照組，這是因為，雖然實驗組和對照組的人工地震動都以RRS匹配為目標擬合，但實驗組的人工波同時包絡(luò)了PSD，因此，在PGA接近的情況下，實驗組在頻域中包含的能量實際大于對照組，這一點從圖9和圖2的PSD對比中就能夠看到，而振動臺試驗的結(jié)果也驗證了這個理論?？紤]到同時匹配RRS及包絡(luò)PSD的人工地震動在頻域中比僅匹配RRS的人工地震動在頻域具有更高的幅值，核設(shè)備在其激勵下有更大的響應(yīng)，因此，可以認為“包絡(luò)PSD”對于核設(shè)備動力分析輸入地震動的擬合是一項有意義的、更保守的設(shè)計條件。

圖11 實驗組及對照組地震動作用下模型內(nèi)件頂層響應(yīng)加速度PSD平均值對比圖Fig.11 Comparison of mean value of response accelerations on the top layer of inner part between experimental group and control group

3 結(jié)論

針對核設(shè)備的設(shè)計地震動需要同時匹配RRS及包絡(luò)標準PSD的問題，在2014版SRP建議方法的基礎(chǔ)上做了提升迭代效率的修改，生成了滿足要求的標準PSD；同時，在傳統(tǒng)的頻域法擬合人工波的算法中引入了新的迭代判定條件，使得最終擬合的人工地震波不僅匹配了RRS，還在規(guī)定的頻域范圍內(nèi)包絡(luò)標準PSD。

通過振動臺模型試驗，驗證了同時匹配RRS及包絡(luò)標準PSD的人工地震動引起的設(shè)備響應(yīng)高于僅匹配RRS的人工地震動這一結(jié)論，說明包絡(luò)標準PSD對于核設(shè)備的抗震設(shè)計輸入地震動的計算有一定的意義，是更保守的設(shè)計條件。但從本例看，PSD包絡(luò)對實際結(jié)構(gòu)地震的地震響應(yīng)影響有限。