葛新峰， 張 敬， 祝雙桔， 錢巨林， 趙 雷， 鄭圣義， 鄭 源， 徐 旭， BINAMA Maxime

(1.河海大學(xué) 能源與電氣學(xué)院，南京 210098；2.重慶航運(yùn)建設(shè)發(fā)展有限公司，重慶 401121；3．國網(wǎng)新源水電有限公司富春江水力發(fā)電廠，浙江 桐廬 311504)

在國外對一般連接螺栓的強(qiáng)度問題研究較早，王玨等[1-2]對承受軸向載荷情況下的螺紋連接進(jìn)行了應(yīng)力分析；Chen[3]分析了螺栓的應(yīng)力松弛特性，并分析了允許應(yīng)力的計算方法。Waszczak等[4]對比了扭曲能準(zhǔn)則、最大應(yīng)力準(zhǔn)則、最大應(yīng)變準(zhǔn)則這三個失效準(zhǔn)則預(yù)測下的受壓螺栓試樣的失效模式和強(qiáng)度；Dimitri等[5]以連接螺栓的尺寸及載荷類型(靜載荷及動載荷)作為隨機(jī)變量研究了螺栓的拉伸及剪切強(qiáng)度并進(jìn)行了優(yōu)化分析；Naruse等[6-7]通過VDI2230及有限元手段對風(fēng)力發(fā)電機(jī)組連接螺栓進(jìn)行了理論及仿真分析，對比計算結(jié)果較為貼近；Onur[8]通過hypermesh-nastran對每個組件進(jìn)行了有限元分析(FEA);Yamanaka等[9]為了深入研究螺栓應(yīng)力，變形和相互作用機(jī)理,基于有限的拉格朗日定律,通過有限元,對螺栓進(jìn)行大變形計算和分析。Kato等[10]基于彈塑性理論進(jìn)一步給出了螺栓強(qiáng)度達(dá)到屈服極限后的模型，是對連接螺栓超出彈性范圍后的補(bǔ)充。Chen等[11]發(fā)現(xiàn)殘余應(yīng)力在對于螺栓強(qiáng)度和可靠性方面起著重要作用。

國內(nèi)學(xué)者針對螺栓靜強(qiáng)度相關(guān)也進(jìn)行了大量的研究，工程上對于螺栓強(qiáng)度的設(shè)計及復(fù)核常采用安全系數(shù)法。陳真等[12-16]通過VDI2230及有限元仿真對風(fēng)力發(fā)電機(jī)塔筒法蘭連接螺栓進(jìn)行分析，計算了實(shí)際工況中螺栓的等效應(yīng)力及連接強(qiáng)度；熊欣等[17]針對水泵水輪機(jī)頂蓋座環(huán)連接螺栓，對比了《機(jī)械設(shè)計手冊》、VDI2230準(zhǔn)則及有限元三種方法進(jìn)行螺栓強(qiáng)度分析的差異，指出機(jī)械設(shè)計手冊中僅考慮同心夾緊、同心加載，VDI2230準(zhǔn)則考慮了偏心夾緊及加載，而有限元基本可以真實(shí)反映螺栓的應(yīng)力集中；何少潤等[18]針對大型抽水蓄能機(jī)組頂蓋螺栓探究了螺栓預(yù)緊力、剩余預(yù)緊力和工作載荷間的關(guān)系，初文華等[19]基于SPH-FEM耦合算法，計算爆炸螺栓與其連接結(jié)構(gòu)間的相互作用,王軍評等[20]研究了預(yù)緊力對爆炸螺栓結(jié)構(gòu)響應(yīng)的影響。楊風(fēng)利等[21]研究了溫度對螺栓應(yīng)力、螺栓預(yù)緊力以及螺栓松動性的影響。葛新峰等[22-23]對不同工況下的水輪機(jī)頂蓋螺栓剛強(qiáng)度進(jìn)行計算與分析,并對某混流式水輪機(jī)部分螺栓在發(fā)生斷裂后剩余螺栓的強(qiáng)度進(jìn)行分析。

水輪機(jī)頂蓋螺栓的結(jié)構(gòu)與風(fēng)機(jī)塔筒螺栓及壓力容器法蘭螺栓等結(jié)構(gòu)極為相似，均采用法蘭螺栓連接形式，因此將上述文獻(xiàn)中對螺栓靜強(qiáng)度的研究方法應(yīng)用于水輪機(jī)頂蓋螺栓依然可行。前人已對螺栓靜強(qiáng)度理論進(jìn)行了的研究，而針對軸流式水輪機(jī)頂蓋螺栓，建立了水輪機(jī)全三維過流部件模型流體模型和水輪機(jī)頂蓋和螺栓結(jié)構(gòu)模型，開展流固體耦合計算還沒有類似研究，本文的研究可以為水電站的安全穩(wěn)定運(yùn)行及檢修提供理論依據(jù)。

