靜電懸浮加速度計(jì)的地面重力傾角標(biāo)定方法

薛大同

(中國空間技術(shù)研究院蘭州物理研究所，蘭州730000)

0 引言

靜電懸浮加速度計(jì)是重力場(chǎng)測(cè)量衛(wèi)星上的主要載荷之一，用于測(cè)量衛(wèi)星受到的非重力，主要有:軌道高度處由于熱層大氣密度、水平中性風(fēng)、太陽和地球的輻射壓力引起的拖曳［1］，以及衛(wèi)星自身產(chǎn)生的干擾力。為了準(zhǔn)確恢復(fù)地球重力場(chǎng)，衛(wèi)星在軌飛行時(shí)，需要以高的精度水平測(cè)定加速度計(jì)的標(biāo)定參數(shù)，這需要靠在每天的幾個(gè)弧段中處理GPS和加速度計(jì)數(shù)據(jù)，特別是使用可論證為最準(zhǔn)確的重力場(chǎng)模型EIGEN－GL04C，和海洋潮汐模型 FES2004來實(shí)現(xiàn)［2］。然而，這絕不意味著加速度計(jì)飛行件在正式交付之前不需要作性能測(cè)試。為此，Marque等人指出，加速度計(jì)檢驗(yàn)質(zhì)量需要在地面重力下懸浮，以便為功能試驗(yàn)提供廣泛的可能性。例如，加速度計(jì)傳感頭飛行件在環(huán)境試驗(yàn)前后需要開展懸浮狀況、傳遞函數(shù)、標(biāo)稱靜電剛度、線性、標(biāo)度因數(shù)等綜合試驗(yàn)，以檢查完整的加速度計(jì)傳感頭狀態(tài)和確認(rèn)環(huán)境試驗(yàn)被通過，可以驗(yàn)收。這些試驗(yàn)是在專用的擺臺(tái)上完成的，該擺臺(tái)可以把懸吊平臺(tái)精密控制在水平狀態(tài)和實(shí)施傾斜擺動(dòng)，并限制由地面重力投影引起的臺(tái)面方向加速度在科學(xué)模式范圍內(nèi)［3］。法國ONERA早在上世紀(jì)90年代初就研制出雙級(jí)擺臺(tái)并用于靜電懸浮加速度計(jì)噪聲檢測(cè)［4］，在世紀(jì)之交又采用主動(dòng)控制技術(shù)抑制擺臺(tái)低于0.1 Hz的干擾，這靠擺臺(tái)上同時(shí)放置的另一臺(tái)靜電懸浮加速度計(jì)提供反饋控制信號(hào)來實(shí)現(xiàn)［5］。進(jìn)入本世紀(jì)后推崇用擺臺(tái)上成對(duì)靜電懸浮加速度計(jì)噪聲頻譜差模來抵消共模環(huán)境噪聲［6］，但存在原理上的錯(cuò)誤［7］。Marque 等人近年來主要用擺臺(tái)來做差分標(biāo)度因數(shù)和二次因數(shù)驗(yàn)證［3］。本文在此基礎(chǔ)上以SuperSTAR加速度計(jì)的性能參數(shù)為例，分析了在擺臺(tái)上用傳統(tǒng)的靜態(tài)標(biāo)定法確定模型方程各參數(shù)的困難，指出在擺臺(tái)上無論是采用動(dòng)態(tài)標(biāo)定法大致判斷標(biāo)度因數(shù)和對(duì)成對(duì)加速度計(jì)進(jìn)行模型方程各參數(shù)的一致性檢驗(yàn)，還是采用電模擬法獲取二階非線性系數(shù)，都必須采取技術(shù)措施使得三階非線性系數(shù)可以忽略。對(duì)動(dòng)態(tài)標(biāo)定法，提出了防止頻譜泄漏和幅度譜中壓低噪聲干擾幅度的措施。對(duì)電模擬法，用實(shí)例給出了具體實(shí)施方案和效果。以上分析結(jié)果對(duì)GRADIO等更高指標(biāo)的加速度計(jì)具有參考意義。

