尹 釗，侯旭濤

（1. 中國空間技術(shù)研究院 載人航天總體部，北京 100094；2. 航天恒星科技有限公司，天津 300000）

?

載人航天器在軌力學環(huán)境測量系統(tǒng)試驗研究

尹 釗1，侯旭濤2

（1. 中國空間技術(shù)研究院 載人航天總體部，北京 100094；2. 航天恒星科技有限公司，天津 300000）

摘要：為獲取載人航天器上升段與在軌期間的力學環(huán)境，設計了在軌力學環(huán)境測量系統(tǒng)，并在載人航天器振動試驗期間對其性能進行了測試，獲取了航天器典型結(jié)構(gòu)的力學環(huán)境數(shù)據(jù)并進行了頻譜分析。將該系統(tǒng)測量數(shù)據(jù)與試驗中載人航天器平臺測點的測量數(shù)據(jù)進行了比對，結(jié)果驗證了該系統(tǒng)的有效性，為載人航天器載荷條件的制修訂提供了參考。

關鍵詞：載人航天器；在軌力學環(huán)境；測量系統(tǒng)；試驗研究

0 引言

大型載人航天器在上升段將經(jīng)歷全新器箭耦合特性下復雜的載荷條件。若能夠獲取上升段載人航天器艙體結(jié)構(gòu)的力學響應，將對后續(xù)載人航天器載荷條件優(yōu)化、結(jié)構(gòu)長壽命設計、力學仿真和試驗提供有利數(shù)據(jù)[1-4]。

在某大型載人航天器地面振動試驗期間，對所設計的在軌力學環(huán)境測量系統(tǒng)進行測試。本文將力學測量系統(tǒng)所采集得到的數(shù)據(jù)與型號平臺力學試驗數(shù)據(jù)進行比對，以驗證該系統(tǒng)采集數(shù)據(jù)的有效性。

1 在軌力學環(huán)境測量系統(tǒng)

在軌力學環(huán)境測量系統(tǒng)是通過布置在航天器艙壁與艙內(nèi)的傳感器，來監(jiān)測發(fā)射段與在軌段艙體典型結(jié)構(gòu)受到的振動及沖擊響應。通過對采集數(shù)據(jù)的分析處理，可獲得載人航天器艙體結(jié)構(gòu)的力學響應特性，形成密封艙體結(jié)構(gòu)的載荷譜。

系統(tǒng)由傳感器、采編單元與電纜組成（見圖1）。傳感器用于采集振動、沖擊加速度信號。采編單元對采集到的信號進行編碼、存儲和通過以太網(wǎng)傳輸至交換機，并下傳至地面，供地面人員進行數(shù)據(jù)分析。采編單元在發(fā)射前、發(fā)射段以及在軌期間的定時開機、關機通過數(shù)據(jù)注入來控制；該單元由器上電源供電，并利用DC/DC變換器進行電壓變換。電纜用于各單元之間的信號及供電傳輸。

圖1　力學環(huán)境測量系統(tǒng)組成框圖Fig. 1 Configuration of load environment measurement system

2 振動試驗及測試結(jié)果分析

2.1試驗目的

在載人航天器振動試驗期間，對力學環(huán)境測量系統(tǒng)進行了測試，目的有：

1）獲取試驗期間載人航天器艙體主結(jié)構(gòu)的力學響應特性；

2）將該系統(tǒng)的試驗測量數(shù)據(jù)和與載人航天器平臺位置相近的測點數(shù)據(jù)進行比對，驗證二者數(shù)據(jù)的一致性；

3）將分析處理后的測量數(shù)據(jù)與平臺試驗數(shù)據(jù)進行比對，驗證該系統(tǒng)的有效性。

2.2測點布局

力學環(huán)境測量系統(tǒng)共有60個傳感器，180路模擬通道。為驗證該系統(tǒng)的有效性，傳感器在載人航天器主要結(jié)構(gòu)部位的布局位置與平臺的試驗測點布局位置基本一致。該系統(tǒng)在載人航天器平臺力學試驗期間進行加電測試，采集結(jié)構(gòu)的力學響應。表1給出了部分力學環(huán)境測量系統(tǒng)傳感器與平臺測點的布局位置及對應關系。

