大口徑望遠(yuǎn)鏡風(fēng)載分析綜述*

安其昌，張景旭，楊 飛，趙宏超

(中國科學(xué)院 長春光學(xué)精密機械與物理研究所，吉林 長春 130033)

0 引言

隨著大口徑望遠(yuǎn)鏡的發(fā)展，風(fēng)載的影響也得到了越來越多的關(guān)注。在上個世紀(jì)八十年代之前，對于風(fēng)載的認(rèn)識僅限于如何使之與望遠(yuǎn)鏡不發(fā)生直接作用，即圓頂?shù)拈_口大小與數(shù)量都控制得盡量低。但是隨后的研究表明，一定程度的空氣流動對于改善視寧度，所以圓頂?shù)念愋鸵灿蛇^去的封閉型逐漸演化為敞開與半敞開型［1-5］。

近年來對于風(fēng)載的影響主要集中在3個方面:①望遠(yuǎn)鏡周圍地形對于望遠(yuǎn)鏡所受風(fēng)載的影響;②望遠(yuǎn)鏡圓頂內(nèi)部的布局與結(jié)構(gòu)對于作用到望遠(yuǎn)鏡上的載荷的影響;③不同的風(fēng)載對于望遠(yuǎn)鏡的機械結(jié)構(gòu)與光學(xué)性能產(chǎn)生的不同影響。值得一提的是，第一方面的研究，主要是在望遠(yuǎn)鏡選址的時候進行參考，但是由于其采樣時間與傳感器排布的限制，造成對站址的評價帶有一定的局限性，往往對于最后的選擇參考的意義有限，同時站址的選擇是一個集合交通、地質(zhì)，經(jīng)濟、文化的復(fù)雜過程，不能僅考慮風(fēng)載這一項。對于第二方面，不同的圓頂結(jié)構(gòu)對于外界風(fēng)載的主要作用在于改變風(fēng)載中不同頻率分量所攜帶的能量，一般來說，當(dāng)氣流通過阻擋物時，其低頻分量會明顯削弱，相應(yīng)地，其中高頻分量會提高。對于第三方面的研究，是近年來研究的熱點，也是本研究論述的重點［6-9］。

研究風(fēng)載對于望遠(yuǎn)鏡結(jié)構(gòu)與成像能力的影響，主要的目的一方面是對于之前已由的理論的驗證與發(fā)展，另一方面是對于下一代大口徑望遠(yuǎn)鏡的建設(shè)進行指導(dǎo)。目前對于大口徑望遠(yuǎn)鏡風(fēng)載分析的主要方法有功率譜法、流體有限元軟件分析法以及實測標(biāo)定的方法。功率譜法已經(jīng)在大量的論述中進行了論述，但是實際上，由于所給出的結(jié)果只能具有統(tǒng)計學(xué)上的意義，針對特殊的工況并未給予足夠的說明，在本研究的第1 節(jié)中，對于功率譜實際的應(yīng)用與局限做出了一定的說明。另一方面，對于流體軟件的分析方法，其主要的問題在于建模的精度問題，該問題的具體論述會在第2 節(jié)展開;最后，對于實際的望遠(yuǎn)鏡建設(shè)來說，現(xiàn)場的測量的意義毋庸置疑［10-17］。雙子星望遠(yuǎn)鏡(GEMINI)是加拿大、智利、巴西、阿根廷和澳大利亞共同出資建造的8 m 級大口徑望遠(yuǎn)鏡，在美國與南美各有一臺，其中GEMINI South 在組裝的過程中，Dr.cho 及其團隊對于安裝了模擬鏡的望遠(yuǎn)鏡進行了動力學(xué)響應(yīng)測試。系統(tǒng)的激勵源為外界風(fēng)載，風(fēng)載的測量使用的是位于模擬主鏡的風(fēng)壓計;同時，系統(tǒng)的響應(yīng)被若干位于望遠(yuǎn)鏡上的加速度計拾取。在測量過程中，Dr.cho 對于實際測量的結(jié)果與傳統(tǒng)的分析結(jié)果進行了對比，得出了很多有價值的結(jié)論［18-19］，該部分會在第3 節(jié)進行論述。

本研究分別對3 種主流的風(fēng)載研究進行闡述，對于大口徑望遠(yuǎn)鏡進行風(fēng)載分析。

1 基于功率譜的風(fēng)載分析

由經(jīng)典的譜分析理論認(rèn)為，風(fēng)載作用主要分為靜風(fēng)載與動態(tài)風(fēng)載，靜風(fēng)載又稱為靜風(fēng)壓。兩者的表達(dá)式如下所示:

