王維，宣兆龍，段志強，程澤，張曉藝，黃鈺明

(1.軍械工程學(xué)院，石家莊050003;2.72478部隊，濟南250310; 3.68065部隊，甘肅武威733200;4.73833部隊，福州350003)

彈藥儲運方艙是針對我軍彈藥技術(shù)特點及作戰(zhàn)運用實際設(shè)計的新型集裝單元［1］，為了集裝、固定藥筒和引信，防止在裝卸和運輸過程中發(fā)生碰撞，同時為了存取方便，采用支撐件放置藥筒和引信。由于藥筒的筒壁較薄，只有0.7～2 mm［2］，受到壓力后很容易變形，影響其作戰(zhàn)效能，同樣引信也不能受到壓力，因此，支撐件不僅要限制藥筒、引信的運動，還要支撐其重量，避免底層藥筒、引信支撐上部藥筒和引信的重量。由于其結(jié)構(gòu)和受力情況較為復(fù)雜，用傳統(tǒng)的力學(xué)分析方法或工程板梁理論來分析支撐件的受力狀態(tài)較為困難［3］，所以需要應(yīng)用有限元分析軟件進行力學(xué)分析、驗證支撐件的受力及變形情況，以判斷該結(jié)構(gòu)設(shè)計是否滿足強度和剛度的要求。

1 支撐件構(gòu)成

文中以某儲運方艙為例，支撐件固定在該彈藥儲運方艙的中間部位。艙體為“工”字型結(jié)構(gòu)，有利于提高方艙艙體的強度［4］，其結(jié)構(gòu)關(guān)系如圖1所示。藥筒和引信在支撐件中采用的是橫裝方式，分為4層，每層放置3個藥筒、3個引信，為了方便取放彈藥，方艙設(shè)計成前后兩側(cè)開門，其中6個藥筒、6個引信從前側(cè)取放，其余的從后側(cè)取放。為了固定藥筒和引信，在支撐件中設(shè)置藥筒和引信限位孔。由于支撐件的體積較大，為了減輕質(zhì)量，在設(shè)計支撐件材料時，選擇密度較小的硬質(zhì)聚氨酯泡沫，而且硬質(zhì)聚氨酯泡沫還具有保溫性能好、尺寸穩(wěn)定性好、易注塑成形和一定的防火能力等特點［5］。

圖1 彈藥儲運方艙示意Fig.1 Sketch figure for the ammunition storage and transportation shelter

2 有限元模型的建立

2.1 建立幾何模型

有限元模型的建模精度對于分析結(jié)果的影響很大，模型的好壞甚至決定了分析結(jié)果的應(yīng)用價值［6］。可以采用 ANSYS Workbench中的 Design-Modeler工具直接建模，也可以采用從外部導(dǎo)入幾何體的方式［7］，本文選擇后一種建模方式。首先在軟件SolidWorks設(shè)計三維模型，然后導(dǎo)入到ANSYS Workbench中，由于支撐件采用整體注塑成形，結(jié)構(gòu)較為規(guī)則，不需要進行簡化處理。

2.2 有限元網(wǎng)格劃分

有限元網(wǎng)格數(shù)目過少，容易產(chǎn)生畸變，并影響計算精度。細劃網(wǎng)格可以使結(jié)果更精確，但是會增加CPU計算時間和需要更大的存儲空間［8］。根據(jù)實體模型的大小及考慮網(wǎng)格規(guī)模，取單元大小為20 mm進行劃分。

六面體網(wǎng)格可以減少單元數(shù)量，加快求解收斂，可以提高分析精度，減少數(shù)值錯誤［9］，對質(zhì)量好的幾何模型應(yīng)首選六面體網(wǎng)格劃分。本文采用六面體網(wǎng)絡(luò)劃分法劃分網(wǎng)格后的有限元模型如圖2所示。劃分網(wǎng)格后生成220 765個節(jié)點，57 663個單元。

圖2 支撐件模型網(wǎng)格劃分Fig.2 Meshing of the supportmodel

2.3 材料參數(shù)

硬質(zhì)聚氨酯的主要參數(shù):密度為40 kg/m3;彈性模量為10 MPa;泊松比為0.3;壓縮強度極限為0.4 MPa;拉伸強度極限為0.6 MPa［10］。

3 靜剛度、強度分析

3.1 載荷及約束處理

