壓電引信大著角發(fā)火性數(shù)值仿真方法

周 浪，王雨時，聞 泉，王光宇，張志彪

(南京理工大學(xué)機械工程學(xué)院，江蘇 南京 210094)

0 引言

壓電引信是現(xiàn)代彈藥中廣泛使用的觸發(fā)式引信之一，其特點是具有極高的瞬發(fā)度，能夠減少彈藥因著速不同造成的炸高散布[1]。為了確保壓電引信在實戰(zhàn)時能夠可靠作用，研制時需要對其作用可靠性進行考核。對于壓電引信工作性能的考核，傳統(tǒng)的方法是采用靶場試驗的方法，但是該方法需要消耗大量的人力和財力，并且難以深入分析試驗過程出現(xiàn)的故障原因和引信的具體失效過程。對于壓電引信的研究，現(xiàn)有文獻都是基于理論和試驗兩方面，通過數(shù)值仿真手段對其進行研究的文獻較少，僅有文獻[2]應(yīng)用Autodyn仿真軟件對某空氣炮彈碰目標過程進行模擬，再現(xiàn)了其彈頭受力、變形過程，獲得了彈頭中壓力、應(yīng)力、速度等各種力學(xué)參數(shù)的分布和時程圖，從而間接獲知壓電陶瓷輸出電能性能，但未從彈頭結(jié)構(gòu)差異方面對彈頭中壓力、應(yīng)力、速度等各種力學(xué)參數(shù)的影響進行深入研究。文獻[3]根據(jù)應(yīng)力波的基本理論近似地分析計算了彈丸碰擊鋼甲時，頭部碰擊式壓電引信壓電機構(gòu)中應(yīng)力波的傳播過程，得出了不同結(jié)構(gòu)中壓電陶瓷軸向應(yīng)力的近似值，在計算過程中將引信頭部結(jié)構(gòu)簡化成平頭，與實際引信頭結(jié)構(gòu)差別較大，計算結(jié)果精度不高。文獻[4]以壓電引信電-1和電-2為例，討論了壓電引信發(fā)火和瞬發(fā)度的影響因素，并以大量的試驗數(shù)據(jù)為支撐，得出影響壓電引信瞬發(fā)度的因素主要包括彈速、目標、著角、彈頭部即引信頭部結(jié)構(gòu)等，同時得出壓電晶體老化不會影響引信發(fā)火可靠性，但是對電雷管的起爆時間會產(chǎn)生1～2 μs的延遲。文獻[5]經(jīng)推導(dǎo)壓電陶瓷壓電方程，得出影響壓電陶瓷壓電效應(yīng)的應(yīng)力包括三種正應(yīng)力和兩種剪應(yīng)力，其中與壓電陶瓷壓電效應(yīng)直接相關(guān)的應(yīng)力是與壓電晶體極化方向一致的正應(yīng)力。文獻[6]對裝有壓電引信的低速破甲彈破甲不穩(wěn)定的問題進行研究，分析得出影響破甲彈破甲不穩(wěn)定的主要原因是目標(或靶板)接地電阻不同，降低了壓電引信瞬發(fā)度。本文針對試驗過程中壓電引信大著角對靶射擊試驗過程中瞎火原因難以分析的問題，采用ANSYS/LS-DYNA仿真軟件對引信頭采用不同結(jié)構(gòu)和不同材料以不同著角碰擊均質(zhì)鋼靶板時引信頭部變形對壓電晶體所產(chǎn)生的應(yīng)力進行數(shù)值模擬，得出不同引信結(jié)構(gòu)以及同種結(jié)構(gòu)不同引信頭材料以不同著角碰擊均質(zhì)鋼板時壓電晶體上的應(yīng)力變化規(guī)律，藉此對引信頭部結(jié)構(gòu)和參數(shù)設(shè)計進行優(yōu)化，從而為提高引信大著角發(fā)火正確性探討可行途徑。

1 碰擊式壓電機構(gòu)結(jié)構(gòu)形式和壓電陶瓷壓電方程

