趙洪波，任順清，李 巍

(哈爾濱工業(yè)大學空間控制與慣性技術研究中心，哈爾濱150001，h84b@163.com)

半球諧振陀螺儀是一種哥氏振動陀螺儀，已在衛(wèi)星和宇宙飛船中得到廣泛的應用［1－2］，其高精度、長壽命特性使其在空間技術方面有廣闊的應用前景，國內(nèi)從十一·五期間已開展了半球諧振陀螺儀的研制工作.在理論與實際應用中半球諧振子的固有頻率是一個很重要的參數(shù)，特別是在由于諧振子密度、厚度等的不均勻性造成的半球諧振子頻率裂解的分析中有重要的作用.文獻［3］通過有限元法計算了半球諧振陀螺儀諧振子的固有頻率，文獻［4］給出了半球諧振陀螺儀諧振子的固有頻率的計算公式，但根據(jù)此公式編程計算的值與實測的固有頻率相差太大.因此本文經(jīng)過一系列的推導，建立了與文獻［4］不同的動力學方程，給出了固有頻率的計算方法，并結合某半球諧振陀螺儀的具體參數(shù)給出了計算的固有頻率，這個頻率與實測結果相近.再取文獻［3］的一組半球諧振陀螺儀的基本參數(shù)，按照本文的方法進行計算并和文獻［3］的結果比較，證明了本文方法的正確性.

1 半球諧振子動力學方程

1.1 半球諧振子任意點的應變方程

如圖1所示，半球殼諧振子中曲面一點的矢徑為R，R在坐標系O－XYZ下的分量為Rsin αcos β，Rsin αsin β，Rcos α(其中R為矢量R的模).把經(jīng)線和緯線作為坐標曲線α，β，它們的切線單位矢量e1，e2和法線單位矢量e3組成一個右手局部坐標系.假設半球殼的中曲面任一點在局部坐標系的位移為M=ue1+ve2+we3，法線方向上和此點距離為Z的一點的位移為MZ= uze1+vze2+wze3，假設變形后坐標系變?yōu)閑'1，e'2，e'3，以θ表示在e1和e3確定的平面內(nèi)，e3向e1的轉角;φ表示e2和e3確定的平面內(nèi)，e3向e2的轉角;w1和w2分別表示在切平面內(nèi)，e1向e2，e2向e1的轉角.R為中曲面內(nèi)一點變形前的矢徑，P為同一點變形后的矢徑，在局部坐標系下表示為基于基希霍夫－李亞夫假設，即薄殼變形前任何垂直于殼體中間表面的直線在變形后依然垂直于這一表面;沿殼體厚度的法線段的長度在變形過程中保持常值;在相鄰的平行于中間表面的薄殼層表面間產(chǎn)生的法向應力與應力張量的其它分量相比是小量，可以忽略不計.再根據(jù)板殼理論、曲面論的知識和幾何關系，得到半球諧振子中距離中曲面為Z的任一點的應變方程［5］為

圖1 半球薄殼一點的位移，變形及局部坐標系

1.2 半球薄殼單元平衡方程

薄殼單元的應力分析如圖2所示，可求出沿中曲面為單位寬度的面積上的內(nèi)力.

在以α，β為法線的截面上，正應力分別為σα，σβ，切應力為ταβ，ταz，τβα，τβz.

圖2 薄殼單元應力圖

圖3為薄殼微元受力分析圖，其中N1，S1，Q1，N2，S2，Q2為沿中曲面為單位寬度的面積上的內(nèi)力，M1，M2，M12，M21為沿中曲面為單位寬度的面積上的力矩，q1，q2，q3分別是外力沿坐標系方向的投影，χ為扭率.根據(jù)面素的平衡條件建立力平衡方程，再求出其力矩平衡方程，從而得到半球薄殼單元平衡方程，如式(2)～(4)所示.

圖3 薄殼微元受力圖

其中:

以下力平衡和力矩平衡方程的推導，可參照文獻［6］.

1.3 半球諧振子二階固有振型的動力學方程

在前面提到的半球薄殼單元平衡方程的推導過程中，在諧振的外力中必須考慮諧振子在慣性空間以角速度 [ΩxΩyΩz]T運動時的慣性力.根據(jù)哥氏定理計算得到諧振子的慣性力［X Y Z］T，并在時得

其中:ρ為諧振子的材料密度

得到以角速率旋轉并處于外分布載荷作用下的半球諧振子邊緣的運動方程為

將諧振子各點的位移按不可拉伸薄殼的二階固有振型展開得

式中:

U(α)，V(α)，W(α)為半球薄殼二階固有振型振動的瑞利函數(shù)，則半球諧振子邊緣某點的位移方程為

用Ax，Ay，Az分別表示式(5)～(7)的左邊，則半球諧振子邊緣的運動方程L可表示為