采用橫裝交叉方式放置藥筒和引信，作用在支撐件上的載荷主要由藥筒、引信和支撐件的重力產(chǎn)生。藥筒作用在支撐件上的載荷并不是均布載荷，重心靠近藥筒底部，為了便于計算，簡化為均布載荷作用于藥筒限位孔內(nèi)。由于藥筒是圓柱形，加載時選用軸承載荷(Bearing Load)，軟件中采用的計算方法是根據(jù)徑向分量的投影面積來分布壓力載荷的。根據(jù)強度設(shè)計準(zhǔn)則，需要的安全系數(shù)為 1.7～2.0［11］，藥筒重量G1=m1g=156.96 N。為安全起見，取2.0作為安全系數(shù)，則加載在藥筒限位孔的載荷F1=2.0×G1=313.92 N。引信重量G2=m2g= 9.81 N，同理加載在引信限位孔的載荷F2=2.0×G2= 19.62 N。支撐件的重量G3=m3g=128.51 N。支撐件加載重力時，選擇施加豎直向下的重力加速度。

為保持彈藥儲運方艙內(nèi)支撐件的位置固定，需要將兩側(cè)面及底面與蒙皮相連接。為了簡化約束條件，假設(shè)與蒙皮相連接的部分固定，則約束條件為加載兩側(cè)面及底面固定約束(Fixed Support)。

3.2 有限元分析結(jié)果

在設(shè)計樹中加入所要求的結(jié)果，要求應(yīng)力(Stress)云圖和形變(Deformation)云圖，應(yīng)用ANSYSWorkbench進行仿真計算。

3.2.1 應(yīng)力分析

支撐件的應(yīng)力分析結(jié)果如圖3所示，最大應(yīng)力發(fā)生在底側(cè)藥筒支撐部分，最大應(yīng)力值為2633.7 Pa，遠遠小于材料強度極限，滿足材料的強度極限。

3.2.2 變形分析

變形分析結(jié)果如圖4所示，最大變形發(fā)生在最上側(cè)中間部分，最大變形量為5.6847×10-5m。在設(shè)計時為了便于藥筒、引信放置，圓孔直徑比藥筒、引信直徑大2 mm，所以最大變形遠遠小于2 mm，滿足強度要求。

圖3 支撐件應(yīng)力Fig.3 Stress nephogram of the support

圖4 支撐件形變Fig.4 Deformation nephogram of the support

4 模態(tài)分析

模態(tài)是機械結(jié)構(gòu)的固有振動特性，每一個模態(tài)具有特定的固有頻率、阻尼比和模態(tài)振型［12］。模態(tài)分析主要用于確定結(jié)構(gòu)或機器部件的振動特性(即固有頻率和振型)，固有頻率和振型是承受動載荷機構(gòu)設(shè)計中的重要參數(shù)［13］。通過仿真可以得到支撐件固有的振動特性即固有頻率和振型［14］，在此基礎(chǔ)上，可以進一步分析支撐件在各種動態(tài)激勵下的響應(yīng)。在結(jié)構(gòu)的動態(tài)分析中，各階模態(tài)所具有的權(quán)因子大小與該模態(tài)頻率的倒數(shù)成反比，即頻率越低，權(quán)重越大，就是說低階模態(tài)特性基本決定了產(chǎn)品的動態(tài)性能［15］。

該支撐件涉及到的動力學(xué)主要是在運輸及裝卸過程中，所以需要分析在放置藥筒、引信情況下支撐件的模態(tài)，需要計算時加載預(yù)應(yīng)力及約束條件。文中計算了前6階固有模態(tài)，前6階模態(tài)振型如圖5所示，在預(yù)應(yīng)力作用下的固有頻率和振型見表1。

表1 支撐件前6階模態(tài)頻率振型表Table 1 The preceding 6 ordered inherent frequencies and vibration modes

圖5 模態(tài)振型Fig.5 Ordermodal vibrationmodes

5 結(jié)論

文中應(yīng)用有限元軟件ANSYSWorkbench分析建立了支撐件的有限元模型，分別進行了靜態(tài)剛度、強度分析和動態(tài)的模態(tài)分析。結(jié)果表明，在最大載荷工況下，支撐件的最大應(yīng)力值遠小于材料強度極限，最大變形不超過2 mm，滿足結(jié)構(gòu)的剛強度要求。振動形式有彎曲振動和扭曲振動，振幅較大區(qū)域集中于頂側(cè)中部及第1、第2、第3層中部區(qū)域;振動的前6階頻率集中在117.5～286.7 Hz之間，為下一步進行瞬態(tài)動力學(xué)分析、隨機振動分析奠定了基礎(chǔ)。

［1］ 宣兆龍，趙瑾，劉亞超.鋼絲繩隔振器及其在彈藥方艙中的應(yīng)用［J］.裝備環(huán)境工程，2012，9(4):79—81.XUAN Zhao-long，ZHAO Jin，LIU Ya-chao.Application of Wire-rope Vibration Isolator in Ammunition Shelter［J］.Equipment Environmental Engineering，2012，9(4):79—81.

［2］ 趙曉利，王軍波.彈藥學(xué)［M］.北京:解放軍出版社，1998.ZHAO Xiao-li，WANG Jun-bo.Ammunition［M］.Beijing: Publishing House of PLA，1998.

