液體運載火箭一維縱橫扭一體化建模技術(shù)

唐玉花，狄文斌，劉靖華

(上海宇航系統(tǒng)工程研究所，上海201109)

液體運載火箭一維縱橫扭一體化建模技術(shù)

唐玉花，狄文斌，劉靖華

(上海宇航系統(tǒng)工程研究所，上海201109)

提出一種新的液體運載火箭一維縱橫扭一體化模型建模方法，在傳統(tǒng)的橫向和扭轉(zhuǎn)振動一維梁模型基礎(chǔ)上，融入了縱向彈簧－質(zhì)量模型，并引入了多液體質(zhì)量塊縱向彈簧－質(zhì)量模型。將該建模方法應(yīng)用于某型運載火箭液體助推器，通過與試驗值的比較證明:該一維縱橫扭一體化模型可同時計算橫向、縱向和扭轉(zhuǎn)模態(tài)及耦合模態(tài)，橫向、縱向和扭轉(zhuǎn)頻率的計算精度較高，滿足工程應(yīng)用要求。同時，增加液體質(zhì)量塊的數(shù)量可以提高縱向頻率的計算精度，然而當(dāng)貯箱內(nèi)液體較少時，增加液體質(zhì)量塊的數(shù)量并不能改善計算結(jié)果的精度。

液體運載火箭;動特性;一維縱橫扭一體化模型;多液體質(zhì)量塊

0 引 言

液體發(fā)動機(jī)比沖高，能多次啟動，便于調(diào)節(jié)推力，因而大型運載火箭以液體火箭為主。液體運載火箭的動力學(xué)特性是進(jìn)行姿態(tài)控制系統(tǒng)設(shè)計和星箭耦合分析及POGO分析的基礎(chǔ)，同時也是載荷設(shè)計的依據(jù)。

隨著國內(nèi)液體運載火箭技術(shù)的快速發(fā)展，特別是帶捆綁的大型運載火箭，由于捆綁火箭與芯級火箭之間存在橫向、縱向和扭轉(zhuǎn)模態(tài)的耦合問題［1］，以往橫向扭轉(zhuǎn)與縱向分別為兩套模型已無法滿足計算要求。為研究捆綁火箭全箭動力學(xué)特性，王建民等［2］對捆綁火箭建立了三維有限元模型，并給出了捆綁火箭橫向、縱向和扭轉(zhuǎn)模態(tài)相互耦合的特征。王建民等［3］認(rèn)為一維梁模型無法使用虛質(zhì)量等三維液體建模方法，無法反映火箭結(jié)構(gòu)的加筋、開口和錐殼等剛度特征，無法準(zhǔn)確模擬捆綁連接機(jī)構(gòu)的剛度等，因而提出了運載火箭三維建模方法。然而，其建模工作量極大。采用虛質(zhì)量法模擬貯箱內(nèi)液體推進(jìn)劑，雖然更貼近液體的實際動力學(xué)行為，精度較高［4］，但虛質(zhì)量建模方法會形成一個貯箱與液體接觸面自由度規(guī)模的滿質(zhì)量矩陣，計算量過大，工程應(yīng)用困難。

液體推進(jìn)劑在液體運載火箭全箭起飛質(zhì)量中所占比例較大，約為全箭起飛質(zhì)量的90%，對液體推進(jìn)劑的模擬是火箭動力學(xué)模型建模的關(guān)鍵。由于液體推進(jìn)劑對橫向和縱向的作用不同，國內(nèi)在進(jìn)行火箭動力學(xué)分析時采用的一維模型，其中橫向和扭轉(zhuǎn)振動特性分析采用梁模型，即無質(zhì)量的彈性梁和集中質(zhì)量模型，縱向振動特性分析采用彈簧－質(zhì)量模型，目前還沒有較為成熟的將橫向、縱向和扭轉(zhuǎn)合并為一套一維模型的理論。潘忠文［5］在進(jìn)行運載火箭動力學(xué)模型建模計算時，貯箱采用梁模型，貯箱內(nèi)液體推進(jìn)劑用集中質(zhì)量單元模擬，得到的低階橫向模態(tài)預(yù)示精度滿足工程設(shè)計要求。文獻(xiàn)［6］采用的運載火箭縱橫扭一體化梁模型中，對液體推進(jìn)劑同樣用集中質(zhì)量單元模擬。之后，潘忠文等［7］提出了基于梁模型的火箭橫向、縱向和扭轉(zhuǎn)一體化模型建模方法，對液體推進(jìn)劑建立了以附加質(zhì)量表示的耦合質(zhì)量矩陣。文獻(xiàn)［8］詳細(xì)闡述了液體推進(jìn)劑集中質(zhì)量方法和耦合質(zhì)量方法，認(rèn)為耦合質(zhì)量方法反映了液體推進(jìn)劑在橫向、縱向和扭轉(zhuǎn)方向的不同作用。然而文獻(xiàn)［7－8］僅考慮了液體推進(jìn)劑對貯箱橫向和縱向的不同質(zhì)量效應(yīng)。文獻(xiàn)［3］對液體推進(jìn)劑的建模方法進(jìn)行了總結(jié)，并提出了一種液體桿單元建模方法，此方法與耦合質(zhì)量方法具有同等效果。林宏等［9］在進(jìn)行運載火箭動力學(xué)模型修正時，對液體推進(jìn)劑同樣采用耦合質(zhì)量方法模擬，由文獻(xiàn)［9］可以看出縱向頻率計算精度欠佳。文獻(xiàn)［10］對運載火箭結(jié)構(gòu)動力學(xué)模擬技術(shù)研究進(jìn)展進(jìn)行了綜述，認(rèn)為充液貯箱流固耦合模型有待進(jìn)一步完善。

