程帥+李守義+司政+畢晨曦+楊陽(yáng)+陳演

摘要：水電站廠房的自振頻率關(guān)乎到其在動(dòng)力荷載下的安全，在設(shè)計(jì)階段加以預(yù)測(cè)和控制具有重要意義。運(yùn)用ANSYS有限元軟件對(duì)水電站廠房進(jìn)行模態(tài)分析，得到廠房結(jié)構(gòu)在一期、二期混凝土及地基的不同材料參數(shù)（彈模與容重）下的自振頻率，采用基于灰色關(guān)聯(lián)度的系統(tǒng)分析法，確定了水電站廠房自振頻率與各主要影響因子的關(guān)聯(lián)性及密切程度。結(jié)果表明：只考慮材料彈性模量時(shí)，地基彈模對(duì)廠房自振頻率的影響度最大，其次為廠房一期混凝土彈模和二期混凝土彈模；只考慮混凝土容重時(shí)，二期混凝土容重對(duì)廠房自振頻率的影響度略大于一期混凝土容重；同時(shí)考慮廠房混凝土彈模和容重時(shí)，結(jié)構(gòu)容重對(duì)其自振頻率的影響度大于彈模，二者與水電站廠房自振頻率的關(guān)聯(lián)度相差甚微。研究成果對(duì)控制類似結(jié)構(gòu)的自振頻率具有參考意義。

關(guān)鍵詞：水工結(jié)構(gòu)工程；水電站廠房；自振頻率；材料參數(shù)；灰色關(guān)聯(lián)度；量化影響因子

中圖分類號(hào)：TV314 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1672-1683（2017）03-0190-07

Abstract：Natural vibration frequency is related to the safety of hydropower stations under dynamic loads.It is significant to predict and control the frequency during the design phase.We performed modal analysis on the hydropower house using the ANSYS finite element software.Thus，we obtained the natural vibration frequency of the structure under the different material parameters （elasticity modulus and bulk density） of the first- and second-stage concrete and the foundation.We estimated the correlation between the structure's natural vibration frequency and the main material parameters using the systematical analysis method based on grey relational grade.The following conclusions were reached.When only the elasticity modulus was considered，the elasticity modulus of the foundation had the largest influence on the natural vibration frequency，followed by the first-stage concrete and the second-stage concrete of the hydropower house.When only the concrete bulk density was considered，the influence of the second-stage concrete bulk density was slightly larger than that of the first-stage concrete.When both elasticity modulus and concrete bulk density were considered，the concrete bulk density had a greater influence than elasticity modulus on the natural vibration frequency of the hydropower station，and their correlations with the natural vibration frequency showed little difference.The research results can provide reference for controlling the natural vibration frequency of similar structures.

Key words：hydraulic structure engineering；hydropower house；natural vibration frequency；material parameters；grey relational grade；to quantify impact factors

自振特性是水電站廠房本身固有且極為重要的力學(xué)性能，包括自振頻率、振型及阻尼比等，是衡量其剛度合理性的重要指標(biāo)[1]。作為自振特性的重要參數(shù)之一，自振頻率直接影響到廠房結(jié)構(gòu)的動(dòng)力響應(yīng)，是水電站廠房動(dòng)力分析的基礎(chǔ)。近些年，許多專家和學(xué)者針對(duì)水電站廠房自振頻率展開了大量的研究，并取得了豐富的成果。文獻(xiàn)[2]研究了水電站廠房上部不同網(wǎng)架支承方式及模擬方法對(duì)廠房整體自振頻率的影響，分析了機(jī)組振動(dòng)荷載作用下的響應(yīng)。其研究表明，網(wǎng)架支撐及模擬方式對(duì)廠房上部結(jié)構(gòu)的振動(dòng)頻率影響顯著，需按照實(shí)際尺寸建立包含網(wǎng)架結(jié)構(gòu)的計(jì)算模型。文獻(xiàn)[3]以勢(shì)流體單元模擬水電站廠房流道內(nèi)水體，以黏彈性邊界模擬地基效應(yīng)，對(duì)比分析了不考慮水體模型及考慮水體的流固耦合模型下廠房結(jié)構(gòu)自振頻率及地震響應(yīng)，認(rèn)為水體對(duì)其自振頻率及動(dòng)力響應(yīng)有較大的影響，僅模擬庫(kù)水和尾水的附加質(zhì)量會(huì)帶來(lái)較大的誤差。文獻(xiàn)[4]研究了不同混凝土強(qiáng)度和鋼管厚度下，水電站廠房排架結(jié)構(gòu)自振頻率的變化及對(duì)抗震性能的影響。其研究表明，在地震荷載作用下，改變混凝土強(qiáng)度相對(duì)于改變鋼管的壁厚對(duì)排架結(jié)構(gòu)的自振頻率影響更大。文獻(xiàn)[5]以某燈泡貫流式機(jī)組水電站廠房為對(duì)象，研究廠頂溢流新型布置型式下廠房的自振頻率，對(duì)部分振源與動(dòng)力響應(yīng)的相關(guān)性進(jìn)行了分析，進(jìn)行了共振校核并得出相應(yīng)振源的影響系數(shù)，為相應(yīng)結(jié)構(gòu)的設(shè)計(jì)提供理論依據(jù)。文獻(xiàn)[6]分析對(duì)比了不同彈性模量、阻尼比、泊松比等參數(shù)下水電站廠房自振頻率及動(dòng)力反應(yīng)，探究了基巖模擬范圍、基巖邊界條件、廠房?jī)?nèi)外水體、廠房樓板結(jié)構(gòu)形式、排架柱結(jié)構(gòu)形式等對(duì)水電站廠房自振頻率的影響，得出了一些定性的結(jié)論，為水電站廠房結(jié)構(gòu)的抗震優(yōu)化設(shè)計(jì)提供了可靠參考。另有大量國(guó)內(nèi)外科研人員對(duì)水電站廠房的自振頻率進(jìn)行了研究，并取得了相關(guān)成果[7-10]。

