寧全利，蔣濱安，梁希，李俊，3

（1.陸軍軍官學(xué)院高過載彈藥制導(dǎo)控制與信息感知實驗室，安徽合肥230031；

2.燕山大學(xué)建筑工程與力學(xué)學(xué)院，河北秦皇島066004；3.陸軍軍官學(xué)院二系，安徽合肥230031）

火炮發(fā)射環(huán)境下彈載光電器件結(jié)構(gòu)動態(tài)響應(yīng)研究

寧全利1，蔣濱安1，梁希2，李俊1，3

（1.陸軍軍官學(xué)院高過載彈藥制導(dǎo)控制與信息感知實驗室，安徽合肥230031；

2.燕山大學(xué)建筑工程與力學(xué)學(xué)院，河北秦皇島066004；3.陸軍軍官學(xué)院二系，安徽合肥230031）

為了研究彈載光電器件抗高過載機理，針對火炮發(fā)射環(huán)境特點，以某型特種彈所用彈載光電器件為例，建立了加裝減載組件的彈載光電器件有限元分析模型并進行仿真分析。采用應(yīng)力波理論，建立了反映彈載光電器件動力學(xué)過程的理論分析模型，計算求解作用于鏡頭的應(yīng)力、相對位移和相對速度等，得到了彈載光電器件的動力學(xué)響應(yīng)。結(jié)果表明：數(shù)值計算與仿真模擬結(jié)果有較好的一致性，驗證了所建立的理論分析模型的正確性，可用于火炮發(fā)射環(huán)境下不同彈載器件的動態(tài)特性研究和抗高過載設(shè)計的理論分析。

兵器科學(xué)與技術(shù)；彈載光電器件；數(shù)值計算；應(yīng)力波；火炮發(fā)射

0 引言

彈載光電器件是炮兵新型彈藥中一類重要的組成單元，在火炮發(fā)射時，隨同彈丸承受上萬g的瞬時高過載，其安全與否直接關(guān)系到該類彈藥使用的可靠性。對火炮發(fā)射環(huán)境下彈載光電器件動態(tài)響應(yīng)進行研究，是彈載器件抗過載機理研究的一項重要內(nèi)容，也是解決該類器件在高過載彈藥中可靠應(yīng)用的基礎(chǔ)和前提［1］。

由于火炮發(fā)射環(huán)境下，彈丸在火炮膛內(nèi)的運動時間通常在幾毫秒到十幾毫秒之間，而一般沖擊過程的時間量級通常在數(shù)百微秒左右，因此不能簡單移植沖擊動力學(xué)有關(guān)理論來研究火炮發(fā)射環(huán)境下彈載光電器件的動態(tài)響應(yīng)，需要根據(jù)火炮具體的發(fā)射條件和彈載光電器件的裝配情況進行分析。由于實彈試驗組織程序復(fù)雜，且彈丸內(nèi)裝各類器件的動態(tài)響應(yīng)過程難以檢測，目前也無其他有效的測試和驗證手段。本文針對火炮發(fā)射環(huán)境下彈載光電器件的動態(tài)響應(yīng)，以某型特種彈所用彈載光電器件為例，采用LS-DYNA有限元分析軟件對該器件及其減載組件進行建模仿真，并與采用應(yīng)力波理論得到的數(shù)值計算結(jié)果進行了對比，驗證理論分析方法的正確性。

1 彈載光電器件的動態(tài)響應(yīng)數(shù)值仿真

1.1 載荷環(huán)境的確定

不同口徑的身管火炮，采用不同發(fā)射裝藥時，彈丸在膛內(nèi)的加速度隨時間變化曲線均存在著較大不同［2-4］。為研究問題的方便，以某型火炮2號裝藥為例，將彈丸發(fā)射時在膛內(nèi)的加速度時程曲線擬合為折線，如圖1所示，過載持續(xù)時間為10.5 ms、峰值為10 000 g，峰值對應(yīng)的時刻為4.1 ms.

1.2 有限元模型的建立和分析

