孟 微,趙國宏
(1.中國直升機設計研究所 直升機旋翼動力學重點實驗室,江西 景德鎮 333001;2.中國兵器工業試驗測試研究院,陜西 華陰 714200)
隨著復合材料技術的發展,直升機旋翼設計有了很多新的變化。20世紀60年代開始,無鉸旋翼和無軸承旋翼的陸續出現提高了旋翼的性能,進而提升了直升機的綜合性能。這些新型旋翼的誕生在帶來結構簡單、重量輕等諸多優點的同時,也帶來了復雜的動力學問題。傳統的直升機采用鉸接式槳轂,槳葉在揮舞、擺振、扭轉方向的彈性耦合很小,從動力學角度來看,這種設計的優點就在于用現有理論可以很好地計算出它的動力學性能。新型旋翼無鉸、無軸承的特性則會產生復雜的非線性揮舞、擺振、扭轉運動之間的耦合問題。旋翼動力學性能對這些耦合運動都極為敏感,因此在計算其動力學性能時,必須重新考慮這些問題。
利用旋轉槳葉模型準確地測量數據,驗證槳葉模型理論的適用性與可靠性在國外已開展多年并日趨成熟。本文就國外已有的旋翼槳葉彈性變形測量方法作簡單的介紹,綜合對比分析各方法的優劣。
槳葉彈性變形測量方法大致分為兩類,即光學測量方法和傳感器測量法。光學測量方法在測試領域是目前較為熱門的研究,是隨著激光、電子等技術的發展而發展起來的,目前來講是較為先進的技術。光學測量方法的特點是針對槳葉彈性變形的測量,使用光學儀器,不接觸測試件,且測試獲得的結果即為位移變化量,無需進行大量復雜的間接轉換。光學測量方法主要有:光柵照相法、立體圖像識別法和立體圖像識別改進法。
光柵照相法,non-intrusive optical fringe correlation method,簡稱FCM方法,又稱為projected grid method,PGM方法[1],該方法可以得到槳葉揮舞角和扭轉角結果,但無法得到擺振角。首先在旋轉中心與槳轂中心(坐標原點)相重合的地面水平轉臺兩端分別安裝光柵投影儀和照相機,安裝角度要求保證光柵和拍攝范圍覆蓋全部槳葉,且在整個試驗過程中保持不變。試驗時用光柵投影儀對整個槳葉表面投射光柵,然后由照相機拍攝下光柵變化。投影儀與Z軸的夾角為β,照相機與Z軸的夾角為α,由此可以得出揮舞方向位移(deflection)與光柵Y方向變形量(image displacement)的幾何關系(見公式1)。再根據槳葉不同半徑處的位移變化結果可以得出其揮舞角。槳葉扭轉角可直接根據其與光柵X方向夾角的幾何關系確定。
(1)
立體圖像識別法(stereo pattern recognition,簡稱SPR)是通過測量沿槳葉展向分布的各點在揮舞、擺振方向的位移量,從而獲得槳葉揮舞角、擺振角及槳葉彈性變形產生的扭轉角的方法。具體原理是用照相機分別拍攝裝有反光標記的槳葉在靜止狀態和旋轉狀態下的照片,將兩狀態下的照片對比獲得每一個固定反光點在旋轉坐標系X、Y方向的變化量,再根據X、Y方向變化量與槳葉揮舞、擺振位移的幾何關系(公式1)確定沿槳葉展向各點的揮舞、擺振位移,最終獲得揮舞角和擺振角。槳葉彈性扭轉角則是通過某一展向位置處前后緣兩點在Z方向的位移差值求出扭轉角,再減去槳葉根部的扭轉角和槳葉預扭角即為其彈性扭轉角,該方法的最大誤差為0.52°。
在文獻[2][3]HARTⅡ槳葉運動測量的試驗中,試驗臺下方以旋轉坐標原點(槳轂中心)為中心的同一半徑下分布四臺照相機,一對用于拍攝槳盤前行邊,另一對拍攝槳盤后行邊。分別在方位角90°和135°位置處安裝兩臺槳尖偵測照相機(BTD),用來確定槳尖在水平、垂直方向的位移變化。試驗前,將槳葉下表面涂黑并安裝白色反光標記,反光標記沿展向平均分布在槳葉前緣(0.228弦長)、后緣(0.993弦長)兩排,每排18個。試驗時每臺相機在每個狀態拍攝50至100幅照片,最終結果為所得數據的平均值。
照相測量法是為了提高測量精度在原有的SPR的方法上進行了改進,以獲得更高精度的槳葉揮舞角、擺振角及彈性扭轉角的方法。較原方法不僅增加了槳葉上的反光標記,而且在試驗風洞壁上按方陣高密度地安裝了大量的反光標記,并增加了照相機的數量(8臺)。測量結果不再是以槳轂中心為坐標原點的坐標系中的結果,而是在固定方陣坐標系中的,從而避免了因坐標系原點不穩定帶來的誤差。試驗照片的測量結果與揮舞、擺振位移的轉換即為固定方陣坐標系與旋轉坐標系的轉換。所以需要在旋翼槳轂上做反光標記以確定其在固定方陣中的位置。在UH-60旋翼系統全尺寸槳葉測量試驗中,采用的就是該方法[4]。
照相測量法的具體原理是:首先明確旋翼軸在風洞壁上光標方陣中的位置,確定其前傾角;然后將風洞中的光標和槳葉上的光標位置用非線性最小方差擬合的方法估算出在零相位時槳葉在固定方陣坐標系中(XYZ)的位置;槳葉的扭轉角將被投影到繞X軸轉動,揮舞角被投影到繞Y軸轉動,擺振角投影到排除前傾角后的繞Z軸轉動。
傳感器測量方法是比較傳統的測試方法,它的特點是測試設備簡單,成本較低,但測量獲得的信號需通過事先建立好的模型關系計算方可得到所需結果。該方法應用較早,于20世紀70年代由RAE(皇家航空中心)在測量旋轉槳葉的振動模態時建立并發展起來,在旋翼槳葉變形測量中有著不可替代的作用。在長期的試驗應用中,該類方法日趨成熟,已經在許多不同的試驗測量中得到成功的運用。該類方法主要有kinematical方法和傳統的應變測量方法。
kinematical方法[5]是一種理論計算與試驗測量相結合的方法,該方法通過將傳感器測量的線性位移或角位移信號用kinematical方法處理,即利用鉸的幾何矩陣與測量系統得到一種變換關系,計算出所需的槳葉揮舞、擺振、扭轉角。其中擺振、揮舞、扭轉角與傳感器線位移與角位移的關系表達式為kinematical方程。該方法比較成熟,已經在許多不同的飛機上得到成功運用。
該方法是把其他位移變化(傳感器可以簡單測量得到)應用合適的運算法則換算成角位移,從而取代直接測量三個方向的角位移。通過傳感器線位移與角位移求出擺振、揮舞、扭轉角,并將該矩陣求逆得到用該矩陣表達的關系式,即kinematical方程,其最大優點是大大減少了傳感器的安裝數量,有些地方甚至可以達到僅安裝一個傳感器便可完成測試的效果。其關鍵在于槳轂和槳葉上位移傳感器安裝位置的合理性。
應變電測量方法(簡稱SPA技術)是20世紀70、80年代由RAE[1-3]皇家航空研究中心發展的一套方法,是旋翼槳葉變形測量中最早使用的方法,同時也是目前力學試驗中應用最廣泛的、最為簡單的電阻應變測試法。
(2)
如果第二個要求也成立,那么應變響應也可表示為:
(3)

