馬帥，張相炎，游修東，2

(1.南京理工大學(xué)，江蘇 南京 210094;2.中國(guó)人民解放軍91388 部隊(duì)，廣東 湛江 524022)

0 引言

彈帶擠進(jìn)是一個(gè)瞬態(tài)過(guò)程，持續(xù)時(shí)間短，彈丸運(yùn)動(dòng)行程也不大，但其對(duì)內(nèi)彈道的影響卻很大。傳統(tǒng)的實(shí)驗(yàn)方法，只能通過(guò)分析擠進(jìn)結(jié)束后的彈帶或者通過(guò)其他間接物理量來(lái)推測(cè)擠進(jìn)過(guò)程的一些參數(shù)，這要求實(shí)驗(yàn)者有豐富的經(jīng)驗(yàn)。隨著計(jì)算機(jī)技術(shù)的發(fā)展，現(xiàn)在的有限元仿真已經(jīng)能夠模擬微秒甚至更短時(shí)間內(nèi)的動(dòng)力學(xué)過(guò)程，利用有限元軟件，可以獲得擠進(jìn)過(guò)程的完整的動(dòng)力學(xué)響應(yīng)。

1 計(jì)算模型

1.1 有限元模型［1］

現(xiàn)代的火炮膛線為了減小磨損和擠進(jìn)的沖擊，一般采用漸速膛線，而本文研究的是彈帶的擠進(jìn)過(guò)程，該過(guò)程彈丸行程不大，因此忽略膛線纏角的變化，取膛線的纏角為固定值。因該過(guò)程只涉及到身管很有限的一部分，為了節(jié)省計(jì)算量，只對(duì)身管的一部分進(jìn)行建模，擠進(jìn)的有限元模型見(jiàn)圖1。

圖1 擠進(jìn)過(guò)程的有限元模型

整個(gè)模型分為彈丸、彈帶和身管三個(gè)部件，這三部分都是采用八結(jié)點(diǎn)的六面體單元SOLID164 來(lái)劃分的。這種單元生成剛度矩陣時(shí)采用的是單點(diǎn)積分，采用單點(diǎn)積分可能會(huì)發(fā)生沙漏模式，為了防止沙漏模式，需要打開(kāi)沙漏控制。實(shí)踐已經(jīng)證明，單點(diǎn)積分在顯示動(dòng)力學(xué)分析中有很好的精度，同時(shí)大大減少了計(jì)算量。

在身管的后端面加固定約束，這與實(shí)際發(fā)射過(guò)程相符，彈丸和彈帶通過(guò)共結(jié)點(diǎn)的方式連接在一起，這種連接方式是對(duì)實(shí)際的一種很好的近似，而且可以大大減少計(jì)算量。為了簡(jiǎn)化處理，只在彈底加壓力，由于彈丸擠進(jìn)過(guò)程中擺動(dòng)很小，這也是對(duì)實(shí)際的一種很好的近似。由于要模擬彈丸的旋轉(zhuǎn)，載荷曲線跟實(shí)際越相近越好，本文用線性增加的壓力來(lái)近似。

由于彈丸的擺動(dòng)，擠進(jìn)過(guò)程中會(huì)發(fā)生彈丸和身管的接觸以及彈帶和身管的接觸，定義彈丸和身管的接觸為自動(dòng)接觸，定義彈帶和身管的接觸為侵蝕接觸，這種接觸曾普遍用于彈丸的侵徹等。彈丸從0 開(kāi)始加速，且速度變化范圍較大，在擠進(jìn)結(jié)束時(shí)可達(dá)100 m/s 以上，因此不能將摩擦系數(shù)定義為常數(shù)，LS－DYNA 采用的是指數(shù)衰減的摩擦定律，見(jiàn)式(1):

其中:fd為動(dòng)摩擦系數(shù)，fs為靜摩擦系數(shù)，Dc為衰減系數(shù)，vrel為相對(duì)滑動(dòng)速度。

1.2 材料模型

