李振哲 華璟 孫杰 周徐斌 劉興天?

(1.上海衛(wèi)星工程研究所，上海 201109)(2.上海航天技術(shù)研究院，上海 201109)

引言

高分辨率是遙感探測類衛(wèi)星發(fā)展的永恒方向，分辨率的提升使得搭載的相機(jī)口徑以及天線陣面不斷向大尺度，大撓性方向發(fā)展.此外，衛(wèi)星功率的增大也使太陽陣等撓性結(jié)構(gòu)進(jìn)一步加大.撓性結(jié)構(gòu)的增多會導(dǎo)致衛(wèi)星姿態(tài)機(jī)動或受外界擾動后產(chǎn)生頻率低、持續(xù)時間長的振動，進(jìn)而影響姿態(tài)控制效果和衛(wèi)星成像性能[1].為了保證衛(wèi)星撓性結(jié)構(gòu)的振動快速衰減，采取被動阻尼減振是可靠性高、性能優(yōu)良的方法，該方法已在多個衛(wèi)星上在軌應(yīng)用[2].

著名的哈勃望遠(yuǎn)鏡太陽陣就采取了黏性阻尼減振措施，還使用Honeywell公司設(shè)計的流體阻尼器來隔離飛輪產(chǎn)生的擾動[3]；Davis等人在兩參數(shù)流體阻尼器的基礎(chǔ)上設(shè)計了一種固有頻率為1.5 Hz的三參數(shù)被動流體阻尼器來隔離反作用飛輪產(chǎn)生的微振動[4]；美國國家航空航天局研發(fā)的D-strut流體阻尼器用于航天器桁架結(jié)構(gòu)的減振與控制力矩陀螺的隔振上[5].國內(nèi)的陸澤瑜[6]、鄧國紅[7]、彭志科[8]還使用了非線性阻尼進(jìn)行隔振研究.對于流體阻尼器的仿真國內(nèi)學(xué)者也進(jìn)行了大量研究.徐曉東使用Ansys CFX軟件仿真分析不同頻率的激勵對孔隙式流體阻尼器性能產(chǎn)生的影響[9]；陳威使用MATLAB軟件模擬分析了阻尼孔大小、阻尼孔長度、硅油黏度等對流體阻尼器性能的影響[10]；丁行武使用CFD軟件對孔隙式阻尼器進(jìn)行了仿真計算與實驗分析[11].

本文設(shè)計了一種新型流體阻尼器，通過對結(jié)構(gòu)增加阻尼器件，提高撓性結(jié)構(gòu)的阻尼比，在軌時快速衰減撓性結(jié)構(gòu)的低頻振動，減少撓性附件低頻振動對整星姿軌控系統(tǒng)的影響，提升衛(wèi)星的在軌穩(wěn)定度.

本文的創(chuàng)新點通過有限元軟件對不同參數(shù)的阻尼器進(jìn)行流固耦合仿真分析，得到了阻尼孔徑、阻尼孔長度對阻尼力的影響規(guī)律，驗證了最佳阻尼系數(shù)的存在，為阻尼器的性能分析和工程應(yīng)用提供設(shè)計指導(dǎo).

1 流體阻尼器結(jié)構(gòu)與工作原理

設(shè)計的流體阻尼器結(jié)構(gòu)示意圖如圖1所示，主要由上阻尼筒、U形柔性筒、下阻尼筒、波紋管、彈簧等組成.上阻尼筒、U形柔性筒與下阻尼筒按示意圖中順序相連接，通過結(jié)構(gòu)間的空隙構(gòu)成上腔體.下阻尼筒與波紋管通過焊接相連，波紋管下端與下阻尼筒之間安裝預(yù)緊彈簧提供背壓，通過結(jié)構(gòu)間的空隙構(gòu)成下腔體.上下兩個阻尼腔通過下阻尼筒中的阻尼孔連接，構(gòu)成一個完整的阻尼器內(nèi)部流場.

圖1 流體阻尼器結(jié)構(gòu)示意圖

