姚保太，王漢平

(1.鄭州機(jī)電工程研究所，鄭州 450052;2.北京理工大學(xué)宇航學(xué)院，北京 100081)

導(dǎo)彈彈射系統(tǒng)中緩沖制動(dòng)錐的軸向沖擊特性①

姚保太1，王漢平2

(1.鄭州機(jī)電工程研究所，鄭州 450052;2.北京理工大學(xué)宇航學(xué)院，北京 100081)

針對(duì)壁厚線性變化制動(dòng)錐的沖擊力學(xué)特性，基于軸壓力學(xué)模型，提出了等效速度概念，在綜合考慮材料應(yīng)變強(qiáng)化和應(yīng)變率效應(yīng)的條件下，按等效壁厚、等效速度條件下的分段能量等效原理得到了制動(dòng)錐軸向沖擊力學(xué)模型，并采用理論分析、有限元計(jì)算和實(shí)驗(yàn)研究相結(jié)合的方式研究了彈射系統(tǒng)中緩沖制動(dòng)錐的沖擊皺褶形成形態(tài)、沖擊壓縮量和緩沖力學(xué)特性，3種方法得到的結(jié)果基本一致，驗(yàn)證了所提模型的正確。所提出的制動(dòng)錐軸向沖擊力學(xué)模型可為制動(dòng)錐的初步設(shè)計(jì)、實(shí)驗(yàn)規(guī)劃以及武器系統(tǒng)的發(fā)射動(dòng)力學(xué)分析提供參考。

彈射系統(tǒng);制動(dòng)錐;應(yīng)變硬化;應(yīng)變率效應(yīng);等效壁厚;等效速度;沖擊

0 引言

薄壁圓柱殼作為緩沖吸能裝置具有結(jié)構(gòu)簡(jiǎn)單緊湊、制造加工方便、吸能效率較高、沖擊變形規(guī)則且力學(xué)性能優(yōu)良等特點(diǎn)，廣為航空航天、武器工業(yè)和核電系統(tǒng)所使用［1－3］。針對(duì)薄壁圓柱殼的軸向壓縮破壞模式及其力學(xué)特性的研究，早在20世紀(jì)60年代就已展開，Alexander等［4］基于材料的剛塑性假設(shè)利用固定塑性鉸模型導(dǎo)出了軸對(duì)稱漸進(jìn)屈曲時(shí)的褶皺波長(zhǎng)和平均軸壓力公式，為本領(lǐng)域研究的開展奠定了基礎(chǔ);Johnson和 Reid［5］、Abramowicz和 Jones［6］在研究中考慮了幾何效應(yīng)，提出了有效壓縮長(zhǎng)度的概念，對(duì)Alexander模型進(jìn)行了修正，獲得了更合理的結(jié)果，同時(shí)還對(duì)沖擊載荷作用下的材料應(yīng)變強(qiáng)化和應(yīng)變率效應(yīng)進(jìn)行了簡(jiǎn)化計(jì)算，將軸壓特性推廣到了軸向沖擊;Wierzbicki和Bhat提出了移動(dòng)塑性鉸模型［7］，并據(jù)此模擬了載荷-位移曲線;此前的模型均基于管壁內(nèi)翻或外翻，而試驗(yàn)表明，疊縮后的管壁是以一定比率分布在未疊縮管壁的兩側(cè)，據(jù)此，Wierzbicki等［8］修正了他們的模型，提出了偏心率效應(yīng)，并給出了不同偏心率下的載荷-位移曲線，但假設(shè)偏心率對(duì)吸能效率沒有影響;Singace等［9］改進(jìn)了Wierzbicki等［8］的模型，認(rèn)為吸能效率與偏心率有關(guān)，并給出了最優(yōu)偏心率，同時(shí)通過實(shí)驗(yàn)驗(yàn)證了其結(jié)果［10］。而針對(duì)壁厚線性變化薄壁圓柱殼的研究，公開材料比較少，畢世華和王漢平［11］基于塑性鉸理論，采用功能原理和能量守恒定律進(jìn)行了理論探索，并借助有限元軟件對(duì)軸壓漸進(jìn)屈曲特性進(jìn)行了分析;王漢平和王忠峰［12－13］通過對(duì)上述各種模型進(jìn)行對(duì)比驗(yàn)證，發(fā)現(xiàn)Wierzbicki等提出的不含偏心率效應(yīng)的移動(dòng)塑性鉸模型與制動(dòng)錐塑性變形模式更吻合，并在此基礎(chǔ)上提出了等效壁厚概念，將該理論推廣應(yīng)用到了壁厚線性變化的薄壁圓柱殼的失效和軸壓載荷分析，獲得了較好的效果。壁厚線性變化的薄壁圓柱殼制動(dòng)錐是提拉式導(dǎo)彈彈射裝置中(如俄羅斯車載的S-300、S-400、道爾-M1和艦載的利夫和克里偌克等)的必備部件，其沖擊力學(xué)特性是影響產(chǎn)品設(shè)計(jì)和使用安全性的關(guān)鍵，而當(dāng)前對(duì)沖擊特性的研究尚比較缺乏，尤其是在充分考慮材料應(yīng)變硬化和應(yīng)變率效應(yīng)之后。工程上比較普遍的做法是采用有限元分析與實(shí)驗(yàn)測(cè)試相結(jié)合的模式來獲取最佳設(shè)計(jì)，但這種方法計(jì)算量大、周期長(zhǎng)，且數(shù)據(jù)分析和獲取均不夠直觀，這給工程應(yīng)用帶來諸多不便。

