MatIab計(jì)算在水電站機(jī)組盤車及軸線調(diào)整中的應(yīng)用

賈棟林

（中國(guó)水利水電第十四工程局有限公司，云南 昆明650032）

1 概述

厄瓜多爾科卡科多辛克雷（Coca codo Sinclair）水電站是我國(guó)對(duì)外承建的EPC合同中最大的沖擊式水電站建設(shè)合同，其電站地下廠房裝設(shè)8臺(tái)單機(jī)容量205 MVA懸式立軸六噴嘴沖擊式水輪發(fā)電機(jī)組，額定功率187.5 MW，額定轉(zhuǎn)速300 r/min，飛逸轉(zhuǎn)速530 r/min，額定水頭604.1 m，是我國(guó)目前生產(chǎn)及厄瓜多爾國(guó)單機(jī)容量最大的立軸沖擊式水輪機(jī)組。

機(jī)組盤車及軸線調(diào)整部位包括：集電環(huán)、推力頭、鏡板、轉(zhuǎn)子、水機(jī)軸、轉(zhuǎn)輪。機(jī)組盤車及軸線測(cè)量剖面如圖1所示。

2 機(jī)組盤車

機(jī)組盤車采用在轉(zhuǎn)子軸端部安裝人工盤車裝置，按順時(shí)針?lè)较蛲苿?dòng)機(jī)組轉(zhuǎn)動(dòng)部分旋轉(zhuǎn)。盤車前在每個(gè)徑向擺度測(cè)量部位的+X、+Y方向設(shè)置百分表，在每個(gè)軸向跳動(dòng)測(cè)量部位的+Y方向設(shè)置百分表。抱緊一半數(shù)量的上導(dǎo)瓦，瓦與軸頸間隙控制在0.02～0.04 mm。盤車時(shí)啟動(dòng)高壓油頂起裝置，當(dāng)鏡板底面與推力瓦接觸面間形成一層均勻油膜后，先連續(xù)推動(dòng)旋轉(zhuǎn)1圈，然后按標(biāo)記點(diǎn)進(jìn)行8等分逐點(diǎn)盤車，并記錄各盤車點(diǎn)停穩(wěn)后的靜態(tài)數(shù)據(jù)。盤車過(guò)程中，在推力頭上放置框式水平儀，推力頭在整個(gè)盤車過(guò)程中的水平度應(yīng)≤0.02 mm/m，若有超差，應(yīng)微調(diào)推力瓦受力，使其滿足要求。

3 應(yīng)用MatIab軟件分析盤車數(shù)據(jù)

圖1 機(jī)組盤車及軸線測(cè)量剖面圖

本文以4號(hào)機(jī)組為例介紹機(jī)組軸線數(shù)據(jù)分析方法，機(jī)組盤車數(shù)據(jù)如表1所示，表中省略了集電環(huán)、發(fā)電機(jī)法蘭、水輪機(jī)法蘭及轉(zhuǎn)輪軸向盤車數(shù)據(jù)。

表1 4號(hào)機(jī)組盤車數(shù)據(jù)單位：0.01 mm

從表1的數(shù)據(jù)可以得出鏡板最大擺度為0.08 mm，與推力頭同心度為0.04 mm，滿足后續(xù)盤車要求。根據(jù)鏡板相對(duì)上導(dǎo)的理論擺度為正弦曲線這一概念，假設(shè)其曲線為Y=A+B×sin（X+C），根據(jù)表1得出鏡板徑向相對(duì)上導(dǎo)徑向的擺度如表2所示。

表2 鏡板徑向相對(duì)上導(dǎo)徑向的擺度值單位：0.01 mm

然后采用Matlab軟件的lsqcurvefit函數(shù)對(duì)鏡板相對(duì)上導(dǎo)擺度值進(jìn)行最小二乘非線性擬合，其表達(dá)式如下：

clear all；

x=[0,45,90,135,180,225,270,315]；

y=[1,4,10,14,17,12,7,1]；

fun_curve=@(a,t) a(1)+a(2)×sin((t+a(3))× pi/180)；

x0=[0,45,90]；

parameter = lsqcurvefit(fun_curve，x0，x，y)；

A=parameter(1)；

B=parameter(2)；

C=parameter(3)；

subplot(2,1,1)；

plot(x,y,x,fun_curve(parameter,x),'r')；

xlabel(‘角度’)；

title(‘y值’)；

title(‘y=A+B×sin[(x+c)×pi/180]曲線擬合’)；

legend(‘原始曲線','擬合后曲線’)；

subplot(2,1,2)；

fai=0:1:360；

y_value=A+B×sin((fai+C)×pi/180）；

plot(fai,y_value)；

xlabel(‘角度’)；

ylabel(‘Y值’)；

title(‘y=A+B×sin[(x+c)×pi/180]曲線擬合’)；

legend(‘理論曲線’)；

通過(guò)擬合程序在Matlab中運(yùn)行，可以在Matlab工作區(qū)中得到 A=8.25，B=-7.89，C=102°。當(dāng) a=0°、45°……315°時(shí)其理論擺度值為y_value。其測(cè)量擺度值與理論擺度值偏差很小，如圖2和表3所示；其理論曲線是正弦曲線，如圖3所示；綜合以上結(jié)果，可判斷盤車數(shù)據(jù)準(zhǔn)確，可以作為軸線調(diào)整的依據(jù)。

圖2 Matlab工作區(qū)顯示計(jì)算結(jié)果

表3 測(cè)量擺度值與理論擺度值偏差表單位：0.01 mm

圖3 Matlab機(jī)組盤車數(shù)據(jù)擬合曲線

4 機(jī)組軸線調(diào)整

當(dāng)機(jī)組整體盤車出現(xiàn)軸線與其旋轉(zhuǎn)中心線不重合時(shí)，計(jì)算下導(dǎo)、水導(dǎo)靜擺度，其方法是將下導(dǎo)和水導(dǎo)各點(diǎn)靜擺度值分解成X，Y分量：

x=r×cos（angle×PI/180）

y=r×sin（angle×PI/180）

然后根據(jù)擺度值求解下導(dǎo)、水導(dǎo)與機(jī)組軸線的理論中心偏差，通過(guò)分?jǐn)嗝婀文グl(fā)電機(jī)推力頭處的鎖緊卡環(huán)進(jìn)行機(jī)組軸線的調(diào)整，即根據(jù)計(jì)算發(fā)電機(jī)推力頭處卡環(huán)刮磨量，達(dá)到下導(dǎo)和水導(dǎo)處的擺度值盡可能的小。

δ-卡環(huán)最大點(diǎn)刮磨量；

D-推力頭的直徑；

Φba-軸線傾斜方向上的最大靜擺度；

L-推力頭距水導(dǎo)處百分表的距離。

5 結(jié)束語(yǔ)

伴隨著科學(xué)技術(shù)的高速發(fā)展，Matlab計(jì)算已成為工程領(lǐng)域應(yīng)用最為廣泛的軟件之一，其強(qiáng)大的計(jì)算平臺(tái)，可以解決實(shí)際工程中最優(yōu)化問(wèn)題，也可通過(guò)對(duì)算法的編程實(shí)現(xiàn)相應(yīng)的最優(yōu)化計(jì)算與仿真，能夠有效指導(dǎo)和解決水電站施工中遇到的實(shí)際問(wèn)題，幫助技術(shù)人員完成復(fù)雜計(jì)算與仿真等工作。