喬 遷,解亞軍
(西北工業大學航空學院,西安 710072)
?
盒式天平裝配應力對測量結果的影響研究
喬 遷,解亞軍
(西北工業大學航空學院,西安 710072)
天平安裝過程中會產生裝配應力,對元件輸出值存在干擾,影響測量準度。針對此現象,文中在保證天平整體外形尺寸不變前提下,用有限元的方法研究了浮動框螺栓的緊固程度和浮動框厚度對天平測量準度的影響。結果表明:裝配應力不平衡會使測量結果準度下降,浮動框的厚度也影響測量結果。風洞天平在試驗中應盡量平衡裝配應力,選擇合適的浮動框厚度,以提高測量的準確性。
有限元仿真;盒式天平;螺栓預緊力;浮動框厚度
天平是風洞測力試驗中最重要的測量裝置之一,用于測量作用在模型上的空氣動力載荷(力和力矩)的大小、方向、作用點。測力試驗對風洞天平的技術要求很高,其中最重要的是要求天平具有高的精度與準度[1-3]。
文中選取的六分量盒式天平(型號:0904)外形尺寸為長200 mm,寬100 mm,高80 mm。
該天平主要由浮動框、固定框和彈性連桿組成,固定框上分布有八個懸臂梁,浮動框與固定框之間用8根彈性連桿連接起來,其中Y方向4根,X方向2根,Z方向2根。彈性連桿均使用雙弧形彈性鉸鏈的形式用來分解6個方向的氣動力。天平的校準中心在其幾何中心,如圖1。天平材料為OONi18Co8Mo5 TiAl,彈性模量E=187.25 GPa,泊松比μ=0.3,屈服極限1.754 GPa。天平量程見表1。

圖1 天平校準中心
天平浮動框上使用的螺栓為M12,在進行有限元計算時,設置其基本參數彈性模量E=210 GPa,泊松比μ=0.3。根據德標DIN267,螺栓預緊力采用安全級別3.6,M12螺栓預緊力為10 800 N,文中在ABAQUS中施加Bolt load類型螺栓緊固力,如圖2所示。

圖2 螺栓預緊力示意圖
2.1 天平理想工作狀態下有限元模擬
使用有限元計算軟件ABAQUS[5],對盒式天平采用全局網格控制。天平和螺栓均采用邊長3 mm的四面體線性單元C3D4,單元總數為38萬,網格劃分如圖3所示。
文中為得到較精確的結果,在天平校準中心設立參考點并建立其與浮動框上表面力的約束關系,將各力與力矩加載到參考點上。施加單分量載荷對天平進行測試,即在天平中心參考點分別施加法向力(Fy)、側向力(Fz)、俯仰力矩(Mz)、偏航力矩(My)、滾轉力矩(Mx)、軸向力(Fx)。在固定框底部施加3個方向的固定約束。
該天平靜校報告[6]各分量量程及校驗載荷見表1。
材料1:綠色開花植物在傳粉受精后,胚珠發育為種子,子房發育成果實。番茄、葡萄、茄子等植物若在開花期遇連續陰雨天氣不能正常授粉,不久花朵枯萎脫落,得不到果實。生產實踐中,人們遇到這種情況時,常在未授粉的雌蕊柱頭涂抹或噴灑一定濃度的生長素,子房就能發育成果實,避免惡劣氣候對生產的影響。
首先,在天平浮動框和固定框裝配理想情況(即沒有裝配螺栓預緊力)下對天平進行有限元分析。計算懸臂梁元件的應變大小,并與靜校報告天平工作公式所計算出的應變(實驗值)相比較。進而,分析單元加載時,通過有限元計算所得出的各分量之間的干擾。

表1 單分量載荷
2.2 有限元分析結果
2.2.1 實驗應變值與有限元計算應變值比較
在理想情況下,天平的所有應變均由施加的測量載荷產生,實驗和計算均用校驗載荷,見表1。取模型各處相對應懸臂梁粘貼應變片位置的應變平均值(例如測量施加Fy=4 800 N的載荷時,取固定框底部4個懸臂梁根部應變的平均值,應變云圖如圖4),并與校準實驗天平工作公式計算得出的應變值(實驗值)相比較,見表2。有限元結果相較于校準應變都較小,但總體相差不大。總之,可以認為有限元模擬結果可靠。

表2 應變值對比