表1　力學環(huán)境測量系統(tǒng)部分測點的位置及與平臺測點位置的對應關系Table 1 Position of some testing sensors of the load environment measurement system and their corresponding measurement points of the spacecraft

2.3試驗條件

某大型載人航天器振動環(huán)境試驗程序為沿y、z、 x三個方向分別進行正弦振動試驗，同一方向試驗分別進行特征級、驗收級與鑒定級試驗。下面以對鑒定級試驗數(shù)據(jù)進行分析為例。鑒定級試驗條件如下：

2）掃描幅值變化規(guī)律，如表2所示。

表2　正弦振動試驗條件(0～P峰值)Table 2 Sine vibration test condition (0～P peak value)

2.4分析原理和方法

采用兩種方法對力學環(huán)境測量系統(tǒng)所采集的數(shù)據(jù)與載人航天器平臺測點數(shù)據(jù)進行分析：

1）利用相關度分析方法分析系統(tǒng)所采集的數(shù)據(jù)與平臺測點數(shù)據(jù)的相關性；

2）利用時-頻分析方法分析響應信號隨時間的變化關系。

2.4.1相關度分析方法

相關系數(shù)是反映兩個變量之間相關關系密切程度的統(tǒng)計分析指標。

如果有兩個變量X、Y，最終計算出的相關系數(shù)的含義為：

1）當相關系數(shù)為0時，兩變量為非相關；

2）當X值增大（減?。瑒tY值增大（減?。r，兩變量正相關，相關系數(shù)在0.00與1.00之間；

3）當X值增大（減小），則Y值減?。ㄔ龃螅r，兩變量負相關，相關系數(shù)在-1.00與0.00之間。

相關系數(shù)的絕對值越大，則兩變量的相關度越強。通常情況下，可通過相關系數(shù)的取值范圍來判斷兩變量的相關度，如表3所示。

表3　相關系數(shù)和相關度的關系Table 3 Relationship between relative coefficient and degree of correlation

該系統(tǒng)的測量數(shù)據(jù)和平臺測點數(shù)據(jù)是兩組相互獨立的量，兩組變量的觀測值是成對出現(xiàn)的，在這里采用皮爾遜相關系數(shù)計算方法對兩組數(shù)據(jù)的相關度進行計算，

式中：ρX,Y為相關系數(shù)；N為時頻分析圖中所有數(shù)據(jù)點的個數(shù)。

2.4.2時-頻分析法

時-頻分析方法提供了時間域和頻率域的聯(lián)合分布信息，描述信號頻率隨時間變化的關系。

根據(jù)正弦試驗原理，在（fL, fH）頻帶內(nèi)的對數(shù)掃描，掃描頻率f隨時間t變化的關系可表示為[5-7]

式中：R為單位時間掃描頻率倍頻程增量；β為單位時間對數(shù)頻率變化率，β=Rln2 。

通過該方法計算得到信號頻率隨時間變化關系，從而分析信號在不同頻率下的能量分布隨時間的變化情況。

2.5數(shù)據(jù)處理方法

對正弦掃描振動試驗數(shù)據(jù)處理時，首先應得到時域信號的頻率信息，即頻率跟蹤辨識。本文采用短時傅里葉變換（STFT）進行頻率跟蹤辨識，即對輸出信號進行采樣，采集整周期信號，然后用STFT算法計算每個通道響應信號所處周期內(nèi)的各響應點的頻率。具體計算方法如下：

1）加窗

加窗的目的是為了降低自時域窗有限信號非周期、非整周期等原因引起的頻譜泄漏問題。STFT方法也即窗口傅里葉變換，是將信號分成許多很小的時間間隔，用傅里葉變換分析每一段時間間隔，以便確定該時間間隔存在的頻率，其表達式為

式中：g(t)是有緊支的函數(shù)；f(t)是被分析的時域信號。在變換中e-jωt起著頻限作用，g(t)起著時限作用。隨著時間τ的變化，g(t)所確定的“時間窗”在t軸上移動，使f(t)“逐漸”被分析。