1 模型方程

根據(jù)等效原理，慣性加速度與引力加速度等效，因此地面上可以利用重力加速度在加速度計(jì)輸入軸方向上的分量作為加速度測(cè)量標(biāo)準(zhǔn)，復(fù)現(xiàn)所需要的加速度值，實(shí)施地面重力傾角標(biāo)定。此時(shí)靜電懸浮加速度計(jì)的敏感軸如圖1所示［8］，x軸處于負(fù)的鉛垂方向。

理想情況下靜電懸浮加速度計(jì)每個(gè)軸的反饋控制電壓與該軸受到的加速度成正比［3］。以z軸為例，可表示為:

圖1 靜電懸浮加速度計(jì)的敏感軸［8］Fig.1 Sensitive axes of the electrostatically suspended accelerometer［8］

式中，Vz為z軸的反饋控制電壓(V);Gp，z為靜電懸浮加速度計(jì)反饋控制的物理增益(m·s－2/V)，并有［9］:

式中，ε為真空介電常數(shù)，ε =8.8542×10－12C2/(N·m2);Az為z軸的電極面積;m為檢驗(yàn)質(zhì)量塊的質(zhì)量;dz為檢驗(yàn)質(zhì)量塊與z軸電極間的平均間隙;Vp為檢驗(yàn)質(zhì)量上施加的固定偏壓。對(duì)于SuperSTAR加速度計(jì)而言，Az=2.08 ×10－4m2，m=7.2 ×10－2kg，dz=1.75 × 10－4m［10］，Vp=10V［11］，因此 Gp，z=1.67 ×10－5m·s－2/V。

實(shí)際上，靜電懸浮加速度計(jì)與其他類型的加速度計(jì)類似，存在偏值、標(biāo)度偏差、非線性、交叉耦合，因此，式(1)表達(dá)的γel，z并非z軸受到的實(shí)際加速度，而只是靜電懸浮加速度計(jì)指示的加速度值(以下簡稱加速度計(jì)示值，單位m/s2)。

傳統(tǒng)加速度計(jì)在地面重力下測(cè)試時(shí)的模型方程為［12］:

式中，K0，z為 z軸的偏值(V)，K1，z為 z軸的標(biāo)度因數(shù)(V·m－1s2)，az為z軸的實(shí)際加速度 (m/s2)，K2，z為 z軸的二階非線性系數(shù) (V·m－2s4)，K3，z為z軸的三階非線性系數(shù)(V·m－3s6)，K4，z為z軸和x軸的交叉耦合系數(shù)(V·m－2s4)，ax為x軸的實(shí)際加速度(m/s2)，K5，z為z軸和y軸的交叉耦合系數(shù)(V·m－2s4)，ay為y軸的實(shí)際加速度(m/s2)。

交叉耦合效應(yīng)是加速度計(jì)某個(gè)敏感軸的偏值因該軸向及與之正交方向同時(shí)受到加速度作用而發(fā)生輕微改變的現(xiàn)象，變化量與這兩個(gè)加速度的乘積成正比，比例系數(shù)稱為交叉耦合系數(shù)。但是對(duì)于靜電懸浮加速度計(jì)來說，情況有所不同。兩個(gè)敏感軸之間的耦合主要來自于質(zhì)量塊的不平行度和不垂直度。為此，要求不平行度和不垂直度小于1×10－5rad(2″)，以保證在三個(gè)軸之間高解耦［13－14］。在此前提下，地面測(cè)試時(shí)鉛垂方向殘余耦合作用的恒定部分表現(xiàn)為偏值［15］，起伏不定部分表現(xiàn)為噪聲［5］。另一水平方向的殘余耦合作用可以忽略。因此，式(3)右端的最后兩項(xiàng)不存在:

對(duì)照式(1)和式(4)，可知1/Gp，z為理論標(biāo)度因數(shù)。SuperSTAR加速度計(jì)規(guī)定實(shí)際標(biāo)度因數(shù)的偏差不超過2%［15］。即:

