白亞瓊,師 平
(西安航空職業(yè)技術學院 西安閻良 710089)
由于加工及材料特性等原因,引信體表面及結構內部有可能會出現初始裂紋。本文根據有限元法和線彈性斷裂力學的基本理論,依據 SliodWorks和ANSYS軟件,預測在引信體工作載荷下危險區(qū)域的裂紋。
爆竹的火藥捻子是最早的引信。引信的系統(tǒng)結構包括:發(fā)火控制系統(tǒng)、安全系統(tǒng)、傳爆序列和能源裝置 4個基本部分。根據作用原理,引信可以分為觸發(fā)引信、非觸發(fā)引信和時間引信 3種類型,大都由擊針、火帽、雷管、傳爆藥和保險機構等組成。本文研究的某引信體為觸發(fā)式,彈藥的部位在彈身引信,引信觸發(fā)后碰著物體即起爆。
本文以某彈引信體為例,由于炸藥在引信體結構內爆炸情況十分復雜,必須在實驗條件下對引信體的受力進行實測并找出峰值壓力。在相關條件缺乏的情況下,假設從炸藥爆炸瞬間到擊針脫離引信體這一過程中,爆炸未發(fā)生變化,且將子彈從發(fā)射到離開引信體看作勻加速運動,引信體材料為 45號冷拉圓鋼,HRC45-50,泊松比為 0.309。由于引信體是軸對稱結構,其中每個截面所受到的爆破應力應相等,故進行有限元計算時,取二分之一引信體為計算區(qū)域,該結構的尺寸為φ18(彈尾)×φ9.4(彈頭)×3.5(厚度)mm,如圖 1所示。在有限元中采用 Solid45實體單元,通過自由網格劃分模型。劃分后的引信體如圖2所示。
圖1 引信體三維模型Fig.1 3,D model of a fuze
圖2 引信體網格劃分的模型Fig.2 Meshing mode of a fuze
已知擊針變形時的最大壓力為160,kg,擊針的出口速度大約為 450,m/s。擊針在引信體中運動時所受的阻力 f等同于其變形時受到的壓力。擊針材料為45號鋼,HRC40-45,質量為 0.11,g。
擊針在引信體中的加速度:a=[(vt)2-(v0)2]/2,s=(vt)2/2,s=18,925,000,其中初始速度v0為零;
引信體膛內氣體壓力為:
F=ma=10,500,000,N
引信體膛內最大氣體壓強為:
p=F/S=24,MPa
不考慮溫度和氣候的影響,引信體的載荷是對稱的,可以取其結構的二分之一進行分析。由于要對其進行有限元靜態(tài)分析必須限制其自由度,所以引信體所受載荷的工作壓力為 24,MPa,施加于引信體結構的內表面上。
通過對圖3和圖4所示內容的分析,可以準確地確定最大應力的值和出現的位置,這些最大值見表1。由圖3和圖4可以看出,在引信體頸部與放置擊針處之間(圖1標注的薄弱區(qū)域)應力變形最大,但其并不影響擊針的正常發(fā)射;而引信體前端的變形非常小,此處為引信體的重要部分。如果該部分應力超過許用應力值就會產生斷裂,直接影響擊針的正常發(fā)射。
圖3 引信體的最大拉應力分布Fig.3 Distribution of a fuze’s Max. tension stress
圖4中,引信體頸部與放置擊針處之間(圖1標注的薄弱區(qū)域)所受到的等效應力最大,此處為引信體的應力集中地帶。查閱《合金鋼手冊》,45號鋼的強度極限為 608,MPa。很明顯,引信體等效應力大于極限應力,因此引信體在爆破應力作用下十分容易斷裂。要避免引信體的斷裂,就必須避免此處的應力集中情況??梢酝ㄟ^增大引信體壁厚或調整內部結構的圓角大小,減少引信體受到的爆破應力,但這些措施只是初級手段,尚需進一步研究。
圖4 依據第三強度理論的等效應力分布Fig.4 Distribution of equivalence stress in accordance with the 3 rd strength theory
表1 引信體的三維模型應力計算結果Tab.1 Stress calculation results of a fuze’s 3,D model
運用 ANSYS對引信體進行爆破應力分析可知,引信體頸部與放置擊針處之間(圖1標注的薄弱區(qū)域)所受的等效應力最大,斷裂處從這里開始。因此,要避免引信體發(fā)生斷裂,必須避免應力在此處集中,具體措施可以是增大引信體壁厚或調整內部結構的圓角大小。
[1] 尹柏生. 有限元分析系統(tǒng)的發(fā)展現狀與展望[J]. 計算機世界,2000(3):1-2.
[2] 董方晴. 引信設計原理[M]. 北京:國防工業(yè)出版社,1973:10.
[3] 曹萬有,張文柱,王道宏. 高膛壓火炮技術[M]. 北京:國防工業(yè)出版社,1989.
[4] 王世軍. 有限元分析軟件 ANSYS在結構分析中的應用[M]. 西安:西安理工大學出版,1999.
[5] R. D. 庫克. 有限元分析的概念和應用[M]. 北京:清華大學出版社,2004.