上述研究?jī)?nèi)容主要是對(duì)水電站廠房自振頻率的定性分析，宏觀上分析其各個(gè)因素的影響趨勢(shì)，但信息并不完全，不清楚各影響因素的主次及重要性程度，對(duì)各因素的影響權(quán)重定量描述較少。常用定量分析各因子之間關(guān)系的方法主要有主成分分析、回歸分析、主分量分析等數(shù)理統(tǒng)計(jì)法[10]，這些方法雖在諸多問(wèn)題中得到應(yīng)用，但往往需要大樣本，且要求有典型的概率分布，在部分問(wèn)題中難以實(shí)現(xiàn)?；疑P(guān)聯(lián)分析法不受此限制，該方法已經(jīng)成功應(yīng)用于多個(gè)領(lǐng)域[11-14]，但目前還沒有將灰色關(guān)聯(lián)分析引用到廠房自振頻率中，以分析其影響因素的重要度。

因此，本文運(yùn)用ANSYS有限元軟件對(duì)不同材料參數(shù)的水電站廠房進(jìn)行模態(tài)分析，以灰色關(guān)聯(lián)理論為基礎(chǔ)，通過(guò)建立水電站廠房自振頻率影響因素的灰色關(guān)聯(lián)度分析模型，選擇能夠量化的影響因素（彈模及容重）和自振頻率結(jié)果作為分析序列，通過(guò)分析各部分彈模、容重分別變化和同時(shí)變化下各因素與廠房自振頻率的關(guān)聯(lián)度，根據(jù)關(guān)聯(lián)度大小排序得出各因子對(duì)廠房自振頻率影響的重要度大小。

1 灰色關(guān)聯(lián)度原理和方法

1.1 灰色關(guān)聯(lián)分析基本原理

灰色系統(tǒng)理論是由鄧聚龍教授提出的一種研究少數(shù)據(jù)、貧信息、不確定性問(wèn)題的新方法[15-16]，其中灰色關(guān)聯(lián)度分析（Grey Relational Analysis）是灰色系統(tǒng)理論中一個(gè)重要分支，它將系統(tǒng)主要特性的發(fā)展趨勢(shì)以及主要影響因子進(jìn)行比較分析，由此來(lái)確定這些因素的關(guān)聯(lián)程度或者因子對(duì)效應(yīng)量的貢獻(xiàn)測(cè)度[17]。本質(zhì)上是根據(jù)各序列曲線微觀或宏觀上的幾何接近程度來(lái)判斷不同序列之間的關(guān)聯(lián)是否緊密。基本原理是通過(guò)線性插值法將系統(tǒng)因素的離散行為觀測(cè)值轉(zhuǎn)化為分段連續(xù)的折線，進(jìn)而根據(jù)折線的幾何特征（一階或二階斜率差）來(lái)量化分析因子之間或因子對(duì)主行為的影響程度。序列在整體幾何特征上越接近，發(fā)展態(tài)勢(shì)越相似，則關(guān)聯(lián)度越大[18]。

灰色關(guān)聯(lián)分析法對(duì)于初始樣本數(shù)據(jù)的多少及樣本是否具備典型分布規(guī)律無(wú)太高要求，且其計(jì)算過(guò)程簡(jiǎn)單，定量計(jì)算結(jié)果與基于理論的定性分析較為吻合，具有較強(qiáng)的適用性[17-18]。

1.2 原始數(shù)據(jù)處理