式中，c1，z為z軸實(shí)際標(biāo)度因數(shù)與理論標(biāo)度因數(shù)的比值，k1，z為 z軸的標(biāo)度偏差，|k1，z| ＜ 0.02。

令:

習(xí)慣上仍分別稱 c0，z為 z軸的偏值 (m/s2)，c1，z為z軸的標(biāo)度因數(shù)(無量綱)，c2，z為z軸的二階非線性系數(shù) (m－1s2)，c3，z為 z軸的三階非線性系數(shù)(m－2s4)。

將式(1)、式(5)和式(6)代入式(4)，便可得到:

將式(7)簡化表達(dá)為［15］:

需要注意的是，如前所述，地面重力下測(cè)試時(shí)得到的c0主要是鉛垂方向殘余耦合引起的，不同于在軌飛行時(shí)的偏值。

我們有:

式中，g為當(dāng)?shù)氐闹亓铀俣戎?m/s2)，θ為擺臺(tái)的重力傾角(rad)，θ0為敏感軸與安裝面間的夾角，即失準(zhǔn)角(rad)。

實(shí)施地面重力傾角標(biāo)定時(shí)，加速度計(jì)的輸出必須處于量程范圍內(nèi)。以SuperSTAR加速度計(jì)為例，由于量程僅為 5μg［16］，所以 θ+ θ0不超過 5 ×10－6rad(1″)。在此極限條件下:

將式(10)代入式(8)，得到:

式(11)即為靜電懸浮加速度計(jì)地面重力傾角標(biāo)定的模型方程。

值得注意的是，加速度計(jì)敏感結(jié)構(gòu)安裝在經(jīng)過研磨等精細(xì)加工的殷鋼基座上，該基座的垂直度和平行度公差為5 ×10－5rad(相當(dāng)于 10″)［4］，因此 θ0可以按不超過5×10－5rad估計(jì)。也就是說，為了使加速度計(jì)的輸出處于量程范圍內(nèi)，θ必須部分抵償θ0，以保證 θ+ θ0不超過5 ×10－6rad。

SuperSTAR加速度計(jì)射前檢驗(yàn)的極限指標(biāo)為:|c0|=1 ×10－5m/s2，c1=1.02，|c2|=10m－1s2，|c3|=3 × 105m－2s4［15］。

2 靜態(tài)標(biāo)定法

模型方程各參數(shù)傳統(tǒng)采用靜態(tài)標(biāo)定法確定。對(duì)于靜電懸浮加速度計(jì)來說，就應(yīng)固定在專用擺臺(tái)上，實(shí)施靜態(tài)定標(biāo)傾斜，先后調(diào)出若干個(gè)穩(wěn)定不變、且精確知道其數(shù)值的θ，分別測(cè)出加速度計(jì)示值，并據(jù)此確定加速度計(jì)模型方程各參數(shù)。如第1節(jié)所述，θ0可以按不超過5×10－5rad估計(jì);θ必須部分抵償θ0，以保證 θ+ θ0不超過5 ×10－6rad。

由式(11)得知，未知數(shù)共有 5個(gè):c0，c1，c2，c3，θ0，因此必須至少有5 個(gè) θ值，對(duì)應(yīng)5 個(gè) γel值，形成5個(gè)方程。它們對(duì)于c0，c1，c2，c3來說，是線性方程，可以由其中4個(gè)方程用克萊姆法則聯(lián)立求解，分別得到c0，c1，c2，c3對(duì)θ0的函數(shù)關(guān)系。原則上似乎將它們?cè)俅氲?個(gè)方程，就可以解出θ0，再分別代回 c0，c1，c2，c3對(duì) θ0的函數(shù)關(guān)系，即可解出 c0，c1，c2，c3。但是對(duì)于θ0來說，是三次方程，有三個(gè)根，且運(yùn)算極為復(fù)雜，從工程觀點(diǎn)來說，精度也無法保證。所以迄今未見靜電懸浮加速度計(jì)采用靜態(tài)標(biāo)定法的報(bào)導(dǎo)。

3 動(dòng)態(tài)標(biāo)定法