傳統(tǒng)的縱向彈簧－質(zhì)量模型［11－12］為液體運載火箭縱向動力學(xué)特性分析提供了較好的解決途徑，最近狄文斌等［13］開展了縱向多液體質(zhì)量塊彈簧－質(zhì)量模型的應(yīng)用研究，指出多液體質(zhì)量塊簡化模型計算得到的高階縱向振動頻率精度明顯優(yōu)于單液體質(zhì)量塊模型，然而隨著貯箱內(nèi)推進(jìn)劑液位的降低，彈簧－質(zhì)量模型的精度也降低。

本文在液體運載火箭橫向和扭轉(zhuǎn)振動特性分析梁模型基礎(chǔ)上，將貯箱內(nèi)液體推進(jìn)劑橫向效應(yīng)按照集中質(zhì)量的形式建模，而縱向效應(yīng)采用彈簧－質(zhì)量的形式建模，并考慮貯箱底為彈性的，將用于橫向和扭轉(zhuǎn)計算的梁模型與用于縱向計算的彈簧－質(zhì)量模型融合，并且引入多液體質(zhì)量塊縱向彈簧－質(zhì)量模型，提高了縱向頻率的計算精度。將該一維縱橫扭一體化建模方法應(yīng)用于某型運載火箭液體助推器的動特性分析，得到結(jié)果與試驗結(jié)果吻合較好。

1 液體運載火箭梁模型等效剛度

目前提供姿態(tài)控制系統(tǒng)設(shè)計、星箭耦合分析和載荷設(shè)計等使用的火箭彈性振動參數(shù)，計算時采用的火箭動力學(xué)分析模型，一般除衛(wèi)星、適配器和支承艙等為三維有限元模型外，其余艙段均為梁模型，發(fā)動機(jī)可采用集中質(zhì)量模型，或者機(jī)架采用三維有限元模型，機(jī)架以外的部分僅模擬質(zhì)量。

艙段一般分為蒙皮桁條結(jié)構(gòu)和化銑結(jié)構(gòu)。蒙皮桁條結(jié)構(gòu)在進(jìn)行梁模型建模時，根據(jù)桁條與蒙皮彈性模量的關(guān)系得到某一截面上桁條總的等效面積，則該截面上總的等效厚度為

式中:t為該段艙體等效的截面總厚度，tmp為該段艙體蒙皮的厚度，tht為該段艙體桁條總的等效截面厚度，D為該段艙體的平均直徑，Ehti、Ahti分別為第i種桁條的彈性模量及面積和，Emp為蒙皮的彈性模量。

梁單元的截面剛度特性主要包含梁的材料屬性、梁截面積A、抗彎截面慣性矩Jx、抗扭截面慣性矩Jn、截面抗剪系數(shù)K，計算式如下:

式中:Amp=π(D－tmp)tmp，為蒙皮截面積。

化銑結(jié)構(gòu)的截面剛度主要由結(jié)構(gòu)材料屬性、艙體的平均直徑以及截面等效厚度t確定，截面等效厚度t按照不同的化銑網(wǎng)格形式可采用體積或質(zhì)量等效的方式獲得，梁截面積A、抗彎截面慣性矩Jx、抗扭截面慣性矩Jn、截面抗剪系數(shù)K的計算式如下:

按照上述方法確定了火箭結(jié)構(gòu)的等效剛度模型。

2 縱向彈簧－質(zhì)量模型

對于如圖1所示的充液貯箱，貯箱筒段為圓柱形，箱筒段壁厚為h，箱底為長半軸為a，短半軸為b的橢球底，貯箱內(nèi)推進(jìn)劑液面高度為L。假設(shè)液體無黏，且不可壓縮，并假設(shè)作用于貯箱和液體的外力軸對稱。

單液體質(zhì)量塊縱向彈簧－質(zhì)量模型中，自由度1表示箱筒段與箱后底結(jié)合處軸向位移，自由度2表示貯箱內(nèi)液體質(zhì)心處軸向位移，自由度3表示自由液面處貯箱的軸向位移。單液體質(zhì)量塊柔度系數(shù)矩陣［12］為

柔度系數(shù)矩陣d的逆即為剛度矩陣k:

由此得到圖2中所示單液體質(zhì)量塊縱向彈簧－質(zhì)量模型各自由度之間的彈簧剛度為

雙液體質(zhì)量塊貯箱見圖3，雙液體質(zhì)量塊縱向彈簧－質(zhì)量模型見圖4。相比單液體質(zhì)量塊模型，將貯箱內(nèi)液體分為兩個自由度，自由度2為下段液體質(zhì)量塊質(zhì)心處軸向位移，自由度3為上段液體質(zhì)量塊質(zhì)心處軸向位移，自由度1仍為箱筒段與箱后底結(jié)合處軸向位移，自由度4為自由液面處貯箱的軸向位移。兩個液體質(zhì)量塊根據(jù)不同的貯箱壁厚進(jìn)行分配，如果壁厚一致可根據(jù)液面高度平均分配。雙液體質(zhì)量塊柔度系數(shù)矩陣為

雙液體質(zhì)量塊縱向彈簧－質(zhì)量模型剛度矩陣k:

圖4所示雙液體質(zhì)量塊縱向彈簧－質(zhì)量模型各自由度之間的彈簧剛度為

根據(jù)文獻(xiàn)［13］，三液體質(zhì)量塊及以上的多液體質(zhì)量塊縱向彈簧－質(zhì)量模型同樣可擴(kuò)展得到。

3 質(zhì)量模型建模

本文采用商用有限元軟件MD.Nastran進(jìn)行結(jié)構(gòu)動力學(xué)問題的求解［14］。

MD.Nastran軟件提供了兩種集中質(zhì)量單元，即CONM1單元和CONM2單元，卡片分別見圖5和圖6，CONM1單元在節(jié)點上定義6×6的質(zhì)量矩陣，可分別定義局部坐標(biāo)系下三個方向的質(zhì)量，CONM2單元認(rèn)為節(jié)點在局部坐標(biāo)系三個方向上的質(zhì)量是一致的。

為實現(xiàn)貯箱內(nèi)液體推進(jìn)劑橫向質(zhì)量和縱向質(zhì)量的分開，對貯箱箱體節(jié)點和液體質(zhì)量塊節(jié)點均采用CONM1單元模擬。自由液面以上的貯箱箱體節(jié)點僅有結(jié)構(gòu)質(zhì)量，因而三個方向質(zhì)量相等，自由液面以下貯箱箱體節(jié)點橫向兩個方向質(zhì)量為分配到該節(jié)點的結(jié)構(gòu)質(zhì)量與液體質(zhì)量之和，而縱向質(zhì)量僅為結(jié)構(gòu)質(zhì)量，即帶液箱體節(jié)點縱向不含液體質(zhì)量。將貯箱內(nèi)液體推進(jìn)劑的縱向質(zhì)量按照采用的液體質(zhì)量塊數(shù)量和相應(yīng)的分配關(guān)系，分配到液體質(zhì)量塊節(jié)點的縱向質(zhì)量上，液體質(zhì)量塊節(jié)點的橫向質(zhì)量為0，由于不考慮液體的轉(zhuǎn)動慣量，因而液體質(zhì)量塊節(jié)點的CONM1單元無轉(zhuǎn)動慣量數(shù)據(jù)。

無液艙段節(jié)點的質(zhì)量采用CONM2單元模擬。

4 梁模型與縱向彈簧－質(zhì)量模型合并