1 流固耦合基本理論

1.1 流固耦合基本方程

流固耦合基本控制方程主要包括三個方程：流體控制方程、固體控制方程、耦合面控制方程。

1.1.1 流體控制方程

對于一般的可壓縮牛頓流來說，流體流動需要遵循基本的質(zhì)量守恒定律及動量守恒定律。

質(zhì)量守恒方程：

(1)

動量守恒方程：

(2)

式中:t表示時間;ρf表示流體密度;v表示流體速度矢量;ff是體積力矢量;τf表示剪切力張量。τf可表示為

τf=(-p+μ?·v)I+2μe

(3)

1.1.2 固體控制方程

固體部分的守恒方程可由牛頓第二定律推導(dǎo)而得

(4)

1.1.3 耦合面控制方程

同樣，流固耦合遵循最基本的守恒原則，因此在流固耦合的耦合面上，應(yīng)滿足流體與固體應(yīng)力、位移等變量的相等或守恒，即：

(5)

式中:下標(biāo)f表示流體;下標(biāo)s表示固體。

1.2 耦合方程求解算法

流固耦合方程的求解算法主要有兩種：直接耦合式解法(directly coupled solution)和分離解法(partitioned solution)。

分離解法將流體及固體控制方程通過不同求解器分別求解，是目前大多數(shù)商業(yè)軟件所使用的方法。在流體求解器中求解壓力，在固體求解器中求解位移，其關(guān)鍵在于進(jìn)行流固耦合面的數(shù)據(jù)交換。而實(shí)際水體側(cè)網(wǎng)格與固體側(cè)網(wǎng)格并非完全匹配，此時就需要在數(shù)據(jù)傳遞前對數(shù)據(jù)進(jìn)行相應(yīng)的插值運(yùn)算。根據(jù)數(shù)據(jù)傳遞方向的不同，又可分為單向流固耦合及雙向流固耦合。直接耦合式解法聯(lián)合求解流固控制方程，不存在時間滯后問題，但卻存在同步求解的收斂問題及耗時問題；而分離解法幾乎可以獨(dú)立利用流體及固體求解器求解，能夠大大降低對內(nèi)存的需求，但存在時間滯后問題及耦合面能量不完全守恒問題，本文考慮到僅模擬水體分析結(jié)果傳遞給頂蓋結(jié)構(gòu)，采用了分離解法的單向流固耦合。

2 模型建立及網(wǎng)格劃分

2.1 流體及頂蓋模型建立

采用UG對軸流式水輪機(jī)流體計算域及頂蓋螺栓結(jié)構(gòu)計算域進(jìn)行三維建模，如圖1所示。流體部分主要包括蝸殼、固定導(dǎo)葉及活動導(dǎo)葉區(qū)、轉(zhuǎn)輪區(qū)以及尾水管這4個部分；結(jié)構(gòu)部分包括外頂蓋法蘭、內(nèi)頂蓋及64根螺栓。

圖1 軸流式水輪機(jī)全流道及頂蓋三維模型Fig.1 Three-dimensional model of the full flow channel and roof of the axial flow turbine

2.2 流體網(wǎng)格劃分及無關(guān)性驗(yàn)證

采用非結(jié)構(gòu)四面體網(wǎng)格對軸流式水輪機(jī)全流道進(jìn)行網(wǎng)格劃分，結(jié)果如下圖所示，圖2(a)～2(d)分別為蝸殼、導(dǎo)葉、轉(zhuǎn)輪及尾水管的網(wǎng)格圖，對頂蓋過流面區(qū)域及葉輪、導(dǎo)葉部分進(jìn)行了局部加密。

選取四種方案進(jìn)行流體網(wǎng)格無關(guān)性驗(yàn)證，網(wǎng)格數(shù)量為90萬～520萬，驗(yàn)證工況為額定水頭額定開度工況，計算結(jié)果如表1所示。隨著網(wǎng)格數(shù)量的增加，軸流式水輪機(jī)的效率不斷增大，頂蓋軸向水壓力幾乎沒有變化，當(dāng)網(wǎng)格數(shù)量增加到450萬時，效率計算值變化趨勢便不再明顯。考慮到網(wǎng)格數(shù)量的增加會導(dǎo)致對計算機(jī)性能的要求越來越高，最終選定的網(wǎng)格數(shù)量為450萬。