如果阻尼器受到壓縮或拉伸，將會產(chǎn)生軸向的變形，壓縮或拉伸U形柔性筒，引起上阻尼腔內(nèi)阻尼液體積減小或增大，從而使阻尼液流入或流出阻尼孔，產(chǎn)生阻尼力，起到振動抑制的作用.將柔性筒設(shè)計成U形結(jié)構(gòu)，其外徑遠(yuǎn)大于阻尼孔的直徑，當(dāng)柔性筒發(fā)生很小的形變時，阻尼液的體積會發(fā)生較大的變化，從而可以增加流經(jīng)阻尼孔的流量，在較小激勵下產(chǎn)生較大的阻尼.

2 流體阻尼器仿真計算

為了研究流體阻尼器的阻尼力、滯回耗能性能及等效阻尼系數(shù)等參數(shù)，本文使有限元軟件對阻尼器進(jìn)行流固耦合仿真計算，對硅油通過阻尼器內(nèi)阻尼孔后阻尼器內(nèi)部的流場進(jìn)行仿真.首先用三維建模軟件建立流體和固體的三維模型，將建立好的模型導(dǎo)入有限元軟件中，分別對流體和固體部分進(jìn)行網(wǎng)格劃分和計算設(shè)置，然后把劃分完成的流體和固體進(jìn)行耦合分析.同時在仿真分析時，作出如下假設(shè)：(1)流體不可壓縮；(2)忽略溫度影響；(3)流體流動滿足流量的連續(xù)性條件.

2.1 有限元模型

阻尼器的有限元模型分別由流體部分和固體部分組成.固體部分由上阻尼筒，柔性筒以及下阻尼筒組成，計算模型如圖2所示，其中下阻尼筒的下端與波紋管相連接，波紋管的下端受彈簧壓縮，提供背壓，本次采取在出口設(shè)置背壓的方式進(jìn)行彈簧力的模擬.實體模型采用六面體單元，流體部分采用四面體單元進(jìn)行網(wǎng)格劃分

(a) 流體部分有限元模型

2.2 阻尼系數(shù)計算方法

流體阻尼器動力學(xué)模型的阻尼系數(shù)通常使用有限元軟件仿真分析或者實驗測試得到.本文通過仿真分析得到流體阻尼器固定端支反力與流體阻尼軸向位移間接計算出流體阻尼器動力學(xué)模型的阻尼系數(shù).

由于流體阻尼器的阻尼特性，對流體阻尼器施加正弦位移激勵時，輸入的位移與輸出的流體阻尼器支反力的數(shù)值會形成一個封閉的曲線，這個封閉曲線被稱為遲滯回線，該曲線圍成的面積與流體阻尼器由于阻尼作用在正弦位移激勵下一個周期內(nèi)消耗的能量呈正比.本文通過遲滯回線的方法處理仿真得到的流體阻尼器固定端支反力與流體阻尼軸向位移，以此計算流體阻尼器的阻尼系數(shù).

因此，當(dāng)對流體阻尼器施加一個正弦位移激勵x(t)=x0sin(ω0-θ)t時，產(chǎn)生的阻尼力為：

(1)

在一個周期內(nèi)流體阻尼器消耗的能量為：

(2)

此時流體阻尼器阻尼系數(shù)為：

(3)

如圖3所示，本次仿真將A面設(shè)置為固定端并記錄該面的支反力，對B面施加正弦速度激勵，最大位移為0.025 mm，將內(nèi)部與流體接觸的面設(shè)置為耦合面進(jìn)行仿真分析.通過流固耦合的方法仿真得到流體阻尼器固定端支反力與流體阻尼軸向位移，以此得到各工況的遲滯回線，并以此計算流體阻尼器的阻尼系數(shù).

圖3 流固耦合參數(shù)設(shè)置

本文通過仿真得到的阻尼孔孔徑0.6 mm、0.8 mm、1.0 mm、1.2 mm工況下遲滯回線如圖4所示.通過遲滯回線的面積對流體阻尼器各參數(shù)進(jìn)行計算分析.

(a) 0.6 mm孔徑遲滯回線圖

3 仿真結(jié)果分析

3.1 阻尼器內(nèi)流場特性