動(dòng)態(tài)標(biāo)定法的基本構(gòu)想是:將靜電懸浮加速度計(jì)固定在專用擺臺(tái)上，擺臺(tái)的傾角在加速度計(jì)通帶范圍內(nèi)實(shí)施正弦擺動(dòng)，對(duì)加速度計(jì)示值γel，z隨時(shí)間變化曲線作頻譜分析，根據(jù)幅度譜在各階頻率處的幅值確定加速度計(jì)模型方程各參數(shù)。

為了防止頻譜泄漏，即防止在各次諧波峰的兩旁出現(xiàn)旁瓣峰，影響主峰峰高，要保證整周期采樣［7］。考慮到:① SuperSTAR的采樣率為 10 Sps(Sampling per second)［15］，② 快速傅里葉變換的樣本量必須為2的整數(shù)次冪［7］，③ 基頻應(yīng)在0.1 Hz附近［16］等三個(gè)因素，可選擇擺動(dòng)周期為6.4 s，即基頻為0.15625 Hz，而且應(yīng)足夠穩(wěn)定。

此時(shí)，擺臺(tái)的重力傾角為［6］:

式中，θ1為擺臺(tái)正弦擺動(dòng)的最大角位移(rad)，不超過5μrad(1″);ω為擺動(dòng)的角頻率(rad/s);t為時(shí)間(s)。

將式(12)代入式(11)，并利用三角函數(shù)公式:

可以得到:

根據(jù)第1節(jié)所述極限值，我們?nèi)?c1=1，c2=是，從式(15)得到:

(3)在幅度譜的二倍頻處，c2，c3θ0的貢獻(xiàn)分別為1.2 ×10－8m/s2，5.3 ×10－7m/s2，也就是說，c2的貢獻(xiàn)可能會(huì)被c3的影響所掩蓋。

(4)在幅度譜的三倍頻處，c3的貢獻(xiàn)為8.8×10－9m/s2。

實(shí)際上，γel還受到環(huán)境振動(dòng)噪聲的干擾，包括地震噪聲和人工振源產(chǎn)生的振動(dòng)［14］。圖2為法國ONERA實(shí)驗(yàn)室花岡巖臺(tái)子上由靜電懸浮加速度測(cè)到的噪聲頻譜，可以看到在基頻至三倍頻間(0.15Hz至0.47Hz)干擾大約在(3.5×10－7～1 ×10－6)m·s－2/Hz1/2之間［14］。由于實(shí)施正弦擺動(dòng)時(shí)，擺臺(tái)無法對(duì)噪聲干擾實(shí)施主動(dòng)控制;且擺臺(tái)對(duì)此頻段無被動(dòng)隔振作用，所以擺臺(tái)上的噪聲干擾亦大致如此。如取采樣率10 Sps，樣本量213個(gè)［3］，依據(jù)頻譜分析原理［7］，在幅度譜上基頻至三倍頻間噪聲干擾的幅度將達(dá)到(1.73～4.94)×10－8m/s2。因此，二倍頻處 c2的貢獻(xiàn)、三倍頻處c3的貢獻(xiàn)將被噪聲淹沒。但如果樣本量增加到218個(gè)，在幅度譜上基頻至三倍頻間噪聲干擾的幅度將降低到(3.06 ～8.73) ×10－9m/s2，二倍頻處 c2的貢獻(xiàn)、三倍頻處c3的貢獻(xiàn)就可超出噪聲本底，不被淹沒。數(shù)據(jù)量越大越清晰。

從以上分析可以看到，由于專用擺臺(tái)的振動(dòng)噪聲干擾，動(dòng)態(tài)標(biāo)定法很難精確測(cè)定c3。由于無法分離出θ0，受到θ0和c3的影響，加速度計(jì)模型方程的其他各參數(shù)均無法確定。

如果c3可以忽略，則情況要好得多。此時(shí)，零頻處的幅值主要是c0和θ0的貢獻(xiàn)，從基頻處的幅值可以大致判斷c1的大小，從二倍頻處的幅值可以大致判斷c2的大小。也就是說，主要問題只剩下無法分離 c0和 θ0。