本文在壁厚線性變化制動(dòng)錐軸壓理論的基礎(chǔ)上提出了等效速度的概念，在綜合考慮材料應(yīng)變強(qiáng)化和應(yīng)變率效應(yīng)的條件下，按等效壁厚、等效速度條件下的分段能量等效原理得到了制動(dòng)錐軸向沖擊力學(xué)模型，該模型能準(zhǔn)確預(yù)測(cè)沖擊狀態(tài)下制動(dòng)錐皺褶個(gè)數(shù)、壓縮量以及平均壓縮力變化情況，可為此類制動(dòng)錐沖擊性能預(yù)測(cè)、沖擊實(shí)驗(yàn)規(guī)劃以及武器系統(tǒng)的發(fā)射動(dòng)力學(xué)分析提供參考。

1 緩沖制動(dòng)錐的軸向沖擊理論

1.1 等壁厚圓柱殼的等速軸向沖擊理論

制動(dòng)錐在變形過程中受提拉桿的支撐，只能向外皺褶，這種失效模式與Wierzbicki等提出的不考慮偏心率效應(yīng)的移動(dòng)塑性鉸模型更加吻合［7，12－13］，但也僅適用于擬靜態(tài)的軸壓過程，針對(duì)等速軸向沖擊，由于材料應(yīng)變率效應(yīng)的存在，其軸向沖擊載荷將會(huì)出現(xiàn)顯著差異，因此需要對(duì)軸壓模型予以修正，以包含因沖擊而產(chǎn)生的應(yīng)變率效應(yīng)的影響。

根據(jù)Wierzbicki等的理論，圓柱殼的壓皺半波長(zhǎng)和平均壓皺力為

式中 σ0為材料屈服極限應(yīng)力;D、q為材料應(yīng)變率相關(guān)的系數(shù)。

在靜態(tài)σ0(ε0)的基礎(chǔ)上乘上當(dāng)量應(yīng)變率d時(shí)的應(yīng)變率效應(yīng)系數(shù)就可獲取應(yīng)變率效應(yīng)影響下的當(dāng)量屈服應(yīng)力σd，用σd替換式(2)中的σ0就可包含材料應(yīng)變強(qiáng)化效應(yīng)和應(yīng)變率效應(yīng)的影響。

1.2 壁厚線性變化圓柱殼軸向沖擊模型

根據(jù)等效壁厚概念［12－13］，可得到各分段等效壁厚的形式如下:

式中 n表示第n個(gè)褶皺;tn為第n個(gè)褶皺的等效壁厚;tn1為第n個(gè)褶皺薄壁端壁厚;tn2為第n個(gè)褶皺厚壁端壁厚。

在沖擊過程中，提拉桿的速度是一個(gè)變化的量，為擴(kuò)展利用等壁厚、勻速?zèng)_擊模型，故根據(jù)應(yīng)變率效應(yīng)的應(yīng)力影響系數(shù)與勻速?zèng)_擊速度間的冪函數(shù)關(guān)系，利用功能原理和能量守恒定律就可獲得應(yīng)變率效應(yīng)的當(dāng)量應(yīng)力影響系數(shù)即為皺褶形成始末的應(yīng)力影響系數(shù)的算術(shù)平均，由此提出等效速度概念，即在皺褶形成過程中對(duì)應(yīng)于應(yīng)變率效應(yīng)的當(dāng)量應(yīng)力影響系數(shù)的速度即為等效速度，其形式如下:

式中 Vn為形成第n個(gè)褶皺的等效速度;Vn1為第n個(gè)褶皺開始形成時(shí)的沖擊速度;Vn2為第n+1個(gè)褶皺開始形成時(shí)的沖擊速度。

根據(jù)功能原理、有效壓縮量與皺褶波長(zhǎng)的關(guān)系［3］可得

本段圓柱殼可看作壁厚為tn的等壁厚圓柱殼，圓柱殼的褶皺長(zhǎng)可按等壁厚圓柱殼褶皺長(zhǎng)公式求出:

式中 Hn為第n個(gè)褶皺半波長(zhǎng)。

另外，還存在一個(gè)已知的隱含條件，即壁厚變化楔形角α已知，且

聯(lián)立方程式(8)～式(14)，設(shè)定計(jì)算終止條件為Vn2≤0，采用遞推的方式就可逐級(jí)求解每個(gè)褶皺的等效壁厚、等效速度、褶皺半波長(zhǎng)、軸向沖擊平均力、皺褶形成后的速度以及每個(gè)皺褶吸收的能量。

2 有限元仿真、實(shí)驗(yàn)研究和理論分析結(jié)果對(duì)比

為方便制動(dòng)錐的變參數(shù)仿真和特性優(yōu)化，本文使用ABAQUS的參數(shù)化建模方法將制動(dòng)錐薄壁端厚度、厚壁端厚度、法蘭半徑、制動(dòng)錐內(nèi)徑、上下法蘭厚度、制動(dòng)錐有效長(zhǎng)度、單元最小尺度、沖擊塊質(zhì)量以及沖擊速度設(shè)置為參數(shù)，借助ABAQUS的關(guān)鍵字NODE、NGEN、NFILL、ELEMENT、ELGEN構(gòu)建了制動(dòng)錐參數(shù)化的軸對(duì)稱有限元模型，并編寫了PYTHON腳本，可實(shí)施實(shí)驗(yàn)設(shè)計(jì)(DOE)仿真;還對(duì)制動(dòng)錐的材料進(jìn)行了擬靜態(tài)和動(dòng)態(tài)實(shí)驗(yàn)，獲取了材料的彈塑性本構(gòu)模型。為了驗(yàn)證上述理論模型，對(duì)不同結(jié)構(gòu)尺寸的壁厚線性變化圓柱殼、不同的沖擊速度狀況進(jìn)行了沖擊有限元分析［12，14］。同時(shí)還針對(duì)2款制動(dòng)錐進(jìn)行了不同速度的沖擊實(shí)驗(yàn)，并對(duì)仿真計(jì)算、理論分析以及實(shí)驗(yàn)結(jié)果進(jìn)行了對(duì)比分析。圓柱殼長(zhǎng)度分別為250 mm和150 mm，內(nèi)徑40 mm，兩端壁厚不同。為了方便描述，本文采用了圓柱殼長(zhǎng)度和上下端壁厚值對(duì)圓柱殼進(jìn)行了編號(hào)，如:“250-0624”表示上端壁厚為 0.6 mm，下端壁厚2.4 mm，長(zhǎng)度為250 mm的圓柱殼。實(shí)驗(yàn)采用彈簧儲(chǔ)能釋放沖擊臺(tái)，由于受結(jié)構(gòu)限制，對(duì)于25 kg的沖擊物，該沖擊臺(tái)所能提供的最大出射速度為20 m/s。

因此，本文就250-0624件沖擊速度為17、18、19、20 m/s和150-0513件沖擊速度為16 m/s進(jìn)行了實(shí)驗(yàn)。

2.1 褶皺數(shù)分析

圖1是不同沖擊速度條件下制動(dòng)錐最終的變形情況和仿真結(jié)果的對(duì)比。從中不難看出，實(shí)驗(yàn)和仿真中，制動(dòng)錐均在薄壁端形成皺褶，且同結(jié)構(gòu)尺寸的制動(dòng)錐壓縮量和皺褶個(gè)數(shù)與沖擊速度正相關(guān)。

圖1 不同試件在不同沖擊速度情況下的變形仿真與實(shí)驗(yàn)結(jié)果對(duì)比Fig.1 Comparison of computational and experimental deformation under various impact velocities