對于空箱狀態(tài)，即貯箱內(nèi)無液體推進(jìn)劑的情況，梁模型可以同時計算橫向、扭轉(zhuǎn)和縱向模態(tài)，即空箱狀態(tài)的梁模型同時包含空箱狀態(tài)火箭的橫向、扭轉(zhuǎn)及縱向的剛度和質(zhì)量。液體推進(jìn)劑的加入，使得貯箱箱體節(jié)點質(zhì)量形式改變，并且引入了液體質(zhì)量塊節(jié)點，以及液體質(zhì)量塊節(jié)點之間、液體質(zhì)量塊節(jié)點與箱體自由液面處節(jié)點及箱筒段與箱后底結(jié)合處節(jié)點之間的縱向彈簧。這里不考慮液體推進(jìn)劑對液面以下箱筒段結(jié)構(gòu)剛度的影響，即認(rèn)為圖2中的和圖4中的與箱筒段結(jié)構(gòu)剛度一致。

不同的計算秒狀態(tài)，貯箱內(nèi)自由液面位置不同，因而首先需要確定計算秒狀態(tài)下貯箱的自由液面位置，并在梁模型中自由液面位置處給出一個節(jié)點，如果這個節(jié)點已有，梁模型不做修改，如果這個節(jié)點是新增的，由于引入了這個節(jié)點，需對這部分梁模型進(jìn)行適當(dāng)?shù)母隆?/p>

對圖7所示帶液貯箱及上、下各一個無液艙段的結(jié)構(gòu)，進(jìn)行一維縱橫扭一體化模型的建模，貯箱內(nèi)液體自由液面在節(jié)點3處，將貯箱內(nèi)液體分為兩個質(zhì)量塊，合并后的模型見圖8。圖8中，左側(cè)為縱向彈簧－質(zhì)量模型，由于不考慮液體對液面以下箱筒段結(jié)構(gòu)剛度的影響，因而箱體自由液面處節(jié)點和箱筒段與箱后底結(jié)合處節(jié)點之間的剛度即為梁模型中這部分結(jié)構(gòu)的剛度。圖8中彈簧剛度根據(jù)前述方法得到，將貯箱內(nèi)液體推進(jìn)劑按照相應(yīng)比例分配到節(jié)點11和12的CONM1單元縱向質(zhì)量上;右側(cè)為修改后的梁模型，相比原梁模型，m3～m8由CONM2單元改為CONM1單元，并對橫向質(zhì)量賦分到該節(jié)點的結(jié)構(gòu)質(zhì)量與液體質(zhì)量之和，而縱向質(zhì)量僅賦分到該節(jié)點的結(jié)構(gòu)質(zhì)量。

通過上述方法，得到同時包含橫向、縱向和扭轉(zhuǎn)振動計算要素的一維模型，實現(xiàn)液體運載火箭一維縱橫扭一體化動力學(xué)模型的建模。同時，可以看出，本文提出的一維縱橫扭一體化建模方法主要是解決帶液貯箱的橫向、扭轉(zhuǎn)和縱向模型的合并，對于無液艙段采用一維梁模型或者三維模型均可。

5 應(yīng)用實例

圖9所示為某型運載火箭液體助推器，助推器主體為兩個貯箱，滿狀態(tài)液體推進(jìn)劑質(zhì)量占助推器總質(zhì)量的90%以上。為獲取助推器自身的模態(tài)特性，開展了助推器自由－自由邊界狀態(tài)的模態(tài)試驗。助推器采用彈性支承懸吊以模擬自由－自由邊界條件，分別獲取了2 s、56 s、73 s和150 s帶液狀態(tài)的模態(tài)參數(shù)。

對模態(tài)試驗狀態(tài)的助推器采用本文方法進(jìn)行一維縱橫扭一體化模型的建模，建立的模型及坐標(biāo)系定義見圖10，除發(fā)動機(jī)為三維模型外，其余艙段均采用一維模型，各秒狀態(tài)貯箱內(nèi)液體推進(jìn)劑質(zhì)量為模態(tài)試驗中貯箱實際加注量的估算值，并采用集中質(zhì)量點模擬試驗中懸吊工裝的質(zhì)量特性。將自由－自由邊界條件下助推器30 Hz以內(nèi)的頻率計算結(jié)果與試驗結(jié)果比較，見表1。表2給出了采用不同數(shù)量液體質(zhì)量塊助推器縱向頻率計算結(jié)果，其中兩個貯箱采用相同數(shù)量的液體質(zhì)量塊，如三液體質(zhì)量塊表示上箱和下箱均采用三個液體質(zhì)量塊。

表1 液體助推器頻率結(jié)果比較Table 1 Frequency comparison of the liquid-propellant roll booster