4 檢驗(yàn)二階非線性系數(shù)的方法

二階非線性系數(shù)起因于加速度計(jì)籠的最終缺陷，例如在相對(duì)電極表面或在電子元件真實(shí)數(shù)據(jù)中的差異。結(jié)果是，施加的靜電力中心與籠中心不完全相配［3］。

圖2 在地面上測(cè)到的加速度頻譜［14］Fig.2 Measured acceleration spectra on the ground［14］

第3節(jié)指出，采用動(dòng)態(tài)標(biāo)定法時(shí)，如果c3可以忽略，從二倍頻處的幅值可以大致判斷c2的大小。然而，在幅度譜上二倍頻的幅度遠(yuǎn)低于基頻的幅度，這種情況下擺臺(tái)擺動(dòng)的二次諧波不可以忽略［3］，此時(shí)，即使c3可以忽略，從二倍頻處的幅值也無法判斷c2的大小。理由如下:

當(dāng)擺臺(tái)擺動(dòng)的二次諧波不可以忽略時(shí)，式(12)應(yīng)改寫為:

式中，β ? θ1，z。

忽略c3時(shí)，式(11)簡化為:

將式(16)代入式(17)，并利用式(13)及三角函數(shù)公式:

可以得到:

式(19)表明，加速度計(jì)示值幅度譜在二倍頻處的幅度不僅與c2有關(guān)，還與擺臺(tái)擺動(dòng)的二次諧波幅度β有關(guān)。所以，擺臺(tái)擺動(dòng)的二次諧波不可以忽略時(shí)，不能用動(dòng)態(tài)標(biāo)定法得到c2。

為了檢驗(yàn)c2，可采用電模擬的方法:擺臺(tái)保持水平靜止?fàn)顟B(tài)，在加速度計(jì)PID的輸出處發(fā)送一個(gè)特定的二階非線性系數(shù)標(biāo)定信號(hào)［3］。這樣做，既可以顯著增加信號(hào)幅度，又可以保證信號(hào)不含二次諧波，還可以突出c2的貢獻(xiàn)。

特定的二階非線性系數(shù)標(biāo)定信號(hào)為ASK(正弦載波數(shù)字幅度調(diào)制)信號(hào)，調(diào)制信號(hào)1，0等間隔排列。當(dāng)c3可以忽略時(shí)，加速度計(jì)示值幅度譜在調(diào)制信號(hào)重復(fù)頻率處的尖峰包含有c1和c2的作用。為了突出c2的作用，載波頻率和幅度要盡量高，調(diào)制信號(hào)重復(fù)頻率要盡量低。由于載波幅度高，加速度計(jì)噪聲的影響可以忽略。

載波頻率可超出加速度計(jì)回路帶寬。載波幅度應(yīng)保證帶外抑制后的信號(hào)不超出加速度計(jì)測(cè)量范圍。調(diào)制信號(hào)重復(fù)頻率應(yīng)處于加速度計(jì)測(cè)量帶寬內(nèi)。為了防止頻譜泄漏，調(diào)制信號(hào)重復(fù)頻率、采樣率、載波頻率間應(yīng)符合2的整數(shù)次冪關(guān)系。

SuperSTAR加速度計(jì)閉環(huán)傳遞函數(shù)的幅頻特性曲線如圖 3 所示［10］。

以調(diào)制信號(hào)重復(fù)周期 25.6 s，即重復(fù)頻率0.0390625Hz，載波頻率320Hz為例，從圖3可以看到，載波信號(hào)幅度將衰減至不足10%，考慮到SuperSTAR 加速度計(jì)測(cè)量范圍 5 ×10－5m/s2［16］，載波幅度可以取 5×10－4m/s2，由第 1節(jié)所述 Super－STAR加速度計(jì)的物理增益可知相當(dāng)于30 V，以此作為加速度的電模擬信號(hào)，若以ae表示，則加速度計(jì)感受到的加速度為a=ae+gθ0，以此代入式(8)，得到:

圖3 SuperSTAR閉環(huán)幅頻特性曲線［10］Fig.3 Closed loop magnitude－frequency characteristic curve of the SuperSTAR［10］