灰色關(guān)系分析中首先需要確定比較序列和參考序列組成原始數(shù)據(jù)矩陣，其中比較序列是各個(gè)影響因子對(duì)應(yīng)的子序列，參考序列是影響因子序列矩陣對(duì)應(yīng)條件下系統(tǒng)結(jié)果的母序列。根據(jù)本文所研究的問(wèn)題，分別選取水電站廠房一期、二期混凝土、基巖等不同部位的彈模及容重作為影響因子，即比較序列Xi（k）={xi（1），xi（2），…，xi（n）}，（其中k=1，2，…，n；i=1，2，…，m）；選取各情況下廠房的各階自振頻率作為系統(tǒng)結(jié)果，即參考序列X0（k）={x0（1），x0（2），…，x0（n）}。將參考序列（自振頻率）及其對(duì)應(yīng)的比較序列（不同部位材料參數(shù)）組合成含有m+1個(gè)序列的原始數(shù)據(jù)矩陣Y如式（1）所示

由于參考序列和比較序列中各元素的物理意義是不一樣的，數(shù)據(jù)的量綱也就不一樣。因此，為了保證各因素具有等效性和可比性，在進(jìn)行關(guān)聯(lián)度計(jì)算之前，需對(duì)原始數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理，常用的處理方法有初值化、均值化和區(qū)間相對(duì)值化等。本文采用初值化進(jìn)行無(wú)量綱處理。

將原始數(shù)據(jù)矩陣進(jìn)行無(wú)量綱化后可通過(guò)一系列的映射計(jì)算，得到比較序列與參考序列的關(guān)聯(lián)度。

1.3 灰色關(guān)聯(lián)度計(jì)算模型

由上述原始數(shù)據(jù)矩陣，根據(jù)灰色關(guān)聯(lián)度基本原理，選擇水電站廠房自振頻率作為參考序列，以對(duì)應(yīng)的各材料參數(shù)作為比較序列，建立水電站廠房自振頻率灰色關(guān)聯(lián)度計(jì)算模型見式（4），該模型為鄧聚龍教授所提出的鄧氏關(guān)聯(lián)度模型[16]，充分體現(xiàn)了灰色關(guān)聯(lián)四公理[19]中的約束條件，且考慮點(diǎn)與點(diǎn)之間的距離遠(yuǎn)近對(duì)關(guān)聯(lián)度的影響，是目前最常用的灰色關(guān)聯(lián)度計(jì)算模型，在諸多分析中得以應(yīng)用[10-11，18-20]。

式中：ri為各比較序列Xi（k）與參考序列X0（k）的關(guān)聯(lián)度；Di（k）=|Xi（k）-X0（k）|為比較序列和參考序列的絕對(duì)差值；Dmax=maxi（maxkDi（k））與Dmin=mini（minkDi（k））（i=1，2，…，m；K=1，2，…，n）分別為絕對(duì)差值中的兩級(jí)最大差和兩級(jí)最小差；ρ為分辨系數(shù)，用于提高灰色關(guān)聯(lián)系數(shù)之間的差異顯著性，ρ∈（0，1），根據(jù)相關(guān)文獻(xiàn)的取值原則和方法[16，20]，當(dāng)數(shù)據(jù)序列出現(xiàn)奇異值時(shí)，ρ取較小的值，可克服奇異值的支配作用，當(dāng)數(shù)據(jù)序列較為平穩(wěn)時(shí)，ρ應(yīng)取較大的值，以充分體現(xiàn)關(guān)聯(lián)度的整體性。為平衡這兩種情況帶來(lái)的影響并保持分辨系數(shù)的靜態(tài)性，鄧氏關(guān)聯(lián)度將分辨系數(shù)取為定值0.5，后來(lái)的灰色關(guān)聯(lián)度分析中一直沿用這個(gè)經(jīng)驗(yàn)值。分辨系數(shù)取值不同，關(guān)聯(lián)度的分布區(qū)間不同，但不會(huì)影響關(guān)聯(lián)度的排序。目前為止，雖已存在很多關(guān)于分辨系數(shù)的取值研究，但適用性不強(qiáng)，科學(xué)性有待進(jìn)一步討論。故本文取其常用值ρ=0.5。

計(jì)算出的灰色關(guān)聯(lián)度ri在[0，1]區(qū)間內(nèi)變化，其值越接近于1，說(shuō)明參考序列與比較序列的變化態(tài)勢(shì)越一致，該參考序列受比較序列的影響越大；相反，越接近于0則影響越小。關(guān)聯(lián)度大小的排序即為影響因子敏感性大小的排序。

2 有限元模型與假設(shè)

2.1 有限元模型