研究所用的彈載光電器件為某型特種彈所用的光學(xué)鏡頭，鏡頭旋擰在由鋁合金制成的固定座上，減載組件采用多片內(nèi)心中空的圓形橡膠墊疊合組成。鏡頭及其固定座、減載組件置于鋁質(zhì)套筒內(nèi)，并采用環(huán)形扣板固定。為防止反向沖擊，在鏡頭固定座和扣板之間加裝一片橡膠墊。在有限元分析中，鏡片采用光學(xué)玻璃材料模型，鏡頭其他的部分采用工程塑料模型，鏡頭固定座、鋁制套筒和扣板采用彈性體材料模型，橡膠采用Mooney-Rivlin不可壓縮模型，建立的加裝減載組件的彈載光電器件有限元分析模型如圖2、圖3所示。

圖1 火炮發(fā)射環(huán)境下彈丸在膛內(nèi)的加速度時程曲線Fig.1 Acceleration vs.time

圖2 彈載光電鏡頭及其固定座有限元模型Fig.2 Finite element model of projectile-borne device

經(jīng)仿真計算，試件在過載持續(xù)時間10.5 ms、峰值10 000 g的沖擊載荷作用下，結(jié)構(gòu)變形如圖4所示。從圖4中可以看出，減載組件在第4.3 ms時刻壓縮量達(dá)到最大，較載荷峰值有一定的延遲，壓縮量最大值為4.57 mm.之后減載組件開始反彈，直至鏡頭底座上的橡膠墊撞到套筒扣板。橡膠墊反彈后，基本能夠恢復(fù)原貌。

從仿真分析結(jié)果可知:在載荷作用時間為10.5 ms、峰值為10 000 g的火炮發(fā)射環(huán)境下，彈載光電器件和減載組件之間存在相互運動。這是因為在加載力開始作用的時候，試件中的減載組件先被壓縮，彈載光電器件與套筒扣板分離，這時在彈載光電器件端面形成了自由面，應(yīng)力即被卸載。在整個加載過程中，彈載光電器件上應(yīng)力不斷被加載和卸載，從而形成振蕩運動，當(dāng)加載力由大變小并開始衰減時，被壓縮的減載組件逐漸復(fù)原，直至彈載光電器件與套筒扣板發(fā)生撞擊，此時作用于鏡片上的應(yīng)力會突然增大。

圖4 減載組件變化時間圖（10.5 ms，10 000 g）Fig.4 Load-relieving subassembly vs.time（10.5 ms，10 000 g）

2 理論模型的建立和數(shù)值計算

2.1 理論模型的建立

假設(shè)減載組件的截面受力不變形，彈載光電器件以易碎的鏡頭表征其破壞極限，當(dāng)鏡頭所受應(yīng)力達(dá)到破壞極限時即視為失效。建立如圖5所示的彈載光電器件動態(tài)響應(yīng)等效模型，A、B、C、D桿分別表示固體本件、減載組件、彈載光電器件和套筒，A和D固聯(lián)，視為同步運動。設(shè)作用于A桿的外載為σ（t），B桿的長度為L、截面積為S、密度為ρ、應(yīng)力波速為c，C桿為剛體，質(zhì)量為m，C桿與D桿可接觸可分離。

圖5 應(yīng)力波的特征線Fig.5 Characteristic lines of stress wave

2.2 理論模型的求解

由于外載σ（t）是持續(xù)變化的，在彈載光電器件內(nèi)部也會出現(xiàn)持續(xù)變化的應(yīng)力波傳播，而傳播過程應(yīng)滿足應(yīng)力波的相容關(guān)系和結(jié)構(gòu)本身的邊界條件，這樣只需知道初值條件，結(jié)合相容關(guān)系和邊界條件，即可逐一求解任意時刻的應(yīng)力狀態(tài)，從而得到彈載光電器件的動力學(xué)響應(yīng)。

2.2.1 初值條件和邊界條件

如圖5所示，在B桿內(nèi)，其左右兩端滿足兩個不同的邊界條件，在AB截面時，應(yīng)力值為輸入應(yīng)力σ（t），在BC截面時，應(yīng)滿足介質(zhì)連續(xù)條件和受力平衡的邊界條件［5］:

式中:v為質(zhì)點速度，下標(biāo)“+”表示應(yīng)力波前狀態(tài)，下標(biāo)“-”表示應(yīng)力波后狀態(tài)。