表2 不同數(shù)量液體質(zhì)量塊縱向頻率結(jié)果比較Table 2 Longitudinal frequency comparison with varying liquid mass number

從表1～2可以看出:

1)一維縱橫扭一體化模型同時獲取了橫向、縱向和扭轉(zhuǎn)模態(tài)，計算得到的一階彎曲頻率與試驗值差異小于5%，二階和三階彎曲頻率與試驗值差異小于10%，一階扭轉(zhuǎn)頻率與試驗值差異小于10%，一階縱向頻率與試驗值差異小于5%，二階縱向頻率與試驗值差異小于10%，提供了較高的計算精度，滿足工程應(yīng)用要求。

2)當(dāng)貯箱內(nèi)液體推進(jìn)劑較多時，增加液體質(zhì)量塊數(shù)量可以提高縱向頻率的計算精度。然而當(dāng)貯箱內(nèi)液體較少時，增加液體質(zhì)量塊的數(shù)量并不能改善計算結(jié)果的精度，此時，由于貯箱內(nèi)液體較少，縱向振動時不容易形成較為整體的運動。

6 結(jié)論

本文在傳統(tǒng)的液體運載火箭橫向和扭轉(zhuǎn)振動一維梁模型基礎(chǔ)上，融入縱向彈簧－質(zhì)量模型，形成了可同時獲得橫向、縱向和扭轉(zhuǎn)模態(tài)及三者之間耦合模態(tài)的一維縱橫扭一體化模型。將該建模方法應(yīng)用于某型運載火箭液體助推器的模態(tài)計算，得到的橫向、縱向和扭轉(zhuǎn)頻率結(jié)果與試驗值差異較小，計算精度較高，滿足工程應(yīng)用要求。通過不同數(shù)量液體質(zhì)量塊縱向頻率的比較發(fā)現(xiàn)，增加液體質(zhì)量塊的數(shù)量可以提高縱向頻率的計算精度，然而當(dāng)貯箱內(nèi)液體較少時，增加液體質(zhì)量塊的數(shù)量并不能改善計算結(jié)果的精度。

［1］ 李家文，李道奎，劉軍虎，等.大型捆綁火箭模態(tài)縱向分量對姿控系統(tǒng)的影響分析［J］.宇航學(xué)報，2010，31(5):1390－1397.［Li Jia-wen，Li Dao-kui，Liu Jun-hu，et al.Analysis effect of longitudinal component of a large strap-on launch vehicle's modes on attitude control system［J］.Journal of Astronautics，2010，31(5):1390－1397.］

［2］ 王建民，榮克林，馮穎川，等.捆綁火箭全箭動力學(xué)特性研究［J］.宇航學(xué)報，2009，30(3):821－826.［Wang Jian-min，Rong Ke-lin，F(xiàn)eng Ying-chuan，et al.The research of dynamic characteristics for the strap-on launch vehicle［J］.Journal of Astronautics，2009，30(3):821－826.］

［3］ 王建民，吳艷紅，張忠，等.運載火箭全箭動特性三維建模技術(shù)［J］.中國科學(xué):技術(shù)科學(xué)，2014，44(1):50－61.［Wang Jian-min，Wu Yan-hong，Zhang Zhong，et al.Threedimensional modeling technology for dynamic characteristics of the launch vehicle［J］.Sci.Sin.Tech.，2014，44(1):50－61.］

［4］ 李青，韓增堯，馬興瑞.航天器貯箱液固耦合振動特性的仿真與試驗研究［J］.宇航學(xué)報，2014，35(11):1233－1237.［LiQing，Han Zeng-yao，MaXing-rui.Simulation and experimental research on fluid-structure-interaction dynamics of tank in spacecrafts［J］.Journal of Astronautics，2014，35 (11):1233－1237.］

［5］ 潘忠文.運載火箭動力學(xué)建模及振型斜率預(yù)示技術(shù)［J］.中國科學(xué)，E輯:技術(shù)科學(xué)，2009，39(3):469－473.［Pan Zhong-wen.The dynamical modeling and modal slope indicate technique of launch vehicle［J］.Science in China，Series E: Technological Sciences，2009，39(3):469－473.］

［6］ 葉長水，萬小朋，趙美英.采用縱橫扭一體化建模技術(shù)分析火箭牽制釋放過程動響應(yīng)［J］.航空工程進(jìn)展，2012，3(2): 178－182.［Ye Chang-shui，Wan Xiao-peng，Zhao Mei-ying.Analysis for the vehicle's structure transient response in the phases ofhold-down and release using longitudinal-lateraltorsional integrated modeling technique［J］.Advances in Aeronautical Science and Engineering，2012，3(2):178－182.］