表1 流體網(wǎng)格無關(guān)性驗(yàn)證Tab.1 Fluid grid independence verification

2.3 結(jié)構(gòu)網(wǎng)格劃分及無關(guān)性驗(yàn)證

選擇SOLID186六面體20節(jié)點(diǎn)網(wǎng)格單元對頂蓋進(jìn)行網(wǎng)格劃分，網(wǎng)格劃分如圖3所示，為頂蓋及螺栓網(wǎng)格，聯(lián)接螺栓為M64×4 8.8級，屈服極限為640 MPa(1 825 kN)，有效截面積2 850.8 mm2。選取額定水頭額定開度工況為網(wǎng)格無關(guān)性分析，對螺栓頭倒角處加密網(wǎng)格，以倒角處螺栓最大應(yīng)力作為網(wǎng)格尺寸評判的標(biāo)準(zhǔn)，無關(guān)系驗(yàn)證計算結(jié)果如圖4所示。

圖3 頂蓋及螺栓網(wǎng)格劃分Fig.3 Top cover and bolt meshing

圖4 結(jié)構(gòu)網(wǎng)格無關(guān)性驗(yàn)證Fig.4 Structural grid independence verification

從圖4中可以看出，隨著網(wǎng)格尺寸的加密，網(wǎng)格數(shù)量增加，螺栓應(yīng)力值也在不斷增大；當(dāng)網(wǎng)格尺寸為0.4 mm時，螺栓的應(yīng)力值變化趨于平緩。因此，選取網(wǎng)格尺寸為0.4 mm的網(wǎng)格劃分方案，此時網(wǎng)格總數(shù)為18萬左右，節(jié)點(diǎn)總數(shù)為59萬左右，頂蓋整體及螺栓最終的網(wǎng)格劃分結(jié)果如圖3所示。

2.4 流體計算邊界條件

流體仿真中針對湍流的模擬采用RNGk-ε模型，該模型能夠較好的反應(yīng)出帶有彎曲壁面及強(qiáng)旋流的流動。流體計算區(qū)域的邊界條件主要包括：蝸殼進(jìn)口：采用壓力進(jìn)口邊界條件，給定軸流式水輪機(jī)的總工作水頭；壁面：采用無滑移壁面邊界條件；尾水管出口：采用質(zhì)量流量邊界條件；轉(zhuǎn)輪區(qū)：由于轉(zhuǎn)輪是旋轉(zhuǎn)部件，該區(qū)域采用旋轉(zhuǎn)參考坐標(biāo)系方法，給定轉(zhuǎn)輪的轉(zhuǎn)速值，并在動靜交界面處使用Frozen Rotor(凍結(jié)轉(zhuǎn)子)。

2.5 結(jié)構(gòu)計算荷載加載及邊界條件

由于頂蓋螺栓模型是復(fù)雜的裝配體模型，因此各部件需要設(shè)置接觸約束，本文模型主要存在3對接觸：內(nèi)頂蓋法蘭和外頂蓋法蘭，螺栓頭與內(nèi)頂蓋法蘭間的水平接觸、螺栓底部螺紋與內(nèi)頂蓋螺紋孔間的接觸，前兩種接觸由于裝配壓緊是存在摩擦力的，所以設(shè)置為摩擦接觸，摩擦因數(shù)取0.2，后一對接觸設(shè)為綁定接觸，以此來模擬螺紋間的擰緊配合。圖5所示為有限元模型的邊界條件并施加工況載荷情況，其中，A為當(dāng)?shù)刂亓铀俣?；B為座環(huán)底部固定約束；C～J為螺栓預(yù)緊力，其中，螺栓荷載的施加分兩步加載，第一步施加預(yù)緊載荷400 kN，第二步對螺栓預(yù)緊力進(jìn)行鎖緊，并加載頂蓋及水導(dǎo)軸承座水壓，即外部流體計算的水壓力載荷通過流固耦合面?zhèn)鬟f到頂蓋和水導(dǎo)軸承座上，如圖6所示，水壓力包含部分頂蓋接觸面水壓力P1，水導(dǎo)軸承座接觸面水壓力P2。

圖5 邊界條件與載荷施加Fig.5 Boundary conditions and load application

圖6 耦合面壓力載荷的傳遞Fig.6 Coupling surface pressure load transfer