式中:P—鏡面上的靜風(fēng)壓，ρ—空氣密度，c—基于鏡面尺寸以及空氣雷諾數(shù)的系數(shù)，v—風(fēng)速。

由上式可知，在該理論中，是假設(shè)風(fēng)載一直處于單向作用的情況下，故通過平均作用可以得到該表達(dá)式，但是根據(jù)實際的測量的結(jié)果，由于望遠(yuǎn)鏡結(jié)構(gòu)的影響，風(fēng)載的方向會發(fā)生變化，在較長時間的采用情況下，平均的風(fēng)壓趨近于零。這與沒有考慮到風(fēng)載反向的原始假設(shè)相距甚遠(yuǎn)。

實際中分析動風(fēng)載的統(tǒng)計學(xué)模型有許多種，比如Kolmogrov 譜與改進Kaimal 譜。改進Kaimal 譜如下所示［20-22］:

地面摩擦剪切風(fēng)速:

式中:z0—地面粗糙度。

歸一化無量綱頻率:

式中:n—頻率，U(z)—離地高度z 處的平均風(fēng)速。

當(dāng)風(fēng)作用到大型結(jié)構(gòu)上，比如大型天文望遠(yuǎn)鏡，眾多的結(jié)構(gòu)與建筑，風(fēng)速的變化比之前假設(shè)的要復(fù)雜的多。

同樣，由于沒有考慮風(fēng)載的方向變化，在大口徑系統(tǒng)中，主鏡的正面與背部的差異越來越大，如使用譜進行分析的時候，在風(fēng)載施加與鏡面背部的情況下，實際的響應(yīng)譜會出現(xiàn)額外的尖峰，額外的尖峰如圖1 所示。

圖1 理論譜線與實際譜線

即使在不考慮風(fēng)向復(fù)雜變化的情況下，風(fēng)壓譜也不一定可以完美地符合統(tǒng)計學(xué)模型，實踐表明，只有當(dāng)風(fēng)速較快并且圓頂完全打開的情況下，主鏡鏡面上的譜才可以比較好地符合卡爾莫羅夫譜。其他的情況或多或少會產(chǎn)生一定量的偏移?？柲_夫譜譜線與譜線擬合如圖2 所示。

圖2 譜線與譜線擬合

2 CFD 軟件分析

現(xiàn)在普遍的思路為建立粗略的模型，之后使用CFD 軟件來計算望遠(yuǎn)鏡各個重要邊界上的載荷，之后再建立詳細(xì)的有限元模型，并將之前的所得到的邊界條件施加在詳細(xì)的有限元模型之上，之后得到結(jié)構(gòu)的響應(yīng)。流體模型計算如圖3 所示。問題的重點在于施加外部載荷點的選取以及初始簡化模型的簡化程度。這里涉及到了模型簡化的問題，如果模型過于簡化，在流體軟件里的計算相對簡單，但是對于詳細(xì)有限元模型，其計算精度就很難保證;相對應(yīng)地，如果在流體軟件中的模型十分詳細(xì)，將對應(yīng)的關(guān)鍵節(jié)點導(dǎo)入有限元模型的時候，會得到較好的分級結(jié)果，冗余的自由度會大為減少，但是，流體軟件的計算能力就會成為制約一切的瓶頸。

圖3 流體模型計算

值得一提的是，使用風(fēng)洞技術(shù)來獲得施加在望遠(yuǎn)鏡模型上的載荷，之后再進行分析，可以在一定程度上解決這個問題，但是氣流雷諾數(shù)的差別還是會影響最后模型的準(zhǔn)確度，故對于實際站址及望遠(yuǎn)鏡的測量是最準(zhǔn)確、最直觀的。考慮到實際測量也需要在特征點上布置傳感器，這一過程與望遠(yuǎn)鏡動力學(xué)模型減縮的過程極其相似，故對于實際望遠(yuǎn)鏡的測量，又可以反過來推動仿真技術(shù)的發(fā)展。

3 實際望遠(yuǎn)鏡的標(biāo)定

GEMINI 的有限元模型以及傳感器布置:在主鏡上有13個風(fēng)壓傳感器，在獲得了了風(fēng)壓數(shù)據(jù)之后，利用前八階Zernike 多項式來擬合，進而分析風(fēng)載輸入和輸出的關(guān)系。

