石川千周徐斌張永濤，2杜冬李昊

（1上海衛(wèi)星工程研究所，上海 200240）

（2上海利正衛(wèi)星應(yīng)用技術(shù)有限公司，上海 200240）

點(diǎn)式連接航天器等效剛度模態(tài)分析方法

石川千1周徐斌1張永濤1，2杜冬1李昊1

（1上海衛(wèi)星工程研究所，上海 200240）

（2上海利正衛(wèi)星應(yīng)用技術(shù)有限公司，上海 200240）

針對點(diǎn)式連接航天器受接觸影響而產(chǎn)生的非線性特性，提出了利用線性等效剛度代替非線性剛度以計(jì)算點(diǎn)式連接航天器模態(tài)的分析方法。首先采用有限元軟件對點(diǎn)式連接航天器的非線性剛度進(jìn)行分析，并基于非線性振動和結(jié)構(gòu)動力學(xué)的原理，建立了等效剛度的計(jì)算方法。將等效剛度代入廣義彈簧單元，最終獲得點(diǎn)式連接航天器的簡化模型并分析其模態(tài)。模態(tài)分析結(jié)果與試驗(yàn)結(jié)果對比顯示，相比于傳統(tǒng)分析方法，基于等效剛度的簡化模型精確地預(yù)測了點(diǎn)式連接航天器的模態(tài)，有利于提高點(diǎn)式連接航天器的動力學(xué)仿真分析精度。

航天器；點(diǎn)式連接；模態(tài)；等效剛度；非線性特性

1 引言

隨著航天技術(shù)應(yīng)用時(shí)代和外層空間實(shí)用時(shí)代的到來，結(jié)構(gòu)大型化、復(fù)雜化和模塊化已成為航天器的一個(gè)發(fā)展趨勢［1］，點(diǎn)式連接多被應(yīng)用在航天器本體和部件之間、艙段之間、航天器與運(yùn)載火箭之間以及航天器之間。在保證強(qiáng)度和剛度的要求下，點(diǎn)式連接的螺栓數(shù)量（4～12個(gè)）一般少于普通艙段連接的螺栓數(shù)量（48～64個(gè)，參見國內(nèi)外典型衛(wèi)星平臺），而且點(diǎn)式連接結(jié)構(gòu)簡單輕便，滿足了航天器輕量化的指標(biāo)要求，并具有可靠性高、承載能力強(qiáng)的特點(diǎn)。例如某大型高軌衛(wèi)星艙段間連接節(jié)圓直徑1.2m，采用6個(gè)爆炸螺栓連接；又如大型運(yùn)載火箭長征五號（CZ－5）星箭連接節(jié)圓直徑3m，采用12個(gè)螺栓連接。由于連接結(jié)構(gòu)一般具有應(yīng)力集中、邊界非線性的問題，還可能伴隨材料非線性、幾何非線性的情況，并與預(yù)緊力、摩擦力等因素密切耦合［2］，導(dǎo)致此類結(jié)構(gòu)的連接剛度下降，動力學(xué)特性分析較為困難，成為長期關(guān)注的研究熱點(diǎn)。特別是，連接剛度變化會直接影響航天器的模態(tài)特性［3－5］，比如文獻(xiàn)［6］發(fā)現(xiàn)連接剛度的損失對系統(tǒng)的彎曲模態(tài)有嚴(yán)重的影響。而航天器的基頻及其對應(yīng)模態(tài)是航天器結(jié)構(gòu)設(shè)計(jì)的重要特性。因此，準(zhǔn)確描述包含點(diǎn)式連接結(jié)構(gòu)的動力學(xué)行為是一個(gè)極具應(yīng)用價(jià)值，同時(shí)也亟待解決的問題。