根據(jù)配用彈種的不同，為了解決大著角發(fā)火與低靈敏度的矛盾，碰擊式壓電機構(gòu)有各種不同的結(jié)構(gòu)形式[5]。其中較為典型的碰擊式壓電機構(gòu)包括帶環(huán)形槽的壓電機構(gòu)、帶半球形壓電塊和防潮帽的壓電機構(gòu)、整體塑封的壓電機構(gòu)、帶碰擊開關(guān)的壓電機構(gòu)、帶介質(zhì)的壓電機構(gòu)、帶球殼形陶瓷的壓電機構(gòu)、用火帽沖力壓電的壓電機構(gòu)和壓電陶瓷配置在彈體圓柱部的壓電機構(gòu)8種[5]。某單兵火箭彈采用的碰擊式壓電機構(gòu)屬于帶環(huán)形槽的壓電機構(gòu)，其結(jié)構(gòu)如圖1所示。

該壓電機構(gòu)在引信頭上車制有環(huán)形槽，形成厚度為h的薄弱環(huán)節(jié)，以調(diào)節(jié)引信頭變形的難易程度[5]。陶瓷盒與引信頭之間有一間隙e，引信頭的軸向變形大于e時，才能向陶瓷施壓。在大著角時為了避免壓電陶瓷側(cè)壁局部受力使其提前破碎，陶瓷盒側(cè)面與引信頭之間也留有適當(dāng)?shù)膹较蜷g隙，此徑向間隙一般為2 mm左右[5]。也可以在陶瓷盒里支承壓電陶瓷的陶瓷座上部做成錐孔(如圖1所示)，以保證有一定的徑向間隙，這樣對大著角發(fā)火也有利[5]。

在彈種、電雷管和壓電陶瓷參數(shù)確定以后，主要是通過改變引信頭薄弱部厚度h、間隙e、引信頭的材料和頂部形狀等來同時滿足大著角發(fā)火和低靈敏度的要求[5]。

對于壓電引信來說，壓電陶瓷性能和引信頭結(jié)構(gòu)是決定引信發(fā)火正確性的關(guān)鍵。圖2為從壓電陶瓷中分離出來的單元立方體。圖3為壓電陶瓷受力和極化方向示意圖。

圖1 帶環(huán)形槽的引信壓電機構(gòu)Fig.1 The piezoelectric structure with ring-shaped groove within fuze

圖2 壓電陶瓷單元三維應(yīng)力分布示意圖Fig.2 The scheme of three-dimensional stress

圖3 壓電陶瓷受力和極化方向示意圖Fig.3 The scheme of external force and polari-distribution onpiezoelectric ceramics element zation direction of piezoelectric ceramics

從圖2可看出，作用在壓電陶瓷立方體單元各個面上的應(yīng)力共有六個，其中包括三個正應(yīng)力σ1，σ2，σ3和三個剪應(yīng)力σ4，σ5，σ6。從文獻[5]對壓電陶瓷壓電方程的推導(dǎo)可得出：在壓電陶瓷中，壓電效應(yīng)與陶瓷中極化狀態(tài)的改變有關(guān)，而不是在任何力的作用下在任何方向都能產(chǎn)生壓電效應(yīng)，只有在某些力的作用下在某些特定方向上才能產(chǎn)生壓電效應(yīng)。具體地說，就是在正應(yīng)力σ1，σ2，σ3的作用下，只有z方向，即原來的極化方向才有壓電效應(yīng)；而在剪應(yīng)力σ4，σ5的作用下，只在與極化方向垂直的方向上才有壓電效應(yīng)[5]。

2 數(shù)值仿真方法

2.1 數(shù)值模擬研究思路

