一種自然破片戰(zhàn)斗部空爆沖擊波和高速破片先后作用的臨界爆距計(jì)算方法
【專(zhuān)利摘要】本發(fā)明涉及一種自然破片戰(zhàn)斗部空中爆炸下沖擊波與高速破片先后作用的臨界爆距計(jì)算方法，它有七大步驟：步驟1，計(jì)算破片的初始速度；步驟2，計(jì)算等效裸裝藥大??；步驟3，計(jì)算戰(zhàn)斗部的極限膨脹半徑；步驟4，推導(dǎo)破片在空氣中飛行距離與時(shí)間的關(guān)系；步驟5，采用步驟2的計(jì)算結(jié)果得到空氣中沖擊波波陣面的傳播距離與時(shí)間的關(guān)系；步驟6，在時(shí)間相同，沖擊波波陣面?zhèn)鞑ゾ嚯x與破片飛行距離相等的條件下，計(jì)算得到?jīng)_擊波和高速破片相遇時(shí)運(yùn)動(dòng)的距離；步驟7，將步驟3和步驟6的結(jié)果求和，得到自然破片戰(zhàn)斗部空中爆炸下沖擊波和高速破片先后作用的臨界爆距值。本發(fā)明能較準(zhǔn)確地計(jì)算出自然破片戰(zhàn)斗部空爆沖擊波和高速破片先后作用的臨界爆距值，在毀傷及防護(hù)等領(lǐng)域均具有較高的實(shí)用價(jià)值和廣闊的應(yīng)用前景。
【專(zhuān)利說(shuō)明】
一種自然破片戰(zhàn)斗部空爆沖擊波和高速破片先后作用的臨界 爆距計(jì)算方法
技術(shù)領(lǐng)域
[0001] 本發(fā)明涉及一種自然破片戰(zhàn)斗部空中爆炸下毀傷元素形成特性及傳播(運(yùn)動(dòng)）規(guī) 律的計(jì)算方法，尤其涉及自然破片戰(zhàn)斗部空中爆炸下沖擊波和高速破片先后作用臨界爆距 的動(dòng)力學(xué)計(jì)算方法，可應(yīng)用于毀傷和防護(hù)等領(lǐng)域，以合理評(píng)估戰(zhàn)斗部毀傷效果或?yàn)榉雷o(hù)結(jié) 構(gòu)的設(shè)計(jì)提供參考依據(jù)等。
【背景技術(shù)】
[0002] 破片戰(zhàn)斗部空中爆炸下會(huì)同時(shí)產(chǎn)生沖擊波和高速破片兩種毀傷載荷。另外，由于 炸藥本身化學(xué)反應(yīng)，還會(huì)形成一定的爆轟產(chǎn)物。相比較而言，接觸爆炸對(duì)結(jié)構(gòu)造成的毀傷要 更為嚴(yán)重，但是在實(shí)戰(zhàn)中非接觸爆炸出現(xiàn)的概率要遠(yuǎn)大于接觸爆炸。因此，對(duì)于空中非接觸 爆炸的情形，沖擊波和高速破片則是兩種主要的毀傷元素。
[0003] 破片戰(zhàn)斗部空中爆炸下初始沖擊波的傳播速度和高速破片初始飛行速度相差很 大，又由于兩者在空氣中的速度衰減效應(yīng)差異明顯，因而導(dǎo)致沖擊波和破片在傳播(運(yùn)動(dòng)） 過(guò)程中存在相遇的問(wèn)題，此相遇位置即為沖擊波和高速破片先后作用的臨界爆距值。爆炸 初期，形成的初始沖擊波傳播速度很高，要遠(yuǎn)高于形成的破片的初始速度;但沖擊波在空氣 中的速度衰減非常快，而破片則慢得多。因而，對(duì)于戰(zhàn)斗部空爆而言，在臨界爆距之前，沖擊 波運(yùn)動(dòng)在破片之前;而超過(guò)臨界爆距值后，破片始終運(yùn)動(dòng)在沖擊波之前。
[0004] 當(dāng)目標(biāo)離破片戰(zhàn)斗部較近時(shí)，其爆炸產(chǎn)生的沖擊波會(huì)先作用在目標(biāo)上，使目標(biāo)產(chǎn) 生變形或破損，隨后破片對(duì)目標(biāo)形成穿甲破壞；當(dāng)目標(biāo)離破片戰(zhàn)斗部較遠(yuǎn)時(shí)，爆炸產(chǎn)生的破 片會(huì)先作用于目標(biāo)，形成穿甲破壞，隨后沖擊波對(duì)目標(biāo)形成進(jìn)一步的毀傷;而當(dāng)目標(biāo)處于臨 界爆距位置處時(shí)，沖擊波和破片會(huì)同時(shí)作用于目標(biāo)，兩者的能量會(huì)形成聯(lián)合增強(qiáng)效應(yīng)，可能 會(huì)使目標(biāo)產(chǎn)生更為嚴(yán)重的破壞。對(duì)于戰(zhàn)斗部設(shè)計(jì)者來(lái)說(shuō)，如果能準(zhǔn)確地掌握該戰(zhàn)斗部的臨 界爆距值，則對(duì)進(jìn)一步提高戰(zhàn)斗部毀傷威力的設(shè)計(jì)具有指導(dǎo)意義;而在防護(hù)領(lǐng)域，如果明確 了防御目標(biāo)戰(zhàn)斗部的臨界爆距值，則對(duì)于盡可能減小戰(zhàn)斗部爆炸毀傷效果的防護(hù)結(jié)構(gòu)的設(shè) 計(jì)具有參考價(jià)值。因此，如何能夠較為準(zhǔn)確且方便快捷地得到戰(zhàn)斗部爆炸下沖擊波與破片 先后作用的臨界爆距值，成為毀傷及防護(hù)等領(lǐng)域深入研究的前提和關(guān)鍵。
[0005] 然而，由于問(wèn)題的復(fù)雜性，破片戰(zhàn)斗部空中爆炸下的沖擊波和高速破片先后作用 的臨界爆距計(jì)算問(wèn)題是毀傷和防護(hù)等領(lǐng)域計(jì)算的難點(diǎn)，一直以來(lái)沒(méi)有一種明確的、可用于 工程實(shí)際的分析計(jì)算方法。目前對(duì)于破片戰(zhàn)斗部空爆下沖擊波和高速破片先后作用的計(jì)算 主要采用數(shù)值仿真的方法。但由于涉及殼體材料、炸藥和空氣等多種介質(zhì)的流固耦合問(wèn)題， 使得有限元數(shù)值仿真計(jì)算需耗費(fèi)大量的計(jì)算資源和時(shí)間，且由于有限元仿真計(jì)算過(guò)程受網(wǎng) 格大小的影響較大，其精確度和可靠性需要驗(yàn)證;而無(wú)網(wǎng)格數(shù)值仿真計(jì)算同樣存在耗費(fèi)大 量計(jì)算資源和時(shí)間等問(wèn)題。因此，急需提出一種相應(yīng)的計(jì)算方法。

【發(fā)明內(nèi)容】

[0006] 本發(fā)明所要解決的技術(shù)問(wèn)題是提出一種自然破片戰(zhàn)斗部空爆沖擊波和高速破片 先后作用的臨界爆距計(jì)算方法，能較好地符合自然破片戰(zhàn)斗部空中爆炸下沖擊波和高速破 片這兩種毀傷元素的形成特點(diǎn)和傳播(運(yùn)動(dòng))規(guī)律。
[0007] 為了解決上述問(wèn)題，本發(fā)明提出了一種自然破片戰(zhàn)斗部空爆沖擊波和高速破片先 后作用的臨界爆距計(jì)算方法，包括以下步驟：
[0008] 步驟1，根據(jù)自然破片戰(zhàn)斗部爆炸形成高速破片的理論，計(jì)算出高速破片的初始速 度；
[0009] 步驟2,根據(jù)戰(zhàn)斗部裝藥的總內(nèi)能等于爆轟產(chǎn)物的內(nèi)能和動(dòng)能及高速破片動(dòng)能等 能量之和的總能量守恒關(guān)系，計(jì)算得到戰(zhàn)斗部的等效裸裝藥大??；
[0010] 步驟3:忽略爆轟產(chǎn)物的繞流效應(yīng)引起的能量損失，結(jié)合能量守恒原理，計(jì)算得到 戰(zhàn)斗部的極限膨脹半徑；
[0011] 步驟4,根據(jù)沖量守恒原理和牛頓第二定律，由破片在空氣中飛行所受到的阻力公 式，推導(dǎo)得到破片在空氣中飛行距離與所需時(shí)間的關(guān)系式；
[0012] 步驟5,根據(jù)爆炸動(dòng)力學(xué)理論，采用步驟2的計(jì)算結(jié)果得到空氣中沖擊波波陣面的 傳播距離與時(shí)間的關(guān)系式；
[0013] 步驟6,在時(shí)間相同，沖擊波波陣面?zhèn)鞑ゾ嚯x與破片飛行距離相等的條件下，計(jì)算 得到自然破片戰(zhàn)斗部空中爆炸下沖擊波和高速破片在空中相遇的位置；
[0014] 步驟7,在步驟3和步驟6的基礎(chǔ)上，得到自然破片戰(zhàn)斗部空中爆炸下沖擊波和高速 破片先后作用的臨界爆距值。
[0015] 進(jìn)一步地，該步驟1中所述的高速破片初始速度vo計(jì)算公式為：
[0017] 式中：co為戰(zhàn)斗部的裝藥量，ms為戰(zhàn)斗部殼體質(zhì)量，E為戰(zhàn)斗部裝藥的特征能量。
[0018] 進(jìn)一步地，若戰(zhàn)斗部的形狀為球形，則步驟1中所述的高速破片的初始速度VQ計(jì)算 公式為：
[0020]進(jìn)一步地，若已知戰(zhàn)斗部裝藥的爆速D，則戰(zhàn)斗部裝藥的特征能量E可根據(jù)下式計(jì) 算得到：
[0022]進(jìn)一步地，該步驟2中所述的戰(zhàn)斗部等效裸裝藥大小為：
[0024]式中:Efk為破片的總動(dòng)能，Et〇tai為戰(zhàn)斗部裝藥總內(nèi)能。
[0025]進(jìn)一步地，破片的總動(dòng)能Efk為：
[0026] = 0.5msv〇
[0027] 進(jìn)一步地，戰(zhàn)斗部裝藥的總內(nèi)能Etotal為：
[0028] Etotai = CO ? Qv
[0029] 式中：Qv為炸藥的爆熱。
[0030] 進(jìn)一步地，該步驟3所述的戰(zhàn)斗部極限膨脹半徑Ro需通過(guò)求解下式得到：
[0032]式中：pf和hf分別為戰(zhàn)斗部殼體厚度和材料密度，為裝藥密度，ro為戰(zhàn)斗部殼體 的初始半徑。
[0033]進(jìn)一步地，該步驟4中所述的破片在空氣中飛行距離Rf與所需時(shí)間tf的關(guān)系為：
[0035]式中:Kf為阻力系數(shù)，mf為破片的設(shè)計(jì)質(zhì)量。
[0036] 進(jìn)一步地，阻力系數(shù)Kf為：
[0037] K{ = -cfpaA
[0038] 式中：cf為與破片形狀相關(guān)的迎風(fēng)阻力系數(shù)，pa為空氣的質(zhì)量密度，j為破片的平 均迎風(fēng)面積。
[0039]進(jìn)一步地，與破片形狀相關(guān)的迎風(fēng)阻力系數(shù)《為： '0.97， 球形
[0040] cf =<1.24, 立方體或規(guī)則矩形 1.50， 不規(guī)則矩形
[0041] 進(jìn)一步地，破片的平均迎風(fēng)面積]由下式計(jì)算：
[0042] A = a- mf
[0043]式中:a為破片的形狀系數(shù)。
[0044] 進(jìn)一步地，破片形狀系數(shù)a的取值為(單位:l(T3m2/kg2/ 3): '3.07, 球形 3.09, 立方體
[0045] a = < 3.55, 規(guī)則矩形 4.42, 不規(guī)則矩形
[0046] 進(jìn)一步地，破片的設(shè)計(jì)質(zhì)量mf通過(guò)以下公式計(jì)算：
[0047] mf = 2BX ? hshdl t/inner(l + /?sheli /^inner)
[0048]式中：Bx為炸藥的爆炸常數(shù)(單位kgV2/mV6)，h sheii和dinner分別為戰(zhàn)斗部殼體的平 均厚度和平均內(nèi)徑。
[0049]進(jìn)一步地，該步驟5中所述采用步驟2的計(jì)算結(jié)果得到空氣中沖擊波波陣面的傳播 距咼Rs與時(shí)間ts的關(guān)系式為：
[0050] Rs = (〇.22 Xl04 -CO°eq16/3 -tj1116
[0051]式中即為步驟2中戰(zhàn)斗部等效裸裝藥大小的計(jì)算結(jié)果。
[0052]進(jìn)一步地，該步驟6中所述的，在時(shí)間ts和tf相同，沖擊波波陣面?zhèn)鞑ゾ嚯xR s與破片 飛行距離Rf相等的條件下，通過(guò)求解以下方程得到自然破片戰(zhàn)斗部空中爆炸下破片追趕上 沖擊波波陣面的時(shí)間tm:
[0054]進(jìn)一步地，可求得自然破片戰(zhàn)斗部空中爆炸下高速破片追趕上沖擊波波陣面的距 離1，可采用下面兩式中的任意一式計(jì)算：
[0057] 進(jìn)一步地，該步驟7中所述自然破片戰(zhàn)斗部空中爆炸下沖擊波和高速破片先后作 用的臨界爆距值Rct為：
[0058] Rcr = R〇+Rm
[0059] 式中：Ro為戰(zhàn)斗部的極限膨脹半徑，心為沖擊波和高速破片空中相遇時(shí)運(yùn)動(dòng)(傳播） 的距離。
[0060] 與現(xiàn)有技術(shù)相比，本發(fā)明具有以下優(yōu)點(diǎn)：
[0061] 本發(fā)明提出的一種自然破片戰(zhàn)斗部空爆沖擊波和高速破片先后作用的臨界爆距 計(jì)算方法通過(guò)動(dòng)力學(xué)計(jì)算能夠簡(jiǎn)便且較為準(zhǔn)確地得到自然破片戰(zhàn)斗部空中爆炸下沖擊波 和高速破片先后作用的臨界爆距值大小，從而能夠?yàn)橄嚓P(guān)武器戰(zhàn)斗部或防護(hù)結(jié)構(gòu)等設(shè)計(jì)及 優(yōu)化提供方便快捷且可靠的參考依據(jù)。
【附圖說(shuō)明】
[0062] 圖1是本發(fā)明的一種自然破片戰(zhàn)斗部空中爆炸下沖擊波和高速破片先后作用臨界 爆距的動(dòng)力學(xué)計(jì)算方法的流程圖。
[0063] 圖2是本發(fā)明實(shí)施例中沖擊波波陣面?zhèn)鞑ニ俣群透咚倨破w行速度的衰減曲線。 [0064]圖3是本發(fā)明實(shí)施例中沖擊波波陣面?zhèn)鞑ゾ嚯x和高速破片飛行距離的變化曲線。
[0065] 圖4是本發(fā)明實(shí)施例中沖擊波波陣面?zhèn)鞑ニ俣燃捌鋫鞑ゾ嚯x隨時(shí)間的變化曲線。
[0066] 圖5是本發(fā)明實(shí)施例中高速破片飛行速度及其飛行距離隨時(shí)間的變化曲線。
【具體實(shí)施方式】
[0067] 為使本發(fā)明的上述目的、特征和優(yōu)點(diǎn)能夠更加明顯易懂，下面結(jié)合附圖和具體實(shí) 施方式對(duì)本發(fā)明作進(jìn)一步詳細(xì)的說(shuō)明。
[0068] 本發(fā)明目的是要提出可應(yīng)用于工程實(shí)際的一種自然破片戰(zhàn)斗部空中爆炸下沖擊 波與高速破片先后作用的臨界爆距計(jì)算方法。
[0069] 請(qǐng)參照?qǐng)D1，為解決上述問(wèn)題，本發(fā)明基于自然破片戰(zhàn)斗部空中爆炸所產(chǎn)生的沖擊 波和高速破片的形成特點(diǎn)及傳播(運(yùn)動(dòng))規(guī)律，提出了一種自然破片戰(zhàn)斗部空中爆炸下沖擊 波與破片先后作用的臨界爆距動(dòng)力學(xué)計(jì)算方法，包括以下步驟：
[0070] 步驟1，根據(jù)自然破片戰(zhàn)斗部爆炸破片形成理論，計(jì)算爆炸形成的高速破片的初始 速度；
[0071] 步驟2,根據(jù)戰(zhàn)斗部裝藥的總內(nèi)能等于爆轟產(chǎn)物的內(nèi)能和動(dòng)能及高速破片動(dòng)能等 能量之和的總能量守恒關(guān)系，計(jì)算得到戰(zhàn)斗部的等效裸裝藥大??；
[0072] 步驟3:忽略爆轟產(chǎn)物的繞流效應(yīng)引起的能量損失，結(jié)合能量守恒原理，計(jì)算得到 戰(zhàn)斗部的極限膨脹半徑；
[0073] 步驟4,根據(jù)沖量守恒原理和牛頓第二定律，由破片在空氣中飛行所受到的阻力公 式，推導(dǎo)得到破片在空氣中飛行距離與所需時(shí)間的關(guān)系式；
[0074] 步驟5,根據(jù)爆炸動(dòng)力學(xué)理論，采用步驟2的計(jì)算結(jié)果得到空氣中沖擊波波陣面的 傳播距離與時(shí)間的關(guān)系式；
[0075]步驟6,在時(shí)間相同，沖擊波波陣面?zhèn)鞑ゾ嚯x與破片飛行距離相等的條件下，計(jì)算 得到自然破片戰(zhàn)斗部空中爆炸下沖擊波和高速破片在空氣中的相遇距離；
[0076] 步驟7,在步驟3和步驟6的基礎(chǔ)上，得到自然破片戰(zhàn)斗部空中爆炸下沖擊波和高速 破片先后作用的臨界爆距值。
[0077] 下面結(jié)合實(shí)例對(duì)本發(fā)明的一種自然破片戰(zhàn)斗部空中爆炸下沖擊波和高速破片先 后作用的臨界爆距計(jì)算方法進(jìn)行詳細(xì)說(shuō)明。實(shí)施例為國(guó)外某一自然破片戰(zhàn)斗部空中爆炸的 情形。該破片戰(zhàn)斗部的炸藥當(dāng)量為34.93kg，爆炸所產(chǎn)生的破片設(shè)計(jì)質(zhì)量及其飛行參數(shù)為：
[0078] 表1爆炸所產(chǎn)生的破片設(shè)計(jì)質(zhì)量及其飛行參數(shù)
[0080] 下面是具體步驟：
[0081] 步驟1，根據(jù)自然破片戰(zhàn)斗部爆炸破片形成理論，計(jì)算出爆炸形成的高速破片的初 始速度；
[0082] 根據(jù)前述步驟1中的計(jì)算公式可得出破片的初始速度為887m/s。
[0083] 步驟2,根據(jù)戰(zhàn)斗部裝藥的總內(nèi)能等于爆轟產(chǎn)物的內(nèi)能和動(dòng)能及高速破片動(dòng)能等 能量之和的總能量守恒關(guān)系，計(jì)算得到戰(zhàn)斗部的等效裸裝藥大??；
[0084] 將戰(zhàn)斗部裝藥的總內(nèi)能減去戰(zhàn)斗部產(chǎn)生的破片的總動(dòng)能，然后再將剩下的總能量 轉(zhuǎn)化為裸裝藥質(zhì)量，根據(jù)前述步驟2中相應(yīng)的轉(zhuǎn)化公式得到戰(zhàn)斗部的等效裸裝藥為 23.75kg〇
[0085] 步驟3:忽略爆轟產(chǎn)物的繞流效應(yīng)引起的能量損失，結(jié)合能量守恒原理，計(jì)算得到 戰(zhàn)斗部的極限膨脹半徑；
[0086] 根據(jù)前述步驟3中的方程，代入各參數(shù)求解得到戰(zhàn)斗部的極限膨脹半徑值約為 0.28m〇
[0087]步驟4,根據(jù)沖量守恒原理和牛頓第二定律，由破片在空氣中飛行所受到的阻力公 式，推導(dǎo)得到破片在空氣中飛行距離與所需時(shí)間的關(guān)系式；
[0088]根據(jù)前述步驟3,計(jì)算得到破片在空氣中飛行時(shí)的阻力系數(shù)Kf = -0.3 XI. 25 X 19 ? 72 X 10-4 = -7 ? 4 X 10-4(kg/m) 〇
[0089] 進(jìn)一步地，根據(jù)前述步驟4中的關(guān)系式，得到破片在空氣中飛行距離與飛行時(shí)間的 關(guān)系式為=887.-393.4.4+232.6.? f3 (m)。
[0090] 步驟5,根據(jù)爆炸動(dòng)力學(xué)理論，采用步驟2的計(jì)算結(jié)果得到空氣中沖擊波波陣面的 傳播距離與時(shí)間的關(guān)系式；
[0091] 根據(jù)前述步驟4的關(guān)系式，得到空氣中沖擊波波陣面?zhèn)鞑ゾ嚯x與所需時(shí)間的關(guān)系 式為：i?s= 125.1 .，s° 568 (m)〇
[0092]步驟6,在時(shí)間相同，沖擊波波陣面?zhèn)鞑ゾ嚯x與破片飛行距離相等的條件下，計(jì)算 得到破片戰(zhàn)斗部空中爆炸下沖擊波與破片先后作用的臨界爆距值大小。
[0093]根據(jù)前述步驟6,通過(guò)聯(lián)立步驟4和步驟5中飛行(傳播)距離與時(shí)間的關(guān)系式，在時(shí) 間相同距離相等的條件下，求得破片追趕上沖擊波波陣面的時(shí)間為10.8ms。
[0094] 進(jìn)一步地，根據(jù)步驟4或步驟5中的關(guān)系式，求得該戰(zhàn)斗部空中爆炸下沖擊波與高 速破片相遇時(shí)運(yùn)動(dòng)(傳播)的距離大小為9.56m。
[0095] 步驟7,在步驟3和步驟6的基礎(chǔ)上，得到自然破片戰(zhàn)斗部空中爆炸下沖擊波和高速 破片先后作用的臨界爆距值。
[0096] 進(jìn)一步地，根據(jù)步驟7中關(guān)系式，得到該自然破片戰(zhàn)斗部空中爆炸下沖擊波和高速 破片先后作用的臨界爆距值為9.84m。文獻(xiàn)中實(shí)驗(yàn)值為10.2m，計(jì)算結(jié)果與文獻(xiàn)實(shí)測(cè)值間的 相對(duì)誤差約為-3.5%?？紤]到實(shí)驗(yàn)測(cè)試的精度誤差，計(jì)算結(jié)果與文獻(xiàn)中的實(shí)驗(yàn)值吻合較好， 能夠滿(mǎn)足工程應(yīng)用的需要。
【主權(quán)項(xiàng)】
1. 一種自然破片戰(zhàn)斗部空爆沖擊波和高速破片先后作用的臨界爆距計(jì)算方法，其特征 在于該方法由以下步驟實(shí)現(xiàn):步驟1，計(jì)算破片的初始速度;步驟2，計(jì)算等效裸裝藥大小;步 驟3,計(jì)算戰(zhàn)斗部的極限膨脹半徑;步驟4,推導(dǎo)破片在空氣中飛行距離與時(shí)間的關(guān)系；步驟 5,采用步驟2的計(jì)算結(jié)果得到空氣中沖擊波波陣面的傳播距離與時(shí)間的關(guān)系；步驟6,在時(shí) 間相同，沖擊波波陣面?zhèn)鞑ゾ嚯x與破片飛行距離相等的條件下，計(jì)算得到?jīng)_擊波和高速破 片相遇時(shí)運(yùn)動(dòng)的距離;步驟7,將步驟3和步驟6的結(jié)果求和，得到自然破片戰(zhàn)斗部空中爆炸 下沖擊波和高速破片先后作用的臨界爆距值。2. 如權(quán)利要求1所述的計(jì)算方法，其特征是步驟2、3的順序可以任意更換。3. 如權(quán)利要求1所述的計(jì)算方法，其特征是步驟4、5的順序可以任意更換。4. 如權(quán)利要求1所述的計(jì)算方法，其特征是步驟6中所述的先后作用臨界爆距值大小Rct 等于辦與1之和，Ro表示戰(zhàn)斗部的極限膨脹半徑，Rm表示沖擊波和高速破片相遇時(shí)所運(yùn)動(dòng)的 距離。
【文檔編號(hào)】G06F17/50GK105912744SQ201610145789
【公開(kāi)日】2016年8月31日
【申請(qǐng)日】2016年3月15日
【發(fā)明人】陳長(zhǎng)海, 朱錫, 侯海量, 胡年明
【申請(qǐng)人】中國(guó)人民解放軍海軍工程大學(xué)