本文采用ANSYS有限元軟件建立某河床式水電站廠房中間機(jī)組壩段有限元模型，其中：壩高80 m，中間機(jī)組壩段寬26 m，順河床方向長(zhǎng)度78 m，采用軸流式機(jī)組，上游水頭65.0 m，下游尾水水頭22.5 m。在模型建立中詳細(xì)考慮了進(jìn)水口、流道、蝸殼、攔污柵庫(kù)、門庫(kù)與門槽、排架結(jié)構(gòu)和牛腿、屋頂、排沙孔、尾水管及尾水閘墩等細(xì)部結(jié)構(gòu)。

有限元計(jì)算模型包括壩體與基礎(chǔ)巖體，基巖的模擬范圍取矩形區(qū)域，基巖邊界值：深度方向由結(jié)構(gòu)最低點(diǎn)向下延伸100 m，左、右岸方向與壩段同寬，上下游方向由邊界處分別向上、下游延伸50 m。模型整體坐標(biāo)系的原點(diǎn)設(shè)在壩段右側(cè)底部與地基相交處，沿水流方向指向下游為X軸正方向；沿高度方向鉛直向上為Y軸正方向；垂直水流方向指向右岸為Z軸正方向。整體模型網(wǎng)格剖分基本采用8結(jié)點(diǎn)六面體實(shí)體單元SOLID45，廠房結(jié)構(gòu)網(wǎng)格尺寸不大于1.0 m，地基網(wǎng)格尺寸不大于10.0 m，地基與廠房相接部分通過(guò)四面體單元進(jìn)行過(guò)渡。整體有限元計(jì)算模型見圖1。

2.2 計(jì)算假定

根據(jù)水電站廠房結(jié)構(gòu)特點(diǎn)及其相關(guān)研究成果[21]，本文在三維有限元計(jì)算中，做如下假定：

①各壩段間設(shè)置有橫縫，因此各壩段獨(dú)立承受荷載，壩段間無(wú)相互作用，且由于橫縫作用，在計(jì)算模型中廠房?jī)蓚?cè)無(wú)約束，基巖上下游面及左右側(cè)面施加法向約束，底面施加三向約束。

②充分考慮結(jié)構(gòu)材料分區(qū)，計(jì)算中假定混凝土、基礎(chǔ)巖體為均質(zhì)、彈性、各向同性的連續(xù)體。

③壩體與地基連接可靠，不存在相互脫離的情形，滿足連續(xù)性條件，不考慮非線性接觸。

3 自振頻率與各影響因素的關(guān)聯(lián)度分析

3.1 各部分彈模對(duì)自振頻率的影響度

水電站廠房混凝土結(jié)構(gòu)總體上一般分兩期施工，一期混凝土主要包括流道底板、尾水管、下游胸墻等下游及下部基礎(chǔ)結(jié)構(gòu)，二期混凝土主要包括上游胸墻、蝸殼、發(fā)電機(jī)墩、座環(huán)基礎(chǔ)、流道及發(fā)電機(jī)層以上混凝土結(jié)構(gòu)等。根據(jù)各部分結(jié)構(gòu)特點(diǎn)與功能，一期和二期分別采用不同標(biāo)號(hào)的混凝土，其材料參數(shù)不同。

不同標(biāo)號(hào)混凝土的彈性模量不同，現(xiàn)設(shè)置以下三組方案，分析一期混凝土、二期混凝土及地基三者的彈性模量對(duì)水電站廠房自振頻率的影響，各方案取值見表1。其中A組方案中地基彈模保持不變，廠房一、二期混凝土等級(jí)標(biāo)號(hào)分別取為C10、C15、C20、C25、C30、C35及C40，彈模由17.5 GPa變化至32.5 GPa，且二者保持一致（根據(jù)最新水電工程水工建筑物抗震設(shè)計(jì)規(guī)范[22]，混凝土動(dòng)態(tài)彈性模量可較靜態(tài)提高50%，后文中混凝土及地基的彈模取值皆為靜態(tài)標(biāo)準(zhǔn)值E0，在模態(tài)計(jì)算中取為1.5E0）；B組方案中廠房一、二期混凝土標(biāo)號(hào)不變（C30），地基彈模根據(jù)一般工程地質(zhì)資料由7.4 GPa變化至11.0 GPa；為研究一期和二期混凝土對(duì)廠房自振頻率的影響程度大小，C組方案中地基彈模保持不變，廠房一期和二期混凝土分別取不同標(biāo)號(hào)。三組研究方案中一期、二期混凝土容重都取為25 kN/m3，泊松比取為0.167；地基的泊松比取為0.25，不考慮地基的容重。

基于上述有限元計(jì)算模型，由表1中各方案取值，運(yùn)用ANSYS有限元軟件對(duì)水電站廠房進(jìn)行模態(tài)分析，得到各方案下水電站廠房的自振頻率，為說(shuō)明廠房自振頻率隨結(jié)構(gòu)彈模變化的一般性規(guī)律，將A組中各方案下廠房的前20階自振頻率做以比較，見圖2。