目前，由于有限元軟件不具備直接計(jì)算非線性連接結(jié)構(gòu)模態(tài)的能力，傳統(tǒng)的方法重點(diǎn)關(guān)注連接部位的簡化建模方式，將非線性連接結(jié)構(gòu)轉(zhuǎn)化為線性結(jié)構(gòu)，主要有三種方式：第一種是全固連建模方法，忽略連接的存在，將模型建立成不考慮螺栓連接和接觸界面的全固連連續(xù)系統(tǒng)［6］，但是由于機(jī)械連接結(jié)構(gòu)是不連續(xù)的，加之摩擦對接觸界面的影響，使計(jì)算頻率明顯大于試驗(yàn)數(shù)據(jù)，給結(jié)構(gòu)設(shè)計(jì)帶來風(fēng)險(xiǎn)；第二種是節(jié)點(diǎn)約束方法，僅約束螺栓孔附近的區(qū)域［7］，此方法確實(shí)有效降低了固有頻率，但由于結(jié)構(gòu)的接觸面積隨載荷變化，其忽略了結(jié)構(gòu)壓縮狀態(tài)和彎曲狀態(tài)接觸面積的變化情況，使點(diǎn)式連接結(jié)構(gòu)計(jì)算頻率偏低，增加了不必要的設(shè)計(jì)成本；第三種方法用桿單元、彈簧單元等模擬螺栓結(jié)構(gòu)［8］，這種方法同樣沒有考慮接觸面積隨載荷的變化情況；所以這三種簡化建模方法均不能真實(shí)反映連接特性。尤其針對點(diǎn)式連接結(jié)構(gòu)，傳統(tǒng)的計(jì)算結(jié)果往往與試驗(yàn)結(jié)果相去甚遠(yuǎn)，例如，當(dāng)前國內(nèi)在研的某大型衛(wèi)星（上下串聯(lián)的雙星）的發(fā)射質(zhì)量是3.9t，采用長征三號丙運(yùn)載火箭（CZ－3C）發(fā)射，雙星之間采用6個(gè)均布的分離螺母進(jìn)行連接。運(yùn)載火箭要求該衛(wèi)星的橫向一階頻率大于10Hz，基于傳統(tǒng)的模態(tài)分析方法，建立兩種線性模態(tài)分析模型：第一種僅將雙星連接界面的6個(gè)連接點(diǎn)節(jié)點(diǎn)約束，計(jì)算頻率為8Hz；第二種雙星界面全部固連，計(jì)算頻率為15Hz，兩者差別高達(dá)87.5%，根本無法通過傳統(tǒng)方法考量結(jié)構(gòu)的真實(shí)模態(tài)。近十余年來，伴隨著大規(guī)模工程結(jié)構(gòu)仿真預(yù)測，發(fā)展起來一種描述接觸界面間本構(gòu)關(guān)系模型的方法，以模擬連接結(jié)構(gòu)的動力學(xué)行為［9－10］，但是這些本構(gòu)關(guān)系模型對參數(shù)有嚴(yán)格的要求，往往需要從大量的試驗(yàn)中獲得，因而難以推廣。

本文提出一種通過等效剛度計(jì)算航天器模態(tài)的方法。利用有限元模型［11－13］并結(jié)合非線性振動［14］和結(jié)構(gòu)動力學(xué)［15］的原理計(jì)算點(diǎn)式連接的等效剛度。然后將等效剛度賦予點(diǎn)式連接航天器簡化模型進(jìn)行模態(tài)分析，與試驗(yàn)數(shù)據(jù)對比，驗(yàn)證了此方法的有效性。

2 點(diǎn)式連接航天器模態(tài)與等效剛度分析

2.1 模態(tài)分析

為了準(zhǔn)確計(jì)算點(diǎn)式連接航天器的模態(tài)，本研究采用廣義彈簧單元［16］，建立點(diǎn)式連接航天器有限元簡化模型，通過給廣義彈簧單元賦予連接剛度，可以準(zhǔn)確計(jì)算點(diǎn)式連接航天器的模態(tài)。具體步驟為：

（1）靜剛度計(jì)算：根據(jù)航天器結(jié)構(gòu)，采用有限元分析軟件進(jìn)行接觸分析，選擇點(diǎn)式連接局部部位建立有限元接觸分析模型，計(jì)算其3個(gè)方向的力－位移曲線和力矩－轉(zhuǎn)角曲線，進(jìn)一步計(jì)算拉壓、彎曲和扭轉(zhuǎn)靜剛度曲線。

（2）等效剛度計(jì)算：根據(jù)本文2.2節(jié)所述的等效剛度計(jì)算方法，處理步驟1獲得的靜剛度曲線，分別計(jì)算軸向拉壓等效剛度和彎曲等效剛度。

（3）航天器模態(tài)計(jì)算：建立航天器有限元模態(tài)分析模型，將點(diǎn)式連接部位用廣義彈簧單元代替，將步驟2所得的等效剛度結(jié)果賦予廣義彈簧單元，計(jì)算航天器的模態(tài)。

