劉傳龍，張宇文，王亞東，祁曉斌

（西北工業(yè)大學航海學院，陜西西安710072）

考慮適配器彈性的潛射導彈出筒載荷特性研究

劉傳龍，張宇文，王亞東，祁曉斌

（西北工業(yè)大學航海學院，陜西西安710072）

同心筒水下垂直發(fā)射導彈時潛艇是具備一定艇速的，相當于給彈體增加了側向橫流，其影響在彈體出筒階段尤為顯著，帶來發(fā)射安全性和彈體結構載荷相關的諸多問題。為獲得在時間和空間維度上置信度更高的導彈出筒載荷特性，并采用較少的假設，首次提出并建立了考慮適配器彈性約束和彈體三自由度的耦合求解流場和彈體運動的非定常數(shù)值模型。闡述了筒口氣泡彈性對彈體載荷特性的影響，同時獲得了出筒過程中的彈體運動、受載及適配器變形情況。仿真結果表明：由艇速產生的橫流造成的振蕩特性對導彈發(fā)射產生不利影響。結果可為同心筒水下垂直發(fā)射安全性評估和彈性結構設計提供研究方法。

兵器科學與技術；水下發(fā)射；載荷特性；流場與運動耦合計算；適配器彈性

0　引言

潛射導彈具備發(fā)射平臺隱蔽、機動性高，打擊突然、火力猛烈等一系列的優(yōu)勢，同心筒式熱發(fā)射作為潛射導彈發(fā)展的方向，成為各大國科研人員研究的熱點。為得到合理的水下垂直發(fā)射流場和彈道特性參數(shù)，必先得到滿足結構強度的載荷特性。由于同心筒水下垂直發(fā)射流場的非定常性，特別是涉及了不同介質間界面捕捉問題，常規(guī)的理論方法必須依賴大量的假設來實施，故而丟失了太多的細節(jié)信息，即使能預報載荷特性，也是在理想條件或平均意義下給出。水下發(fā)射過程中引入了水介質和海洋環(huán)境，發(fā)射裝置的工作狀態(tài)發(fā)生了本質的改變：高溫燃氣在水下形成的筒口氣泡和水介質的交互作用，會帶來一系列的流場和載荷特性變化［1］；導彈發(fā)動機在發(fā)射筒內點火，而發(fā)射筒的燃氣排導為筒口流場所限，表現(xiàn)出和空氣介質中發(fā)射不同的流場和載荷分布特性［2］；潛艇在發(fā)射導彈時是具備一定艇速的，相當于給垂直發(fā)射出筒的彈體增加了橫流，其影響在彈體出筒階段尤為顯著，帶來了彈體結構載荷相關的諸多問題［3-7］。

本文選擇了適合于水下垂直發(fā)射建模和進一步研究的商用軟件Fluent作為基本平臺，基于基本的流體控制方程并根據(jù)剛體動力學理論進行二次開發(fā)，結合動網格技術建立了耦合求解流場和彈體運動的非定常仿真模型，建模中增加了彈性適配器［8］約束，并考慮了彈體在出筒過程中的微幅運動，采用了較少的假設，能獲得在時間和空間維度上置信度更高的結果。繼而提高導彈垂直發(fā)射出筒過程中的彈體載荷預報能力。

1　模型建立

1.1彈性適配器建模

在發(fā)射過程中，彈體在橫流作用下偏轉，壓縮適配器，適配器變形產生彈性恢復力迫使彈體作彈性振動。由于彈體在發(fā)射筒內的偏轉角度非常小，忽略適配器的切向變形與應力，認為適配器產生的正壓力僅為最大壓縮量的函數(shù)。如圖1所示，適配器的變形可由彈筒間隙來描述，求得彈筒在O0x0y0剖面內的間隙分布，即可通過切片法求出每個適配器切片產生的彈力。

圖1　適配器變形表示Fig.1 Schematic diagram of adapter deformation