對于系統(tǒng)的響應(yīng)，本研究采用了多個加速度計來拾取系統(tǒng)的三維運動。具體來說，在次鏡的底部沿光軸方向安放三向加速度計，這是因為對于系統(tǒng)來說，地基是一個極大質(zhì)量的剛體，故當(dāng)有外部激勵時，靠近底部的響應(yīng)最小，相應(yīng)的，經(jīng)過整個系統(tǒng)之后，次鏡末端的傳感器獲得的響應(yīng)中包含了望遠(yuǎn)鏡所有的特征頻率。在環(huán)梁上有4個傳感器，可以測得除剛體位移之外的模態(tài)。支撐環(huán)梁的桁架上也有若干加速度計來獲得其運動規(guī)律。

響應(yīng)測試試驗示意圖如圖4 所示。

一般來說，地平式望遠(yuǎn)鏡叉臂的搖擺模態(tài)是最容易產(chǎn)生的幾階模態(tài)之一，但是GEMINI 望遠(yuǎn)鏡測量時，并沒有在叉臂上布置傳感器，主要考慮的是，四通的剛度非常好，故兩側(cè)的叉臂被四通連接成為一體。

圖4 響應(yīng)測試試驗

在測試開始階段，本研究首先檢測了所有通道是否通暢，以及是否與結(jié)構(gòu)完全連接。具體檢測的方法是檢查增益，如果物理位置上相差不多的兩個傳感器其增益應(yīng)該相差不大。

最后測量的結(jié)果使用累積功率譜的方法來表達(dá)，即利用功率譜密度另外一個定義方法:

由圖可知在超過15 Hz 之后，系統(tǒng)的響應(yīng)下降的尤其的快，故只需要考慮15 Hz 以內(nèi)的模態(tài)。

根據(jù)實測，第一階模態(tài)為方位軸的搖擺為1.83 Hz。

在測量的過程中，還得到一些對于8 m 級口徑系統(tǒng)的經(jīng)驗性結(jié)論。首先，對于相同的方位角，當(dāng)俯仰角變化時，對于結(jié)果的影響并不大;其次，在微風(fēng)情況下，主鏡邊緣的運動大約為1 μm RMS，次鏡的運動大約為幾十微米。

通過實際測量來獲得風(fēng)載數(shù)據(jù)還有另外一項重要的作用就是確定當(dāng)?shù)仫L(fēng)載的相關(guān)尺度。結(jié)構(gòu)函數(shù)是一種之前已經(jīng)存在的數(shù)學(xué)工具，自從上個世紀(jì)九十年代之后，越來越廣泛的用于了大口徑望遠(yuǎn)鏡的各項評價之中，如在加工大口徑光學(xué)元件的過程中，利用結(jié)構(gòu)函數(shù)就可以分辨出在何種尺度還存在著起伏。將此數(shù)學(xué)工具用于風(fēng)載的評價，可以看出，GEMINI 在該站址的特征尺度大約為2 m 左右。

風(fēng)壓結(jié)構(gòu)函數(shù)如圖5 所示。

圖5 風(fēng)壓結(jié)構(gòu)函數(shù)

4 結(jié)束語

對于大口徑的風(fēng)載分析，主要來源于土木工程中的風(fēng)致結(jié)構(gòu)振動，望遠(yuǎn)鏡主要受到的載荷主要可以分為平均風(fēng)與脈動風(fēng)。與土木工程不同的是，對于可以調(diào)節(jié)方位、俯仰的大口徑望遠(yuǎn)鏡，脈動風(fēng)對于系統(tǒng)觀測的影響更大，同時，由于望遠(yuǎn)鏡主鏡的特殊性，風(fēng)載會引起主鏡變形，這也是大口徑望遠(yuǎn)鏡風(fēng)載需要特別研究的原因。

本研究介紹了大口徑望遠(yuǎn)鏡風(fēng)載分析的幾種方法，其中以實測法為重點，闡述了GEMINI 望遠(yuǎn)鏡在裝配過程中進行的風(fēng)載響應(yīng)測試，并介紹了累計功率譜以及結(jié)構(gòu)函數(shù)這兩種結(jié)果分析的數(shù)學(xué)工具。本研究的綜述，對于大口徑望遠(yuǎn)鏡的建設(shè)工作具有一定的指導(dǎo)意義。

［1］BECKERS J M，WILLIAMS J T.“Performance of the Multiple Mirror Telescope (MMT)III.Seeing experiments with the MMT，”Advanced Technology Optical Telescopes［C］.SPIE 332 Tucson，AZ，1982:16-23.