彈丸采用LS－DYNA 中的剛體材料模型，剛體材料模型需給定密度，彈性模量和泊松比，這是為了計(jì)算彈丸的慣性并準(zhǔn)確地計(jì)算接觸。彈帶材料選擇的是紫銅，采用的是Johson－Cook 材料模型，該材料模型綜合考慮了應(yīng)變、應(yīng)變率和絕熱溫升對(duì)應(yīng)力的影響，非常適合模擬瞬態(tài)的動(dòng)力學(xué)過(guò)程。身管材料選某炮鋼，采用等向強(qiáng)化的彈塑性材料模型。Johson－Cook 材料模型見(jiàn)式2:

其中:A，B，C，n，m 均為材料常數(shù)，εp為有效塑性應(yīng)變，ε*為相對(duì)應(yīng)變，T=(T－ Troom)/(Tmelt－ Troom)為相對(duì)溫度變化。

對(duì)于材料的失效，LS－DYNA 采用累積損傷模型，累積損傷模型定義一個(gè)損傷變量D，初始時(shí)D=0，當(dāng)D=1時(shí)材料失效，單元被刪除，D 的定義見(jiàn)式(3):

其中:Δεp為塑性應(yīng)變?cè)隽浚舊為材料發(fā)生失效時(shí)的應(yīng)變，見(jiàn)式(4):

其中:D1～D5皆為材料常數(shù)，一般通過(guò)實(shí)驗(yàn)測(cè)定，σ*=p/σeff為應(yīng)力三軸度，p 為由狀態(tài)方程求得的壓力，σeff為Von mises 有效應(yīng)力。

2 擠進(jìn)過(guò)程仿真結(jié)果及分析

2.1 擠進(jìn)過(guò)程結(jié)果

圖2－圖7 顯示了擠進(jìn)過(guò)程中彈帶的應(yīng)力云圖。

圖2 擠進(jìn)前:σeff/105 MPa

圖3 坡膛擠壓第一條彈帶，未刻槽:σeff/105 MPa

圖4 第一條彈帶開(kāi)始刻槽:σeff/105 MPa

圖5 坡膛擠壓第二條彈帶:σeff/105 MPa

圖6 第二條彈帶開(kāi)始刻槽:σeff/105 MPa

圖7 擠進(jìn)完成:σeff/105 MPa

圖2－圖7 展示了彈帶擠進(jìn)的詳細(xì)過(guò)程，但所刻槽的傾角不太明顯，這是由于彈丸的慣性比較大，其轉(zhuǎn)動(dòng)有一個(gè)滯后效應(yīng)。由圖4 與圖6 可以看出，彈帶的最大應(yīng)力位于擠進(jìn)所刻槽的前部，這表明膛線和彈帶剛發(fā)生接觸即處于較高的應(yīng)力狀態(tài)，這種應(yīng)力沖擊可能是膛線磨損的一個(gè)原因。因?yàn)椴牧现圃爝^(guò)程中不可避免的會(huì)有一些缺陷，如微裂紋、微孔洞等，在強(qiáng)的剪切應(yīng)力下或受到應(yīng)力沖擊，極易導(dǎo)致裂紋擴(kuò)展，使材料發(fā)生破壞或磨損［2］。

2.2 單元應(yīng)力結(jié)果

圖8 為某些單元的應(yīng)力曲線。由于采用的是單點(diǎn)積分單元，該種單元為常應(yīng)力單元，因此用LS－DYNA 只能得到單元應(yīng)力結(jié)果，而非結(jié)點(diǎn)應(yīng)力。曲線A 和曲線B 可以看出兩條彈帶對(duì)單元應(yīng)力的影響，曲線都有兩個(gè)明顯的峰，由于擠進(jìn)速度隨時(shí)間在增加，第一個(gè)峰在時(shí)間上寬一些。曲線A 和曲線B 所代表的兩個(gè)單元材料相同，但A單元的應(yīng)力卻比B 單元大很多。通常認(rèn)為，應(yīng)力大的地方磨損也較快，A 單元位于膛線起始部，這可以解釋為什么膛線起始部磨損較快［3］。曲線C 為彈帶上的單元，該單元在擠進(jìn)開(kāi)始后應(yīng)力幾乎不變，表明該處的材料處于塑性流動(dòng)狀態(tài)，符合對(duì)銅這種材料的認(rèn)知。以上分析表明，本文的應(yīng)力計(jì)算結(jié)果符合實(shí)際，可以為實(shí)驗(yàn)研究提供理論依據(jù)。