［3］ 劉亞超.彈藥方艙結(jié)構(gòu)設(shè)計與力學(xué)性能分析［D］.石家莊:軍械工程學(xué)院，2012.LIU Ya-chao.The Structural Design and Mechanical Properties Analysis of Ammunition Shelter［D］.Shijiazhuang: Ordnance Engineering College，2012.

［4］ 劉亞超，宣兆龍，樂惠寧，等.彈藥集裝單元儲存靜力學(xué)有限元分析［J］.軍械工程學(xué)院學(xué)報，2011，23(6): 39—42.LIU Ya-chao，XUAN Zhao-long，YUE Hui-ning，et al.Finite Element Analysis of Statics for Ammunition Packaging Unit［J］.Journal of Ordnance Engineering College，2011，23(6):39—42.

［5］ 王康，黃紅軍，萬國順，等.聚氨酯泡沫材料的粘接性能研究［J］.裝備環(huán)境工程，2012，9(5):122—125.WANG Kang，HUANG Hong-jun，WAN Guo-shun，et al.Study of Bonding Performance of Polyurethane Foam［J］.Equipment Environmental Engineering，2012，9(5): 122—125.

［6］ 劉加凱，齊杏林，崔靜，等.基于ANSYS仿真的引信振動強化試驗夾具設(shè)計［J］.裝備環(huán)境工程，2010，7 (2):95—98.LIU Jia-kai，QI Xing-lin，CUI Jing，et al.Design of Fuse Fixture Used in Vibration Enhancement Test Based on ANSYSSimulation［J］.Equipment Environmental Engineering，2010，7(2):95—98.

［7］ 姚道壯，楊建國，呂志軍.基于ANSYSWorkbench的貨架立柱截面優(yōu)化設(shè)計［J］.東華大學(xué)學(xué)報，2001，7(4): 438—441.YAO Dao-zhuang，YANG Jian-guo，LYU Zhi-jun.Optimization of Steel Storage Rack Column Cross-sections Based on ANSYSWorkbench［J］.Journal of Donghua University，2001，7(4):438—441.

［8］ 浦廣益.ANSYSWorkbench 12基礎(chǔ)教程與實例詳解［M］.北京:中國水利水電出版社，2010.PU Guang-yi.The Basis Course and Example Explanation of ANSYSWorkbench 12［M］.Beijing:China Waterpower Press，2010.

［9］ 許京荊.ANSYS 13.0 Workbench數(shù)值模擬技術(shù)［M］.北京:中國水利水電出版社，2012.XU Jing-jing.The Numerical Simulation Technique of ANSYS 13.0 Workbench［M］.Beijing:ChinaWaterpower Press，2012.

［10］TAN Suqin，ABRAHAM Tim，F(xiàn)ERENCE Don，et al.Rigid Polyurethane Foams from a Soybean Oil-based Polyol［J］.Polymer，2011，52(13):2840—2846.

［11］任文敏，陳艷秋，范欽珊.材料力學(xué)［M］.清華大學(xué)出版社，2004.RENWen-min，CHEN Yan-qiu，F(xiàn)AN Qin-shan.Material Mechanics［M］.Tsinghua University Press，2004.

［12］張國成，樊啟要.層合大板方艙車艙體的模態(tài)分析方法對比研究［J］.機械制造與自動化，2007:80—82.ZHANG Guo-cheng，F(xiàn)AN Qi-yao.Comparison Research on Modal Analysis Methods of Laminated Plate Van［J］.Machine Building＆Automation，2007:80—82.

［13］王良模，吳長風(fēng)，王辰至.特種車輛結(jié)構(gòu)有限元模態(tài)分析［J］.南京理工大學(xué)學(xué)報，2008(11):94—95.WANG Liang-mo，WU Chang-feng，WANG Chen-zhi.Finite Element Analysis of Special Vehicle Shelter Structure［J］.Journal of Nanjing University of Science and Technology，2008(11):94—95.

［14］陳勇武.基于ANSYSWorkbench的噴涂機器人的模態(tài)分析［J］.現(xiàn)代機械，2012(2):44—46.CHEN Yong-wu.Model Analysis of a 6-DOF Painting Robot Based on ANSYSWorkbench［J］.Modern Machinery，2012(2):44—46.

［15］池振坤，楊俊智，周強，等.基于ANSYSWorkbench的大板式方艙模態(tài)分析［J］.汽車工程學(xué)報，2011，1(3): 226—229.CHIZhen-kun，YANG Jun-zhi，ZHOU Qiang，et al.Analysis of Modal Characteristics of the Square Cabin Based on ANSYSWorkbench［J］.Chinese Journal of Automotive Engineering，2011，1(3):226—229.