圖4 施加單分量Fy=4 800 N時天平底部懸臂梁應變云圖
2.2.2 分量間的相互干擾
通過有限元單元加載法,取各懸臂梁處應變值,并按照組橋方式計算后,各分量的一次干擾量A[7]見表3。
雖然目前對天平一次總干擾量的大小沒有嚴格的規定,但是在《風洞天平》(參考文獻[1])一書中指出天平分量之間的總干擾量要控制在10%之內。由表3可以看出,在有限元模擬的理想狀態下,分量的一次總干擾量均在10%以內,干擾較小。實際靜校的總干擾量A0在Fy、Fz、Mx、Mz方向上較小,均在10%以下,而在Fx、My方向上的一次總干擾量較大超過10%。在軸向力方向影響較大的為俯仰力矩Mz,對側向力影響較大的為軸向力Fx、滾轉力矩Mx和偏航力矩My。經初步分析,由于天平在使用過程中測量結果良好,故排除校驗產生較大誤差的可能。分析天平結構認為可能是由于側向力梁和軸向力梁對于天平中心軸線的平行度或對稱度不良造成的一次干擾項過大。

表3 各分量的干擾應變
3.1 幾種不同工況的螺栓預緊力
為了研究螺栓預緊力對天平測量應變的影響,對模型中的螺栓施加螺栓預緊力。對預緊力的研究分為3種工況:工況一,將4個連接螺栓都施加10 kN的預緊力,研究平衡的螺栓預緊力對天平測量的影響;工況二,將半邊螺栓預緊力減小到5 kN,研究預緊力不平衡對天平懸臂梁應變的影響,如圖5所示;工況三,只裝配一個螺栓(預緊力為10 kN),研究單個螺栓對天平懸臂梁應變的影響。

圖5 施加不均勻的螺栓預緊力
3.2 螺栓預緊力的影響結果
工況一,加4個螺栓:此時單分量懸臂梁應變基本相同,取應變平均值,見表4。所測的3個分量Fy、Mx、Mz的應變值與實驗值更為接近(表2),這是因為在校驗天平時,天平用螺栓固定在天平加載架上,校驗時已有螺栓預緊力,此種加載情況也與試驗情況更為吻合。Fz、Fx和My依舊變化不大,均小于3%。說明實驗過程中螺栓預緊力對天平的測量有一定影響。
工況二,施加不平衡的螺栓預緊力,如圖3。其有限元計算結果見表4。取施加單分量Fy=4 800 N時的有限元計算應變云圖,如圖6。不平衡的預緊力,導致懸臂梁的應變也不平衡,從而導致了測量結果的不精確。

圖6 單邊四個螺栓預緊力底部懸臂梁應變云圖
工況三,加單個螺栓:在只有單個螺栓預緊力10 kN的情況下,預緊力不平衡,四根懸臂梁底部均出現了應變不均衡情況。例如,施加Fy=4 800 N時,距施加預緊力螺栓較近處的連桿應變明顯增大,由原來平均的217 με增加為241 με,其它連桿應變變化不大,4個懸臂梁的平均值為223 με,應變云圖如圖7。測量各力/力矩分量的應變平均值見表4。其中法向連桿所測的3個分量Fy、Mx、Mz均有較大增大。而軸向X連桿所測的Fz,Z向連桿所測的Fx的變化均不大,均在3%之內。進一步說明螺栓預緊力不平衡時對天平測量結果有較大影響。

圖7 單螺栓預緊力施加單分量載荷Fy底部懸臂梁應變云圖

分量四個螺栓預緊輸出應變值(με)兩個螺栓預緊輸出應變值(με)一個螺栓預緊輸出應變值(με)Fy241223(242)221(241)Fx315315310Fz218218214Mx285265(284)239(255)My117711701167Mz384372(384)359(384)
注:括號內為懸臂梁最大應變值。
4.1 改變浮動框的厚度
研究不同的浮動框厚度的影響。原天平模型浮動框厚度為39 mm,削減其厚度:浮動框厚度為分別為34 mm、29 mm時,計算施加單分量載荷天平懸臂梁的應變。得出的結果與原天平有限元模擬結果相比較,判斷天平浮動框厚度的影響。
4.2 有限元計算結果
將不同厚度的浮動框理想狀態的有限元計算結果和原厚度結果比較,見表5。

