專利名稱:超聲速流道的實現(xiàn)方法及裝置的制作方法
技術(shù)領(lǐng)域:
本發(fā)明空氣動力設(shè)計領(lǐng)域,具體而言,涉及一種超聲速流道的實現(xiàn)方法及裝置。
背景技術(shù):
超聲速流道,如噴管、進(jìn)氣道、隔離段、轉(zhuǎn)彎段等,常用于風(fēng)洞和高速飛行器,一般是二維或軸對稱構(gòu)型,因此其截面一般是圓形、環(huán)形或矩形。隨著現(xiàn)代空氣動力學(xué)的發(fā)展, (此處刪去回車,并為一段)截面為橢圓形、星形、多邊形等形狀的異型流道有著越來越廣泛的用途,如增強(qiáng)兩股氣流間的混合、提高超聲速流道的結(jié)構(gòu)強(qiáng)度等。由于失去了幾何上的對稱性,傳統(tǒng)的流道設(shè)計方法無法應(yīng)用于異型出口截面流道設(shè)計中,因此,目前異形截面流道還沒有較為成熟的設(shè)計技術(shù),與之相近的是二維或軸對稱超聲速流道設(shè)計方法。現(xiàn)有的超聲速流道設(shè)計方法一般步驟如下(I)根據(jù)收縮比和長徑比的要求,設(shè)計亞聲速段曲線;(2)根據(jù)亞聲速段曲線確定喉部的超聲速初值線;(3)根據(jù)超聲速初值線和一定的理論公式確定超聲速段的馬赫數(shù)分布;(4)根據(jù)超聲速段的馬赫數(shù)分布確定超聲速段壁面。在上述超聲速流道的設(shè)計方法中,只能設(shè)計二維或軸對稱流道,無法設(shè)計出異型出口截面流道。此外,在現(xiàn)有的超聲速流道設(shè)計方法,只通過簡單的曲面插值而不考慮氣動特性,所設(shè)計的異型出口截面流道流場品質(zhì)難以保證。針對相關(guān)技術(shù)中的無法設(shè)計出口截面為異型的超聲速流道問題,目前尚未提出有效的解決方案。
發(fā)明內(nèi)容
針對相關(guān)技術(shù)中的無法設(shè)計異型出口截面流道問題,本發(fā)明提供了一種超聲速流道的實現(xiàn)方法及裝置,以至少解決上述問題。根據(jù)本發(fā)明的一個方面,提供了一種超聲速流道的實現(xiàn)方法,包括根據(jù)預(yù)設(shè)的異型出口截面的形狀獲取該形狀的外接圓;計算外接圓的圓心到異型出口截面的形狀的距離;從超聲速流道的流線族中選擇與該距離相對應(yīng)的流線,其中,選擇出的流線構(gòu)成超聲速流道的壁面的曲面形狀。優(yōu)選的,計算外接圓的圓心到異型出口截面的形狀的距離的步驟包括獲取圓心到外接圓的圓周上的每個點的連線;獲取連線與形狀的邊的第一交點;計算圓心到第一交點的距離,其中,圓心到第一交點的距離為外接圓的圓心到異型出口截面的形狀的距離。優(yōu)選的,從超聲速流道的流線族中選擇與距離相對應(yīng)的流線的步驟包括獲取流線族中的各流線與超聲速流道的出口連接線之間的第二交點,其中,出口連接線為超聲速流道的出口點與該出口點相對于超聲速流道的中心線的鏡像點的連接線;獲取中心線與超聲速流道的出口連接線之間的第三交點;計算第二交點與第三交點之間的距離;若第二交點與第三交點之間的距離為外接圓的圓心到異型出口截面的形狀的距離,則選擇與該第二交點對應(yīng)的流線。優(yōu)選的,通過以下步驟獲取中心線根據(jù)預(yù)設(shè)的異型出口截面的進(jìn)口位置和出口位置以及馬赫數(shù)約束確定中心線的坐標(biāo)和長度。優(yōu)選的,通過以下步驟獲取超聲速流道的出口點以及該出口點的鏡像點根據(jù)馬赫數(shù)分布以及出口流場設(shè)計約束確定出口點;根據(jù)中心線、特征線法以及確定出的出口點生成超聲速流道的曲線型面;將曲線型面相對于中心線進(jìn)行鏡像得到出口點的鏡像點。優(yōu)選的,異型出口截面的形狀包括以下之一多邊形截面、曲線截面、星形截面以及多連通截面。根據(jù)本發(fā)明的另一方面,提供了一種超聲速流道的實現(xiàn)裝置,包括獲取單元,用于根據(jù)預(yù)設(shè)的異型出口截面的形狀獲取形狀的外接圓;計算單元,用于計算外接圓的圓心到異型出口截面的形狀的距離;選擇單元,用于從超聲速流道的流線族中選擇與距離相對應(yīng)的流線,其中,選擇出的流線構(gòu)成超聲速流道的壁面的曲面形狀。優(yōu)選的,計算單元包括第一獲取模塊,用于獲取圓心到外接圓的圓周上的每個點的連線;第二獲取模塊,用于獲取連線與形狀的邊的第一交點;計算模塊,用于計算圓心到第一交點的距離,其中,圓心到第一交點的距離為外接圓的圓心到異型出口截面的形狀的距離。優(yōu)選的,選擇單元包括第三獲取模塊,用于獲取流線族中的各流線與超聲速流道的出口連接線之間的第二交點,其中,出口連接線為超聲速流道的出口點與該出口點相對于超聲速流道的中心線的鏡像點的連接線;第四獲取模塊,用于獲取中心線與超聲速流道的出口連接線之間的第三交點;計算模塊,用于計算第二交點與第三交點之間的距離;選擇模塊,用于選擇第二交點與第三交點之間的距離為外接圓的圓心到異型出口截面的形狀的距離的第二交點對應(yīng)的流線。優(yōu)選的,第三獲取模塊包括第一確定子模塊,用于根據(jù)預(yù)設(shè)的異型出口截面的進(jìn)口位置和出口位置以及馬赫數(shù)約束確定中心線的坐標(biāo)和長度;第二確定子模塊,用于根據(jù)馬赫數(shù)分布以及出口流場設(shè)計約束確定出口點;生成子模塊,用于根據(jù)中心線、特征線法以及確定出的出口點生成超聲速流道的曲線型面;鏡像子模塊,用于將曲線型面相對于中心線進(jìn)行鏡像得到出口點的鏡像點。通過本發(fā)明,先獲取預(yù)設(shè)的異型出口截面的形狀獲取形狀的外接圓,通過確定的外接圓的流線族和預(yù)設(shè)的異型出口截面來獲取出口截面為異型的超聲速流道,解決了相關(guān)技術(shù)中的無法設(shè)計出口截面為異型的超聲速流道問題,達(dá)到了實現(xiàn)超聲速流道的出口截面為異型的技術(shù)效果。
此處所說明的附圖用來提供對本發(fā)明的進(jìn)一步理解,構(gòu)成本申請的一部分,本發(fā)明的示意性實施例及其說明用于解釋本發(fā)明,并不構(gòu)成對本發(fā)明的不當(dāng)限定。在附圖中圖I是根據(jù)本發(fā)明實施例的超聲速流道的實現(xiàn)方法的一種優(yōu)選流程圖;圖2是根據(jù)本發(fā)明實施例的超聲速流道的實現(xiàn)裝置的一種優(yōu)選結(jié)構(gòu)框圖;圖3是根據(jù)本發(fā)明實施例的超聲速流道的實現(xiàn)裝置的另一種優(yōu)選結(jié)構(gòu)框圖;圖4是根據(jù)本發(fā)明實施例的超聲速流道的實現(xiàn)方法中截面形狀為多邊形截面的
5示意圖;圖5是根據(jù)本發(fā)明實施例的超聲速流道的實現(xiàn)方法中截面形狀為曲線截面的示 意圖;圖6是根據(jù)本發(fā)明實施例的超聲速流道的實現(xiàn)方法中截面形狀為星形截面的示 意圖;圖7是根據(jù)本發(fā)明實施例的超聲速流道的實現(xiàn)方法中截面形狀為多連通截面的 示意圖;圖8是根據(jù)本發(fā)明實施例的超聲速流道的實現(xiàn)方法中多邊形截面對應(yīng)的外接圓;圖9是根據(jù)本發(fā)明實施例的超聲速流道的實現(xiàn)方法中曲線截面對應(yīng)的外接圓;圖10是根據(jù)本發(fā)明實施例的超聲速流道的實現(xiàn)方法中星形截面對應(yīng)的外接圓;圖11是根據(jù)本發(fā)明實施例的超聲速流道的實現(xiàn)方法中多連通截面的外接圓;圖12是根據(jù)本發(fā)明實施例的超聲速流道的中心線及其馬赫數(shù)分布;圖13是根據(jù)本發(fā)明實施例的根據(jù)特征線法設(shè)計流道的曲線壁面的示意圖;圖14是根據(jù)本發(fā)明實施例的超聲速流道的實現(xiàn)方法中流線族的示意圖;圖15是根據(jù)本發(fā)明實施例的超聲速流道的實現(xiàn)方法中確定多邊形截面的流線;圖16是根據(jù)本發(fā)明實施例的超聲速流道的實現(xiàn)方法中確定曲線截面的流線;圖17是根據(jù)本發(fā)明實施例的超聲速流道的實現(xiàn)方法中確定星形截面的流線;圖18是根據(jù)本發(fā)明實施例的超聲速流道的實現(xiàn)方法中確定多連通截面的流線; 以及圖19是根據(jù)本發(fā)明實施例的超聲速流道的實現(xiàn)方法中特征線法迭代公式的示意 圖。
具體實施例方式下文中將參考附圖并結(jié)合實施例來詳細(xì)說明本發(fā)明。需要說明的是,在不沖突的 情況下,本申請中的實施例及實施例中的特征可以相互組合。實施例1本發(fā)明提供了一種超聲速流道的實現(xiàn)方法,圖1為本發(fā)明超聲速流道的實現(xiàn)方法 的一種優(yōu)選流程圖,如圖1所示,該方法包括S102,根據(jù)預(yù)設(shè)的異型出口截面的形狀獲取形狀的外接圓;S104,計算外接圓的圓心到異型出口截面的形狀的距離;S106,從超聲速流道的流線族中選擇與距離相對應(yīng)的流線,其中,選擇出的流線構(gòu) 成超聲速流道的壁面的曲面形狀。上述實施例中,先獲取預(yù)設(shè)的異型出口截面的形狀獲取形狀的外接圓,通過確定 的外接圓的流線族和預(yù)設(shè)的異型出口截面來獲取出口截面為異型的超聲速流道,解決了相 關(guān)技術(shù)中的無法設(shè)計出口截面為異型的超聲速流道問題,達(dá)到了實現(xiàn)超聲速流道的出口截 面為異型的技術(shù)效果。為了計算外接圓的圓心到異型出口截面的形狀的距離,本發(fā)明還提供了一種優(yōu)選 的計算方法,具體的,計算外接圓的圓心到異型出口截面的形狀的距離的步驟包括獲取圓 心到外接圓的圓周上的每個點的連線;獲取連線與形狀的邊的第一交點;計算圓心到第一交點的距離,其中,圓心到第一交點的距離為外接圓的圓心到異型出口截面的形狀的距離。 上述技術(shù)方案的實施,可以達(dá)到通過外接圓的半徑來確定外接圓的圓心到異型出口截面的形狀的距離的技術(shù)效果。本發(fā)明還提供了一種優(yōu)選的從超聲速流道的流線族中選擇與距離相對應(yīng)的流線的方法,具體的,該方法包括獲取流線族中的各流線與超聲速流道的出口連接線之間的第二交點,其中,出口連接線為超聲速流道的出口點與該出口點相對于超聲速流道的中心線的鏡像點的連接線;獲取中心線與超聲速流道的出口連接線之間的第三交點;計算第二交點與第三交點之間的距離;若第二交點與第三交點之間的距離為外接圓的圓心到異型出口截面的形狀的距離,則選擇與該第二交點對應(yīng)的流線。上述技術(shù)方案的實施,可以達(dá)到通過簡單的運(yùn)算即可從超聲速流道的流線族中選擇與距離相對應(yīng)的流線的技術(shù)效果。本發(fā)明還提供了一種獲取中心線的方案,優(yōu)選的,該方案包括根據(jù)預(yù)設(shè)的異型出口截面的進(jìn)口位置和出口位置以及馬赫數(shù)約束確定中心線的坐標(biāo)和長度。為了獲取超聲速流道的出口點以及該出口點的鏡像點,本發(fā)明提供了一種優(yōu)選的方法,該方法包括根據(jù)馬赫數(shù)分布以及出口流場設(shè)計約束確定出口點;根據(jù)中心線、特征線法以及確定出的出口點生成超聲速流道的曲線型面;將曲線型面相對于中心線進(jìn)行鏡像得到出口點的鏡像點。上述技術(shù)方案的實施,可以達(dá)到僅僅通過簡單的幾何運(yùn)算即可獲取超聲速流道的出口點以及該出口點的鏡像點的技術(shù)效果。優(yōu)選的,上述異型出口截面的形狀包括以下之一多邊形截面、曲線截面、星形截面以及多連通截面。上述技術(shù)方案的實施,可以達(dá)到該超聲速流道的實現(xiàn)方法對多種異型出口截面的超聲速流道的設(shè)計。實施例2本發(fā)明提供了一種超聲速流道的實現(xiàn)裝置,圖2為本發(fā)明超聲速流道的實現(xiàn)裝置的一種優(yōu)選結(jié)構(gòu)框圖,如圖2所示,該裝置包括獲取單元202,用于根據(jù)預(yù)設(shè)的異型出口截面的形狀獲取形狀的外接圓;計算單元204,用于計算外接圓的圓心到異型出口截面的形狀的距離;選擇單元206,用于從超聲速流道的流線族中選擇與距離相對應(yīng)的流線,其中, 選擇出的流線構(gòu)成超聲速流道的壁面的曲面形狀。上述實施例中,先獲取預(yù)設(shè)的異型出口截面的形狀獲取形狀的外接圓,通過確定的外接圓的流線族和預(yù)設(shè)的異型出口截面來獲取出口截面為異型的超聲速流道,解決了相關(guān)技術(shù)中的無法設(shè)計出口截面為異型的超聲速流道問題,達(dá)到了實現(xiàn)超聲速流道的出口截面為異型的技術(shù)效果。本發(fā)明還對上述計算單元204進(jìn)行了改進(jìn),優(yōu)選的,如圖3所示,計算單元204包括第一獲取模塊2042,用于獲取圓心到外接圓的圓周上的每個點的連線;第二獲取模塊 2044,用于獲取連線與形狀的邊的第一交點;計算模塊2046,用于計算圓心到第一交點的距離,其中,圓心到第一交點的距離為外接圓的圓心到異型出口截面的形狀的距離。上述技術(shù)方案的實施,可以達(dá)到通過外接圓的半徑來確定外接圓的圓心到異型出口截面的形狀的距離的技術(shù)效果。本發(fā)明還對上述選擇單元206進(jìn)行了改進(jìn),優(yōu)選的,如圖3所示,選擇單元206包括第三獲取模塊2062,用于獲取流線族中的各流線與超聲速流道的出口連接線之間的第二交點,其中,出口連接線為超聲速流道的出口點與該出口點相對于超聲速流道的中心線的鏡像點的連接線;第四獲取模塊2064,用于獲取中心線與超聲速流道的出口連接線之間的第三交點;計算模塊2066,用于計算第二交點與第三交點之間的距離;選擇模塊2068,用于選擇第二交點與第三交點之間的距離為外接圓的圓心到異型出口截面的形狀的距離的第二交點對應(yīng)的流線。上述技術(shù)方案的實施,可以達(dá)到通過簡單的運(yùn)算即可從超聲速流道的流線族中選擇與距離相對應(yīng)的流線的技術(shù)效果。本發(fā)明還對上述第三獲取模塊2062進(jìn)行了改進(jìn),優(yōu)選的,如圖3所示,第三獲取模塊2062包括第一確定子模塊20622,用于根據(jù)預(yù)設(shè)的異型出口截面的進(jìn)口位置和出口位置以及馬赫數(shù)約束確定中心線的坐標(biāo)和長度。本發(fā)明還對上述第三獲取模塊2062進(jìn)行了改進(jìn),優(yōu)選的,如圖3所示,第三獲取模塊2062還包括第二確定子模塊20624,用于根據(jù)馬赫數(shù)分布以及出口流場設(shè)計約束確定出口點;生成子模塊20626,用于根據(jù)中心線、特征線法以及確定出的出口點生成超聲速流道的曲線型面;鏡像子模塊20628,用于將曲線型面相對于中心線進(jìn)行鏡像得到出口點的鏡像點。優(yōu)選的,上述異型出口截面的形狀包括以下之一多邊形截面、曲線截面、星形截面以及多連通截面。上述技術(shù)方案的實施,可以達(dá)到該超聲速流道的實現(xiàn)裝置完成對多種異型出口截面的超聲速流道的設(shè)計。實施例3本發(fā)明還提供了另一種優(yōu)選的超聲速流道的實現(xiàn)方法,具體的,該方法包括SI,根據(jù)結(jié)構(gòu)設(shè)計要求確定出口截面形狀;優(yōu)選的,如圖4-圖7所示,出口截面形狀包括但不限于以下形狀之一多邊形截面、曲線截面、星形截面以及多連通截面。S2,確定出口截面的外接圓;圖8-圖11分別對應(yīng)于上述圖4-圖7中出口截面的外接圓,在圖8-圖11中,虛線表示根據(jù)結(jié)構(gòu)設(shè)計要求確定出口截面形狀,實線表示該出口截面的外接圓。S3,采用特征線法設(shè)計所確定的圓形流道,其中,包括a)根據(jù)進(jìn)出口位置和馬赫數(shù)約束,確定中心線長度及其馬赫數(shù)分布;圖12示出所確定的中心線長度及其馬赫數(shù)分布,其中,AB為中心線。b)根據(jù)中心線長度及其馬赫數(shù)分布,和出口流場設(shè)計約束,確定出口點D,如圖13 所示,根據(jù)特征線法設(shè)計流道的曲線型面D-C,根據(jù)對稱性,將D-C關(guān)于A-B鏡像,得到對邊壁面D’ -C’。優(yōu)選的,此處附上特征線法迭代公式已知兩點=CKx1,!T1, M1, 0 ) , P(x2, r2, M2, 0 2)求解第三個點Q(x3,r3,M3,0 3)(四個方程,四個未知數(shù)),如圖19所示,其中,x為 X軸坐標(biāo),在二維情況下,r分別為y軸坐標(biāo),軸對稱情況下,r為柱坐標(biāo)下的徑向坐標(biāo),M為馬赫數(shù),0為流動方向角,下標(biāo)1、2、3分別對應(yīng)各特征點。預(yù)估步驟先求解(x3,r3)u I = SirT1 (I/M),為點 的馬赫角,u 2 = SirT1 (!/M2),為點 P 的馬赫角,
h: = tan[ 0 ^ U J,為點Q發(fā)出的左行特征線的斜率h2 = tan[ 0 2_ ii 2],為點P發(fā)出的右行特征線的斜率根據(jù)差分方程有T^r1 = Ii1 (X3-X1)r3-r2 = h2 (X3-X2)兩式相減可得T1-T2 = {h2-hj XjXih1-X2Ii2由此得到點Q的坐標(biāo)
權(quán)利要求
1.一種超聲速流道的實現(xiàn)方法,其特征在于,包括根據(jù)預(yù)設(shè)的異型出口截面的形狀獲取所述形狀的外接圓;計算所述外接圓的圓心到所述異型出口截面的形狀的距離;從超聲速流道的流線族中選擇與所述距離相對應(yīng)的流線,其中,選擇出的流線構(gòu)成所述超聲速流道的壁面的曲面形狀。
2.根據(jù)權(quán)利要求I所述的方法,其特征在于,所述計算所述外接圓的圓心到所述異型出口截面的形狀的距離的步驟包括獲取所述圓心到所述外接圓的圓周上的每個點的連線;獲取所述連線與所述形狀的邊的第一交點;計算所述圓心到所述第一交點的距離,其中,所述圓心到所述第一交點的距離為所述外接圓的圓心到所述異型出口截面的形狀的距離。
3.根據(jù)權(quán)利要求I所述的方法,其特征在于,所述從超聲速流道的流線族中選擇與所述距離相對應(yīng)的流線的步驟包括獲取所述流線族中的各流線與所述超聲速流道的出口連接線之間的第二交點,其中, 所述出口連接線為所述超聲速流道的出口點與該出口點相對于所述超聲速流道的中心線的鏡像點的連接線;獲取所述中心線與所述超聲速流道的出口連接線之間的第三交點;計算所述第二交點與所述第三交點之間的距離;若所述第二交點與所述第三交點之間的距離為所述外接圓的圓心到所述異型出口截面的形狀的距離,則選擇與該第二交點對應(yīng)的流線。
4.根據(jù)權(quán)利要求3所述的方法,其特征在于,通過以下步驟獲取所述中心線根據(jù)預(yù)設(shè)的異型出口截面的進(jìn)口位置和出口位置以及馬赫數(shù)約束確定所述中心線的坐標(biāo)和長度。
5.根據(jù)權(quán)利要求3所述的方法,其特征在于,通過以下步驟獲取所述超聲速流道的出口點以及該出口點的鏡像點根據(jù)馬赫數(shù)分布以及出口流場設(shè)計約束確定所述出口點;根據(jù)所述中心線、所述特征線法以及確定出的出口點生成所述超聲速流道的曲線型面;將所述曲線型面相對于所述中心線進(jìn)行鏡像得到所述出口點的鏡像點。
6.根據(jù)權(quán)利要求I所述的方法,其特征在于,所述異型出口截面的形狀包括以下之一 多邊形截面、曲線截面、星形截面以及多連通截面。
7.一種超聲速流道的實現(xiàn)裝置,其特征在于,包括獲取單元,用于根據(jù)預(yù)設(shè)的異型出口截面的形狀獲取所述形狀的外接圓;計算單元,用于計算所述外接圓的圓心到所述異型出口截面的形狀的距離;選擇單元,用于從超聲速流道的流線族中選擇與所述距離相對應(yīng)的流線,其中,選擇出的流線構(gòu)成所述超聲速流道的壁面的曲面形狀。
8.根據(jù)權(quán)利要求7所述的裝置,其特征在于,所述計算單元包括第一獲取模塊,用于獲取所述圓心到所述外接圓的圓周上的每個點的連線;第二獲取模塊,用于獲取所述連線與所述形狀的邊的第一交點;計算模塊,用于計算所述圓心到所述第一交點的距離,其中,所述圓心到所述第一交點的距離為所述外接圓的圓心到所述異型出口截面的形狀的距離。
9.根據(jù)權(quán)利要求7所述的裝置,其特征在于,所述選擇單元包括第三獲取模塊,用于獲取所述流線族中的各流線與所述超聲速流道的出口連接線之間的第二交點,其中,所述出口連接線為所述超聲速流道的出口點與該出口點相對于所述超聲速流道的中心線的鏡像點的連接線;第四獲取模塊,用于獲取所述中心線與所述超聲速流道的出口連接線之間的第三交占.計算模塊,用于計算所述第二交點與所述第三交點之間的距離;選擇模塊,用于選擇所述第二交點與所述第三交點之間的距離為所述外接圓的圓心到所述異型出口截面的形狀的距離的所述第二交點對應(yīng)的流線。
10.根據(jù)權(quán)利要求9所述的裝置,其特征在于,所述第三獲取模塊包括第一確定子模塊,用于根據(jù)預(yù)設(shè)的異型出口截面的進(jìn)口位置和出口位置以及馬赫數(shù)約束確定所述中心線的坐標(biāo)和長度;第二確定子模塊,用于根據(jù)馬赫數(shù)分布以及出口流場設(shè)計約束確定所述出口點;生成子模塊,用于根據(jù)所述中心線、所述特征線法以及確定出的出口點生成所述超聲速流道的曲線型面;鏡像子模塊,用于將所述曲線型面相對于所述中心線進(jìn)行鏡像得到所述出口點的鏡像
全文摘要
本發(fā)明公開了一種超聲速流道的實現(xiàn)方法及裝置,其中,該方法包括根據(jù)預(yù)設(shè)的異型出口截面的形狀獲取該形狀的外接圓;計算外接圓的圓心到異型出口截面的形狀的距離;從超聲速流道的流線族中選擇與該距離相對應(yīng)的流線,其中,選擇出的流線構(gòu)成超聲速流道的壁面的曲面形狀。本發(fā)明解決了相關(guān)技術(shù)中的無法設(shè)計出口截面為異型的超聲速流道問題,達(dá)到了實現(xiàn)超聲速流道的出口截面為異型的技術(shù)效果。
文檔編號F15D1/02GK102606565SQ20121010905
公開日2012年7月25日 申請日期2012年4月13日 優(yōu)先權(quán)日2012年4月13日
發(fā)明者王振國, 趙玉新 申請人:中國人民解放軍國防科學(xué)技術(shù)大學(xué)