圖5為激勵速度最大時阻尼器流體區(qū)域的壓力場分布圖，此時流場壓力達(dá)到最大值.從圖中可以看出左右兩側(cè)腔體內(nèi)壓力基本保持恒定，而阻尼孔中壓力從右至左呈線性增大的規(guī)律.

圖5 流體壓力場分布圖

提取了阻尼孔內(nèi)A,B,C三點處流體壓力隨時間的變化規(guī)律，如圖6所示，阻尼孔內(nèi)各點的流體壓力隨時間成正弦變化，當(dāng)激勵速度達(dá)到最大值時，流體壓力也達(dá)到最大值.壓力變化的幅值由A點到C點逐漸減小.

(a) 點A處流場壓力隨時間變化曲線圖

3.2 阻尼孔徑對阻尼系數(shù)的影響

流體阻尼器阻尼孔徑是影響其隔振性能的重要因素之一，阻尼孔的直徑會影響阻尼器的動剛度和阻尼系數(shù).本文分析了改變阻尼孔徑對阻尼器隔振性能的影響，計算參數(shù)如表1所示.

表1 計算參數(shù)

圖7是在1Hz激勵下不同孔徑時阻尼系數(shù)隨孔徑大小的變化規(guī)律.由圖可知，隨著阻尼孔徑的增加，在阻尼活塞中流經(jīng)阻尼孔的流量先增大后減小，因而阻尼系數(shù)隨阻尼孔徑的增大而先增加后減小.

圖7 不同孔徑阻尼系數(shù)變化規(guī)律

3.3 阻尼孔長對阻尼系數(shù)的影響

流體阻尼器阻尼孔長也是影響其隔振性能的重要因素之一.本文分析了改變阻尼孔長對阻尼器隔振性能的影響，計算參數(shù)如表2所示.

表2 計算參數(shù)

圖8是在1 Hz激勵下不同孔長時阻尼系數(shù)隨孔長度的變化規(guī)律.由圖可知，阻尼系數(shù)隨著孔長先增大后減小，因為孔長較小時，主要是阻尼活塞內(nèi)流經(jīng)阻尼孔的流體產(chǎn)生阻尼，隨著阻尼孔的增長，流經(jīng)阻尼孔的流體摩阻耗能增多，阻尼系數(shù)增大，當(dāng)阻尼孔繼續(xù)增長時，液阻產(chǎn)生的影響變大，阻尼活塞內(nèi)被壓縮的流體增多，因此阻尼孔內(nèi)的流量減少，導(dǎo)致阻尼系數(shù)減小.

圖8 不同孔長阻尼系數(shù)變化規(guī)律

3.4 激勵頻率對阻尼系數(shù)的影響

流體阻尼器的阻尼系數(shù)隨著激勵頻率不斷改變，因此，本文分析了激勵頻率對阻尼系數(shù)的影響，計算參數(shù)如表3所示.

表3 計算參數(shù)

圖9是阻尼系數(shù)隨激勵頻率的變化曲線.由圖中可知，隨著激勵頻率的增大，阻尼系數(shù)減小，其原因在于激勵頻率增大時，液阻迅速增大，阻尼活塞內(nèi)被壓縮的流體也增大，因此通過阻尼孔的流量減小，阻尼系數(shù)也減小.

圖9 不同激勵頻率阻尼系數(shù)變化規(guī)律

4 結(jié)論

本文設(shè)計了一種流體阻尼器，并采用流體力學(xué)仿真軟件對其進(jìn)行了流固耦合仿真分析.仿真結(jié)果表明，在低頻激勵下，隨著孔徑與孔長的增加，阻尼系數(shù)先增大后減小，以最大阻尼系數(shù)為目標(biāo)，存在對應(yīng)的最佳阻尼孔徑與孔長，其大于或小于最佳孔徑與孔長，阻尼系數(shù)均減小.隨著激勵頻率的增加，阻尼系數(shù)減小.在實際應(yīng)用與撓性附件減振時，應(yīng)根據(jù)結(jié)構(gòu)和基頻信息合理設(shè)計參數(shù)以達(dá)到最佳阻尼比和減振效果.本文對此類低頻振動控制用阻尼器的設(shè)計具有一定的指導(dǎo)作用.