隨時間變化的σ（t）會使應(yīng)力波在傳播中持續(xù)變化，為了求解應(yīng)力波的傳播過程，按圖5所示，將計算區(qū)域離散分化成若干個子區(qū)域，每一個子區(qū)域內(nèi)部的應(yīng)力狀態(tài)皆相同，相鄰子區(qū)域滿足相容關(guān)系。為得到σ，將（1）式展開，區(qū)域0為初始狀態(tài)，將每一個時間周期劃分為k等分，以區(qū)域編號ij為例，i∈｛1，2，3，…｝，j∈｛1，2，3，…，k+1｝，子區(qū)域ij滿足如下的相容關(guān)系:

可得彈載光電器件滿足的邊界條件差分格式為

2.2.2 數(shù)值計算方法

理論上，選取的k值越大，代表劃分的越細(xì)致，數(shù)值計算也就越精確，對于線彈性體，求解（3）式、（4）式，即可以求出整個計算區(qū)域的應(yīng)力狀態(tài)。而對于粘彈性材料，應(yīng)變率效應(yīng)不能忽略，具體來說，對于粘彈性材料的動力學(xué)響應(yīng)，應(yīng)滿足連續(xù)方程、運動方程和材料的本構(gòu)方程，構(gòu)建其差分格式如下:

2.2.3 計算程序設(shè)計

基于上述離散求解方法可知:目標(biāo)子區(qū)域的應(yīng)力狀態(tài)，都依賴于其周圍子區(qū)域的應(yīng)力狀態(tài)，而任意子區(qū)域的應(yīng)力狀態(tài)，也只會影響其相鄰的4個子區(qū)域，介于這樣的力學(xué)性質(zhì)，通過構(gòu)建差分方程組（5）式～（7）式，并使用FORTRAN90編寫基于Lagrange坐標(biāo)下的數(shù)值計算程序來求解，即可得到作用于彈載光電器件的應(yīng)力σ、應(yīng)變ε和質(zhì)點速度v，也就知道了火炮發(fā)射環(huán)境下彈載光電器件的動力學(xué)響應(yīng)。

3 數(shù)值計算與仿真結(jié)果的對比

按照彈丸實際結(jié)構(gòu)和載荷環(huán)境，將加載力、彈載光電器件及其減載器件尺寸參數(shù)等代入計算方程進行求解，采用數(shù)值計算與仿真模擬兩種方法，分別得到作用于鏡頭上的最大應(yīng)力和減載組件最大壓縮量結(jié)果見表1.數(shù)值計算時，外載σ（t）將采用圖1簡化加載條件，將實際加載環(huán)境簡化為3條線段的組合，4個端點坐標(biāo)為（0 ms，0 g）、（4.3 ms，10 000 g）、（10 ms，3 700 g）以及（10.5 ms，0 g），以此4個端點組成的3條線段，即為外載σ（t）的函數(shù)值曲線。

表1 計算結(jié)果對比Tab.1 Simulated and calculated results

圖6是數(shù)值計算和仿真模擬得到的鏡頭相對位移、相對速度曲線對比，從圖中可以看出，兩種不同求解方法得到的計算結(jié)果一致性較好，但也存在一定誤差，主要來源于參數(shù)的設(shè)置和算法的區(qū)別。此外，從圖6（b）中可以發(fā)現(xiàn)，在加載力開始作用時，橡膠墊振蕩壓縮，鏡頭相對速度為負(fù)值，在4.3 ms時刻，即橡膠墊壓縮量達(dá)到最大后，橡膠墊振蕩復(fù)原，鏡頭相對速度由負(fù)值變?yōu)檎?，振幅與加載力的變化率直接相關(guān)，而鏡頭在運動過程中不斷振蕩，是由于鏡頭受到持續(xù)不斷的加載波和卸載波作用。當(dāng)加載力結(jié)束后，鏡頭的相對速度出現(xiàn)劇烈變化，是由于鏡頭反彈使鏡頭底座上的橡膠墊與套筒扣板碰撞所致。

圖6 數(shù)值計算和仿真模擬結(jié)果對比Fig.6 Calculated and simulated relative velocities and displacements of camera

4 結(jié)論

