【正文】 參數(shù)符號數(shù)值(cm)參數(shù)符號數(shù)值(cm)dadbdcd2l1l219232122l3l4l5Llalb1923lclpabc計算簡圖如下：圖47 導(dǎo)葉校核計算簡圖導(dǎo)葉的強度是以三點支梁的型式，應(yīng)用馬克士威—馬爾公式計算的，在本次設(shè)計中采用編程計算，其計算程序見附錄中的附錄1。復(fù)核若干斷面的程序見附錄2程序運行結(jié)果如下： m a t b 由結(jié)果看出復(fù)核的各個斷面的應(yīng)力都小于許用應(yīng)力。致 謝本文是在羅興锜教授的悉心指導(dǎo)下完成的，在此期間，廖偉麗老師、郭鵬程老師和鄭小波老師給予了極大的幫助并提出許多寶貴的意見，在此表示衷心的感謝。參考文獻(xiàn)［1］哈爾濱大電機研究所編著.《水輪機設(shè)計手冊》.北京：機械工業(yè)出版社，1976［2］劉大愷編著.《水輪機》.北京：中國水利水電出版社，［3］水電站機電設(shè)計手冊編寫組.《水電站機電設(shè)計手冊水力機械分冊》.北京：水利電力出版社，1983［4］武漢水利電力學(xué)院主編.《水輪機》.上冊：北京：水利電力出版社，1980［5］哈爾濱大電機研究所譯．《水輪機設(shè)計者手冊》．機械工業(yè)出版社，1976［6］高建銘等編著．《水輪機及葉片泵結(jié)構(gòu)》．清華大學(xué)出版社，1992［7］陳新方等編著．《水輪機結(jié)構(gòu)分析》．水利電力出版社，1994［8］天津電氣傳動設(shè)計研究所編．《水輪機結(jié)構(gòu)圖冊》．科學(xué)出版社，1978附 錄附錄1：導(dǎo)葉校核程序 program main parameter(m=10,n=5,pi=) real ra(m+1),rb(m+1),rc(m+1),a(n,m+1),b(n,m+1),qa(m+1),qb(m+1) $,qc(m+1),yc(m+1),k real jd,ja,jb,jc,j11,m11,m22,m33 h1= d0= b0=180 z0= jd= wd= dcc= dmm= p= e= dy= dc= an= r= da1= db1= dc1= d2= z1= z2= z3= z4= z5= z= za= zb= zc= zp= ps=*pi*d0*b0*h1/z0 ps0=ps/b0 pc=pi*p*(2*dcc**2dmm**2)/4 pp=pc*dy*(sin(pi*r/))/(z0*dc*(cos(pi*an/))) ja=pi*da1**4/ wa=pi*da1**3/ jb=pi*db1**4/ wb=pi*db1**3/ jc=pi*dc1**4/ wc=pi*dc1**3/ j11=pi*d2**4/ w11=pi*d2**3/ a1=jd/ja a2=jd/jb q1=z1/z q2=z2/z q3=q1+q2 q4=(z3+*z2)/z q5=z4/z q6=z5/z k=z**3/(e*jd) yc1=k*ps*q5*(q1**3*q4*a1/+(q4*q3**3/**3*q4/**3/ $q1*q2**2/)+(q4/**3*q4/**3/+q4**3/**2/2+ $q2**2*q3/)*a2) yc2=k*pp*q5*(q1**3*q6*a1+q6*(q3**3q1**3)+(q6+*q5*q6q5**2/2 $q3**3*q6)*a2)/ y=yc1+yc2 yc(1)= do 10 i=1,m+1 rc(i)=(yyc(i))*(k*q5**2*(q1**3*a1+q3**3q1**3+(**3+ $q5)*a2)) ra(i)=ps*q4+pp*q6rc(i)*q5 rb(i)=ps*(q1+*q2)pp*(1+q6)+rc(i)*(1+q5) qa(i)=ra(i)/(za*da1) qb(i)=rb(i)/(zb*db1) qc(i)=rc(i)/(zc*dc1) if() then yc(i+1)=yc(i)+*y endif if() then yc(i+1)=y endif10 continue tb=pp*zp/(2*wb) tc=pp*zp/(2*wc) do 20 i=1,n do 20 j=1,m+1 if() then a(i,j)=ra(j)*z1/wa elseif() then a(i,j)=(ra(j)*(z1+*z2)ps0*z2**2/)/wd elseif() then a(i,j)=(ra(j)*(z1+z2)ps0*z2**2/)/wb b(i,j)=sqrt(a(i,j)**2+4*tb**2) elseif() then a(i,j)=(pp*z5rc(j)*z4)/wb b(i,j)=sqrt(a(i,j)**2+4*tb**2) else a(i,j)=pp*(z5z4)/wc b(i,j)=sqrt(a(i,j)**2+4*tc**2) endif20 continue do 30 i=1,n if() then print*,39。,i do 40 j=1,m+1 print39。,yc(j),a(i,j)40 continue else print*,39。,i do 50 k=1,m+1 print39。,yc(k),b(i,k)50 continue endif30 continue do 70 j=1,m+1 print39。,yc(j),ra(j),rb(j) $,rc(j),qa(j),qb(j),qc(j)70 continue end附錄2：導(dǎo)葉部分?jǐn)嗝鎻?fù)核程序 program main parameter(m=10,n=5,pi=) real ra(m+1),rb(m+1),rc(m+1),a(n,m+1),b(n,m+1),qa(m+1),qb(m+1) $,qc(m+1),yc(m+1),k,rcc real jd,ja,jb,jc,j11,m11,m22,m33 h1= d0= b0= z0= jd= wd= dcc= dmm= p= e= dy= dc= an=6 r= da1= db1= dc1= d2= z1= z2= z3= z4= z5= z= za= zb= zc= zp=94 ps=*pi*d0*b0*h1/z0 ps0=ps/b0 pc=pi*p*(2*dcc**2dmm**2)/4 pp=pc*dy*(sin(pi*r/))/(z0*dc*(cos(pi*an/))) ja=pi*da1**4/ wa=pi*da1**3/ jb=pi*db1**4/ wb=pi*db1**3/ jc=pi*dc1**4/ wc=pi*dc1**3/ j11=pi*d2**4/ w11=pi*d2**3/ a1=jd/ja a2=jd/jb q1=z1/z q2=z2/z q3=q1+q2 q4=(z3+*z2)/z q5=z4/z q6=z5/z k=z**3/(e*jd) yc1=k*ps*q5*(q1**3*q4*a1/+(q4*q3**3/**3*q4/**3/ $q1*q2**2/)+(q4/**3*q4/**3/+q4**3/**2/2+ $q2**2*q3/)*a2) yc2=k*pp*q5*(q1**3*q6*a1+q6*(q3**3q1**3)+(q6+*q5*q6q5**2/2 $q3**3*q6)*a2)/ y=yc1+yc2 ycc= zaa= zbb= zcc= rcc=(yycc)*3/(k*q5**2*(q1**3*a1+q3**3q1**3+(1q3**3+q5)*a2)) m11=pp*(z5zaa)rcc*(z4zaa) a11=m11/w11 t11=pp*zp/(2*w11) b11=sqrt(a11**2+4*t11**2) m22=pp*(z5z4+zbb)rcc*zbb a22=m22/wc t22=pp*zp/(2*wc) b22=sqrt(a22*2+4*t22**2) m33=pp*(z5z4zcc) a33=m33/wc t33=pp*zp/(2*wc) b33=sqrt(a33**2+4*t33**2) write(*,100) write(*,200)m11,a11,t11,b11 write(*,200)m22,a22,t22,b22 write(*,200)m33,a33,t33,b33200 format(1x,)100 format(5x,39。,9x,39。,8x,39。,8x,39。) end附錄3：程序中各符號所代表的意義實際符號程序中代號實際符號程序中代號H’D0b0Z0JDWDdcdmpEDyDCβγdadbh1d0b0z0jdwddccdmmpedydcanrda1db1dcdel1l2l3l4l5Llalblclpabcdc1d2z1z2z3z4z5zzazbzczpzaazbbzcc