圖8 單元的von mises 應(yīng)力－時(shí)間曲線:σeff/105Mpa

3 擠進(jìn)阻力的計(jì)算及分析

圖9 和圖10 分別為彈丸加速度和位移曲線。由達(dá)朗貝爾原理可知，物體所受的力和它的慣性力是平衡力系。對(duì)彈丸做軸線方向的受力分析，彈丸受到彈底壓力，擠進(jìn)阻力和慣性力，利用matlab 對(duì)數(shù)據(jù)進(jìn)行處理，便可求出擠進(jìn)阻力隨時(shí)間的變化，再結(jié)合圖10 位移隨時(shí)間的變化曲線，即可得擠進(jìn)阻力隨擠進(jìn)行程的曲線，見(jiàn)圖11。

圖9 彈丸加速度曲線

圖10 彈丸位移曲線

圖12 為文獻(xiàn)［4］求得的擠進(jìn)阻力曲線，張浩等通過(guò)對(duì)實(shí)際物理過(guò)程所做的一系列假設(shè)，對(duì)擠進(jìn)過(guò)程進(jìn)行了理論分析，分階段得出了擠進(jìn)阻力的公式。

圖11 擠進(jìn)阻力隨擠進(jìn)行程的曲線

圖12 文獻(xiàn)［4］得到的理論解

將圖11 中第一條彈帶的擠進(jìn)阻力曲線(擠進(jìn)行程在0.014 m 以前的部分)與圖12 進(jìn)行對(duì)比，可以發(fā)現(xiàn)兩者的變化趨勢(shì)基本相同，兩者都是在彈帶與坡膛發(fā)生接觸到開(kāi)始刻槽前的這一段位移中阻力急劇上升，在外側(cè)刻槽完畢后阻力開(kāi)始急劇下降，兩者之間是一個(gè)阻力緩慢變化區(qū)間。但文獻(xiàn)［4］是一小口徑火炮的結(jié)果，而本文是對(duì)某一大口徑火炮求得結(jié)果，兩者數(shù)值顯然不同。通過(guò)兩條曲線的對(duì)比，表明本文的曲線可以體現(xiàn)擠進(jìn)過(guò)程中擠進(jìn)阻力的變化規(guī)律，而且有限元分析可以得到雙彈帶對(duì)擠進(jìn)阻力的影響。

4 結(jié)語(yǔ)

基于ANSYS/LS－ DYNA 對(duì)彈帶擠進(jìn)過(guò)程進(jìn)行了仿真，通過(guò)輸出關(guān)鍵時(shí)刻的圖形，展示了整個(gè)擠進(jìn)過(guò)程中彈帶的變形規(guī)律，然后通過(guò)分析模型中某些單元的應(yīng)力曲線，對(duì)膛線起始部的磨損提出了看法。通過(guò)仿真獲得的加速度曲線和速度曲線，基于達(dá)朗貝爾原理，利用matlab 對(duì)數(shù)據(jù)進(jìn)行處理，得到了擠進(jìn)阻力曲線，通過(guò)與文獻(xiàn)［4］阻力曲線的對(duì)比，驗(yàn)證了文中計(jì)算結(jié)果的正確性，表明ANSYS/LS－DYNA 是研究擠進(jìn)過(guò)程的一個(gè)很好的工具。

［1］談樂(lè)斌，張相炎，管紅根，等.火炮概論［M］.北京:北京理工大學(xué)出版社，2005:49-57.

［2］孫河洋，馬吉?jiǎng)?，李偉，?坡膛結(jié)構(gòu)變化對(duì)彈帶擠進(jìn)過(guò)程影響的研究［J］.振動(dòng)與沖擊，2011，30(3):30-33.

［3］米糧川，高樹(shù)滋，賴長(zhǎng)纓.內(nèi)彈道彈帶擠進(jìn)過(guò)程仿真研究［J］.彈箭與制導(dǎo)學(xué)報(bào)，2012(2)，32(1):133-135.

［4］張浩，周彥煌.埋頭彈火炮擠進(jìn)過(guò)程研究［J］.彈道學(xué)報(bào)，2006，18(1):76-79.