采用ANSYS/LS-DYNA軟件，對該單兵火箭彈以不同著角碰擊均質(zhì)鋼靶板時引信頭部變形對壓電晶體作用過程進行數(shù)值仿真，得到不同結(jié)構(gòu)、不同材料狀態(tài)下的引信頭以不同著角碰擊均質(zhì)鋼靶板時引信壓電晶體表面法向最大應(yīng)力峰值平均值和等效綜合最大應(yīng)力峰值平均值。通過比較兩種應(yīng)力峰值在不同碰擊著角情況下平均應(yīng)力值，得出火箭彈以不同著角碰擊均質(zhì)鋼靶板時壓電晶體表面應(yīng)力的變化規(guī)律。由于不同引信結(jié)構(gòu)、不同結(jié)構(gòu)材料在碰擊過程中所產(chǎn)生的應(yīng)力不同，傳到壓電晶體中的應(yīng)力也就不同，所以對所測得的壓電晶體表面應(yīng)力值進行對比，就可以對引信頭部結(jié)構(gòu)進行優(yōu)化并對最佳結(jié)構(gòu)材料進行優(yōu)選，使得引信頭部結(jié)構(gòu)所產(chǎn)生的應(yīng)力能夠更好地傳遞到壓電晶體上。

2.2 仿真模型建立

已知某單兵火箭彈參數(shù)：彈重2.25 kg；初速120±10 m/s，發(fā)動機工作后，在70 m處著靶速度295 m/s；彈丸最大直徑Dmax=80 mm；彈丸長度924.43～929.5 mm，飛行時彈丸長度920.57～929.5 mm；彈丸重心距彈底409 mm。

該火箭彈主要由引信彈頭壓電部分、藥柱、彈體、尾翼和引信彈底起爆部分組成。由于主要研究對象為引信彈頭壓電部分，所以為了三維建模和重心位置調(diào)整方便，把藥柱、彈體、尾翼和引信彈底起爆部分簡化成一個整體——彈體，而引信彈頭壓電部分則按產(chǎn)品具體結(jié)構(gòu)簡化為引信頭、壓電陶瓷、陶瓷盒、絕緣座、絕緣墊、加強板、電極板，共七部分。仿真火箭彈垂直碰擊靶板過程簡化后的彈丸整體模型和引信彈頭壓電部分結(jié)構(gòu)模型如圖4所示。

圖4 火箭彈垂直碰擊靶板三維簡化模型Fig.4 The simplified three-dimensional finite element model of rocket hitting target at verticality

為了簡化有限元模型，提高計算效率，作如下假設(shè)[2]：

1)對引信零件做簡化處理，保留其外形特征，去掉倒角。絕緣墊簡化為一形狀規(guī)則的階梯形空心柱體；

2)對簡化后的整體彈丸模型進行配重處理，使其與實體質(zhì)量一致；

3)彈丸和靶板均為連續(xù)介質(zhì)；

4)彈丸碰擊靶板的過程為絕熱過程，不考慮熱效應(yīng)；

5)忽略重力加速度、空氣阻力以及靶板的整體運動；

6)忽略尾翼碰擊靶板時產(chǎn)生的沖擊；

7)彈丸和靶板初始應(yīng)力均為零。

建模時，為了節(jié)省計算時間，對彈丸和靶板取1/2模型，靶板厚度按實際靶板厚度(180 mm)進行建模，在模型對稱面施加反對稱約束，在靶板周圍定義無反射邊界，模型單元類型選用八節(jié)點六面體單元，網(wǎng)格劃分采用映射網(wǎng)格法，模型計算采用拉格朗日法，引信頭與靶板之間采用面-面侵蝕算法，其余各零件之間的接觸采用面-面自動接觸，模型計算單位為g·cm·μs?；鸺龔棿笾桥鰮舭邪宓挠邢拊Ｐ腿鐖D5所示。

圖5 火箭彈大著角碰擊靶板的有限元模型Fig.5 The finite element model of rocket hitting target at large impact angle

2.3 材料模型