1）對火炮發(fā)射環(huán)境下彈載光電器件動態(tài)響應(yīng)進行研究，是彈載器件抗過載機理研究的一項重要內(nèi)容。本文以某型特種彈為例，采用應(yīng)力波理論，建立了彈載光電器件的動力學(xué)分析模型，數(shù)值計算與仿真模擬結(jié)果具有較好的一致性，驗證了理論分析方法的正確性，為火炮發(fā)射環(huán)境下彈載器件的動態(tài)特性研究提供了一種有效的理論分析和計算方法。

2）通過仿真模擬和數(shù)值計算可知，加裝減載組件的彈載光電器件在彈丸發(fā)射過程中，在應(yīng)力波作用下做持續(xù)振蕩運動，振蕩幅值會隨著加載力的變化率增大而增大，而且在加載力結(jié)束后，鏡頭會出現(xiàn)較強反彈。此外，通過本文建立的理論分析模型，可在已知減載組件材料屬性、幾何尺寸和彈載器件質(zhì)量等參數(shù)情況下，即可通過求解方程，得到彈載器件在火炮發(fā)射環(huán)境下的運動和受力情況，為不同類型彈載器件的抗高過載設(shè)計提供參考和借鑒。

（

）

［1］ 錢立志.信息化炮彈技術(shù)［M］.北京:解放軍出版社，2006: 24-25. QIAN Li-zhi.Informatization artillery technology［M］.Beijing: Chinese People's Liberation Army Publishing House，2006:24-25.（in Chinese）

［2］ 金志明.槍炮內(nèi)彈道學(xué)［M］.北京:北京理工大學(xué)出版社，2004:35-36. JIN Zhi-ming.Gun ballistics［M］.Beijing:Beijing Institute of Technology Press，2004:35-36.（in Chinese）

［3］ 寧全利.高過載環(huán)境下彈載光電器件結(jié)構(gòu)動態(tài)響應(yīng)研究［D］.合肥:陸軍軍官學(xué)院，2014:35-36. NING Quan-li.Research on structural dynamical response of projectile-based equipment under high overload［D］.Hefei:Army Officer Academy of PLA，2014:35-36.（in Chinese）

［4］ Jayasooriya R，Thambiratnam D P，Perera N J，et al.Blast and residual capacity analysis of reinforced concrete framed buildings［J］.Engineering Structures，2011，33（12）:3483-3495.

［5］ 王禮立，任輝啟，虞吉林，等.非線性應(yīng)力波傳播理論的發(fā)展及應(yīng)用［J］.固體力學(xué)學(xué)報，2013，34（3）:217-240. WANG Li-li，REN Hui-qi，YU Ji-lin，et al.Development and application of the theory of nonlinear stress wave propagation［J］. Chinese Journal of Solid Mechanics，2013，34（3）:217-240.（in Chinese）

The Structural Dynamic Response of Projectile-borne Optoelectronic Devices during Launching

NING Quan-li1，JIANG Bin-an1，LIANG Xi2，LI Jun1，3
（1.Laboratory of Guidance Control and Information Perception Technology of High Overload Projectiles，Army Officer Academy of PLA，Hefei 230031，Anhui，China；2.College of Civil Engineering&Mechanics，Yanshan University，Qinhuangdao 066004，Hebei，China；3.Department 2，Army Officer Academy of PLA，Hefei 230031，Anhui，China）

In order to study the anti-high-overload mechanism of projectile-borne optoelectronic devices，a numerical simulation model of projectile-borne opto-electronic devices is established according to the artillery launching environmental characteristics.A theoretical analysis model of projectile-borne optoelectronic devices is also established to solve the stress，strain and relative velocity of projectile-borne optoelectronic devices.The result shows that the numerically calculated results are in good agreement with the simulated results to verify the correctness of the theoretical analysis model.The proposed model can be used for the reseasch on the dynamic characteristics of the projectile-borne devicess and the theoretical analysis of anti-high-overload design.

ordnance science and technology；projectile-borne optoelectronic device；numerical simulation；stress wave；artillery launching

TJ410.3

A

1000-1093（2015）11-2185-05

10.3969/j.issn.1000-1093.2015.11.025

2015-01-18

國家自然科學(xué)基金項目（11172331）

寧全利（1978—），男，副教授。E-mail:522736949@qq.com