2）傅里葉變換

對于離散時間序列信號，X(t)=X(n?t), (n= 0,1,2,…,N-1)，其離散傅里葉變換為

式中：k為傅里葉變換中所取的實部個數(shù)，k=0,1,2,…,N/2-1；為子樣均值為0的第m幀離散時間序列；Xk,m,R為Xk,m的實部分量；Xk,m,I為Xk,m的虛部分量。

本文中N取211=2048。

2.6試驗結(jié)果及分析

2.6.1試驗結(jié)果

按照表1所示取4組測量數(shù)據(jù)與平臺數(shù)據(jù)進行比對。對數(shù)據(jù)分析和處理后，獲得4組數(shù)據(jù)的頻域曲線以及時-頻分析結(jié)果，如圖2～圖5所示。

圖2　測點D1和P1的x向時頻分析Fig. 2 Time-frequency analysis at measurement points D1 and P1 in x direction

圖3　測點D2和P2的y向時-頻分析Fig. 3 Time-frequency analysis at measurement points D2 and P2 in y direction

圖4　測點D3和P3的y向時-頻分析Fig. 4 Time-frequency analysis at measurement points D3 and P3 in y direction

(a) D4與P4的z向頻域曲線

圖5　測點D4和P4的z向時-頻分析Fig. 5 Time-frequency analysis at measurement points D4 and P4 in z direction

2.6.2結(jié)果分析

1）數(shù)據(jù)相關度分析

4組數(shù)據(jù)的相關度分析結(jié)果如表4所示。由表4可見，鑒定級試驗中，力學環(huán)境測量系統(tǒng)所采集的數(shù)據(jù)與載人航天器平臺數(shù)據(jù)的相關系數(shù)均在0.98以上，相關度均為極強。這說明力學環(huán)境測量系統(tǒng)測得的數(shù)據(jù)是可信而有效的。

表4　頻域曲線相關性分析結(jié)果Table 4 Correlations of frequency responses

2）時-頻分析

以測點D1和P1的數(shù)據(jù)為例，其時-頻分析結(jié)果見圖2(b)和圖2(c)，其中：橫坐標是時間，縱坐標是頻率，反映了信號的能量分布隨時間的變化情況；時-頻分析曲線的顏色反映出對應某一時刻的響應在某一頻率附近的積聚程度，且與顏色深淺成比例。振動試驗過程中，激勵頻率由接近5Hz掃描至100Hz，因此結(jié)構(gòu)對應某一具體時刻的響應能量便聚集在該時刻所對應的激勵頻率附近，這一關系由時-頻分析曲線中最下端的紅線體現(xiàn)。圖2中其他位置上出現(xiàn)了暗紅色線條，表明了測量系統(tǒng)中的不同頻率諧振信號的影響（測量信號未經(jīng)過濾波）；與結(jié)構(gòu)的實際響應信號相比，諧振信號的能量低很多。

從時-頻分析結(jié)果可看出：最下面的一條紅線的變化趨勢符合式(2)的指數(shù)變化關系；在接近第140s時，D1的響應信號的頻率達到100Hz，而P1的響應信號的頻率同樣在接近140s時達到100Hz。這說明二者的數(shù)據(jù)一致性較好。

3 結(jié)論

本文對力學環(huán)境測量設備采集到的相關數(shù)據(jù)進行了處理，分別利用相關度分析法、時-頻分析法對結(jié)果進行了分析，并與平臺相應測點的數(shù)據(jù)進行比對，驗證了力學測量系統(tǒng)采集數(shù)據(jù)的有效性。試驗得出了以下結(jié)論：

1）力學傳感器測量系統(tǒng)測量數(shù)據(jù)與試驗中平臺測點測得數(shù)據(jù)的相關度分析結(jié)果表明，該系統(tǒng)可有效采集載人航天器力學試驗的振動響應數(shù)據(jù)；

2）從力學傳感器測量系統(tǒng)與平臺測點所采集信號的時-頻分析圖可以看出，二者的時-頻分析圖中均含有諧振信號，這說明振動過程中存在結(jié)構(gòu)的非線性因素。