將適配器沿發(fā)射筒軸線方向切片，獲得每一切片中心在地面坐標系下的坐標xi0、寬度wi和厚度di.在導彈發(fā)射出筒過程中的任一時刻，第i切片處彈筒間隙計算方法如下：

1）將切片中心坐標轉換到彈體坐標系，得到切片相對于彈體的位置

式中：xi為切片在彈體坐標系下的坐標；x0為彈體質心坐標。

2）如果切片位于圓柱段之外，則彈體不會壓縮切片，切片變形量Δdi=0.

3）如果切片位于圓柱段上，則計算切片處彈體中心偏離筒軸的距離，作為該切片的變形量Δdi= xisin θ+yi，y0為導彈在地面坐標系下對應坐標。

適配器的徑向變形與作用力之間的關系通常利用多項式擬合：F=k1Δd+k2Δd2+k3Δd3+….由于不同的變形量范圍，作用力與變形量之間的對應規(guī)律可能不同，這種情況需要分段擬合。現(xiàn)以兩段擬合：第1段為線性關系、第2段為二次方關系為例，說明擬合關系式的建立與應用。

設適配器徑向變形量在0≤Δd＜Δd1范圍內，作用力與變形量呈線性關系；變形量在Δd1≤Δd≤Δd2范圍內，作用力與變形量呈二次方關系。

第1段利用關系式F=k1Δd擬合，求出k1；第2段利用關系式F=F1+k2（Δd-Δd1）+k3（Δd-Δd1）2擬合，求得k2和k3.

根據(jù)上述的擬合公式可以得到適配器任一截面（單位高度）的徑向變形量Δdr與作用力的關系：

式中：ha為適配器高度；a1為常數(shù)項，

a2為系數(shù)，

所有適配器的作用力通過對所有切片作用力累加得到

所有適配器的作用力矩通過對所有切片力矩累加得到

1.2彈體數(shù)值模型

同心筒發(fā)射裝置以垂直方式布置于發(fā)射潛艇中，通過導彈助推器自力熱發(fā)射出筒，本文研究中不考慮艇體表面附著的小部件和艇體外形，僅建模同心筒內部主體結構和導彈外形，外圍邊界僅考慮艇體表面約束，其余邊界均取為自由流邊界，計算域如圖2所示。

圖2　計算區(qū)域Fig.2 Compute area

導彈發(fā)射非定常模型主要分為流場求解、運動求解和邊界運動3個模塊。其中流場求解由Fluent求解器完成，提供運動求解的基本素材；運動求解由UDF完成，從流場求解結果中獲取彈體受力和力矩，根據(jù)運動和動力學關系求解運動；邊界運動由UDF控制，F(xiàn)luent求解器完成。非定常仿真實現(xiàn)的基本流程如圖3所示。

圖3　非定常模型求解流程Fig.3 Solving process of unsteady flow model

由于彈體運動不再約束為一自由度運動，增加了側向運動和繞z0軸的轉動，變?yōu)槿杂啥冗\動，且發(fā)生側向運動和轉動后壓縮適配器會產生適配器的約束反力、力矩，故而導彈運動和動力學方程為

當運動維度增加到三自由度時，由于彈體存在側向運動和轉動，發(fā)生了相對于網格的位移，相應地要求網格發(fā)生變化，甚至是重構以適應相對位移，此時包裹彈體區(qū)域采用結構化網格劃分已經無法滿足要求?？紤]到在整個出筒過程中，彈體側向運動和轉動幅度很小，拉扯網格的力度不會特別大，故而在包裹彈體的局部區(qū)域使用四面體網格，其他區(qū)域則使用結構化網格，網格劃分如圖4所示。

圖4　網格劃分Fig.4 Mesh partitioning

根據(jù)上述網格劃分，賦予彈體三自由度的運動，流體運動區(qū)域僅賦予軸向運動。對于軸線上的流體區(qū)域采用Layering方法更新，包裹彈體的區(qū)域則采用Smoothing和Remeshing方法更新。