彈體材料采用30CrMnSiA鋼(按質(zhì)量等效原則確定等效密度ρ)；引信頭、加強板和電極板材料為2A12鋁合金；陶瓷盒材料為T2；壓電陶瓷材料為GTQ-3；絕緣座材料為FX505壓塑料；絕緣墊材料為ABS樹脂；靶板材料為均質(zhì)鋼。絕緣座和絕緣墊采用PLASTIC_KINEMATIC材料模型；壓電晶體采用JOHNSON_HOLMQUIST_CERAMICS材料模型；彈體、引信頭、電極板、加強板、晶體盒和靶板采用Johnson-Cook材料模型。各材料模型主要仿真參數(shù)如表1、表2和表3所列。

表1 Johnson-Cook材料模型主要參數(shù)[2, 7-9]Tab.1 The parameters of Johnson-Cook constitutive model

表2 PLASTIC_KINEMATIC材料模型主要參數(shù)[8]Tab.2 The parameters of PLASTIC_KINEMATIC material model

表3 JOHNSON_HOLMQUIST_CERAMICS材料模型主要參數(shù)[10]Tab.3 The parameters of JOHNSON_HOLMQUIST_CERAMICS material model

3 仿真結(jié)果及分析

3.1 引信頭結(jié)構(gòu)對壓電晶體應(yīng)力的影響分析

火箭彈以100 m/s，200 m/s，295 m/s三種不同的初速分別以55°，63°，68°，73°著角碰擊180 mm厚均質(zhì)鋼靶板。為獲得不同引信頭結(jié)構(gòu)作用于壓電晶體上壓電晶體應(yīng)力的變化情況，擬定了圖6所示的四種不同引信頭結(jié)構(gòu)方案，四種不同引信頭結(jié)構(gòu)方案在仿真時所用材料都為2A12鋁合金。

圖6 四種不同結(jié)構(gòu)的引信頭有限元結(jié)構(gòu)模型Fig.6 The finite element models of fuze head with different structures

方案一為目前制式單兵火箭彈所采用的結(jié)構(gòu)，其引信頭頭部錐角為100°，見圖6(a)。

方案二在方案一基礎(chǔ)上修改，減小位置1處環(huán)形槽的寬度(由原3.2 mm減小為3.0 mm)，在位置2處加厚(由原0.6 mm增大到0.8 mm)，在位置3處減薄壁厚(由原0.95 mm減小到0.90 mm)，在位置4處加錐角，在位置5處減小凸臺直徑(由原14.4 mm減小到14.0 mm)并對位置5處凸臺圓弧加倒角(C0.5)，見圖6(b)。

方案三在方案二的基礎(chǔ)上修改，在圖中標識部位加圓弧過渡處理(過渡圓弧R2)，見圖6(c)。

方案四在方案一的基礎(chǔ)上修改，在圖中標識部位加圓角或倒角(R0.5 mm或C0.5)，見圖6(d)。

為了避免仿真過程中計測壓電晶體應(yīng)力時所選壓電晶體單元數(shù)量對仿真結(jié)果的影響，以圖6(a)所示的方案一引信頭結(jié)構(gòu)為研究對象、引信頭材料選用2A12鋁合金，彈丸以295 m/s的碰靶速度分別以55°，63°，68°，73°著角碰擊180 mm厚均質(zhì)鋼靶板。仿真計測應(yīng)力時所選壓電晶體單元方案分別如圖7(a)、圖7(b)和圖7(c)所示，所選壓電晶體單元為緊靠引信頭部壓電晶體表面上的單元，對應(yīng)的壓電晶體單元數(shù)量分別為11個、25個和38個。仿真所得三種方案中計測應(yīng)力所選壓電晶體單元上法向應(yīng)力峰值平均值和等效綜合應(yīng)力峰值平均值如表4所列。

圖7 仿真計測應(yīng)力時所選壓電晶體單元Fig.7 The chosen piezoelectric ceramics elements for stress measurement

表4 仿真計測應(yīng)力時所選不同壓電晶體單元上應(yīng)力峰值平均值Tab.4 The average of stress peak calculated for different piezoelectric ceramics elements