假設(shè)加速度計(jì)無通帶限制，采樣率5120 Sps，樣本量220個(gè)，采用第3節(jié)的模型方程各參數(shù)，得到式(20)中的 c1(ae+gθ0) 項(xiàng)、c2(ae+gθ0)2項(xiàng)、c3(ae+gθ0)3項(xiàng)幅度譜分別如圖4、圖5、圖6所示。可以看到，調(diào)制信號(hào)重復(fù)頻率處的尖峰幅度在圖4中最低，僅3.8 ×10－8m/s2;在圖5 中其次，為8.0 ×10－7m/s2;在圖6中最高，達(dá)到 3.5×10－5m/s2。也就是說，c2對(duì)調(diào)制信號(hào)重復(fù)頻率處尖峰的貢獻(xiàn)可能會(huì)被θ0和c3的影響所掩蓋。而c1的影響不嚴(yán)重。所以，為了突出c2的貢獻(xiàn)，必須設(shè)法使c3可以忽略。此時(shí)，式(20)變?yōu)?

為了了解θ0對(duì)調(diào)制信號(hào)重復(fù)頻率處尖峰幅度的影響，我們改設(shè)θ0=0 rad，而其他條件不變。此時(shí)，式(20)變?yōu)?

圖7、圖8、圖9分別為式(22)中的c1ae項(xiàng)項(xiàng)項(xiàng)幅度譜。可以看到，調(diào)制信號(hào)重復(fù)頻率處的尖峰幅度在圖7中為3.8×10－8m/s2，與圖4相同;在圖8 中為8.0 ×10－7m/s2，與圖5 相同;在圖9中為1.9 ×10－9m/s2，僅為圖 6 的 0.0054%。這說明θ0對(duì)調(diào)制信號(hào)重復(fù)頻率處尖峰幅度的影響，僅限于c3的貢獻(xiàn)。也就是說，當(dāng)c3可以忽略時(shí)，不需要知道θ0的具體數(shù)值，就可以根據(jù)調(diào)制信號(hào)重復(fù)頻率處的尖峰，確定c2的實(shí)際值。

圖4 在指定參數(shù)下的c1(ae+gθ0)項(xiàng)幅度譜Fig.4 Amplitude spectrum of the c1(ae+gθ0)item with the appointed parameters

圖5 在指定參數(shù)下的c2(ae+gθ0)2項(xiàng)幅度譜Fig.5 Amplitude spectrum of the c2(ae+gθ0)2 item with the appointed parameters

式(21)所示加速度計(jì)示值幅度譜如圖10所示。可以看到，調(diào)制信號(hào)重復(fù)頻率處的尖峰幅度為8.0 ×10－7m/s2，與圖 5、圖8 相同，證明它確實(shí)主要取決于二階非線性系數(shù)。由于峰高與c2值成正比，由實(shí)際峰高便可確定c2的實(shí)際值。需要說明的是，實(shí)際檢測(cè)加速度計(jì)示值時(shí)，由于受到測(cè)量通帶的限制，得到的幅度譜在測(cè)量通帶以外會(huì)受到抑制。至于實(shí)際采樣率較低，樣本量較少，只要有與之相應(yīng)低通，從而符合采樣定理，且樣本量保證清晰識(shí)別調(diào)制信號(hào)重復(fù)頻率，就不會(huì)影響判讀。但如果實(shí)際采用的ASK信號(hào)與此不同，則應(yīng)根據(jù)實(shí)際采用的ASK信號(hào)重新作頻譜分析。

圖6 在指定參數(shù)下的c3(ae+gθ0)3項(xiàng)幅度譜Fig.6 Amplitude spectrum of the c3(ae+gθ0)3 item with the appointed parameters

圖7 在指定參數(shù)下的c1ae項(xiàng)幅度譜Fig.7 Amplitude spectrum of the c1ae item with the appointed parameters

圖8 在指定參數(shù)下的項(xiàng)幅度譜Fig.8 Amplitude spectrum of theitem with the appointed parameters