［2］WOOLF N J，ULICH B L.“Gone with the wind，or sailing and seeing with a giant telescope，”Site Testing for Future Large Telescopes［C］.ESO Proc.18，La Silla，Chile，1983:163-183.

［3］GILLINGHAM P R.Seeing Measurements at the Anglo-Australian Telescope［C］//Advanced Technology Optical Telescopes II，SPIE 444.London，England:［s.n.］，1983:165-174.

［4］CAO Xiao-tao，SUN Tian-yu，ZHAO Yun-long et al.Current status and development tendency of image stabilization system of large aperture space telescope［J］.Chinese Optics，2014，7(5):739-748.

［5］BENISTON M，HERTIG J A.Analysis and design of the wind-shield for the ESO Very Large Telescope［C］//Advanced Technology Optical Telescopes Ⅲ，SPIE 628.Tucson，Arizona:［s.n.］，1986:360-368.

［6］SCHNEERMANN M.VLT Enclosures Wind Tunnel Tests and Fluid Dynamic Analyses［C］//Advanced Technology Optical Telescopes V，SPIE 2199，Kona，Hawaii:［s.n.］，1994:465-476.

［7］WANG Fu-guo，YANG Fei，ZHAO Hong-chao，et al.Progress in TMT M3 System［J］.Chinese Optics，2013，10(6):643-651.

［8］ITOH N，MIKAMI I，NOGUCHI T，et al.“Mechanical structure of JNLT—Analysis of mirror deflection due to wind loading，”［C］//Advanced Technology Optical Telescopes IV，SPIE 1236，Tucson，Arizona:［s.n.］，1990:866-877.

［9］ZHAO Hong-chao，ZHANG Jin-xu，YANG Fei，et al.Preloading eight-van spider for supporting structure of secondary mirror［J］.Opt.Precision Engineering，2013，21(5):1199-1204.

［10］David Smith.“Implications of the Nobeyama Wind Testing Results for the GSMT，”NIO-RPT-0005［M］.AURA New Initiatives Office，2002.

［11］WU Xiao-xia，LI Jian-feng，SONG Shu-mei et al.Active support system for 4 m SiC lightweight primary mirror［J］.Optics and Precision Engineering，2014，22(9):2451-2457.

［12］CHO M K，STEPP L，KIM S.Wind buffeting effects on the Gemini 8m primary mirrors［C］//Opto mechanical Design and Engineering 2001.Hathaway，SPIE 4444，San Diego，California:［s.n.］，2001:302-314.

［13］DOUGLAS R N，SEBAG J，WARNER M，et al.Wind Induced Image Degradation(Jitter)of the LSST Telescope［J］.SPIE Ground based and Airborne Telescopes II，2008，7012(1):1-11.

［14］SOBEK R D.Mitigating wind induced telescope jitter［J］.SPIE Acquisition，Tracking and Pointing，2005，5810(1):1-10.

［15］IRARRAZAVAL B，BULERI C，JOHNS M.Wind responses of Giant Magellan telescope［J］.SPIE Modeling System Engineering and Project Management for Astronomy VI，2014，9150(22):1-18.

［16］LIAO Hong-bo，F(xiàn)AN Shi-xun，Hei Mo，et al.Modeling and parameter identification for electro-optical stabilized platform servo systems［J］.Optics and Precision Engineering，2015，23(2):477-484.

［17］KAN F W，EGGERS D W.Wind Vibration Analyses of Giant Magellan Telescope［J］.SPIE Modeling System Engineering and Project Management for Astronomy II，2006，6271(1):1-13.

［18］CHO M K，STEPP L，KIM S.“Wind buffeting effects on the Gemini 8m primary mirrors，”Proc .SPIE Optomechanical Engineering，2001.

［19］CHO M K，STEPP L M，SMITH D R.Wind Loading of Large Telescopes［C］//SPIE symposium on Astro-nomical Telescopes and Instrumentation.Waikoloa，HI:［s.n.］，August，2002:352-367.

［20］DOUGLAS E，ADAMS，MARK A.Rumsey“Theoretical analysis of acceleration measurements in a model of an operating wind turbine”［J］.Active and Passive Smart Structures and Integrated Systems，2010(12):76431F.

［21］FORBES F，GABOR G.“Wind Loading of Large Astronomical Telescopes”［C］//Proc.SPIE Advanced Technology Telescopes，332，1982:198-205.

［22］M Mynowski，D G，VOGIATZIS，K ANGELI，G Z，et al.Wind Loads on Ground-Based Telescopes［J］.Applied Optics，2006，45(30):7912-7923.