1.3數(shù)值模型設置

模型設置要注意以下3點：

1）采用SST k-ω湍流模型，選用Mixture多相流模型，且不考慮相之間的滑移；

2）導彈水下垂直發(fā)射出筒時間較短，水汽相變效應不顯著，不考慮水汽相變效應；

3）壓力項離散方法選用PRESTO，副相體積分數(shù)離散方法采用一階迎風格式；

4）壓力速度耦合算法選用PISO.

2　載荷特性

本文研究中僅選用縮比的模型彈和發(fā)射裝置進行，設彈性適配器變形和約束力呈線性關系，以抗壓剛度來表征彈性系數(shù)。同心筒發(fā)射裝置中的參數(shù)如表1所示。

表1　發(fā)射裝置參數(shù)Tab.1 Parameters of launcher

圖5為8組彈性適配器在同心筒發(fā)射裝置中的布置示意圖，給出了適配器的安裝位置、外形、間隔尺寸及編號規(guī)則。

圖5　適配器布置示意圖Fig.5 Schematic diagram of adapter arrangement

設定發(fā)射艇速為1.0 m/s，發(fā)射深度為9 m，得到了考慮彈體三自由度運動和適配器彈性約束的仿真計算結果，圖6為不同時刻的筒口氣泡形態(tài)演變過程。從圖6可得到與筒口氣泡演變相關的載荷產生和變化的原因：

1）筒口氣泡受橫流影響偏移造成的不對稱附著為橫向載荷提供了條件；

2）氣泡彈性效應所致的壓力周期脈動是橫向載荷的直接成因，氣泡在過壓縮時泡內壓力升高，導致彈體背流面壓力大于迎流面，產生y軸正向力，和橫向流所致力疊加生成正峰值，反之氣泡過膨脹時生成負峰值；

3）側向力矩與側向力大小及壓心位置相關，彈體出筒過程中，筒口氣泡形態(tài)和附著位置均在不斷變化，從而側向力和力矩也隨之變化；

4）圖7為彈體在發(fā)射過程中運動特性和受流體載荷特性，圖8為特性適配器產生的約束反力、反力矩和變形特性。

由圖7和圖8可得以下結論：

1）在艇速引起的橫流、彈性適配器約束綜合作用下，彈體在出筒過程中運動參數(shù)出現(xiàn)了明顯的振蕩過程，振蕩幅度隨出筒長度增加而增大，彈尾接近離筒時，彈身偏向橫流指向方向；

2）離筒過程中，適配器約束力和約束力矩振蕩幅值隨出筒距離的增大而增大，最大值和相應的流體載荷數(shù)值相當甚至比其大一個量級；

3）彈體出筒距離增大后，迎流面積增加，流體動力增大，故筒口附近的幾個適配器在彈體離筒末階段承受較大的結構載荷，載荷方向由彈體側向位移和俯仰姿態(tài)綜合決定。

圖6　筒口氣泡形態(tài)演變Fig.6 Evolution of outlet bubble shape

圖7　彈體特性Fig.7 Missile characteristics

圖8　彈性適配器特性Fig.8 Elastic adapter characteristics

3　結論

1）本文建立了考慮適配器彈性的耦合求解流場和彈體運動的非定常數(shù)值模型，對由艇速引起的橫流作用下導彈的出筒運動和載荷進行了仿真計算。結果表明所建模型是可行的，為分析研究同心筒水下垂直發(fā)射研究提供了一條途徑。

2）由于艇速引起的橫流所造成的影響導致導彈運動特性出現(xiàn)振蕩，筒口附近的適配器承受較大的結構載荷，會對導彈發(fā)射產生較大的不利影響。

（

）