圖9 在指定參數(shù)下的項(xiàng)幅度譜Fig.9 Amplitude spectrum of theitem with the appointed parameters

5 模型方程各參數(shù)的一致性檢驗(yàn)方法

動(dòng)態(tài)標(biāo)定法可以用來直接檢驗(yàn)兩臺(tái)加速度計(jì)模型方程各參數(shù)的一致性［3］。GRACE和GOCE均需要成對(duì)使用靜電懸浮加速度計(jì)，特別關(guān)注每對(duì)加速度計(jì)的一致性，直接檢驗(yàn)?zāi)Ｐ头匠谈鲄?shù)的一致性是非常必要的。

成對(duì)加速度計(jì)同時(shí)固定在一個(gè)專用擺臺(tái)上，擺臺(tái)按照第3節(jié)的要求實(shí)施正弦擺動(dòng)。將式(5)代入式(15)，得到兩臺(tái)加速度計(jì)示值幅度譜在基頻處的幅度分別為:

式中，對(duì)于編號(hào) i和編號(hào) j的加速度計(jì)，A1，i，A1，j分別為各自示值在基頻處的幅度，k1，i，k1，j分別為各自的標(biāo)度偏差，c2，i，c2，j分別為各自的二階非線性系數(shù)，θ0，i，θ0，j分別為各自的失準(zhǔn)角，c3，i，c3，j分別為各自的三階非線性系數(shù)。

稱兩臺(tái)加速度計(jì)示值的平均值為共模，半差分值為差模，則差模和共模幅度譜在基頻處幅度之比為:

圖10 在指定參數(shù)下的加速度計(jì)示值幅度譜Fig.10 Amplitude spectrum of the accelerometers indication with the appointed parameters

式中，A1，D為基頻處的差模幅度，A1，C為基頻處的共模幅度。將式(23)代入式(24)，并考慮到第1節(jié)所有參數(shù)的極限值，得到:

在第1節(jié)給出的參數(shù)極限值范圍內(nèi)，|k1|＜

(1)在第1節(jié)給出的參數(shù)極限值范圍內(nèi)，由式(25)得到的比值有可能主要反映兩臺(tái)加速度計(jì)c3，θ0的不一致性。只有當(dāng)c3可以忽略時(shí)，c1的不一致性才占主要地位。

(2)式(25)得到的比值越小，兩臺(tái)加速度計(jì)的c3，θ0，進(jìn)而 c1，c2越相近。

6 結(jié)論

(1)靜電懸浮加速度計(jì)檢驗(yàn)質(zhì)量需要在地面重力下懸浮，以便為功能試驗(yàn)提供廣泛的可能性。

(2)由于敏感軸之間的耦合機(jī)理不同于傳統(tǒng)加速度計(jì)，所用的模型方程亦有所差別。

(3)地面重力下測(cè)試時(shí)，偏值主要是鉛垂方向殘余耦合引起的，不同于在軌飛行時(shí)的偏值。

(4)由于失準(zhǔn)角遠(yuǎn)大于量程范圍內(nèi)的重力傾角，無法采用傳統(tǒng)的靜態(tài)標(biāo)定法確定模型方程各參數(shù)。

(5)必須采取技術(shù)措施使得三階非線性系數(shù)可以忽略，才能在專用擺臺(tái)上用動(dòng)態(tài)標(biāo)定法大致判斷標(biāo)度因數(shù)和檢驗(yàn)成對(duì)加速度計(jì)檢驗(yàn)?zāi)Ｐ头匠谈鲄?shù)的一致性，用電模擬法得到二階非線性系數(shù)。

(6)動(dòng)態(tài)標(biāo)定作幅度譜分析時(shí)為防止頻譜泄漏，即防止在各次諧波峰的兩旁出現(xiàn)旁瓣峰，影響主峰峰高，要保證整周期采樣。

(7)動(dòng)態(tài)標(biāo)定作幅度譜分析時(shí)可采用顯著加大樣本量的辦法降低噪聲干擾幅度。