由圖2可知，廠房的自振頻率隨著混凝土彈性模量的增大而增大，且這種增大效應(yīng)在高階頻率上更為明顯，在低頻增幅較??；其次，彈模的變化幅度與頻率變動(dòng)幅度基本一致，符合頻率增量與彈模增量呈線性關(guān)系的理論規(guī)律。B組六個(gè)方案中，由于地基的彈性模量變化范圍較小，其各階頻率偏差不大，與方案A5相似。C組六個(gè)方案中，廠房一期和二期混凝土標(biāo)號(hào)不同，但都各自由小變大，其各階頻率的規(guī)律性與A組七個(gè)方案相類似，在此不做贅述。

采用上述三組共17個(gè)方案的計(jì)算結(jié)果，以廠房的各階自振頻率為參考序列，以一期混凝土彈模、二期混凝土彈模及地基彈模三個(gè)影響因子作為比較序列，根據(jù)（1）～（4）式編寫Matlab程序計(jì)算出影響因子與參考序列的關(guān)聯(lián)度見表2。

由表2可知，總體上，水電站廠房的自振頻率對(duì)地基的彈性模量最為敏感，其次是一期混凝土，二期混凝土彈模對(duì)廠房自振頻率的影響最小。其中一、二期混凝土彈模與自振頻率的關(guān)聯(lián)系度較為接近，大部分在0.45至0.5之間，變化趨勢(shì)一致且波動(dòng)不大，地基彈模與各階自振頻率的關(guān)聯(lián)度浮動(dòng)較大，但都大于其它兩因素的關(guān)聯(lián)度，大部分在0.65以上，相關(guān)性較好，與水電站廠房自振頻率的關(guān)聯(lián)最為密切。前20階自振頻率與一期彈模、二期彈模及地基彈模的關(guān)聯(lián)度平均值分別為0.486 0、0.450 0、0.658 0。

3.2 各部分容重對(duì)自振頻率的影響度

為研究一期、二期混凝土容重對(duì)水電站廠房自振頻率的影響，根據(jù)實(shí)際工程中的取值范圍，將混凝土容重由23.0 kN/m3變化至26.0 kN/m3，設(shè)置以下兩組方案見表3。其中，D組中兩期混凝土容重變化一致；E組中為得出各部分混凝土容重對(duì)自振頻率的影響度，根據(jù)關(guān)聯(lián)度分析中原始數(shù)據(jù)的需要，兩期混凝土容重不同且有各自的變化規(guī)律。各方案中一期、二期混凝土彈性模量取為28 GPa，泊松比取為0.167；地基彈性模量取為8 GPa，泊松比取為0.25。

根據(jù)表3各方案，對(duì)水電站廠房進(jìn)行模態(tài)分析，得到水電站廠房的各階自振頻率。為分析廠房自振頻率隨混凝土容重變化的一般性規(guī)律性，將D組中各方案下水電站廠房的前20階自振頻率做以比較，見圖3。可知，與彈模對(duì)自振頻率的影響不同，水電站廠房的自振頻率隨著混凝土容重的增大而減小，但其減小程度與彈模的影響一致，在高階頻率上更為明顯，在低頻減幅較小。

采用D、E兩組方案的模態(tài)分析結(jié)果，以廠房各階自振頻率為參考序列，以一期、二期混凝土容重兩個(gè)影響因子作為比較序列，計(jì)算出各期混凝土容重與各階自振頻率的關(guān)聯(lián)度如表4所示。

由表4可知，一期混凝土容重與廠房各階自振頻率的性關(guān)度略小于二期混凝土，二者極為接近，差值小于0.038，前20階自振頻率與一期混凝土容重、二期混凝土容重關(guān)聯(lián)度平均值分別為0.509 5、0.544 7?？梢?，二期混凝土容重對(duì)廠房自振頻率影響略大于一期混凝土，其重要性程度二者相差不大。由于結(jié)構(gòu)混凝土容重的實(shí)際變化范圍相對(duì)較小，上述結(jié)論只針對(duì)僅考慮容重變化時(shí)的規(guī)律性，當(dāng)同時(shí)考慮容重和其他因素的變化時(shí)，容重相對(duì)于其他因素對(duì)水電站廠房自振頻率的影響度如何，在下文會(huì)做進(jìn)一步研究。

3.3 彈模和容重影響度綜合分析