表5 不同厚度懸臂梁應變值
可以看出,當浮動框厚度為34 mm時,計算結果與原厚度略有減小。單分量載荷Fy、Mx和Mz應變值與原厚度應變值相差不大。而Fx、Fz和My分量,計算結果卻相差較大,最大差值為My分量,應變計算結果減小了5.8%。
當浮動框厚度為29 mm時,計算結果與原厚度計算結果相差較大。Fx、Fz分量的計算結果比原天平減小百分之十幾。My分量計算應變值比原天平減小了14.1%。由有限元計算的應力云圖分析,產生此結果是因為浮動框厚度減小使得浮動框與橫向彈性連桿連接處剛度減小,出現較大的應力和形變,從而分解了懸臂梁的應力,導致懸臂梁應變減小,測量值出現較大誤差。取誤差較大Fz分量3個厚度浮動框與彈性連桿相接處的應力云圖進行比較,如圖8所示:原厚度時,如圖8(a),連接處并沒有出現應力集中情況,如圖9所示的P1、P2兩點的相對位移ΔU為0.000 7 mm;當厚度為34 mm時,如圖8(b)所示的彈性連桿連接處應力增大,其大小為6.5 MPa,ΔU=0.001 32 mm,剛度減小;當厚度為29 mm時圖8(c),連接處應力明顯增大,其大小為9.2 MPa,ΔU=0.002 95 mm,剛度減小更多。
其它3個分量Fy、Mx和Mz的計算的應變值的減小主要是由浮動框厚度減小導致其剛度的減小變形增大造成的。

圖8 單分量載荷Fz作用時不同厚度浮動框與彈性連桿連接處應力云圖
文中利用ABAQUS有限元軟件分析實驗室某型號盒式天平,研究螺栓預緊力和浮動框厚度對天平測量的影響,得出以下結論:

圖9 P1、P2的位置
1)在天平測量過程中,浮動框螺栓預緊力對固定框底部4個懸臂梁的測量結果影響較大,使測得的Fy、Mx和Mz三個分量明顯增大。而對Fx、Fz和My的測量結果影響很小。
2)施加平衡的螺栓預緊力使得模擬結果與實驗結果吻合,不平衡時會有較大誤差。這就要求在盒式天平使用過程中盡可能施加均勻螺栓預緊力,甚至必須使用帶數字顯示的扭力扳手進行操作。
3)浮動框厚度的減小會影響天平的剛度,從而影響天平的測量精準度。當浮動框較薄時,浮動框出現應力集中,使得測量的應變誤差變大。并且,隨著浮動框的厚度減小,螺栓預緊力對測量結果的影響也越來越大。所以,在設計天平時,注意厚度的影響,保證天平有足夠的剛度。
[1] 賀德馨. 風洞天平 [M]. 北京: 國防工業出版社, 2001: 379-396.
[2] 戰培國, 譚飛程. 國外風洞天平校準技術研究進展 [J]. 航空工程進展, 2012(2): 18-20.
[3] ZHANG Zhaoming, HAN Buzhang. Development for high precision six component strain gauge balance [J]. Transactions of Nanjing University of Aeronautics&Astronautics, 2004(6): 5-11.
[4] 茍文選. 材料力學 [M]. 2版. 北京: 科學出版社, 2010: 48-53.
[5] 石亦平, 周玉蓉. ABAQUS有限元分析實例詳解 [M]. 北京: 機械工業出版社, 2006: 393-420.
[6] 解亞軍, 白靜. 500型盒式天平靜校報告 [R]. 西安: 西北工業大學, 2011.
[7] 趙亮亮, 胡國風, 薛偉. 基于虛擬校準技術的應變天平設計準度分析 [C]∥中國空氣動力學會測控技術專業委員會第六屆四次學術交流會論文集. 中國空氣動力學會, 2013.
The Analysis of Assembly Stress of Wind Cassette Balance
QIAO Qian,XIE Yajun
(School of Aeronautics, Northwestern Polytechnical University, Xi’an 710072, China)
Assembly stress from balance installation disturbs output value of component and affects measurement accuracy. In view of this, effect of float box bolt’s fastening degree and thickness of float box on balance measurement accuracy were studied using finite element method without changing overall balance dimension. The results show that imbalanced assembly stress reduces measurement accuracy, floating box thickness also affects the measurement result. In test, wind tunnel balance should try to balance assembly stress, appropriate floating box thickness should be selected to improve accuracy of the measurement.
cassette balance; bolt pre-tightening force; thickness of float box; FEM
2015-09-14
喬遷(1989-),女,河北石家莊人,碩士研究生,研究方向:風洞試驗,風洞天平。
v211.72
A