3 數(shù)值模擬結(jié)果分析

3.1 頂蓋過流面水壓力分布

對額定水頭Hr=14.5 m及最大水頭Hmax=18.5 m不同導(dǎo)葉開度下頂蓋過流面的水壓力分布進(jìn)行分析，如圖7為兩個水頭下不同導(dǎo)葉開度的水輪機(jī)頂蓋及水導(dǎo)軸承座過流面的整體壓力分布圖，從 圖7(a)將過流面的整體壓力分布圖分為三個部分：活動導(dǎo)葉處水壓力P0，內(nèi)頂蓋處水壓力P1，水導(dǎo)軸承座水壓力P2。

從圖7中可以看出，活動導(dǎo)葉分布圓周圍壓力P0沿圓周方向呈現(xiàn)很明顯的不均勻性，活動導(dǎo)葉進(jìn)口前的壓力大于出口處的壓力，根據(jù)伯努利方程，導(dǎo)葉出口處的水流流速便大于導(dǎo)葉進(jìn)口處，這也說明了流體經(jīng)過導(dǎo)葉，經(jīng)導(dǎo)葉引流產(chǎn)生了明顯的加速，同時也增加了水流進(jìn)入轉(zhuǎn)輪之前的環(huán)量；活動導(dǎo)葉之后內(nèi)頂蓋整體壓力P1分布較為均勻，呈軸對稱分布；隨著半徑的減小，水輪機(jī)頂蓋及水導(dǎo)軸承座整體過流面的壓力P2逐漸減小。根據(jù)不同導(dǎo)葉開度下的壓力云圖可以發(fā)現(xiàn)，隨著導(dǎo)葉開度的增大，頂蓋整體的高壓區(qū)沿著徑向向內(nèi)擴(kuò)展，這將導(dǎo)致大開度下頂蓋過流面的壓力增大，從而便增大了傳遞到頂蓋螺栓上的工作載荷。

從圖7中可見內(nèi)頂蓋耦合面壓力是基本是軸對稱分布的，如圖8(a)以A特征線為例，研究不同導(dǎo)葉開度下壓力分布的不同規(guī)律。圖8(b)和圖8(c)分為額定水頭及最大水頭不同導(dǎo)葉開度頂蓋耦合面的壓力隨半徑變化的曲線，從圖中可知，工作水頭的增大，耦合面壓力也隨之增大，如Hr工況最大壓力為135 kPa，Hmax工況最大壓力為173 kPa；隨著導(dǎo)葉開度的增大，內(nèi)頂蓋耦合面上的總壓力增加，如Hr工況下40%開度平均壓力為100 kPa；100%開度下平均壓力為130 kPa；隨著導(dǎo)葉開度的增大，耦合面上的壓力降減小，100%開度工況耦合面上的壓力值幾乎相同，額定水頭和最大水頭的壓力降都僅有8 kPa左右。

(a) 內(nèi)頂蓋耦合面壓力特征線

(b) 額定水頭Hr

(c) 最大水頭Hmax圖8 內(nèi)頂蓋耦合面壓力隨半徑變化曲線Fig.8 Change curve of inner top cover coupling surface pressure with radius

3.2 頂蓋螺栓整體剛強(qiáng)度分析

圖9 頂蓋及螺栓整體變形云圖Fig.9 Overall deformation cloud map of the top cover and bolts

圖10 頂蓋及螺栓整體應(yīng)力云圖Fig.10 Overall stress cloud diagram of top cover and bolts

3.3 螺栓截面應(yīng)力分析

為了研究螺栓不同位置截面等效應(yīng)力的危險值，以最危險工況，即最大水頭大開度工況為例，選取了螺栓的三個截面，如圖11(a)所示，M1代表螺栓中間截面；W1代表靠近螺栓頭處的危險截面；W2代表靠近螺栓光桿與螺紋交接處的危險截面，從圖11(b)～11(d)中可以看出，螺栓三個截面的等效應(yīng)力分布都不均勻，應(yīng)力最大值點(diǎn)與應(yīng)力最小值點(diǎn)幾乎沿圓截面中性軸對稱分布，因此螺栓在承受軸向應(yīng)力的同時還存在著偏心夾緊及偏心載荷所引起的彎曲應(yīng)力，螺栓三個截面處的等效應(yīng)力最大值并不相同，W1和W2截面的最大等效應(yīng)力(182.62 MPa、209.87 MPa)都大于M1截面的最大等效應(yīng)力值169.34 MPa，這也表明螺栓頭部及螺栓光桿與螺紋交接處相比螺栓光桿區(qū)域更易先發(fā)生破壞，以螺栓M1截面所有計算節(jié)點(diǎn)的平均值作為螺栓整體的平均應(yīng)力，即有限元計算的螺栓工作應(yīng)力，各個工況螺栓工作應(yīng)力有限元統(tǒng)計如表2中“中截面等效應(yīng)力平均值”所示，從表中可以看出，有限元及理論計算值誤差很小并且各個計算值隨工況載荷變化的趨勢也一樣。