上述以單因子變量法分別分析了水電站廠房各部分彈性模量和容重與其自振頻率的相關(guān)度。得到了只考慮彈模和只考慮容重時(shí)，各部分對(duì)應(yīng)參數(shù)與廠房自振頻率的普遍規(guī)律及其影響程度。然而，廠房自振頻率與這兩個(gè)因素都有著較為密切的聯(lián)系，實(shí)際工程中往往需要掌握最為敏感的因素，以便更加有效地對(duì)其自振頻率加以控制?，F(xiàn)同時(shí)考慮彈模和容重兩個(gè)因素，分析廠房的自振頻率與此二者的關(guān)聯(lián)性，以量化彈模和容重對(duì)廠房自振頻率的敏感度。

根據(jù)實(shí)際工程，在混凝土彈性模量和容重的合理取值范圍內(nèi)，設(shè)置12種不同取值組合見表5（為與前文區(qū)別，此12種方案設(shè)為F組）。各方案中一期、二期混凝土泊松比都取為0.167；地基的泊松比取為0.25，彈性模量取為8 000 MPa，不考慮地基的容重。

采用F組12種方案取值設(shè)置混凝土彈模和容重，對(duì)水電站廠房進(jìn)行模態(tài)分析，得到個(gè)方案下廠房的前20階自振頻率。以廠房的各階自振頻率為參考序列，分別以混凝土彈性模量和容重作為比較序列，計(jì)算出彈模和容重與水電站廠房自振頻率的關(guān)聯(lián)度見表6。

由表6可知，混凝土彈模和容重對(duì)水電站廠房各階自振頻率的影響效果相當(dāng)，大部分關(guān)聯(lián)度在0.6以上，相關(guān)關(guān)系較好，各階頻率與彈模和容重的關(guān)聯(lián)度大小無(wú)固定排序，對(duì)基頻的影響，容重因素大于彈模因素，其它大部分階次中，自振頻率與容重的關(guān)聯(lián)度亦大于與彈模的關(guān)聯(lián)度，前20階自振頻率與混凝土彈模和容重的關(guān)聯(lián)度平均值分別為0.603 3、0.625 1。可知，總體上，混凝土容重對(duì)水電站廠房自振頻率的影響大于其彈性模量的影響。

至此，以三種情況（只考慮彈模、只考慮容重、同時(shí)考慮彈模和容重）分析了部分材料參數(shù)對(duì)水電站廠房自振頻率的影響，每種情況都分別求出了各考慮因素與廠房自振頻率的關(guān)聯(lián)度。圖4為三種情況下對(duì)應(yīng)影響因素與水電站廠房前20階自振頻率關(guān)聯(lián)度的平均值對(duì)比。4 結(jié)論

本文用ANSYS有限元軟件對(duì)水電站廠房結(jié)構(gòu)進(jìn)行模態(tài)分析，以灰色關(guān)聯(lián)分析法為基礎(chǔ)編寫Matlab計(jì)算程序，對(duì)廠房自振頻率與其材料影響參數(shù)的相關(guān)度進(jìn)行了計(jì)算分析。主要結(jié)論如下：（1）水電站廠房自振頻率與其各部分彈性模量和容重密切相關(guān)。其中，自振頻率會(huì)隨著彈性模量的增大而增大，隨容重的增大而減小，且頻率變幅與彈?；蛉葜氐淖兎境示€性關(guān)系；（2）只考慮彈性模量時(shí)，地基彈模對(duì)廠房自振頻率的影響最大，關(guān)聯(lián)度達(dá)到0.65以上，明顯大于結(jié)構(gòu)混凝土與其自振頻率的關(guān)聯(lián)度，說(shuō)明在廠房選址、地基勘探試驗(yàn)及數(shù)值模擬中合理確定地基參數(shù)的重要性，當(dāng)水電站廠房結(jié)構(gòu)形式和材料分區(qū)與文中所研究的廠房相類似時(shí)，一期混凝土彈模與廠房自振頻率的關(guān)聯(lián)度略大于二期混凝土；（3）只考慮混凝土容重時(shí)，對(duì)于類似結(jié)構(gòu)的廠房，二期混凝土容重與其自振頻率的關(guān)聯(lián)度大于一期混凝土，二者相差不大，關(guān)聯(lián)度值都在0.5以上，在選擇各部分混凝土標(biāo)號(hào)及材料參數(shù)時(shí)，可適當(dāng)考慮此因素的影響；（4）同時(shí)考慮彈模和容重因素時(shí)，容重對(duì)水電站廠房自振頻率的影響程度略大于彈性模量，二者關(guān)聯(lián)度大小分別為0.6251、0.6033，可見在控制廠房結(jié)構(gòu)自振頻率過(guò)程中應(yīng)優(yōu)選影響權(quán)重高的容重作為主要控制因素。