(8)電模擬法只需擺臺(tái)保持水平靜止?fàn)顟B(tài)，并在加速度計(jì)PID的輸出處發(fā)送一個(gè)特定的ASK信號(hào)，就可以根據(jù)調(diào)制信號(hào)重復(fù)頻率處的尖峰，確定c2的實(shí)際值，而不需要知道失準(zhǔn)角θ0的具體數(shù)值。

［1］ Flury J，Bettadpur S，Tapley B D.Precise accelerometry onboard the GRACE gravity field satellite mission［J］.Advances in Space Research，2008，42:1414－1423.

［2］ Bruinsma S，Biancale R，Perosanz F.Calibration parameters of the CHAMPand GRACE accelerometers［OL］.http://www.dgfi.badw.de/typo3_mt/fileadmin/20071015 － 17 － Potsdam/di_1800_05_bruinsma.pdf

［3］ Marque JP，Christophe B，Liorzou F，et al.The ultra sensitive accelerometers of the ESA GOCE mission［C］.The59th International Astronautical Congress，Glasgow，UK，September 28 －October 2，2008，IAC －08 －B1.3.7.

［4］ Bernard A，Touboul P.The GRADIO accelerometer:design and development status［C］.Proceedings of the Workshop on Solid －Earth Mission ARISTOTELES，Anacapri，Italy，1991，ESA SP－329，61－67.

［5］ Touboul P.μSCOPE［C］.ESA －CERN Workshop，F(xiàn)undamental Physics in Space and Related Topics，CERN，April 5 －7，2000.

［6］ Touboul P.MICROSCOPE status［C］.GREX Workshop，Pisa，Italy，October 7－10，2002.

［7］ 薛大同.靜電懸浮加速度計(jì)噪聲測(cè)試數(shù)據(jù)的頻譜分析方法［J］.空間科學(xué)學(xué)報(bào)，2008，28(1):55 －63.［Xue Da－tong.Method for spectral analysis of electrostatically suspended accelerometer noise data［J］.Chinese Journal of Space Science，2008，28(1):55 －63.］

［8］ Nati M，Bernard A，F(xiàn)oulon B，et al.ASTRE － a highly performant accelerometer for the low frequency range of the microgravity environment［C］.The 5th European Symposium on Space Environmental Control Systems and the 24th International Conference on Environmental Systems，F(xiàn)riedrichshafen，Germany，June 20－23，1994.

［9］ 薛大同.靜電懸浮加速度計(jì)伺服控制分析［J］.空間科學(xué)學(xué)報(bào)，2009，29(1):102 － 106.［Xue Da－tong.Servo－control analysis of electrostatically suspended accelerometer［J］.Chinese Journal of Space Science，2009，29(1):102 －106.］

［10］ Frommknecht B. Simulation des sensorverhaltens bei der GRACE－mission ［D］. Diplomarbeit TUM (Technische Universit?t München)1874366，August 31，2001.

［11］ Davis E S，Dunn C E，Stanton R H，et al.The GRACE mission:Meeting the technical challenges［C］.The50th International Astronautical Congress，Amsterdam，Netherlands，October 4－8，1999，IAF －99 －B.2.05.

［12］ 何鐵春，周世勤.慣性導(dǎo)航加速度計(jì)［M］.北京:國防工業(yè)出版社，1983:335.

［13］ Touboul P，F(xiàn)oulon B，Willemenot E.Electrostatic space accelerometers for present and future missions［C］.The 47th International Astronautical Congress，Beijing，China，October 7 － 11，1996，IAF －96 － J.1.02.

［14］ Touboul P，F(xiàn)oulon B.Space accelerometer developments and drop tower experiments［J］.Space Forum，1998，4:145-165.

［15］ Jet Propulsion Laboratory，University of texas center for space reserch，GeoForschungZentrum Potsdam.GRACE Science＆ Mission Requirements Document 327-200，Revision D［R］.March 6，2002，JPL D －15928.

［16］ Touboul P， Foulon B， LISA senseur gravitationnel ＆accéléromètres electrostatiques［C］，1ères journées LISA －France，Paris，Jan 20－21，2005.