朱世清,連蔭俊,郝永波,3,李 博,臧立彬
(1.北京控制工程研究所,北京 100190; 2.天津航天機電設(shè)備研究所,天津 300301;3.天津市微低重力環(huán)境模擬技術(shù)重點實驗室,天津 300301)
空間飛行器在軌服務(wù)與維護是未來航天發(fā)展的重要方向,也是航天領(lǐng)域研究的熱點之一[1-9]。應(yīng)用在軌服務(wù)飛行器執(zhí)行任務(wù)之前,必須對其開展半物理仿真和全物理仿真,以驗證其控制策略及執(zhí)行效能。為在地面微重力模擬環(huán)境下對空間飛行器進行有效驗證,需建立一個大范圍運動模擬半物理仿真平臺[10],以模擬服務(wù)目標(中小型)在復(fù)雜初始運動條件下被操作時的運動學(xué)過程。但在地面上進行運動模擬時,機械傳動部件之間的摩擦、運動系統(tǒng)控制精度和反饋控制中的動態(tài)響應(yīng)誤差等會使運動模擬系統(tǒng)產(chǎn)生振動[11-12];尤其對于大范圍運動模擬系統(tǒng),縱向運動的二維平動系統(tǒng)主梁和橫向運動的小車在急?;騿訒r產(chǎn)生的慣性力容易誘發(fā)系統(tǒng)的振動[13]。此外,由于質(zhì)量受控,大跨度系統(tǒng)剛性降低,系統(tǒng)動態(tài)運行時含有復(fù)雜的激勵源,如二維平動系統(tǒng)起停和加減速,模擬目標星周期自旋運動等,都可能引起整個運動模擬系統(tǒng)的振動。
有文獻提出通過減小加速度、增強系統(tǒng)實時性的方法來抑制運動系統(tǒng)的振動,但這種方法會影響系統(tǒng)的響應(yīng)速度[14],因此不適用于大范圍三維平動系統(tǒng)。通過優(yōu)化電機系統(tǒng)的PID 參數(shù)也可以起到一定的抑振效果,但對于大范圍運動模擬系統(tǒng)來說,抑振效果有限[15-16]。
因此,本文通過簡化三維平動系統(tǒng)的力學(xué)模型,分析系統(tǒng)對慣性激勵的響應(yīng)。探究試驗小車在縱向二維平動系統(tǒng)主梁不同位置下的系統(tǒng)響應(yīng)、固有頻率,以及伸縮立柱長度的影響,確定系統(tǒng)振動的主要來源;并提出一種由彈簧、質(zhì)量塊及阻尼器組成的調(diào)諧質(zhì)量阻尼器(tuned mass damper, TMD)抑振方法,并通過仿真分析和實驗驗證所采用的動能吸振裝置的減振效果,旨在減小系統(tǒng)結(jié)構(gòu)破壞的風(fēng)險,提高系統(tǒng)的運行可靠性。
大范圍運動模擬系統(tǒng)為龍門懸吊式三維平動裝置,主要由橫向二維平動系統(tǒng)(包括x向主梁和y向試驗小車)和豎向平動系統(tǒng)(z向伸縮立柱)等組成,具體結(jié)構(gòu)參見圖1。該運動模擬系統(tǒng)的x向運動行程可達60 m,y向運動行程可達40 m。本文通過有限元方法建立該運動模擬系統(tǒng)的精確模型,用來分析在軌服務(wù)飛行器運動過程中所產(chǎn)生的振動特性。從圖1 還可看到整個龍門懸吊式大范圍三維平動裝置在GNC 廠房中的布局關(guān)系,實驗室原二維平動系統(tǒng)(天車)位于南側(cè),龍門懸吊式三維平動裝置位于北側(cè)。大范圍三維平動裝置用于安裝和承載博弈目標運動模擬器(即并聯(lián)機構(gòu)和模擬目標星),并實現(xiàn)其在實驗室三維空間的大范圍平動運動。
圖1 大范圍運動模擬系統(tǒng)結(jié)構(gòu)及布局示意Fig.1 Structure and layout of the large-scale kinematic simulation system
三維平動系統(tǒng)的主梁采用雙梁桁架結(jié)構(gòu),單根桁架梁由4 段拼接而成,拼接處用銷軸和高強度鉸接螺栓連接固定,梁跨距45.6 m、高度2.1 m、端梁寬度8.0 m,如圖2 所示。單段桁架梁采用型材焊接而成,為了減重并提高剛性,桁架梁截面形狀設(shè)計為三角形:下弦梁選用2 根135 mm×135 mm×12 mm 的方管,上弦梁采用250H 型鋼,斜腹桿采用φ76 mm×6 mm 圓管,材料均為Q355B。伸縮立柱同為桁架梁焊接結(jié)構(gòu),采用100 mm×100 mm×5 mm 的方管焊接而成,材料為Q345 鋼。
圖2 桁架梁結(jié)構(gòu)及有限元模型Fig.2 Structure and finite element model of the truss beam
為分析桁架梁剛度,邊界條件定義如圖3 所示,梁結(jié)構(gòu)作簡支約束,分析時計算帶載和不帶載2 種工況,帶載時單根主梁的承載為10 t/2=5 t,載荷與桁架主梁的z向距離為5.6 m。
圖3 桁架梁邊界條件設(shè)置Fig.3 Boundary condition setting for the truss beam
在建立模型前,為便于仿真分析,對系統(tǒng)做如下簡化:
1)去除縱向大車行走平臺的4 個車輪,在車輪與水平滑臺的4 個接觸圓柱面處添加固定約束;
2)將水平滑臺和縱向二維平動系統(tǒng)主梁之間的輪-軌接觸改為面-面接觸;
3)將立柱和博弈目標運動模擬器(即并聯(lián)機構(gòu)和模擬目標星)用等效質(zhì)量塊替代。
為簡化雙梁計算,采用荷載減半的形式分析單梁變形情況,二維平動系統(tǒng)中段承載豎直向下5 t載荷(加載方式見圖4(b))和空載時的桁架梁變形如圖4 所示。仿真結(jié)果顯示二維平動系統(tǒng)最大變形發(fā)生在中段,桁架梁整體因載荷作用呈向下彎曲的形態(tài),空載和5 t 載荷下的最大變形分別為46.1 mm和78.1 mm。
圖4 桁架梁靜態(tài)載荷及形變分析Fig.4 Static load and deformation analysis of the truss beam
由于該運動模擬系統(tǒng)具有多個自由度運動,振動來源比較復(fù)雜,主要有:
1)縱向運動的二維平動系統(tǒng)主梁和橫向運動的小車在啟動或剎車/急停時產(chǎn)生的慣性力;
2)模擬目標星高速旋轉(zhuǎn)過程中由于不平衡產(chǎn)生的簡諧振動;
3)并聯(lián)機構(gòu)作動時產(chǎn)生的簡諧激勵。
試驗小車在主梁上運動到不同位置時,主梁的變形是不同的;伸縮立柱伸長量不同時,立柱結(jié)構(gòu)的尺寸和剛度均發(fā)生變化。可見,整個結(jié)構(gòu)的基頻是隨著運動工況發(fā)生變化的,系統(tǒng)為一個變結(jié)構(gòu)系統(tǒng)。因此,運動系統(tǒng)的模態(tài)分析應(yīng)針對不同工況分別進行。
2.2.1 工況劃分和編號
按試驗小車在主梁上的位置和伸縮立柱長度對系統(tǒng)的工況進行劃分和編號,詳見表1。
表1 工況編號Table 1 Numbering of working conditions
2.2.2 模態(tài)分析結(jié)果
計算不同工況下運動系統(tǒng)前6 階固有頻率,如表2 所示。分析表中數(shù)據(jù)可知:當伸縮立柱長度保持不變時,試驗小車偏離主梁中心的距離越大,運動系統(tǒng)的后4 階固有頻率就越小,而對前2 階固有頻率影響不大;當試驗小車在主梁上的位置保持不變時,伸縮立柱長度對后4 階固有頻率影響不大。
表2 系統(tǒng)模態(tài)頻率Table 2 System modal frequencies
2.3.1 材料及邊界條件設(shè)置
進行有限元仿真分析時,首先要設(shè)置模型各部分材料屬性,添加邊界條件,進行網(wǎng)格剖分處理等。
1)將二維平動系統(tǒng)主梁、試驗小車、伸縮立柱和并聯(lián)機構(gòu)的材料設(shè)置為結(jié)構(gòu)鋼,將模擬目標星等效為質(zhì)量塊;
2)將模型中的全部裝配接觸關(guān)系(即所有接觸面)均設(shè)為Bonded 類型;
3)對于約束條件,在承載輪與試驗小車的4 個接觸圓柱面處添加固定約束。
2.3.2 載荷設(shè)置
如2.1 節(jié)所述,引起振動的載荷主要有慣性力、模擬目標星偏心旋轉(zhuǎn)激勵和并聯(lián)機構(gòu)簡諧激勵。
1)慣性力
根據(jù)實際工況,給予系統(tǒng)x向和y向大小為0.2 m/s2的加速度載荷,載荷為階躍函數(shù)。
2)模擬目標星偏心旋轉(zhuǎn)激勵
模擬目標星質(zhì)量為200 kg,質(zhì)心偏心100 mm,自旋角速度為65 (°)/s。在模擬目標星自旋過程中,由于偏心產(chǎn)生的激勵可以分解為沿y向和沿z向的簡諧激勵fy和fz,均為正弦激勵函數(shù)。
3)并聯(lián)機構(gòu)簡諧激勵
根據(jù)并聯(lián)機構(gòu)的作動能力,最大可使平臺z向產(chǎn)生0.1 mm 的位移,頻率為10 Hz。因此在仿真分析時給予并聯(lián)機構(gòu)一個幅值0.1 mm、頻率10 Hz 的簡諧位移激勵。
2.3.3 結(jié)果分析
不同工況在3 種激勵分別作用下,伸縮立柱末端的響應(yīng)如表3 所示。分析表中數(shù)據(jù)可知,主振動方向是由主梁x向作動引起的。因此減振器的功能應(yīng)為減少衛(wèi)星等效質(zhì)量塊沿x方向、頻率為2 Hz 的振動,這符合弱剛性模型的假設(shè),可以采用TMD來抑制振動。
表3 伸縮立柱末端在3 種激勵下的響應(yīng)Table 3 Response of telescopic rod end under three excitations
通過分析得知,模擬目標星的偏心、并聯(lián)機構(gòu)以及慣性力均會對系統(tǒng)產(chǎn)生振動激勵,且均會導(dǎo)致主梁發(fā)生x、y、z三個方向的形變,其中x向激勵產(chǎn)生的振動幅值最大。若減振系統(tǒng)可以實現(xiàn)x方向的振幅減小,則可以在其余兩方向也加裝減振系統(tǒng)以實現(xiàn)振動隔離。因此本文將作用在主梁上x方向的激勵作為振動的主要激振力,采用TMD 來吸收主梁所受到的振動激勵。
在待減振結(jié)構(gòu)上安裝TMD 實現(xiàn)振動隔離的系統(tǒng)可以等效為二自由度系統(tǒng),其理論模型見圖5。圖中m1、k1、c1、x1分別為系統(tǒng)主結(jié)構(gòu)的質(zhì)量、彈簧剛度、阻尼系數(shù)、位移響應(yīng);m2、k2、c2、x2分別為TMD 系統(tǒng)的質(zhì)量、彈簧剛度、阻尼系數(shù)、位移響應(yīng);F0表示系統(tǒng)施加的外部激勵。
圖5 TMD 理論簡化模型Fig.5 Theoretically simplified model of TMD
將TMD 安裝在系統(tǒng)主結(jié)構(gòu)上,主要是將彈簧和阻尼器連接到待減振結(jié)構(gòu)上,則在外部激勵作用下,待減振結(jié)構(gòu)與 TMD 同時振動,TMD 振動會吸收待減振結(jié)構(gòu)的振動能;此外,由于 TMD 的調(diào)諧減振使原有系統(tǒng)在共振頻率處的振幅減小,故可降低共振對原有系統(tǒng)造成的破壞。
建立圖5 所示二自由度振動系統(tǒng)動力學(xué)模型[17]
根據(jù)文獻[17],獲得主系統(tǒng)位移振幅比量綱為1 的表達式
其中:μ為等效質(zhì)量比,μ=m2/m1;s為激振力與主系統(tǒng)固有頻率之比;α為TMD 與主系統(tǒng)固有頻率之比;ζ為阻尼比。
根據(jù)3.1 節(jié)的理論分析,對TMD 進行設(shè)計,取質(zhì)量比μ=0.02,阻尼比ζ=0.1,頻率比s=0.98。TMD主要由1 個質(zhì)量塊、4 根彈簧以及2 個阻尼器組成,外形尺寸為550 mm×500 mm×500 mm,結(jié)構(gòu)如圖6 所示。
圖6 TMD 結(jié)構(gòu)三維模型Fig.6 Three-dimensional model of the TMD structure
為探究TMD 對主梁的減振效果,將所設(shè)計TMD 安裝在伸縮立柱根部,對系統(tǒng)在慣性載荷激勵下x方向上的振動情況進行仿真,結(jié)果如圖7 所示??梢钥吹?,TMD 能夠在固有頻率范圍內(nèi)實現(xiàn)對系統(tǒng)(x方向)的有效減振,大幅降低振動幅值。
圖7 TMD 對主梁x 方向振動的抑制效果Fig.7 Suppresing effect of TMD on x direction vibration of the main beam
本文針對在軌服務(wù)目標星大范圍運動模擬系統(tǒng)低頻振動問題,通過仿真方法研究了運動模擬系統(tǒng)振動機理,提出振動抑制方法。首先分析了不同運行工況下的系統(tǒng)響應(yīng)、固有頻率以及激勵分別作用時模擬目標星在各個方向上的最大位移,確定了系統(tǒng)振動的主要來源;然后設(shè)計TMD 來降低系統(tǒng)振動幅值以提高系統(tǒng)穩(wěn)定性。研究結(jié)果表明:在不同運行工況下,系統(tǒng)的固有頻率不同,為變結(jié)構(gòu)系統(tǒng);采用TMD 可以有效抑制系統(tǒng)振動,降低系統(tǒng)發(fā)生共振的風(fēng)險,進而提高系統(tǒng)的運行可靠性。