但是連接剛度隨外激勵(lì)的變化具有非線性特性，通過商用有限元軟件對點(diǎn)式連接部位建立精細(xì)化有限元模型，能夠獲得點(diǎn)式連接的3條平動剛度曲線和3條轉(zhuǎn)動剛度曲線，大量研究表明，除了扭轉(zhuǎn)剛度曲線基本是線性的，其他連接剛度曲線都具有明顯的非線性特性。但廣義彈簧單元需要被賦予準(zhǔn)確的剛度值才能計(jì)算點(diǎn)式連接的模態(tài)，如何通過非線性剛度曲線計(jì)算廣義彈簧單元的剛度值（本文稱為等效剛度），成為點(diǎn)式連接航天器模態(tài)分析的關(guān)鍵。因此，下文將圍繞等效剛度計(jì)算進(jìn)行深入探討。

2.2 等效剛度分析

為了獲得等效剛度，這里采用非線性振動［14］和結(jié)構(gòu)動力學(xué)［15］的原理，計(jì)算非線性連接剛度曲線的等效剛度。圖1是某航天器的器箭艙段連接結(jié)構(gòu)示意圖，它采用螺栓－法蘭點(diǎn)式連接結(jié)構(gòu)。為了便于理解，以圖1坐標(biāo)系為例介紹連接等效剛度的計(jì)算方法。因?yàn)楹教炱饕话憔哂袑ΨQ性，剪切剛度Ky＝Kz（Ky是y方向的剪切剛度，Kz是z方向的剪切剛度），彎曲剛度Kθy＝Kθz（Kθy是y方向的彎曲剛度，Kθz是z方向的彎曲剛度，θ為轉(zhuǎn)角）。同時(shí)，點(diǎn)式連接航天器一般為長細(xì)結(jié)構(gòu)，剪切剛度Ky（Kz）對結(jié)構(gòu)模態(tài)的影響較小［17］，因此，本文將重點(diǎn)計(jì)算軸向拉壓等效剛度Kx、彎曲等效剛度Kθy（Kθz）和扭轉(zhuǎn)等效剛度Kθx，并假設(shè)剪切剛度為軸向拉壓剛度的10倍，即Ky（Kz）＝10 Kx。

圖1 點(diǎn)式器箭艙段連接結(jié)構(gòu)示意圖Fig.1 Separation nuts and joined parts

由于扭轉(zhuǎn)剛度曲線是線性的，可以直接讀取剛度值，而拉壓剛度曲線和彎曲剛度曲線具有明顯的非線性特性，須要計(jì)算等效剛度。

令K為點(diǎn)式連接非線性系統(tǒng)（軸向拉壓或彎曲）的等效剛度，由結(jié)構(gòu)動力學(xué)［15］可知

式中：T是非線性系統(tǒng)（軸向拉壓或彎曲）自由振蕩的運(yùn)動周期，m是外力作用點(diǎn)處的集中質(zhì)量。假設(shè)連接結(jié)構(gòu)為保守系統(tǒng)，采用非線性振動［14］中計(jì)算分段線性保守系統(tǒng)的方法計(jì)算等效剛度。

2.2.1 軸向拉壓等效剛度

大量研究表明，點(diǎn)式連接的拉伸和壓縮剛度通常存在較大差異，因此整體顯示出非線性特性［3］（參見3.1節(jié)）。忽略零點(diǎn)過渡段，點(diǎn)式連接在僅受拉或者受壓時(shí)均可簡化為獨(dú)立的線性保守系統(tǒng)。保守系統(tǒng)的動力學(xué)方程一般形式為

艙段連接系統(tǒng)的拉壓恢復(fù)力f（x）是位移x的分段函數(shù)，

式中：KxT和KxC分別是拉伸段和壓縮段的剛度，對應(yīng)的勢能和相軌跡方程分別為

式中：V（x）是保守系統(tǒng)的勢能，Em1和Em2分別是系統(tǒng)的拉伸段和壓縮段的機(jī)械能。由相軌跡方程（6）可知，當(dāng)速度y＝0時(shí)，得拉伸段和壓縮段的幅值分別為

2.2.2 彎曲等效剛度

前期試驗(yàn)顯示，點(diǎn)式連接的彎曲剛度具有分段線性特征（參見3.1節(jié)）。當(dāng)外載荷較小時(shí)，彎曲剛度較小；隨著載荷的增加，彎曲剛度增大。忽略過渡段，將彎矩作用下的點(diǎn)式連接簡化為分段線性系統(tǒng)，則對應(yīng)的彎曲恢復(fù)力f（θ）是轉(zhuǎn)角θ的分段線性函數(shù)：