文中利用灰色關(guān)聯(lián)分析法直觀量化了各材料參數(shù)對(duì)水電站廠房自振頻率的影響程度，但只是局限于部分可以量化的因素，對(duì)其他因素及其他自振特性參數(shù)未做討論。其次，廠房的自振頻率與其結(jié)構(gòu)形式有很大關(guān)系，文中所得到的部分結(jié)論僅適用于類似結(jié)構(gòu)的河床式廠房，但對(duì)于其他結(jié)構(gòu)形式的水電站廠房，亦可用本文的研究方法分析其影響因素的重要性程度，優(yōu)選控制參數(shù)。此外，現(xiàn)行灰色關(guān)聯(lián)度理論及算法，其模型本身具有一定缺陷，很難同時(shí)滿足保序性和規(guī)范性，其模型檢驗(yàn)準(zhǔn)則和具體量化標(biāo)準(zhǔn)還不夠成熟。這些問(wèn)題還有待于在今后的實(shí)踐中進(jìn)一步深化研究。

參考文獻(xiàn)（References）：

[1] 練繼建，王海軍，秦亮.水電站廠房結(jié)構(gòu)研究[M].北京：中國(guó)水利水電出版社，2007.（LIAN Ji-jian，WANG Hai-jun，QIN Liang.Study on the structure of hydropower station [M].Beijing：China Water & Power Press，2007.（in Chinese））

[2] 黑燦，伍鶴皋，傅丹，等.基于不同網(wǎng)架模擬方式的水電站廠房振動(dòng)特性研究[J].三峽大學(xué)學(xué)報(bào)（自然科學(xué)版），2016，38（1）：26-30.（HEI Can，WU He-gao，F(xiàn)U Dan，et al.Study of vibration characteristics of hydropower house based on different simulation methods of grid roof [J].Journal of China Three Gorges University（Natural Sciences），2016，38（1）：26-30.（in Chinese）） DOI：10.13393/j.cnki.issn.1672-948X.2016.01.006

[3] 孫偉，何蘊(yùn)龍，苗君，等.水體對(duì)河床式水電站廠房動(dòng)力特性和地震動(dòng)力響應(yīng)的影響分析[J].水力發(fā)電學(xué)報(bào)，2015，34（9）：119-127.（SUN Wei，HE Yun-long，Miao Jun，et al.Effects of water body on dynamic characteristics and seismic responses of run-of-river hydropower house [J].Journal of Hydroelectric Engineering，2015，34（9）：119-127.（in Chinese）） DOI：10.11660/slfdxb.20150916

[4] 陳娟，張燎軍.混凝土強(qiáng)度和鋼管厚度對(duì)鋼管混凝土排架結(jié)構(gòu)抗震性能的影響[J].水電能源科學(xué)，2011（11）：90-93.（CHEN Juan，ZHANG Liao-jun.Influence of concrete and steel pipe thickness on aseismic performance of concrete-filled steel tubular bent structure [J].Water Resources and Power，2011（11）：90-93.（in Chinese））

[5] 張龑，練繼建，劉昉，等.基于原型觀測(cè)的廠頂溢流式水電站廠房結(jié)構(gòu)振動(dòng)特性研究[J].天津大學(xué)學(xué)報(bào)：自然科學(xué)與工程技術(shù)版，2015（7）：584-590.（ZHANG Yan，LIAN Ji-jian，LIU Fang，et al.Vibration characteristics of powerhouse structure of roof overflow hydropower station based on prototype observation [J].Journal of Tianjin University （Science and Technology），2015（7）：584-590.（in Chinese））DOI：10.11784/tdxbz201312046

[6] 霍學(xué)平.大型水電站地面廠房動(dòng)力特性與抗振分析[D].大連：大連理工大學(xué)，2015.（HUO Xue-ping.Dynamic characteristics and vibration analysis for ground hydropower house [D].Dalian：Dalian University of Technology，2015.（in Chinese））

[7] Amorosi A，Boldini D，di Lernia A.Seismic ground response at Lotung：Hysteretic elasto-plastic-based 3D analyses [J].Soil Dynamics and Earthquake Engineering，2016，85：44-61.DOI：10.1016/j.soildyn.2016.03.001

[8] Wei S，Zhang L.Vibration analysis of hydropower house based on fluid-structure coupling numerical method [J].Water Science and Engineering，2010，3（1）：75-84.DOI：10.3882/j.issn.1674-2370.2010.01.008

[9] 陳玲玲，戴湘和，彭定.三峽電站廠房鋼襯蝸殼不同埋設(shè)方式抗震分析[J].長(zhǎng)江科學(xué)院院報(bào)，2012，29（2）：51-54.（CHEN Ling-ling，DAI Xiang-he，PENG Ding.Anti-seismic study on different embedment options for steel-lined spiral case in the Three Gorges Hydropower Plant [J].Journal of Yangtze River Scientific Research Institute，2012，29（2）：51-54.（in Chinese））