(4)


(5)
其中D是δ1,δ2,δ3,……δp的(nxp)矩陣。
總的來說,應變片的數量和分布是由所要測量的位移量的復雜程度決定的,比如測量旋翼槳葉的低階模態時所需要的應變片數量就比測量高階時少??紤]到噪聲限制,應變響應的數量r和所需測量位移點的數量n要大致相同,應變片的分布也要大致沿著槳葉結構分布。
通過前面對各種方法的詳細介紹,綜合分析兩大類槳葉彈性變形的測量方法,文中所闡述的光柵照相法、立體圖像識別法和立體圖像識別改進法等光學方法,如果沒有足夠復雜的光學設備,往往不能保證測量精度,也不適合于確定相位關系,而且試驗成本較高。
雖然SPA技術是早期發展起來的,但相比傳感器測量方法里各種不同的位移傳感器,應變片更容易安裝在槳葉上。其最大優點是質量小,低成本,不需要相位特征;缺點是由于槳葉的揮舞、扭轉耦合問題及應變片在槳葉上布置粘貼問題,很難保證復雜外形槳葉的測量精度,例如有后掠尖削的ERATO[1]槳葉。各種方法的優、缺點對比見表1。

表1 各種方法的綜合比較
總而言之,目前的旋翼槳葉彈性變形測量方法中,并沒有哪一種方法可以取代其他方法,各有優缺點。針對較為頻繁和常規的試驗來說,考慮到應變片成本和微小加速計及其他方法成本的比較,在直升機旋翼槳葉的懸停及風洞振動模型測量中,SPA技術仍是目前最為經濟實用且最常采用的方法。
[1] Muller R H G, Pengel K, Van der Wall B G. Blade Deflection Measurement at the Low Noise ERATO Rotor[C].26th EUROPEAN RORORCRAFT FORUM, The Hafue, Netherlands, September 2000:104.
[2] Schneider O, Van der Wall B G, Pengel K. Final Analysis of HART-II Blade Deflection Measurement[C].29th European Rotorcraft Forum, Friedrichshafen, Gemary, September 2003.
[3] Schneider O, Van der Wall B G, Pengel K. HART-II Blade Motion Measured by Stereo Pattern Recognition (SPR)[C]. 59th Annual Forum of the American Helicopter Society, Phoenix, USA, May 6-8 2003.
[4] Olson L E, Abrego A I, Barrows D A, et al. Blade Deflection Measurements of a Full-Scale UH-60A Rotor System[C]. The Ameri--can Helicopter Society Aeromechanics Specialist’ Conference, January 20-22 2010.
[5] Attilio C, Alessandro L. Measuring Blade Angular Motions: a Kinematical Approach[C].30th European Rotorcraft Forum,Augusta 2004.
[6] Gaukroger D R, Hassal G J W. Measurement of Vibratory Displacement of a Rotating Blade[J].Vertica , 1978(2):111-120.
[7] Walker A R. Further Application and Development of Strain Pattern Analysis[C].8th European Rotorcraft and Powered Lift Forum, 31 August-3 September,1982.
[8] Tourjansky N, Szechenyi E. The Measurement of Blade Deflections A New Implementation of the Strain Pattern Analysis[C].Associantion Aeronautique et Astronautique de France,1992.