式中：KI和KII分別是小剛度段和大剛度段線性系統(tǒng)的剛度，θ1是彎曲剛度改變時(shí)對應(yīng)的角度值。

一般地，在力學(xué)模態(tài)試驗(yàn)中，小剛度段范圍較窄，彎曲剛度很快從小剛度段過渡到大剛度段，所以彎曲等效剛度包括小剛度KI和大剛度KII兩部分，推導(dǎo)過程參考軸向拉壓等效周期的計(jì)算方法，彎曲方向運(yùn)動的周期為

將式（12）代入式（1）獲得彎曲等效剛度為

3 典型衛(wèi)星承力筒結(jié)構(gòu)模態(tài)分析算例

考慮到承力筒是國內(nèi)外航天器極為常見的主承載結(jié)構(gòu)，如DFH－3、DFH－4、SAST3000、SAST5000、EuroStar3000、AlphaBus等衛(wèi)星平臺。因此參照航天器的典型艙段連接結(jié)構(gòu)，設(shè)計(jì)了一個(gè)衛(wèi)星承力筒模型（見圖2），由上承力筒1、下承力筒2和正方形配重塊組成，通過若干個(gè)M6的螺栓將兩筒連為一體（螺栓數(shù)目分別取4，8，16，24，32，48個(gè)），螺栓數(shù)目的遞增代表了從點(diǎn)式連接過渡到普通艙段連接。模型尺寸如圖2所示，材料參數(shù)見表1。

圖2 衛(wèi)星承力筒幾何模型Fig.2 Assembled structure with bolted joints

表1 衛(wèi)星承力筒模型材料特性Table 1 Material properties of bolted joints

為了計(jì)算連接部位的等效剛度，取螺栓法蘭區(qū)域建立三維有限元分析模型，尺寸如圖3所示。外力作用點(diǎn)位于承力筒1上表面中點(diǎn)，與承力筒1上表面耦合約束，承力筒2的下表面固支約束。采用ABAQUS有限元軟件對模型進(jìn)行分析。有限元模型中，上下承力筒、法蘭和螺栓均采用8節(jié)點(diǎn)實(shí)體單元。螺栓與法蘭界面、上下法蘭界面之間建立接觸（如圖3的粗實(shí)線所示），接觸特性為法向硬接觸、切向庫倫摩擦（摩擦因數(shù)根據(jù)機(jī)械設(shè)計(jì)手冊取0.05）。在每組工況下，參考緊固件擰緊力矩標(biāo)準(zhǔn)，螺栓的預(yù)緊力取4000N。模型約有293 472個(gè)單元和364 928個(gè)節(jié)點(diǎn)，可以滿足計(jì)算收斂性的需求。

圖3 螺栓法蘭連接結(jié)構(gòu)有限元模型Fig.3 FEM bolted flange joints

3.1 衛(wèi)星承力筒連接剛度曲線

如圖4所示，拉力FxT施加在外力作用點(diǎn)處，通過ABAQUS軟件計(jì)算得到作用點(diǎn)的拉伸位移uxT，拉伸剛度KxT是拉力－拉伸位移曲線（圖4（a）虛線右側(cè)拉伸段）的斜率，即

同理，壓縮剛度KxC是壓力FxC和壓縮位移uxC曲線（圖4（a）虛線左側(cè)壓縮段）的斜率，即

圖4（b）是拉壓剛度曲線，虛線右端是拉伸剛度，虛線左端是壓縮剛度?？梢钥闯?，拉壓剛度具有明顯的非線性特性，即拉伸剛度遠(yuǎn)小于壓縮剛度［3］。同時(shí)，壓縮剛度與螺栓的數(shù)目基本無關(guān)，而拉伸剛度隨螺栓數(shù)目的增加而增加。