表1數(shù)據(jù)表明，軸向沖擊特性條件下將壁厚線性變化的緩沖制動(dòng)錐按等效壁厚、等效速度分段等效簡(jiǎn)化后計(jì)算所得的褶皺數(shù)與有限元仿真以及實(shí)驗(yàn)結(jié)果一致。理論、仿真及實(shí)驗(yàn)結(jié)果均表明壁厚線性變化圓柱制動(dòng)錐在受到軸向沖擊時(shí)的最終軸向褶皺數(shù)隨著沖擊速度的增加而增多，且壓縮量也相應(yīng)增大，由于沖擊速度差異不大，最終有些皺褶并未完全成形。

表1 不同沖擊速度下試件的褶皺數(shù)和壓縮量對(duì)比Table 1 Comparison of buckling waves and compression length under various impact velocities

2.2 軸向沖擊平均力分析

沖擊載荷作用下緩沖制動(dòng)錐的緩沖制動(dòng)時(shí)間較短，由于選用的傳感器采樣頻率不夠，同時(shí)傳感器與有源放大器之間參數(shù)存在匹配問題，實(shí)驗(yàn)時(shí)未能很好記錄到完整皺褶形成過程的緩沖載荷歷程，但所記錄到的最大載荷也具有一定的參考價(jià)值，表2是實(shí)驗(yàn)記錄最大緩沖載荷與計(jì)算緩沖載荷的對(duì)比。

表2 不同沖擊速度下250-0624件最大緩沖力的仿真與實(shí)驗(yàn)對(duì)比Table 2 Comparison of maximum impact loads under various impact velocities

考慮到仿真計(jì)算的皺褶波數(shù)以及沖擊后制動(dòng)錐的壓縮量與實(shí)驗(yàn)結(jié)果吻合較好，本文在進(jìn)行軸向沖擊平均力分析時(shí)直接對(duì)仿真緩沖載荷-行程歷程采用目標(biāo)函數(shù)f(x)=A+B·xC進(jìn)行了曲線擬合，同時(shí)將擬合曲線、仿真曲線和理論計(jì)算的軸向沖擊平均力進(jìn)行了對(duì)比，不同沖擊速度條件下的對(duì)比情況見圖2。

圖2 不同試件在不同沖擊速度條件下軸向平均沖擊力理論計(jì)算與仿真結(jié)果對(duì)比Fig.2 Comparison of theoretical and FEA computational curves for linear thickness braking cylindrical shells under various impact velocity

由圖2可看出，不同沖擊速度條件下有限元計(jì)算 緩沖載荷擬合曲線與理論計(jì)算曲線基本吻合，但總體上理論計(jì)算值偏大。另外，高速攝像(盡管因采樣頻率不足，無法完整記錄皺褶發(fā)展過程)和有限元仿真動(dòng)畫表明，沖擊過程中制動(dòng)錐的總體變形效果是薄壁端先出現(xiàn)褶皺，然后依次向厚壁端推進(jìn)。仿真數(shù)據(jù)表明，在緩沖器剛開始承受沖擊時(shí)，載荷存在較大的脈沖峰值，而隨后軸向沖擊載荷隨褶皺向厚壁端推進(jìn)而逐漸增大，且載荷在增大過程中還呈現(xiàn)出一定程度的波動(dòng)特性，波動(dòng)幅度也隨壁厚增大而增大。

3 結(jié)論

(1)針對(duì)壁厚線性變化緩沖制動(dòng)錐的沖擊力學(xué)特性，提出了等效速度概念，并聯(lián)立制動(dòng)錐等效壁厚和有效壓縮量概念，利用圓柱殼分段能量等效辦法建立了圓柱殼的沖擊力學(xué)模型;模型中綜合考慮了材料應(yīng)變強(qiáng)化和應(yīng)變率效應(yīng)的影響，使得模型具有更好的精度。

(2)通過材料動(dòng)態(tài)試驗(yàn)獲得的本構(gòu)模型進(jìn)行了理論分析、有限元計(jì)算和沖擊實(shí)驗(yàn)的對(duì)比研究表明，制動(dòng)錐皺褶數(shù)一致、沖擊壓縮量吻合、最大沖擊載荷基本相符，這驗(yàn)證了等效速度和等效壁厚概念聯(lián)立使用的制動(dòng)錐沖擊力學(xué)模型的可信性，該模型可用于制動(dòng)錐的工程設(shè)計(jì)及其軸向沖擊力特性的預(yù)測(cè)，也可用于提拉式彈射裝置的彈射發(fā)射動(dòng)力學(xué)的建模。

