對轉螺旋槳升力面設計方法
第10卷第2期2006年4月
文章編號:1007-7294(2006)02-0040-07
船舶力學
JournalofShipMechanics
Vol.10No.2Apr.2006
對轉螺旋槳升力面設計方法
辛公正,丁恩寶,唐登海
(中國船舶科學研究中心,江蘇無錫214082)
摘要:本文針對對轉螺旋槳這類組合式推進器的設計問題,采用了升力面設計方法進行設計,在設計過程中通過前后槳的相互迭代求出彼此間的誘導速度,以此來考慮它們之間的相互作用,并采用此方法作了一對對轉螺旋槳的設計,同時應用面元法進行了水動力預報。關鍵詞:對轉螺旋槳;升力面方法;面元法;誘導速度中圖分類號:U661.31+3
文獻標識碼:A
Adesignmethodforcontra-rotatingpropellerby
lifting-surfacemethod
XINGong-zheng,DINGEn-bao,TANGDeng-hai
(ChinaShipScientificResearchCenter,Wuxi214082,China)
Abstract:Acontra-rotatingpropellerdesignmethodbylifting-surfacetheoryispresented.Itistoconsid-ertheinducedvelocitythroughtheiterationbetweentheforwardpropellerandtheafterpropeller.Usingthismethod,acontra-rotatingpropellerwasdesigned,thenthehydrodynamicwithpanelmethodwaspre-dicted.Theresultsshowthatthismethodcanmeetengineeringrequirements.
Keywords:contra-rotatingpropeller;lifting-surfacemethod;panelmethod;boundarycondition;
inducedvelocity
1前
言
對轉螺旋槳通常具有高效率,易于改善航行體操縱性的優點,其設計方法雖然已經發展了很多年,但是隨著市場需求的發展,例如水下航行器導管對轉槳、AUV對轉槳、Z推以及對轉—吊艙推進器等等,它們的螺旋槳運轉方式雖然也是相對旋轉的,但其推進器組成結構中還多了其它的附屬結構—導管、支柱,這就需要我們在設計過程中更加充分地考慮到各組成部件的相互作用。因此,我們首先需要解決的是如何把前后槳之間的相互干擾引入到設計過程中,而升力面設計方法在螺旋槳的設計過程中仍被廣泛采用。對于常規的對轉螺旋槳是在同一軸上由前后兩個旋轉方向相反的螺旋槳組成。這樣可起到兩個作用,一是推力負荷由兩個槳分擔,故每個槳的推力負荷減輕了,另一個作用是前槳尾流的旋轉動能損失可被后槳所利用而得到回收,致使二槳合成的尾流旋轉能量損失減小,從而可以提高推進器的效率,同時前后槳的扭矩也可以相互平衡增加航行體的穩定性;但由于軸系結構復雜,在大船上難以采用,大多數用于水下航行器上。收稿日期:2005-11-21
作者簡介:辛公正(1978-),男,中國船舶科學研究中心工程師。
第2期辛公正等:對轉螺旋槳升力面設計方法41
對轉螺旋槳是組合式推進器的一種。所謂組合式推進器是由幾個推進部件所組成。對轉螺旋槳即由前、后槳兩個推進部件所組成。此外還有串列式螺旋槳,導管螺旋槳等組合式推進器。組合式推進器的水動力學問題的復雜性在于每一部件對其它部件有擾動速度,也稱為誘導速度。如果部件之間有相對運動,如對轉螺旋槳的前后槳、帶前置或后置導葉的轉定子之間則干擾速度是非定常的,只可能是個非定常問題。對于這類問題有二個處理方法。一種是在平均意義下求解,即把干擾速度進行周向平均,或者說把在一個周期內的平均值作為干擾速度,即把它假定為與周向坐標或者說與時間無關。如此在周向不變的來流條件下,則問題作為定常問題來處理。另一種方法是真正按瞬時的干擾速度條件作為非定常問題處理。早在1976年董世湯和陳立強發表了對轉螺旋槳的升力線理論計算加升力面修正的設計方法[1],并應用于產品型號的設計,滿足了當時的工程需要。本文的對轉螺旋槳設計方法是以升力面理論為基礎建立的設計方法,針對前槳和后槳,分別用各自的程序進行設計,在二者之間通過迭代來考慮相互的影響。
2升力面設計方法
2.1進流場的分析
設計過程中前、后槳采用適伴流設計,伴流場分布只考慮為半徑r的函數,軸向位置取槳盤面位置,周向采用平均。前、后槳之間的干擾速度也是取各自槳盤面位置進行周向平均。由此將各個槳的伴流分布加上另一個槳在槳盤面的誘導速度可得到實效伴流場的分布,在此速度場下進行螺旋槳的升力面設計,以此來考慮二者之間的相互影響。給定設計條件及設計要求2.2坐標系、
在設計計算中,前槳和后槳均定義為右旋槳,這里建立了兩個坐標系:一個是柱坐標系(x,r,!),與螺旋槳固定在一起,以槳葉參考線與槳軸交點為坐標原點,x仍以向下游為正,以1號葉片的槳葉參考線定義!=0,!方向按右手定則定義正方向,r向外為正;另一個是直角坐標系(x,y,z),x方向與柱坐標系相同,y方向是豎直向上為正,z方向是與x,y方向滿足右手定則。
在進入升力面設計之前,給定的設計條件通常有:伴流分數分布,直徑D,葉片數Z,航速VS,轉速
n,初始螺距角"p0,葉片的弦長分布c(r),葉片的厚度分布t(t,s),縱斜分布ZR(r),側斜角分布!s(r),!(s),設計的水動力要求。徑向環量分布形式F(r)和弦向環量分布形式#m
通過升力面設計計算,希望達到的設計要求為:設計得到前后槳葉片各葉剖面的螺距角及拱弧線形狀,使前后槳產生所需要的總推力T或者用盡主機功率。通過迭代計算使前后槳的扭矩達到平衡,最后疊加厚度分布得出型值。
2.3升力面設計方法的數值模型
螺旋槳的'升力面設計計算采用渦格法,升力面建立在拱弧面上,邊界條件采用了能考慮徑向速度和側斜、縱傾分布的比較完善的邊界條件。由于在設計計算的迭代過程中,拱弧面的幾何形狀在改變,故每次迭代需要修改升力面的幾何條件,把網格重新劃分。對于升力面設計的渦格法這里不再細
[4]
述,具體方法可參見以及關于邊界條件分析比較[5]。“與粘流耦合的導管螺旋槳升力面設計方法”
3對轉槳設計
3.1表達式定義
對轉槳與常規螺旋槳在數據表達上有一些區別,為方便計算和分析,下面給出本文的約定及采用的對轉槳推扭力系數、效率等定義:
42
下標1———前槳;——前槳直徑(m);D1—
——螺旋槳進速(m/s);V—
——前槳推力(N);T1—
——前槳扭矩(N?Q1—m);———水的重量密度(N/m3);ρ
——前槳敞水效率;η1—
船舶力學
下標2———后槳;——后槳直徑(m);D2—
——總推力(N);T—
——后槳推力(N);T2—
——后槳扭矩(N?Q2—m);——總的敞水效率;η—
——后槳敞水效率;η2—
第10卷第2
期
——螺旋槳轉速(rps),前后槳轉速方向相反;Q———總扭矩(N?n—m);
KT1=
T1!nD1
Q1!nD1
KT2=KQ2=
T2!nD1Q2
KQ1=
!nD1
KT=KT1+KT2K"1=JT1
Q1K
%=JT
Q
3.2設計過程
KQ=KQ1+KQ2K
"2=JT2
Q2J=V1
采用上述對轉槳設計程序進行設計,首先輸入前后槳的主參數,先進行前槳的設計,此時后槳在前槳槳盤面處的誘導速度u21(r)設為零,前槳的進流條件即為前槳的伴流分數w1(r),通過對前槳的設計計算可求出前槳在后槳槳盤面處的誘導速度u12(r),將其作為前槳對后槳的干擾;然后進行后槳的設計,此時后槳的進流條件為后槳的伴流分數w2(r)加上前槳在后槳槳盤面處對后槳的誘導速度u12(r),同樣通過對后槳的設計計算可求出后槳在前槳槳盤面處的誘導速度u21(r),將其作為后槳對前槳的干擾;如此前槳的進流條件發生了變化,變為前槳的伴流分數w1(r)加上后槳在前槳槳盤面處對前槳的誘導速度u21(r),因此前槳需要重新進行升力面設計,計算出新的對后槳的誘導速度u12(r),由此可見這個過程需要相互迭代直到前槳與后槳間的相互作用即誘導速度收斂后為止,這樣可認為在前
圖1對轉螺旋槳升力面設計流程圖
Fig.1Flowchartofcontra-propellerliftingsurfacedesign
后槳的一體化設計中比較充分地考慮了相互間的干擾。具體的迭代流程圖如圖1所示,其中對于每一次前后槳的升力面設計都是在各自進流條件下迭代到滿足水動力要求結束的,因此整個設計過程實際上有兩層迭代,一層是升力面設計自身的迭代,另一層是前后槳相互干擾的迭代。
第2期辛公正等:對轉螺旋槳升力面設計方法43
4面元法水動力預報
面元法作為螺旋槳的水動力性能預報工具已經發展得比較成熟,其計算精度已經能滿足各類工在程問題的要求。因此升力面設計結束后,我們采用CSSRC面元法[6,7]對設計結果進行了水動力預報。面元法的預報過程中采用與升力面設計方法類似的迭代過程來考慮前后槳之間的相互作用。
5實例計算
本文采用前面所述的設計預報方法,分別對某水下航行器設計了一對對轉螺旋槳。
5.1算例設計參數和設計條件
前后槳的主參數和設計條件,分別見表1、表2和表3。
表1算例前槳主參數
表2算例后槳主參數
Tab.1Forwardpropellerparameter
ofthescheme
r/R0.2920.30.40.50.60.70.80.90.951
C/D0.13260.13400.15360.17760.19880.21440.22300.20440.16400.0500
T/C0.25150.18220.12830.09020.06260.04290.03080.02820.06000.2515
T/D0.03340.03300.02800.02280.01800.01340.00960.00630.00460.0030
F(r)0.35880.38470.64180.87291.00001.02940.95240.71650.47470.0000
wx(r)0.51000.50160.46540.36350.27210.19270.12730.07740.04480.0355
r/R0.2350.30.40.50.60.70.80.90.951
Tab.2Afterpropellerparameter
ofthescheme
C/D0.11800.14120.17840.21050.24180.26210.26610.24420.19570.0500
T/C0.30870.23010.14840.09920.06630.04540.03280.02410.02250.0596
T/D0.03640.03250.02650.02090.01600.01190.00870.00590.00440.0030
F(r)-0.2394-0.4680-0.6920-0.8731-0.9766-1.0000-0.9251-0.6560-0.36570.0000
wx(r)0.46020.39860.31370.23370.16190.10150.05570.02780.02150.0209
5.2算例設計和預報結果
螺旋槳設計采用最佳環量分布,以獲得較高的效率。前槳平均伴流分數為0.19,后槳平均伴流分數為0.17。設計時以滿足扭矩KQ要求為目標,并使得前后槳的扭矩達到比較好的平衡,經過迭代計算,最終獲得前后槳設計方案,如圖2;前后槳螺距比、拱度比的分布如圖3和4;設計迭代過程中前后槳相互間的干擾如圖5~10,從圖中可見隨著迭代的進行前后槳之間的相互作用的誘導速度很快收斂,同時由于后槳是工作在前槳的尾流場中,前槳對后槳的誘導速度更大且干擾更強,而后槳的前槳的干擾主要體現在軸向誘導速度上,
葉數直徑(m)轉速(rpm)
表3算例設計條件
Tab.3Designconditionsofthescheme
前
槳
后
槳
50.4165000.436
右
40.4085000.410
左
Ao/Ae
旋向航速Vs(kn)
50.74190.02207
50.7419-0.02207
Js
設計KQ剖面形式
NACA66moda=0.8
切向誘導速度幾乎為零。設計結束后利用面元法分別對該方案進行均流狀態的水動力預報,預報時相對于設計狀態下進速系數為0.601,預報時考慮前后槳的相互影響,計算網格如圖11,預報結果見表
4。預報結果均按前槳進行無量綱化,從預報結果可以看到,由于設計過程中是在考慮前后槳相互作用
通過迭代后設計出來的,前后槳的KT、KQ預報值與設計值相比偏差在3%左右,其中KQ值比較接近所要求的結果,后槳的扭矩偏大,前槳扭矩偏小,二者的扭矩不平衡度為5%,能夠滿足工程的需要。
44船舶力學第10卷第2期
圖2算例對轉槳設計方案圖3前后槳螺距比分布
Fig.2Designschemeofthecontra-rotatingpropellerFig.3Thepitchratiodistributionofforwardand
afterpropeller
圖4前后槳拱度比分布圖5設計迭代過程中前槳對后槳的軸向誘導速度
Fig.4Thecamberratiodistributionofforward
andafterpropeller
Fig.5Theaxialinducedvelocityofforwardpropellerto
afterpropelleriniteratedprocess
圖6設計迭代過程中前槳對后槳的徑向誘導速度圖7設計迭代過程中前槳對后槳的切向誘導速度
Fig.6Theradialinducedvelocityofforwardpropeller
toafterpropelleriniteratedprocessFig.7Thetangentialinducedvelocityofforwardpropeller
toafterpropelleriniteratedprocess
第2期辛公正等:對轉螺旋槳升力面設計方法45
圖8設計迭代過程中后槳對前槳的軸向誘導速度圖9設計迭代過程中后槳對前槳的徑向誘導速度
Fig.8Theaxialinducedvelocityofafterpropellerto
forwardpropelleriniteratedprocess
Fig.9Theradialinducedvelocityofafterpropellerto
forwardpropelleriniteratedprocess
圖10設計迭代過程中后槳對前槳的切向誘導速度圖11面元法預報網格的劃分
Fig.10Thetangentialinducedvelocityofafterpropellerto
forwardpropelleriniteratedprocess
Fig.11Meshdistributionofthepanelmethod
表4對轉槳水動力設計和預報結果
Tab.4Thedesignandpredictionresultsofthecontra-rotatingpropellerhydrodynamic
設計值
預報值
偏差
扭矩值(N.m)
———
扭矩不平衡度
前槳
KT10KQ
0.14100.22070.1297-0.2206
0.13560.21520.1280-0.2261
-3.81%-2.51%-1.30%2.50%
19.0802
-
5.09%
———
后槳
KT10KQ
-20.0520
6結論
我們針對某水下航行器,采用對轉螺旋槳升力面設計程序進行了方案設計,并用面元法進行了水動力預報,在設計和預報過程中通過前后槳的相互迭代來考慮它們之間的影響,通過對整個過程的分析,我們有以下建議和結論:
(1)根據上述的理論分析和實例計算,本文建立的對轉螺旋槳升力面設計方法和面元法預報方
46
法能夠滿足工程問題的需要。
船舶力學第10卷第2期
(2)對于對轉螺旋槳升力面設計方法和面元法預報方法,由于在理論上的數值模擬與實際螺旋槳工作環境有顯著差別,因此其計算精度還有待改善,這只有通過更多的設計實例和相應的試驗分析比較才能實現。
(3)CFD方法在近幾年發展得比較成熟,將其應用在對轉螺旋槳這類組合式推進器的預報分析上應該更能捕捉到前后槳之間的這種強干擾的粘性流動,因此用CFD方法來計算對轉螺旋槳問題應該在以后得到更多的研究。參考文獻:
[1]董世湯,陳立強.對轉螺旋槳環流理論計算方法[J].艦船性能研究,1976(2).
[2]YangChenjun,TamashimaM,WangGuoqiang,YamazakiR.Predictionofthesteadyperformanceofcontra-rotating
prpellersbyliftingsurfacetheory[J].TransactionsoftheWest-JapanSocietyofNavalArchitects,1991(82).
[3]YangChenjun,TamashimaM,WangGuoqiang,YamazakiR,KoizukaH.Predictionoftheunsteadyperformanceof
contra-rotatingprpellersbyliftingsurfacetheory[J].TransactionsoftheWest-JapanSocietyofNavalArchitects,1992(83).
[4]辛公正.與粘流耦合的導管螺旋槳升力面設計方法[A].船舶水動力學學術會議論文集[C].武漢,2004.[5]辛公正,唐登海,董世湯.螺旋槳升力面設計邊界條件的處理分析[J].船舶力學,2004,8(2):16-24.
[6]TangDenghai,DongShitang,ChenJiadong,ZhouWeixin.ComparativecalculationsofpropellerperformancebyRANS/
PANELmethod[A].Proceedingsof22ndITTCPropulsionCommitteePropellerRANS/PanelMethodWorkshop[C].France,1998.
[7]董世湯,唐登海,周偉新.CSSRC螺旋槳定常面元法[J].船舶力學,2005,9(5):46-60.
【對轉螺旋槳升力面設計方法】相關文章:
飛艇螺旋槳氣動優化設計07-19
升力是怎樣產生的06-15
高升力翼型的氣動優化設計和實驗研究08-02
基于Kirchhoff方法的跨音速螺旋槳的氣動聲學計算07-07
飛機的升力是怎樣產生的 -工程01-01
升力系數不變時跨音速機翼阻力優化設計08-15
大型螺旋槳靜壓支承靜平衡儀設計研究11-30
機場瀝青道面結構設計方法與指標探討07-15
健身應該怎么提升力量02-17