在外力作用點(diǎn)施加±y（±z）方向的彎矩和±x方向的扭矩，獲得相應(yīng)的轉(zhuǎn)角，參考式（14）或式（15）計(jì)算彎曲剛度Kθy（Kθz）和扭轉(zhuǎn)剛度Kθx，計(jì)算結(jié)果見圖5和圖6所示。圖5（a）是彎矩My和轉(zhuǎn)角θy的曲線，圖5（b）是彎曲剛度Kθy曲線，虛線左右分別是－θy轉(zhuǎn)角和＋θy轉(zhuǎn)角，彎曲剛度具有非線性特性。圖6（a）是扭矩Mx和轉(zhuǎn)角θx的曲線，圖6（b）是扭轉(zhuǎn)剛度Kθx曲線，虛線左右分別是－θx轉(zhuǎn)角和＋θx轉(zhuǎn)角，扭轉(zhuǎn)剛度是線性的。與拉伸剛度類似，彎曲和扭轉(zhuǎn)剛度隨螺栓數(shù)目的增加而增加。

圖4 軸向拉壓靜剛度曲線Fig.4 Tension and compression stiffness curves

圖5 彎曲靜剛度曲線Fig.5 Bending stiffness curves

圖6 扭轉(zhuǎn)靜剛度曲線Fig.6 Torsional stiffness curves

3.2 衛(wèi)星承力筒等效剛度計(jì)算

圖4是拉壓剛度曲線，忽略零點(diǎn)的過渡區(qū)域，拉伸剛度和壓縮剛度近似看作恒定的，可采用式（10）計(jì)算軸向等效剛度。圖5中彎曲剛度曲線具有兩條線性段，當(dāng)外載荷較小時(shí)，彎曲剛度幅值較小且近似看作恒定的；隨著載荷增加，彎曲剛度幅值成階梯狀增加，穩(wěn)定后的結(jié)果也近似看作線性的；忽略過渡區(qū)域，將彎矩作用下的點(diǎn)式連接簡化為分段線性系統(tǒng)，采用式（13）計(jì)算彎曲等效剛度。切向等效剛度取軸向等效剛度的10倍，扭轉(zhuǎn)剛度直接從圖6所示扭轉(zhuǎn)剛度曲線讀取，等效剛度具體數(shù)值見表2。

表2 等效剛度Table 2 Equivalent joint stiffness

從圖4、圖5和圖6可以看出，連接結(jié)構(gòu)的連接剛度隨螺栓數(shù)目的增加而增加，表2也列出不同螺栓數(shù)目的等效剛度取值，為了進(jìn)一步分析螺栓數(shù)目對連接剛度的影響，對等效剛度的取值歸一化得

圖7 剛度變化趨勢曲線Fig.7 Stiffness trend curves

特別是，點(diǎn)式連接呈現(xiàn)明顯的剛度衰減現(xiàn)象，相比于48個(gè)螺栓的等效剛度，4個(gè)螺栓的衰減幅值高達(dá)85%。顯然，若不在建模過程中加以考慮點(diǎn)式連接的剛度衰減現(xiàn)象，將直接影響動力學(xué)分析結(jié)果的正確性。

3.3 衛(wèi)星承力筒模態(tài)分析與試驗(yàn)驗(yàn)證

圖8是衛(wèi)星承力筒的有限元簡化模型，采用MSC／Nastran有限元軟件建模并分析模態(tài)。與圖2中幾何模型相同，結(jié)構(gòu)由承力筒1、承力筒2和質(zhì)量為22.92kg的配重組成。其中螺栓法蘭連接結(jié)構(gòu)（見圖3）由廣義彈簧單元（Nastran中即為Bush單元，其剛度矩陣是由6個(gè)自由度的對角陣組成［16］）代替。Bush單元的兩端與承力筒1和2通過Mpc約束連接，Bush單元的剛度取表2中的等效剛度（根據(jù)螺栓數(shù)目的不同計(jì)算6組工況）。承力筒采用4節(jié)點(diǎn)殼單元，配重采用8節(jié)點(diǎn)實(shí)體單元，殼單元與實(shí)體單元通過Mpc約束連接［18］。承力筒2下表面48個(gè)螺栓孔附近的區(qū)域固支約束，材料參數(shù)與表1相同。