[10] 尹吉娜，楊杰，任杰，等.基于灰色關(guān)聯(lián)分析的邊坡穩(wěn)定影響因素敏感性評(píng)估[J].合肥工業(yè)大學(xué)學(xué)報(bào)：自然科學(xué)版，2015（1）：75-78.（YIN Ji-na，YANG Jie，REN Jie，et al.Study on impact factors sensitivity of slope stability based on grey correlation analysis [J].Journal of Hefei University of Technology：Natural Science Edition，2015（1）：75-78.（in Chinese）） DOI：10.3969/j.issn.1003-5060.2015.01.016

[11] 姬宏，田盛.基于灰色關(guān)聯(lián)度的干旱區(qū)社會(huì)經(jīng)濟(jì)子系統(tǒng)需水量協(xié)調(diào)研究[J].南水北調(diào)與水利科技，2009，7（3）：42-44.（JI Hong，TIAN Sheng.Coordination study on water demand of socioeconomic subsystem in arid area based on grey relational grade [J].South-to-North Water Transfers and Water Science & Technology，2009，7（3）：42-44.（in Chinese）） DOI：10.3969/j.issn.1672-1683.2009.03.012

[12] Kuram E，Ozcelik B.Multi-objective optimization using Taguchi based grey relational analysis for micro-milling of Al 7075 material with ball nose end mill [J].Measurement，2013，46（6）：1849-1864.DOI：10.1016/j.measurement.2013.02.002

[13] Thanigaivelan R，Arunachalam R.Optimization of process parameters on machining rate and overcut in electrochemical micromachining using grey relational analysis [J].Journal of Scientific and Industrial Research，2013，72：36-42.

[14] 曹洪洋，王禹，滿兵.基于改進(jìn)灰色關(guān)聯(lián)分析的泥石流危險(xiǎn)性評(píng)價(jià)[J].南水北調(diào)與水利科技，2015，13（1）：91-94.（CAO Hong-yang，WANG Yu，MAN Bing.Risk evaluation of potential debris flow based on the improved grey correlation method [J].South-to-North Water Transfers and Water Science & Technology，2015，13（1）：91-94.（in Chinese）） DOI：10.13476/j.cnki.nsbdpk.2015.01.021

[15] Deng J L.Control problems of grey systems [J].Systems & Control Letters，1982，5（5）：288-294.

[16] 鄧聚龍.灰理論基礎(chǔ)[M].武漢：華中科技出版社，2002.（DENG Ju-long.The basis of grey theory [M].Wuhan：Huazhong University of Science & Technology Press，2002.（in Chinese））

[17] 劉思峰，蔡華，楊英杰，等.灰色關(guān)聯(lián)分析模型研究進(jìn)展[J].系統(tǒng)工程理論與實(shí)踐，2013，33（8）：2041-2046.（LIU Si-feng，CAI Hua，YANG Ying-jie，et al.Advance in Grey incidence analysis modeling [J].Systems Engineering Theory & Practice，2013，33（8）：2041-2046.（in Chinese））

[18] 張莎.灰色關(guān)聯(lián)分析新算法研究及其意義[D].長(zhǎng)春：東北師范大學(xué)，2012.（ZHANG Sha.The study on the new algorithms of grey relational analysis and its significances [D].Changchun：Northeast Normal University，2012.（in Chinese））

[19] 劉思峰，謝乃明.灰色系統(tǒng)理論及其應(yīng)用[M].北京：科學(xué)出版社，2008.（LIU Si-feng，XIE Nai-ming.Grey system theory and application [M].Beijing：Science Press，2008.（in Chinese））

[20] 東亞斌，段志善.灰色關(guān)聯(lián)度分辨系數(shù)的一種新的確定方法[J].西安建筑科技大學(xué)學(xué)報(bào)（自然科學(xué)版），2008，40（4）：589-592.（DONG Ya-bin，DUAN Zhi-shan.A new determination method for identification coefficient of grey relational grade [J].Journal of Xi′an University of Architecture & Technology （Natural Science Edition），2008，40（4）：589-592.（in Chinese））

[21] 李守義，肖靜.河床式水電站廠房橫縫止水布置形式研究[J].水力發(fā)電學(xué)報(bào)，2014，33（5）：165-168.（LI Shou-yi，XIAO Jing.Study on waterstop arrangements of transverse joints of hydropower house in river channel [J].Journal of Hydroelectric Engineering，2014，33（5）：165-168.（in Chinese））

[22] NB 35047-2015.水電工程水工建筑物抗震設(shè)計(jì)規(guī)范[S].（NB 35047-2015.Code for seismic design of hydraulic structures of hydropower projects [S].（in Chinese））