將式(8)代入式(5)～(7)并利用布勃諾夫－加廖爾金法［4］得

其中V為半球殼區(qū)域，化為球坐標

化簡得

進行積分運算，可得動力學方程為

其中:

其中c0與文獻［4］中的計算結果不同，增加了很多項，使諧振子的動力學方程更加精確，為進行半球諧振子運動機理分析及其質(zhì)量，密度，厚度，楊氏模量等的不均勻性的誤差分析打下基礎.

下面與中電集團二十六所研制的半球諧振陀螺儀的實際參數(shù)進行計算比較.

2 固有頻率的計算及比較

根據(jù)式(10)可知固有頻率ω0的計算公式為

由U(α)=V(α)=sinαtan2(α/2)，W(α)=－(2+cosα)tan2(α/2)，對它們求取關于α的各階導數(shù)U'(α)，U″(α)，U?(α)及W'(α)，W″(α)，W?(α)，W(4)(α)，然后代入c0，d，e和m0的公式中得

目前國內(nèi)生產(chǎn)諧振子的材料為熔融石英，其密度ρ=2 200(kg·m－3)，楊氏模量E=7.67× 1010(N·m－2)，泊松比γ=0.17，它的幾何參數(shù)為中曲面半徑R=0.015 m，厚度h=0.85× 10－3m，將以上各值帶入式(11)得:

這與實際測量結果是一致的.而按照文獻［4］給定的公式，計算結果為ω0=10 950 Hz，與實際不符.

文獻［3］給定的參數(shù)為ρ=2 500(kg·m－3)，γ=0.17，h=10－3m，R=0.014 5 m，E=7.67× 1010(N·m－2)，按照本文的方法計算得 ω0= 2 623.6 Hz，半球殼支柱的長度和半徑對固有頻率影響較大［3］，故考慮支柱的影響后文獻［3］計算結果為ω0=3 640 Hz，與本文的按完整半球殼計算的固有頻率有一定的差別.按照文獻［4］計算得ω0=10 622 Hz，此結果與實際嚴重不符，是不正確的.

3 結論

1)在建立半球諧振子的應變方程，力平衡方程的基礎上得到了半球諧振子的動力學方程，比文獻［4］中的動力學方程更加準確;

2)通過與有限元法計算的固有頻率比較，證明了本文固有頻率的計算公式的正確性;

3)通過實際的半球諧振陀螺儀的參數(shù)計算并與實際結果比較，證明了本文固有頻率的計算公式及結果的正確性;

4)建立的半球諧振子動力學方程是今后半球諧振子的質(zhì)量，密度，厚度，楊氏模量等的不均勻性的誤差分析的基礎.

［1］呂志清.半球諧振陀螺在宇宙飛船中的應用［J］.壓電與聲光，1999，21(5):349－353.

［2］EMILY L B，ALLAN Y L.In－flight characterization of cassini inertial reference units［C］//AIAA Guidance，Navigation and Control Conference and Exhibit.Hilton Head，South Carolina:AIAA，2007.

［3］楊倩.基于半球諧振陀螺的慣性導航系統(tǒng)研究［D］.哈爾濱:哈爾濱工業(yè)大學，2007:16－22.

［4］МАТВЕЕВ В А，ЛИПАТНИКОВ В И，АЛЕХИН А В. Проектирование.Волнового твердотелъного гироскопа［M］.［S.l.］:Издателъство МГТУ имени Н.Э.Баумана，1998:1－43.

［5］劉鴻文，林建興，曹曼玲.板殼理論［M］.杭州:浙江大學出版社，1987.

［6］QI Jiayi，CHEN Xue，BEN Shunqing，et al.Dynamic error mechanism analysis of HRG［C］//ISSCAA2008.Shenzhen:［s.n.］，2008.

［7］梁嵬.半球陀螺儀諧振子振動特性及其結構研究［D］.長春:長春理工大學，2008:17－18.

［8］雷霆.半球諧振陀螺控制技術研究［D］.重慶:重慶大學，2006:9－10.

［9］ZHBANOV Y K.Vibration of a hemispherical shell gyro excited by an electrostatic field［J］.Mechanics of Solids，2008，43(3):328－332.

［10］ZHBANOV Y K，ZHURAVLEV V P.Effect of movability of the resonator center on the operation of a hemispherical resonator gyro［J］.Mechanics of Solids，2007，44(3):851－859.

［11］高勝利，吳簡彤.半球諧振陀螺的漂移機理及其控制［J］.彈箭與制導學報，2008，28(3):61－64.

［12］FAN Shangchun，LIU Guangyu，WANG Zhenjun.On Flexural Vibration Of Hemispherical Shell［J］.Applied Mathematics and Mechanics，1991，12(10):1023－1030.