圖9是螺栓法蘭連接結(jié)構(gòu)的模態(tài)試驗(yàn)示意圖，試驗(yàn)系統(tǒng)包含試驗(yàn)件、加速度傳感器（靈敏度是40mV／g）、力錘、LMS公司數(shù)據(jù)采集設(shè)備、LMS公司軟件分析系統(tǒng)等。承力筒與配重的安裝方式與圖2相同。增加加速度傳感器數(shù)量可以提高測試精度，但由于傳感器質(zhì)量增加對模態(tài)有影響，這里采用13個(gè)傳感器，1個(gè)布置在配重頂部中心位置，另12個(gè)布置在承力筒1的上法蘭根部、腹部和下法蘭根部，均沿周向呈90°分布。力錘的錘擊點(diǎn)分為A、B兩處，點(diǎn)A在配重上表面中點(diǎn)，點(diǎn)B在配重側(cè)面中點(diǎn)，每個(gè)測點(diǎn)用力錘敲擊5次取樣，通過LMS軟件求均值，根據(jù)螺栓數(shù)目的不同測試6組工況。

圖8 有限元模型示意圖Fig.8 FEM mode analytical specimen

圖9 模態(tài)試驗(yàn)示意圖Fig.9 Modal test specimen

表3列出了6組工況下試驗(yàn)和仿真的模態(tài)頻率對比結(jié)果，航天器模態(tài)分析重點(diǎn)關(guān)注橫向一階彎曲模態(tài)、縱向一階模態(tài)、橫向二階彎曲模態(tài)和扭轉(zhuǎn)一階模態(tài)，試驗(yàn)和仿真結(jié)果顯示，頻率誤差基本在15%以內(nèi)，仿真結(jié)果與試驗(yàn)數(shù)據(jù)吻合很好，等效剛度法適用于艙段連接結(jié)構(gòu)的模態(tài)分析。尤其對于點(diǎn)式連接航天器，因?yàn)椴捎脗鹘y(tǒng)建模方法無法獲得精確的模態(tài)，而等效剛度模態(tài)分析方法解決了連接剛度非線性對模態(tài)分析精度的影響。

表3 仿真與試驗(yàn)?zāi)B(tài)頻率對比Table 3 Frequency obtained from simulation and test Hz

從表3可以看出模態(tài)頻率隨螺栓數(shù)目的增加而增加，但增加的幅度隨螺栓數(shù)目的增加而減小。對比點(diǎn)式連接（4個(gè)螺栓）和普通艙段連接（48個(gè)螺栓），頻率下降高達(dá)30%，說明點(diǎn)式連接的螺栓數(shù)量對結(jié)構(gòu)模態(tài)的影響很大，考慮螺栓數(shù)目變化引起連接剛度的變化，準(zhǔn)確分析點(diǎn)式連接航天器的模態(tài)是很有必要的。

雖然螺栓數(shù)目有所不同，但不同工況下的振型基本相同，圖10列出了8個(gè)螺栓仿真和試驗(yàn)振型的對比，證明仿真結(jié)果的正確性。

圖10 仿真與試驗(yàn)振型對比Fig.10 Mode of vibration obtained from simulation and test

3.4 等效剛度模態(tài)法與傳統(tǒng)分析方法對比

以4個(gè)螺栓的衛(wèi)星承力筒點(diǎn)式連接計(jì)算模型為例，分別采用全固連建模方法和連接節(jié)點(diǎn)約束方法計(jì)算點(diǎn)式連接模型的模態(tài)，并與等效剛度模態(tài)分析方法和模態(tài)試驗(yàn)結(jié)果比較。全固連方法將上下承力筒連接部分節(jié)點(diǎn)融合，使其成為一個(gè)線性結(jié)構(gòu)；連接節(jié)點(diǎn)約束方法僅將上下承力筒對應(yīng)節(jié)點(diǎn)通過MPC約束，這里僅建立4個(gè)約束（見表4）。

表4 模態(tài)分析結(jié)果Table 4 Results of modal analyses Hz

從表4模態(tài)分析結(jié)果的比較來看，全固連的模態(tài)計(jì)算結(jié)果與試驗(yàn)結(jié)果誤差高達(dá)64%，連接節(jié)點(diǎn)約束方法無法計(jì)算出真實(shí)的模態(tài)頻率，誤差高達(dá)99%。因此，等效剛度模態(tài)分析方法對于連接結(jié)構(gòu)的計(jì)算精度與其他方法的相比具有顯著的優(yōu)勢，同時(shí)也說明傳統(tǒng)的分析方法過于粗放和保守，不能滿足現(xiàn)代航天器精細(xì)化設(shè)計(jì)計(jì)算的需求。

4 結(jié)論