表2 頂蓋螺栓工作應(yīng)力有限元計算結(jié)果與現(xiàn)場試驗(yàn)對比Tab.2 Comparison of finite element calculation results of roof bolt working stress with field test

M1截面遠(yuǎn)離上下兩個接觸面，能夠忽略圣維南邊界效應(yīng)的影響，可以將其截面的平均應(yīng)力作為螺栓整體的平均應(yīng)力,而W1和W2位置是工程實(shí)踐中螺栓最常發(fā)生破壞的兩個位置。

4 現(xiàn)場試驗(yàn)及數(shù)模驗(yàn)證

中空空心圓柱筒型受力傳感器(以下簡稱：中空測力傳感器)由內(nèi)置彈性元件、電阻應(yīng)變片、溫度補(bǔ)償電路和轉(zhuǎn)換電路構(gòu)成，其中彈性元件對于測量來說具有極其重要的作用，當(dāng)被測軸向力作用于筒型彈性敏感元件上時，筒型彈性元件會發(fā)生形變，此時內(nèi)置在傳感器內(nèi)部的彈性元件上的應(yīng)變片就會感知該力作用的形變，并且引起應(yīng)變片電阻阻值的變化，通過電路轉(zhuǎn)換，輸出一個與被測信號成線性關(guān)系的電壓信號。彈性元件的受力原理框圖如圖12(a)，中空空心圓柱筒型受力傳感器外形如圖12(b)、(c)，受力傳感器材料為不銹鋼1Cr18Ni9Ti，精度1%FS，量程2 000 kN，可滿足現(xiàn)場強(qiáng)度及精度的要求的需要，傳感器實(shí)際測得的螺栓總載荷為不考慮應(yīng)力集中效應(yīng)的平均應(yīng)力，如圖12(d)所示為單個中空測力傳感器安裝于頂蓋螺栓和內(nèi)頂蓋之間，如圖12(e)為現(xiàn)場均勻安裝了6個傳感器?，F(xiàn)場試驗(yàn)工況為高水頭下的機(jī)組從停機(jī)狀態(tài)到開機(jī)，然后改變開度，最后停機(jī)的整個機(jī)組啟停過程，活動導(dǎo)葉開度每次調(diào)整10%，由于在30%開度以下機(jī)組不能穩(wěn)定運(yùn)行，數(shù)據(jù)采集時間較短，如圖13為Hmax水頭下6個頂蓋螺栓受力曲線圖，每個測點(diǎn)是30秒時間段內(nèi)的平均值，試驗(yàn)時間約為40分鐘，由圖13可知：

(a) 彈性元件受力原理圖

(b) 中空測力傳感器-俯視

(c) 中空測力傳感器-測試

(d) 單個傳感器現(xiàn)場安裝圖

(e) 6個中空壓力傳感器安裝位置圖圖12 中空測力傳感器基本原理及現(xiàn)場布置Fig.12 Basic principle and site layout of hollow force sensor

(1) 可以通過停機(jī)的螺栓受力得到，6個螺栓測點(diǎn)p1～p6的預(yù)緊力不同，分別為352、420、372、405、406、410 kN，因此6個螺栓的受力曲線并未重合，但整體變化趨勢一致；

(2) 螺栓在開機(jī)和停機(jī)兩個階段螺栓受力有明顯區(qū)別，開機(jī)正常運(yùn)行時螺栓受力比停機(jī)時增加了20～30 kN；

(3) 在較低負(fù)荷0～30%開度時段為振動區(qū)，頂蓋螺栓受力最大較正常運(yùn)行高5 kN；

(4) 在30%～100%導(dǎo)葉開度不同負(fù)荷時螺栓受力保持穩(wěn)定。