從表4可看出：三種方案的計測壓電晶體單元上的應(yīng)力平均值雖然在數(shù)值上存在一定差異，但都隨著角的增大而減小，仿真結(jié)果規(guī)律性一致。所選計測單元數(shù)量增多，雖然更能體現(xiàn)整個壓電晶體表面的應(yīng)力值，但是數(shù)據(jù)處理所耗時間長，因此在對結(jié)論不會產(chǎn)生影響的前提下，應(yīng)該選擇合適的單元數(shù)量對仿真結(jié)果進行處理。

影響計測應(yīng)力所選壓電晶體單元上應(yīng)力峰值準確性的主要因素是數(shù)據(jù)采樣時間間隔，而與仿真數(shù)據(jù)采樣時間間隔密切相關(guān)的是仿真計算結(jié)果輸出文件d3plot的輸出時間間隔，對應(yīng)的關(guān)鍵字為“DATABASE_BINARY_D3PLOT”。為了避免所設(shè)置的仿真計算結(jié)果輸出文件d3plot的輸出時間間隔參數(shù)影響仿真結(jié)果，分別取時間間隔0.01 μs，0.02 μs，0.05 μs，0.2 μs，0.4 μs，0.6 μs，0.8 μs和1.0 μs進行仿真。仿真模型采用圖6(a)所示的方案一引信頭結(jié)構(gòu)、材料選用2A12鋁合金，彈丸以295 m/s速度分別以55°，63°，68°，73°著角碰擊180 mm厚均質(zhì)鋼靶板。仿真中計測應(yīng)力所選壓電晶體單元選擇圖7(a)中方案，仿真所得計測應(yīng)力所選壓電晶體單元上法向應(yīng)力峰值平均值和等效綜合應(yīng)力峰值平均值如表5所列。

表5 不同時間間隔下仿真所得計測應(yīng)力時所選不同壓電晶體單元上應(yīng)力峰值平均值Tab.5 The average of stress peak of piezoelectric ceramics elements using different output frequency

從表5可看出：在相同速度、不同著角情況下，隨著數(shù)據(jù)采樣時間間隔的縮小，計測應(yīng)力時所選壓電晶體單元上的法向應(yīng)力平均值和等效綜合應(yīng)力平均值總體上呈現(xiàn)增大趨勢，當(dāng)所取時間間隔精度達到一定時，時間間隔精度不再影響仿真結(jié)果。在數(shù)據(jù)采樣時間間隔相同情況下，隨著著角的增大，計測應(yīng)力時所選壓電晶體單元上的法向應(yīng)力平均值和等效綜合應(yīng)力平均值逐漸減小。當(dāng)數(shù)據(jù)采樣時間間隔從0.05 μs，0.03 μs，0.01 μs變化過程中，所得仿真結(jié)果已經(jīng)趨于穩(wěn)定，其中0.05 μs與0.01 μs計算結(jié)果誤差在4%以內(nèi)，說明所取數(shù)據(jù)采樣時間間隔基本不會影響仿真結(jié)果，所取時間間隔為0.05 μs時能夠滿足計算精度要求。因此以下仿真中計測應(yīng)力時所選壓電晶體單元采用圖7(a)所示方案、數(shù)據(jù)采樣時間間隔為0.05 μs。

表6為上述四種不同引信頭結(jié)構(gòu)分別以55°，63°，68°，73°著角碰擊180 mm厚均質(zhì)鋼靶板時仿真得到的壓電晶體上各壓電晶體單元(對應(yīng)圖7(a)所示壓電晶體上所選的壓電晶體單元)的應(yīng)力峰值平均值。其中法向應(yīng)力峰值對應(yīng)壓電晶體單元上的正應(yīng)力σ3，而等效綜合應(yīng)力即Von Mises應(yīng)力是基于剪切應(yīng)變能的一種等效應(yīng)力。

表6 不同結(jié)構(gòu)方案以55°，63°，68°，73°著角碰擊均質(zhì)鋼時壓電晶體單元上的應(yīng)力峰值平均值Tab.6 The average of stress peak value of piezoelectric ceramics elements in different structural designs with various impact angles