［1］馬慶鵬.潛射航行體出筒過程載荷特性的數(shù)值研究［D］.哈爾濱：哈爾濱工業(yè)大學，2011. MA Qing-peng.Numerical simulation on the tube-exit load characteristics of submarine-launched vehicle［D］.Harbin：Harbin Institute of Technology，2011.（in Chinese）

［2］Weiland C J，Vlachos P P，Yagla J J.Concept analysis and laboratory observations on a water piercing missile launcher［J］.Ocean Engineering，2010，37（11）：959-965.

［3］殷崇一.潛射導彈發(fā)射與出水載荷研究［D］.西安：西北工業(yè)大學，2004. YIN Chong-yi.Research ofsubmarine-launched missile and effluent load［D］.Xi'an：Northwestern Polytechnical University，2004.（in Chinese）

［4］裴譞，張宇文，袁緒龍.潛射導彈垂直發(fā)射橫向振動特性仿真分析［J］.兵工學報，2009，30（8）：1056-1060. PEI Xuan，ZHANG Yu-wen，YUAN Xu-long.A simulation andanalysis on the lateral vibration characteristics of vertical launching for submarine-borne missile［J］.Acta Armamentarii，2009，30（8）：1056-1060.（in Chinese）

［5］裴譞，張宇文，李代金.導彈水下動機座垂直發(fā)射載荷建模與仿真［J］.計算機仿真，2010，27（2）：104-107. PEI Xuan，ZHANG Yu-wen，LI Dai-jin.Missile vertical motive sub-marine launching payload modeling and simulation［J］.Computer Simulation，2010，27（2）：104-107.（in Chinese）

［6］王英華.軸向載荷對潛射導彈橫向振動特性影響分析［J］.導彈與航天運載技術，2002（6）：9-15. WANG Ying-hua.Analysis of effect ofaxial load on lateral vibration of underwater-to-surface missile［J］.Missile and Space Vehicles，2002（6）：9-15.（in Chinese）

［7］覃東升，方勇.水下垂直發(fā)射出筒載荷響應分析［J］.計算機仿真，2013，30（2）：112-119. QIN Dong-sheng，F(xiàn)ANG Yong.Analysis of tube-departing load response of SLBS vertical launching system［J］.Computer Simulation，2013，30（2）：112-119.（in Chinese）

［8］尚書聰，孫建中.兩種彈筒適配方式對導彈出筒姿態(tài)的影響［J］.彈道學報，2012，24（2）：96-100. SHANG Shu-cong，SUN Jian-zhong.Effect of two adaptation methods of missile and tube on outlet attitude of missile［J］.Journal of Ballistics，2012，24（2）：96-100.（in Chinese）

Investigation into Load Characteristics of Submarine-launched Missile Being Ejected from Launch Tube Considering the Adapter Elasticity

LIU Chuan-long，ZHANG Yu-wen，WANG Ya-dong，QI Xiao-bin
（School of Marine Science and Technology，Northwestern Polytechnical University，Xi'an 710072，Shaanxi，China）

A cross flow is applied on a missile，which has a great influence on the stage of out of tube，when a missile is launched from the concentric canister launcher during the movement of submarine.This leads to a series of launching safety and load problems.In order to obtain a higher degree of confidence of load characteristics out of tube both in time and spatial dimensions and use the less hypotheses，an unsteady simulation model which can solve the flow field and missile motion in consideration of the adapter elasticity and 3 DOF motion of missile is proposed and established by the first time.The effect of bubble elasticity on the load characteristics is analyzed.The motion parameters，load characteristics and adapter deformation during the launching of missile from launch tube are obtained.The results show that the oscillation resulting from the cross flow has an adverse effect on missile launching.

ordnance science and technology；underwater launching；load characteristics；coupling calculation of flow field and motion；adapter elasticity

TJ760.1

A

1000-1093（2015）02-0379-06

10.3969/j.issn.1000-1093.2015.02.027

2014-05-08

劉傳龍（1988—），男，博士研究生。E-mail：liu2011100322@163.com；張宇文（1946—），男，教授，博士生導師。