為驗(yàn)證有限元計算結(jié)果，選取安裝預(yù)緊力最為接近400 kN的螺栓應(yīng)力曲線，p4，p5兩個測點(diǎn)，螺栓截面面積為2 850 mm2，則僅有預(yù)緊力時螺栓中截面的等效應(yīng)力為400 kN/2 850 mm2=142.5 MPa，將圖13中p4，p5兩個測點(diǎn)預(yù)緊力轉(zhuǎn)換為預(yù)緊應(yīng)力如圖14，同時將表2中Hmax工況下的停機(jī)、40%、70%、100%時機(jī)組的有限元計算螺栓中截面的平均值畫在圖14中，從圖14和表2中可知在停機(jī)及三個穩(wěn)定運(yùn)行工況下數(shù)模計算的數(shù)值和現(xiàn)場試驗(yàn)比較吻合，誤差不大于5%，驗(yàn)證了有限元計算的準(zhǔn)確性，說明中空測力傳感器能夠比較準(zhǔn)確的反應(yīng)螺栓的軸向受力。

圖13 現(xiàn)場螺栓受力測試結(jié)果(Hmax)Fig.13 On-site bolt force test results(Hmax)

圖14 中空壓力傳感器及現(xiàn)場布置情況Fig.14 Hollow pressure sensor and site layout

5 結(jié)論

本文基于流固耦合的方法對軸流式水輪機(jī)全流道及頂蓋螺栓進(jìn)行了計算分析，并開展現(xiàn)場試驗(yàn)研究，研究對比了不同水頭及導(dǎo)葉開度工況下頂蓋過流面水壓力分布及螺栓受力特性，主要結(jié)論如下：

(1) 頂蓋耦合面的水壓力為軸對稱分布，大開度工況下內(nèi)頂蓋耦合面的壓力降較小，幾乎是均勻分布的；頂蓋水壓力隨著水頭和導(dǎo)葉開度的增大而增大，最大水頭大開度工況頂蓋總軸向壓力值達(dá)到最大。

(2) 通過流固耦合數(shù)值計算可以得到頂蓋及螺栓的變形及應(yīng)力，驗(yàn)證是否超標(biāo)，螺栓不同位置截面應(yīng)力分布不均，螺栓同時承受著軸向應(yīng)力與彎曲應(yīng)力，螺栓頭部、螺栓與螺紋交接處相比螺栓比螺栓中部截面受力大。

(3) 現(xiàn)場采用中空測力傳感器進(jìn)行測試，由于6個螺栓測點(diǎn)預(yù)緊力不同，受力曲線并未重合，但整體變化趨勢一致，這驗(yàn)證了頂蓋耦合面的水壓力為軸對稱分布的特效；螺栓在開機(jī)和停機(jī)兩個階段螺栓受力有明顯區(qū)別，開機(jī)正常運(yùn)行時螺栓受力比停機(jī)時增加了20～30 kN，這是由于開機(jī)時水壓力的加載左右；在較低負(fù)荷振動區(qū)，頂蓋螺栓受力最大較正常運(yùn)行高5 kN，在穩(wěn)定運(yùn)行區(qū)螺栓受力保持穩(wěn)定。

(4) 對比Hmax工況下的停機(jī)、40%、70%、100%導(dǎo)葉開度情況下螺栓中截面的平均值數(shù)值模擬和現(xiàn)場試驗(yàn)數(shù)值，停機(jī)及三個穩(wěn)定運(yùn)行工況下數(shù)模計算的數(shù)值和現(xiàn)場試驗(yàn)比較吻合，誤差不大于5%，驗(yàn)證了有限元計算的準(zhǔn)確性，說明中空測力傳感器能夠比較準(zhǔn)確的反應(yīng)螺栓的軸向受力。

本文開展了頂蓋螺栓數(shù)值及試驗(yàn)研究，驗(yàn)證了螺栓中截面的受力分析，但是螺栓受力的分析和現(xiàn)場測試還存在一些值得深入研究的問題：①中空測力傳感器能夠比較好的螺栓軸向受力，并且只能是平均值；有限元計算可以獲得螺栓危險位位置、不同截面應(yīng)力分布，如何在現(xiàn)場能測得螺栓的受力分布，這需要開發(fā)更高功能的受力傳感器;②現(xiàn)場測試在小導(dǎo)葉開度機(jī)組振動區(qū)螺栓受力比較大，通過有限元計算應(yīng)該在小開度計算的時候水壓力比較小，這存在不一致性，這主要是由于在不穩(wěn)定區(qū)存在壓力脈動，同時現(xiàn)場存在較大的機(jī)械振動，此時的數(shù)值模擬應(yīng)考慮到機(jī)械振動及壓力脈動等因素，如采用全數(shù)值模擬的方法研究涉及多場耦合的非定常計算問題。