由表6可看出：

1)在上述四種方案中，同種結(jié)構(gòu)的引信頭在相同速度下分別以55°，63°，68°，73°著角撞擊均質(zhì)鋼靶板時，壓電晶體單元上的法向應(yīng)力值和等效綜合應(yīng)力值都隨著角的增大而減小，與文獻[5]中結(jié)論一致，說明仿真是可信的。

2)同種結(jié)構(gòu)的引信頭結(jié)構(gòu)在相同著角下分別以100 m/s，200 m/s，295 m/s初速撞擊均質(zhì)鋼靶板時，壓電晶體單元上的法向應(yīng)力值和等效綜合應(yīng)力值都隨著角的增大而減小。

3)對比三種改進方案即方案二、方案三和方案四中的法向應(yīng)力值和等效綜合應(yīng)力值可以得出：三種改進方案在碰擊著角大于55°時壓電晶體單元法向應(yīng)力和等效綜合應(yīng)力值均比原結(jié)構(gòu)的要小，也就是說，通過改變引信頭部結(jié)構(gòu)難以有效改善在大著角情況下引信頭作用于壓電晶體上的應(yīng)力值，原結(jié)構(gòu)及其參數(shù)設(shè)計已經(jīng)處于優(yōu)化狀態(tài)。

4)在數(shù)值仿真中，壓電晶體單元法向應(yīng)力值和壓電晶體等效綜合應(yīng)力值都可以反映出壓電晶體所受應(yīng)力的變化規(guī)律。

3.2 引信頭錐角對壓電晶體應(yīng)力的影響分析

為驗證引信頭部錐角變化對壓電晶體單元應(yīng)力的影響，本次仿真在方案一結(jié)構(gòu)基礎(chǔ)上改變引信頭頭部錐角大小，分別以80°，90°，98°，99°，100°，101°，102°和110°共八種不同錐角進行數(shù)值仿真，仿真結(jié)果如表7所列。

表7 不同頭部錐角以55°，63°，68°，73°著角碰擊時壓電晶體單元上的應(yīng)力峰值平均值Tab.7 The average of stress peak value of piezoelectric ceramics elements in different head cones with various impact angles

從表7可看出，在大著角工況下，彈丸以相同著角、相同速度撞擊均質(zhì)鋼靶板時壓電晶體單元法向應(yīng)力以及等效綜合應(yīng)力均隨引信頭部錐角的增大而出現(xiàn)先增大后減小的趨勢；彈丸以相同頭部錐角、相同初速撞擊均質(zhì)鋼靶板時，壓電晶體單元法向應(yīng)力以及等效綜合應(yīng)力隨碰擊角度的增大而減?。辉陬^部錐角為100°時壓電晶體單元法向應(yīng)力及等效綜合應(yīng)力相對于其他錐角時的要大，也就是說引信頭部錐角存在最優(yōu)錐角值，使得彈以不同碰擊角撞擊目標時壓電晶體單元的應(yīng)力值達到最大。該頭部錐角即為原有結(jié)構(gòu)設(shè)計參數(shù)，即100°?；鸺龔椧?95 m/s初速大傾角撞擊均質(zhì)鋼靶板時壓電晶體的應(yīng)力隨撞擊角度的變化曲線如圖8和圖9所示。

圖8 碰擊著角與壓電晶體單元法向應(yīng)力和等效綜合應(yīng)力值曲線Fig.8 The relationship between impact angle and normal stress and equivalent stress of piezoelectric ceramics elements

圖9 引信頭部錐角與壓電晶體單元法向應(yīng)力和等效綜合應(yīng)力值曲線Fig.9 The relationship between head cone angle and normal stress and equivalent stress of piezoelectric ceramics elements

3.3 引信頭材料對壓電晶體應(yīng)力的影響分析