［7］ 潘忠文，王旭，刑譽峰，等.基于梁模型的火箭縱橫扭一體化建模技術(shù)［J］.宇航學(xué)報，2010，31(5):1310－1316.［Pan Zhong-wen，Wang Xu，Xing Yu-feng，et al.A beam model based longitudinal-lateral-torsional integrated modeling technique for launch vehicle［J］.Journal of Astronautics，2010，31(5):1310－1316.］

［8］ 潘忠文，刑譽峰，朱禮文，等.運載火箭動力學(xué)建模中液體推進(jìn)劑模擬技術(shù)［J］.中國科學(xué):技術(shù)科學(xué)，2010，40(8): 920－928.［Pan Zhong-wen，Xing Yu-feng，Zhu Li-wen，et al.Liquid propellant analogy technique in dynamic modeling of launch vehicle［J］.Sci.China Tech.Sci.，2010，40(8):920－928.］

［9］ 林宏，羅恒，潘忠文，等.運載火箭動特性有限元模型修正技術(shù)研究［J］.載人航天，2011，(6):30－34.［Lin Hong，Luo Heng，Pan Zhong-wen，et al.Research on dynamic FEM model updating technique of launch vehicle［J］.Manned Spaceflight，2011，(6):30－34.］

［10］ 潘忠文，曾耀祥，廉永正，等.運載火箭結(jié)構(gòu)動力學(xué)模擬技術(shù)研究進(jìn)展［J］.力學(xué)進(jìn)展，2012，42(4):406－415.［Pan Zhong-wen，Zeng Yao-xiang，Lian Yong-zheng，et al.Reviews in structural dynamics analogy technique of launch vehicle［J］.Advances in Mechanics，2012，42(4):406－415.］

［11］ Pinson L D.Longitudinal spring constants for liquid-propellant tanks with ellipsoidal ends［R］.NASA TN D－2220，1964.

［12］ Staley J A.Dynamicstability ofspace vehicles， Vol. II， determination of longitudinal vibration modes［R］.NASA CR－936，1967.

［13］ 狄文斌，唐玉花.液體運載火箭縱向振動結(jié)構(gòu)動力學(xué)模型應(yīng)用研究［J］.強度與環(huán)境，2015，42(6):1－10.［Di Wenbin，Tang Yu-hua.Study on application of dynamical models for the longitudinal vibration of liquid-propellant launch vehicles［J］.Structure＆Environment Engineering，2015，42(6):1－10.］

［14］ MSC.Software Corporation.MSC Nastran 2012.2 quick reference guide［M］.Santa Ana:MSC.Software Corporation，2012.

通信地址:上海市閔行區(qū)元江路3888號(201109)

電話:(021)24186773

E-mail:yh_tang@163.com

(編輯:牛苗苗)

A One-Dimension Longitudinal-Lateral-Torsional Integrated Modeling Technique for Liquid-Propellant Launch Vehicle

TANG Yu-hua，DI Wen-bin，LIU Jing-hua
(Aerospace System Engineering Shanghai，Shanghai 201109，China)

A new one-dimension longitudinal-lateral-torsional integrated modeling method for liquid-propellant launch vehicle is proposed in this paper，which is based on the traditional beam model and combined with the longitudinal massspring model.Multi-mass of liquid for the longitudinal mass-spring model is also adopted.This modeling method is applied to a liquid-propellant roll booster of a launch vehicle.By comparing with the test results，it is proved that the lateral mode，longitudinal mode，torsional mode and coupling mode of them can be achieved by the one-dimension longitudinal-lateraltorsional integrated model，and the lateral frequency，longitudinal frequency and torsional frequency from this model have high precision.Thus，this modeling method can be applied to engineering.Increasing liquid mass number improves the precision of the longitudinal frequency，while it provides no effect on the precision of the longitudinal frequency when there is less liquid in the tank.

Liquid-propellant launch vehicle;Dynamic characteristic;One-dimension longitudinal-lateral-torsional integrated model;Multi-mass of liquid

V475.1;V414

A

1000-1328(2017)01-0089-08

10.3873/j.issn.1000-1328.2017.01.012

唐玉花(1984－)，女，碩士，主要從事運載火箭結(jié)構(gòu)動力學(xué)研究。

2016-06-20;

2016-10-09