點(diǎn)式連接已廣泛應(yīng)用于大型航天器艙段及部件級結(jié)構(gòu)，本文針對模態(tài)分析中存在的技術(shù)難題，提出了一種采用等效剛度計(jì)算點(diǎn)式連接航天器模態(tài)的方法。運(yùn)用非線性振動和結(jié)構(gòu)動力學(xué)相關(guān)理論，計(jì)算出非線性連接結(jié)構(gòu)的等效剛度，結(jié)合廣義彈簧單元計(jì)算點(diǎn)式連接航天器的模態(tài)。主要結(jié)論有：

（1）采用廣義彈簧單元計(jì)算的模態(tài)結(jié)果與模態(tài)試驗(yàn)數(shù)據(jù)對比，航天器的橫向一階、二階彎曲，縱向一階，一階扭轉(zhuǎn)模態(tài)誤差小于15%，基本上在10%左右，證明了該方法計(jì)算點(diǎn)式連接航天器的有效性。

（2）本文定量地證明了點(diǎn)式連接的連接剛度隨著螺栓數(shù)量的減少將呈現(xiàn)明顯衰減，衰減幅度高達(dá)85%。

（3）仿真與試驗(yàn)結(jié)果均顯示，點(diǎn)式連接與普通艙段連接相比，航天器整體的固有頻率下降幅度高達(dá)30%，說明航天器設(shè)計(jì)中點(diǎn)式連接的靜／動力學(xué)分析應(yīng)予以高度重視。

（4）等效剛度模態(tài)分析方法可以準(zhǔn)確分析航天器的模態(tài)，尤其針對連接點(diǎn)數(shù)目較少的點(diǎn)式連接航天器，解決了傳統(tǒng)的全固連建模方法或連接節(jié)點(diǎn)約束方法計(jì)算不準(zhǔn)確的問題。

本文所提出的等效剛度方法對于提高點(diǎn)式連接的航天器模態(tài)計(jì)算精度具有參考意義。

（References）

［1］趙超，周軍，周鳳岐.大型復(fù)合航天器的建模與分散控制技術(shù)［J］.飛行力學(xué)1998，16（3）：22－27 Zhao Chao，Zhou Jun，Zhou Fengqi.Modeling and decentralized control techniques for large compound spacecraft［J］.Flight Dynamics，1998，16（3）：22－27（in Chinese）

［2］鄭曉亞，徐超，王燾，等.螺栓－法蘭連接結(jié)構(gòu)非線性優(yōu)化設(shè)計(jì)方法研究綜述［J］.強(qiáng)度與環(huán)境，2008，35（3）：7－13 Zheng Xiaoya，Xu Chao，Wang Tao，et al.Study on nonlinear optimization design for bolted flanged connection［J］.Structure ＆Environment Engineering，2008，35（3）：7－13（in Chinese）

［3］欒宇，劉松，關(guān)振群，等.小變形下螺栓法蘭連接結(jié)構(gòu)的靜剛度［J］.強(qiáng)度與環(huán)境，2011，38（3）：29－35 Luan Yu，Liu Song，Guan Zhenqun，et al.The nonlinear character of static stiffness in the structure with bolted flanged connection under small deformation［J］.Structure ＆Environment Engineering，2011，38（3）：29－35（in Chinese）

［4］鄭曉亞，張鐸，姜晉慶.連接剛度對導(dǎo)彈固有特性的影響［J］.彈箭與制導(dǎo)學(xué)報(bào)，2005，25（4）：667－669 Zheng Xiaoya，Zhang Duo，Jiang Jinqing.The joint stiffness effect upon missile’s normal modes［J］.Journal of Projectiles，Rockets，Missiles and Guidance，2005，25（4）：667－669（in Chinese）

［5］Zhang O，Poirier J A.New analytical model of bolted joints［J］.Journal of Mechanical Design，2004，126（4）：721－728

［6］Wang C，Zhang D Y，Zhu X B，et al.Study on the stiffness loss and the dynamic influence on rotor system of the bolted flange joint［C］／／ASME Turbo Expo 2014：Turbine Technical Conference and Exposition.New York：ASME，2014

［7］程來，宋言明，李賀，等.基于等效接觸區(qū)域的螺栓連接結(jié)構(gòu)模態(tài)分析［J］.機(jī)械設(shè)計(jì)與研究，2012，28（2）：44－47 Cheng Lai，Song Yanming，Li He，et al.Modal analysis based on the equivalent contact area of the structure with bolt connections［J］.Machine Design ＆Research，2012，28（2）：44－47（in Chinese）