為驗證引信頭部材料對壓電晶體應(yīng)力的影響，在方案一所示引信頭結(jié)構(gòu)基礎(chǔ)上(引信頭部錐角為100°)，對頭部材料選用7種不同材料進行數(shù)值仿真，其中2A12鋁合金材料模型參數(shù)如表1所列，其余材料模型仿真參數(shù)如表8所列，仿真結(jié)果如表9所列。

表8 Johnson-Cook材料模型主要參數(shù)[2, 7, 8, 11-15]Tab.8 The parameters of Johnson-Cook constitutive model

表9 不同引信頭材料時壓電晶體單元應(yīng)力峰值平均值Tab.9 The average of stress peak of piezoelectric ceramics elements on fuze head with different materials

仿真結(jié)果表明，引信頭部材料為7A04時，在火箭彈初速為295 m/s并以不小于63°的碰擊著角碰擊均質(zhì)鋼靶板時，引信壓電晶體上的應(yīng)力值相比于其他幾種材料時的都要顯著偏大；同種材料的引信頭，碰擊速度對壓電晶體表面上的應(yīng)力衰減影響較大，碰靶初速越大，在大著角工況下，壓電晶體表面上應(yīng)力衰減量越??；不同材料的引信頭在相同速度下，2A11，2A12和7075材料引信頭對壓電晶體作用效果接近，但2A11和7075不如2A12；45鋼、Q235材料引信頭對壓電晶體的作用效果優(yōu)于2A12但不如7A04，且45鋼及Q235密度較大，不太適用于輕武器；鉛黃銅H69對壓電晶體的作用效果優(yōu)于紫銅T2，但兩者作用效果都不如2A12。從表9可看出，當(dāng)7A04引信頭以63°，68°，73°法向著角碰擊均質(zhì)鋼靶板時，壓電晶體上的應(yīng)力值衰減較小。而對于其他幾種材料，隨著著角逐漸增大壓電晶體上應(yīng)力值衰減較快，顯然對于大著角碰擊情形而言，7A04是壓電引信頭材料的最佳選擇。與原2A12相比，壓電晶體上的應(yīng)力至少增大了30%。

4 結(jié)論

本文提出了壓電引信大著角發(fā)火性數(shù)值仿真方法。該方法通過比較仿真得到的壓電晶體上不同位置處壓電晶體單元的法向應(yīng)力峰值平均值和等效綜合應(yīng)力峰值平均值來對各種狀態(tài)、各種工況下的大著角發(fā)火正確性進行比較和判斷。仿真結(jié)果表明，隨著著角的增大，壓電晶體表面應(yīng)力值呈現(xiàn)減小趨勢。在相同著角時，壓電晶體單元的應(yīng)力值隨著碰擊初速的增大而增大，與理論分析結(jié)果一致，說明該方法可行。另外，改變引信頭部結(jié)構(gòu)可以改變壓電晶體表面上的應(yīng)力值，其中頭部錐角的變化對壓電晶體應(yīng)力的影響較大；相同材料的引信頭，其頭部錐角在100°時壓電晶體單元應(yīng)力值相對于其他角度錐角的壓電晶體上的應(yīng)力值要大，也就是說，錐形壓電引信頭的頭部錐角最優(yōu)值是100°(既有設(shè)計即取此值，說明已經(jīng)優(yōu)化，同時也說明本文仿真工作基本可信)，使得壓電引信在以大著角碰擊目標時，壓電晶體上的應(yīng)力值達到最大；在相同速度下，引信頭材料2A11，2A12，7075和H69對壓電晶體作用效果接近，但2A12略優(yōu)于2A11，7075和H69；Q235，45鋼對壓電晶體作用效果優(yōu)于2A12，但不如7A04。但Q235，45鋼密度較大，不太適用于輕武器；采用7A04鋁合金代替原有的2A12作為引信頭材料可以使大著角碰擊目標時壓電晶體表面的碰靶應(yīng)力至少增大30%以上，這對引信保證大著角發(fā)火作用正確性、降低瞎火率具有一定的參考作用。