本文的研究成果為該系統(tǒng)在載人航天器的在軌應用奠定了基礎，可用于載人航天器在軌力學環(huán)境的測量，監(jiān)測航天器發(fā)射、在軌運行及返回期間的動力學響應特性，為后續(xù)航天器結(jié)構(gòu)設計與優(yōu)化、載荷條件的修訂以及結(jié)構(gòu)的損傷容限與壽命分析提供參考。

參考文獻（References）

[1] MIL-STD-810G Environment engineering considerations and laboratory tests[S], 2008

[2] 謝久林, 楊松, 張俊剛, 等. 航天器聲振動力學環(huán)境響應分析[J]. 航天器環(huán)境工程, 2006, 23(2): 83-89 Xie Jiulin, Yang Song, Zhang Jungang, et al. The response prediction of the spacecraft under acoustic vibration environment[J]. Spacecraft Environment Engineering, 2006, 23(2): 83-89

[3] 張正平. 航天運載器力學環(huán)境工程技術(shù)發(fā)展回顧及展望[J]. 航天器環(huán)境工程, 2008, 25(3): 233-236 Zhang Zhengping. A review of studies on dynamic environments of launch vehicles[J]. Spacecraft Environment Engineering, 2008, 25(3): 233-236

[4] 趙煜, 周東強. “海洋二號”衛(wèi)星主動段、自由飛行段力學環(huán)境測量與分析[J]. 航天器環(huán)境工程, 2012, 29(4): 458-463 Zhao Yu, Zhou Dongqiang. In-situ measurement and analysis of mechanical environment on HY-2 satellite[J]. Spacecraft Environment Engineering, 2012, 29(4): 458-463

[5] 沈鳳霞. 正弦振動試驗及結(jié)構(gòu)響應數(shù)據(jù)處理方法[J].航天器環(huán)境工程, 2004, 21(4): 7-13 Shen Fengxia. Sine vibration test techniques and data processing method of structure response[J]. Spacecraft Environment Engineering, 2004, 21(4): 7-13

[6] 金恂叔. 航天器力學環(huán)境試驗技術(shù)的發(fā)展概況和趨勢[J]. 航天器環(huán)境工程, 2003, 20(6): 15-21 Jin Xunshu. The development status and trends of spacecraft dynamic environment testing[J]. Spacecraft Environment Engineering, 2003, 20(6): 15-21

[7] 閆凱, 李艷琴, 李丹丹. 振動試驗及其相關參數(shù)的分析[J]. 環(huán)境技術(shù), 2014(4): 23-24 Yan Kai, Li Yanqin, Li Dandan. Vibration test and its relevant parameters analysis[J]. Environmental Testing, 2014(4): 23-24

（編輯：許京媛）

Experimental study of orbital environment measurement system for manned spacecraft

Yin Zhao1, Hou Xutao2
(1. Institute of Manned Space System Engineering, China Academy of Space Technology, Beijing 100094, China; 2. Space Star Technology Co. Ltd., Tianjin 300000, China)

Abstract:In order to acquire the mechanical environment data in launching and on orbit section of a manned spacecraft, an orbital environment measurement system for the spacecraft is designed and tested in the vibration experiment, and the mechanical environmental data of the representative structure of the manned spacecraft is acquired and analyzed in the frequency domain. A comparison is made between the data of the measurement system and the data of the manned spacecraft in the vibration experiment. The result verifies the availability of the system. The research provides a reference for the formulation and revision of the orbit condition of the manned spacecraft.

Key words:manned spacecraft; in-orbit mechanical environment; measurement system; experimental study

作者簡介：尹 釗（1986—），男，碩士學位，主要從事載人航天器總體設計工作。E-mail: yinzhao08@126.com。

基金項目：國家重大科技專項工程

收稿日期：2015-08-17；修回日期：2016-01-20

DOI:10.3969/j.issn.1673-1379.2016.01.017

中圖分類號：V416.8

文獻標志碼：A

文章編號：1673-1379(2016)01-0095-05