［1］諾曼·瓊斯.結(jié)構(gòu)沖擊［M］.蔣平，譯.成都:四川教育出版社，1994.

［2］余同希，盧國(guó)興.材料與結(jié)構(gòu)的能量吸收-耐撞性 包裝安全防護(hù)［M］.華云龍，譯.北京:化學(xué)工業(yè)出版社，2006.

［3］杜星文，宋宏偉.圓柱殼沖擊動(dòng)力學(xué)及耐撞性設(shè)計(jì)［M］.北京:科學(xué)出版社，2004.

［4］Alexander J M.An approximate analysis of the collapse of thin cylindrical shells under axial loading［J］.Q.J.Mech.，1960，13(1):10-15.

［5］Johnson W，Reid S.Metallic energy dissipating system［J］.Applied Mechanics Reviews，1978(31):277-288.

［6］Abramowicz W，Jones N.Dynamic progressive buckling of circular and square tubes［J］.Int.J.Impact Engng，1986，4(4):243-270.

［7］Wierzbicki T，Bhat S U.A moving hinge solution for axisymmetric crushing of tubes［J］.International Journal of Mechanic and Science，1985，28(3):135-151.

［8］Wierzbicki T，Bhat S U，Abramowicz W，et al.Alexander revisited—a two folding element model of progressive crushing of tubes［J］.International Journal of Solids and Structures，1992(29):3269-3288.

［9］Singace A，Elsobky H，T Y R.On the eccentricity factor in the progressive crushing of tubes［J］.International Journal of Solids and Structures，1995(32):3589-3602.

［10］Singace A，Elsobky H.Further experimental investigation on the eccentricity factor in the progressive crushing of tubes［J］.International Journal of Solids and Structures，1996(33):3517-3538.

［11］畢世華，王漢平，梁征.導(dǎo)彈垂直彈射過程中制動(dòng)錐的動(dòng)力學(xué)特性研究［J］.北京理工大學(xué)學(xué)報(bào)，2004，24(9):762-765.

［12］王忠峰.提拉桿式彈射系統(tǒng)緩沖制動(dòng)錐力學(xué)特性研究［D］.北京:北京理工大學(xué)，2006.

［13］王漢平，王忠峰.導(dǎo)彈彈射系統(tǒng)中緩沖制動(dòng)錐的軸壓特性［J］.北京理工大學(xué)學(xué)報(bào)，2007，27(2):99-102.

［14］Karagiozova D，Alves M，Jones N.Inertia effects in axisymmetrically deformed cylindrical shells under axial impact［J］.International Journal of Impact Engineering，2000，24:1083-1115.

(編輯:呂耀輝)

Mechanical characteristics of the braking cylindrical shells of missile ejector under axial impact

YAO Bao-tai1，WANG Han-ping2
(1.Zhengzhou Mechano-Electronic Engineering Research Institute，Zhengzhou 450052，China;2.School of Aerospace Engineering，Beijing Institute of Technology，Beijing 100081，China)

For impact mechanical characteristics of the linear thickness braking cylindrical shells，a concept of equivalent velocity was proposed.Based on the energy equivalent principle of subsections under equivalent thickness and eguivalent velocity condition，an impact mechanical model of linear thickness braking cylindrical shells，which takes into account the strain hardening and strain rate effect of structural material，was built.The wrinkle creation forms，compress stroke and mechanical characteristics of the braking cylindrical shells of missile ejector under axial impact has been studied by means of theory analysis，F(xiàn)EM and experiment.Three methods get identical results，which verifies the correctness of the model.The theory model of linear thickness braking cylindrical shells can be used in the initial engineering design and experiment planning of the braking cylindrical shells and in the launching dynamics analysis of the weapon system.

missile ejector;braking cylindrical shell;strain hardening;strain rate effect;equivalent thickness;equivalent velocity;impact

TJ768.2

A

1006-2793(2014)06-0863-05

10.7673/j.issn.1006-2793.2014.06.023

2013-11-05;

2014-05-26。

姚保太(1976—)，男，高工，主要研究方向?yàn)楸靼l(fā)射理論與技術(shù)。E-mail:yaobaotai@163.com

王漢平(1971—)，男，副教授。E-mail:whp161@bit.edu.cn