［8］薛闖，賈建軍，舒嶸，等.裝配體結(jié)構(gòu)有限元分析中的螺釘連接模型［J］.科學(xué)技術(shù)與工程，2006，6（7）：825－828 Xue Chuang，Jia Jianjun，Shu Rong，et al.Bolted joint models of finite element assemblies［J］.Science Technology and Engineering，2006，6（7）：825－828（in Chinese）

［9］Segalman D J.A four－parameter Iwan model for laptype joints［J］.Journal of Applied Mechanics，2005，72（5）：752－760

［10］Mayer M H，Gaul L.Segment－to－segment contact elements for modelling joint interfaces in finite element analysis［J］，Mechanical System Signal Processing，2007，21（2）：724－774

［11］石亦平，周玉蓉.ABAQUS有限元分析實(shí)例詳解［M］.北京：機(jī)械工業(yè)出版社，2006 Shi Yiping，Zhou Yurong.Detailed examples of ABAQUS finite element analysis［M］.Beijing：China Machine Press，2006（in Chinese）

［12］Kim J，Yoon J C，Kang B S.Finite element analysis and modeling of structure with bolted joints［J］.Applied Mathematical Modelling，2007，31（5）：895－911

［13］Yorgun C，Dalcl S，Altay G A.Finite element modeling of bolted steel connections designed by double channel［J］.Computers ＆Structures，2004，82（29）：2563－2571

［14］劉延柱，陳立群.非線性振動［M］.北京：高等教育出版社，2001：19－27 Liu Yanzhu，Chen Liqun.Nonlinear vibrations［M］.Beijing：Higher Education Press，2001：19－27

［15］余旭東.飛行器結(jié)構(gòu)動力學(xué)［M］.西安：西北工業(yè)大學(xué)出版社，2012 Yu Xudong.Aircraft structure dynamics［M］.Xi’an：Northwestern Polytechnical University Press，2012（in Chinese）

［16］鄒元杰.Bush有限單元原理及其在航天器結(jié)構(gòu)建模中的應(yīng)用［J］.航天器工程，2010，19（1）：99－105 Zou Yuanjie.Bush finite element principle and its application in building structural models for spacecraft［J］.Spacecraft Engineering，2010，19（1）：99－105（in Chinese）

［17］茍文選.材料力學(xué)［M］.西安：西北工業(yè)大學(xué)出版社，2000：183－185 Gou Wenxuan.Mechanics of materials［M］.Xi’an：Northwestern Polytechnical University Press，2012（in Chinese）

［18］龍凱.Patran2010與Nastran2010有限元分析從入門到精通［M］.北京：機(jī)械工業(yè)出版社，2010 Long Kai.Rudiments and mastery of Patran2010and Nastran2010FEA［M］.Beijing：China Machine Press，2010（in Chinese）

（編輯：張小琳）

Modal Analysis of Spacecraft with Dotted Joints Based on Equivalent Stiffness

SHI Chuanqian1ZHOU Xubin1ZHANG Yongtao1，2DU Dong1LI Hao1
（1Shanghai Institute of Satellite and Engineering，Shanghai 200240，China）（2Shanghai Lizheng Satellite Application Technology Co.，Ltd.，Shanghai 200240，China）

To include the interaction effect in modal analysis，a novel model analysis method employing the equivalent stiffness of dotted joints is proposed in this paper.Firstly，a full－size finite element model of a dotted joint for a spacecraft structure is established.The nonlinear connection stiffness of dotted joint corresponding to different bolt quantities and loads are predicted.A novel method based on nonlinear dynamic theory is proposed to calculate the equivalent stiffness of dotted joints.Then，a simplified FE model of the spacecraft structure is established by employing the generalized spring element and the calculated equivalent stiffness.The modes of the spacecraft structure are predicted by the simplified FE model and are also measured in experiments.The predicted modes coincide quite well with the experimental results.The proposed modal analysis method is proven to be accurate and will facilitate the dynamic analysis of spacecraft structure with dotted joints.

spacecraft；dotted joints；mode；equivalent stiffness；nonlinearity

V414.1

：ADOI：10.3969／j.issn.1673－8748.2016.01.005

2015－10－30；

：2015－12－28

國家重大航天工程

石川千，女，碩士研究生，從事衛(wèi)星結(jié)構(gòu)設(shè)計(jì)和結(jié)構(gòu)非線性動力學(xué)研究工作。Email：shichuanqian＠163.com。