• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    槍管內(nèi)膛損傷對(duì)彈頭外彈道過程的影響研究

    2020-07-15 07:42:12沈超周克棟赫雷陸野李峻松
    關(guān)鍵詞:內(nèi)膛彈孔彈頭

    沈超,周克棟,赫雷,陸野,李峻松

    (1.南京理工大學(xué) 機(jī)械工程學(xué)院,江蘇,南京 210094;2.中國兵器工業(yè)第208研究所,北京 102202)

    槍彈發(fā)射過程中,膛線與彈頭殼之間緊密接觸以減少火藥氣體的泄露,并且陽線會(huì)在彈頭上形成若干道刻槽,具有一定纏度的陽線與刻槽之間的接觸力,可以分解為沿槍管軸向的擠進(jìn)阻力和沿周向的導(dǎo)轉(zhuǎn)力,后者使得彈頭在沿槍管軸向運(yùn)動(dòng)的同時(shí)產(chǎn)生高速自轉(zhuǎn),以保證彈頭出膛后的飛行穩(wěn)定性[1-2].線膛的結(jié)構(gòu)參數(shù)(如陰、陽線寬度,陽線高度,纏度等)直接決定了彈頭的表面形貌.隨著射彈量的增加,槍管內(nèi)膛磨損導(dǎo)致的陽線高度、寬度等的降低會(huì)使得彈頭殼上的刻槽高度、寬度等隨之降低.此外,內(nèi)膛表面還分布著裂紋,燒蝕坑和鉻層剝落等損傷形式,在彈頭高速運(yùn)動(dòng)過程中會(huì)劃傷彈頭殼,破壞彈頭殼的表面完整性[3].因此,不同壽命階段槍管所發(fā)射彈頭的表面形貌是有差異的,彈頭的氣動(dòng)特性必然會(huì)隨之發(fā)生改變.研究槍管內(nèi)膛損傷對(duì)彈頭氣動(dòng)特性的影響規(guī)律,是獲得射擊精度、橫彈孔率等(槍管壽終與否的直接判據(jù))隨射彈量增加的變化規(guī)律的基礎(chǔ),對(duì)于加深對(duì)槍管壽終機(jī)理的認(rèn)識(shí)有重要意義.

    國內(nèi)外學(xué)者們對(duì)旋轉(zhuǎn)彈頭或彈丸氣動(dòng)特性和外彈道過程的實(shí)驗(yàn)測(cè)試與數(shù)值計(jì)算方面進(jìn)行了較為深入的研究,并取得了一定的成果.De Spirito等[4-5]使用雷諾平均和大渦模擬(RANS/LES)混合模型獲得了穩(wěn)定旋轉(zhuǎn)的M910彈丸在亞聲速、跨聲速及超聲速來流下的氣動(dòng)特性參數(shù),并通過與實(shí)驗(yàn)結(jié)果的對(duì)比,得到了不同的數(shù)值模擬方法對(duì)彈丸各氣動(dòng)特性參數(shù)預(yù)測(cè)的準(zhǔn)確度.Silton等[6]采用實(shí)驗(yàn)測(cè)試與數(shù)值模擬相結(jié)合的方法研究了膛線在小口徑彈丸上留下的刻槽對(duì)彈丸氣動(dòng)特性的影響,發(fā)現(xiàn)刻槽對(duì)不同氣動(dòng)參數(shù)的影響程度也存在較大差異.Jacobson 等[7]采用有限差分法計(jì)算了旋成體在超聲速流中的馬格努斯效應(yīng),研究了各種邊界層狀態(tài)下馬赫數(shù)、攻角和轉(zhuǎn)速等變化對(duì)馬格努斯效應(yīng)的影響.孟鵬等[8]基于SSTk-ω湍流模型對(duì)3維Navier-Stokes方程進(jìn)行求解,通過數(shù)值模擬研究了彈帶對(duì)高速旋轉(zhuǎn)彈丸氣動(dòng)特性的影響規(guī)律;陳亮等[9]結(jié)合數(shù)值仿真及近似理論方法分析了彈箭高速滾轉(zhuǎn)運(yùn)動(dòng)對(duì)其氣動(dòng)參數(shù)的影響,并獲得了不同馬赫數(shù)下的平衡轉(zhuǎn)速隨尾翼結(jié)構(gòu)參數(shù)的變化規(guī)律.趙清杰等[10]提出了一種新的大迎角非線性非定常氣動(dòng)力及氣動(dòng)力矩建模方法,降低了氣動(dòng)參數(shù)的預(yù)測(cè)結(jié)果與真實(shí)數(shù)據(jù)的誤差.

    對(duì)于旋轉(zhuǎn)彈頭或彈丸的氣動(dòng)特性,現(xiàn)有的試驗(yàn)分析及數(shù)值模擬大多研究的是無膛線刻槽、無損傷的彈頭或彈丸的氣動(dòng)特性,也有部分學(xué)者考慮了膛線在彈頭表面形成的刻槽對(duì)彈頭氣動(dòng)特性的影響,但均未對(duì)不同壽命階段的含內(nèi)膛損傷槍管所發(fā)射的彈頭的氣動(dòng)特性及外彈道過程受內(nèi)膛損傷的影響規(guī)律進(jìn)行深入分析,而是用初始狀態(tài)下彈頭的氣動(dòng)力參數(shù)去計(jì)算損傷槍管所發(fā)射彈頭的外彈道過程參數(shù),這種處理方法的準(zhǔn)確度及可信度均較低.

    為了更深入地研究槍管內(nèi)膛損傷導(dǎo)致槍管壽終的機(jī)理,本文基于某12.7 mm機(jī)槍槍管壽命試驗(yàn),研究對(duì)象為12.7 mm重機(jī)槍槍管在1個(gè)壽命周期中的4個(gè)不同的壽命階段所發(fā)射的具有不同表面形貌的彈頭,通過對(duì)4個(gè)階段槍管的完整內(nèi)彈道過程進(jìn)行仿真計(jì)算,獲得了各階段槍管所發(fā)射彈頭出膛時(shí)的表面形貌狀態(tài).基于SSTk-ω湍流模型對(duì)3維Navier-Stokes方程進(jìn)行雷諾時(shí)間平均,采用滑移網(wǎng)格技術(shù)處理彈頭的高速旋轉(zhuǎn)所引起的動(dòng)邊界問題,對(duì)4種不同表面形貌的彈頭的氣動(dòng)特性進(jìn)行了數(shù)值模擬;采用均勻設(shè)計(jì)方法安排隨機(jī)因素影響下的彈頭內(nèi)彈道計(jì)算過程,獲得了各壽命階段槍管發(fā)射彈頭的膛口擾動(dòng)隨機(jī)響應(yīng)狀態(tài);基于彈頭的6自由度剛體外彈道方程對(duì)其外彈道過程進(jìn)行了求解,結(jié)合槍管壽命試驗(yàn)結(jié)果對(duì)內(nèi)膛損傷導(dǎo)致彈頭橢圓彈孔率超標(biāo)進(jìn)而導(dǎo)致槍管壽終的過程進(jìn)行了深入研究.

    1 數(shù)值計(jì)算方法

    本文以3維Navier-Stokes方程(簡(jiǎn)稱N-S方程)為基礎(chǔ),采用剪應(yīng)力輸運(yùn)湍流模型(shear stress transportk-ω,簡(jiǎn)稱SSTk-ω),針對(duì)槍管壽命試驗(yàn)中彈頭初速變化較小,而運(yùn)動(dòng)姿態(tài)變化較大的情況[11],對(duì)4個(gè)壽命階段的槍管所發(fā)射的具有不同表面形貌的彈頭的氣動(dòng)特性及外部流場(chǎng)分布情況進(jìn)行數(shù)值模擬.由于膛線對(duì)彈頭的導(dǎo)轉(zhuǎn)力作用,彈頭出膛時(shí)轉(zhuǎn)速可以達(dá)到甚至超過10 000 rad/s,為準(zhǔn)確模擬彈頭旋轉(zhuǎn)導(dǎo)致的非定常流對(duì)其氣動(dòng)特性的影響,本文采用滑移網(wǎng)格方法處理彈頭高速旋轉(zhuǎn)所引起的動(dòng)邊界問題.

    采用滑移網(wǎng)格技術(shù)需將計(jì)算域分為兩個(gè)區(qū)域,即包圍彈頭并隨彈頭以相同角速度旋轉(zhuǎn)的內(nèi)部旋轉(zhuǎn)區(qū)和該內(nèi)部旋轉(zhuǎn)區(qū)之外的外部固定區(qū)(內(nèi)部旋轉(zhuǎn)區(qū)相對(duì)于外部固定區(qū)產(chǎn)生滑移運(yùn)動(dòng)),兩區(qū)域之間有一對(duì)交界面,交界面上的網(wǎng)格節(jié)點(diǎn)不需要相互重合,在兩區(qū)域的滑移交界面上采用數(shù)值插值以保證區(qū)域間的通量守恒.這樣處理既便于在流動(dòng)核心區(qū)域建立更密的網(wǎng)格,又可以降低遠(yuǎn)離核心區(qū)域的網(wǎng)格數(shù)量,且內(nèi)部旋轉(zhuǎn)區(qū)網(wǎng)格在運(yùn)動(dòng)過程中不發(fā)生變形,求解精度高,計(jì)算速度快.

    1.1 控制方程

    采用滑移網(wǎng)格技術(shù)時(shí),內(nèi)部旋轉(zhuǎn)區(qū)相對(duì)外部固定區(qū)做定常運(yùn)動(dòng)以模擬彈頭的旋轉(zhuǎn),內(nèi)部旋轉(zhuǎn)區(qū)網(wǎng)格轉(zhuǎn)速與彈頭轉(zhuǎn)速相等,外部固定區(qū)網(wǎng)格保持不動(dòng),且由于旋轉(zhuǎn)彈頭飛行速度可達(dá)數(shù)倍聲速,空氣的壓縮性不可忽略,因此N-S方程應(yīng)用于求解滑移網(wǎng)格流場(chǎng)域時(shí)的基本控制方程組可表示為

    (1)

    式中:V為控制體體積;ρ為空氣密度;?V為控制體的邊界面;u,ug分別為流速矢量和內(nèi)部旋轉(zhuǎn)區(qū)旋轉(zhuǎn)速度矢量;A為控制體的表面區(qū)域矢量;S為比焓源項(xiàng),由黏性耗散、可壓縮性和輻射等引起.

    1.2 湍流模型

    本文采用剪應(yīng)力輸運(yùn)SSTk-ω湍流模型,該模型是標(biāo)準(zhǔn)k-ω模型修正后的兩方程模型,綜合了標(biāo)準(zhǔn)k-ω模型能夠很好地模擬邊界層內(nèi)的低雷諾數(shù)流動(dòng),以及標(biāo)準(zhǔn)k-ε模型能夠很好地模擬邊界層外的充分發(fā)展的湍流流動(dòng)的優(yōu)勢(shì),能適應(yīng)壓力梯度變化的各種復(fù)雜物理現(xiàn)象.該模型對(duì)湍流黏性系數(shù)進(jìn)行了修正,且考慮了湍流剪切應(yīng)力從而不會(huì)對(duì)渦流黏度進(jìn)行過度預(yù)測(cè),因而能更好地適用于跨聲速及超聲速來流下的彈頭擾流場(chǎng)的模擬.

    SSTk-ω湍流模型的數(shù)學(xué)表達(dá)式為

    (2)

    (3)

    1.3 計(jì)算模型及邊界條件

    1.3.1計(jì)算模型

    計(jì)算模型為帶膛線刻槽的某12.7 mm穿甲燃燒彈,該彈頭主要包括弧形部,圓柱部和尾錐部3個(gè)部分,如圖1所示,圓柱部直徑D=13 mm.如圖2所示,將計(jì)算域分為3部分,從內(nèi)到外分別為邊界層、內(nèi)部旋轉(zhuǎn)區(qū)和外部固定區(qū).

    由于彈頭表面刻槽增加了其幾何形狀的復(fù)雜性,為增強(qiáng)解的收斂性與穩(wěn)定性,全彈面網(wǎng)格及全計(jì)算域空間網(wǎng)格分別采用正交性較好的四邊形網(wǎng)格和六面體網(wǎng)格,并對(duì)有刻槽的彈頭圓柱部網(wǎng)格進(jìn)行了局部加密.綜合多次計(jì)算并收斂后的結(jié)果,對(duì)邊界層劃分了20層網(wǎng)格,第一層網(wǎng)格厚度滿足y+≤0.5,相鄰網(wǎng)格厚度沿徑向的增加率小于1.15,邊界層、內(nèi)部旋轉(zhuǎn)區(qū)和外部固定區(qū)網(wǎng)格數(shù)分別為100萬、140萬和90萬.

    1.3.2邊界條件

    根據(jù)槍管壽命試驗(yàn)及文獻(xiàn)[11],一個(gè)槍管壽命周期內(nèi),彈頭的初速變化很小,試驗(yàn)時(shí)采用的100 m靶道內(nèi)彈頭的速度下降值也很小,而彈頭出膛時(shí)的運(yùn)動(dòng)姿態(tài)改變更為明顯,后者是導(dǎo)致100 m靶處橢圓彈孔率超標(biāo)(長(zhǎng)軸與短軸比大于1.25)進(jìn)而判斷槍管壽終的主要原因.彈孔的長(zhǎng)短軸之比由入射角、攻角和彈頭的外形尺寸決定,自動(dòng)武器多采用低伸彈道,100 m范圍內(nèi)幾乎相當(dāng)于是平射,在入射角很小時(shí),根據(jù)幾何關(guān)系可得,當(dāng)彈頭攻角大于3°時(shí),就已滿足彈孔長(zhǎng)軸與短軸之比大于1.25.因此,本文采用的計(jì)算條件為:來流馬赫數(shù)Ma為2.2,攻角α為0°~3°,量綱一化轉(zhuǎn)速ω*=0.215(ω*=ωD/v∞),其中:ω為彈頭旋轉(zhuǎn)角速度;D為彈頭圓柱部直徑;v∞為來流速度.

    彈頭附近網(wǎng)格劃分如圖3所示,彈頭表面采用無滑移壁面邊界條件,內(nèi)部旋轉(zhuǎn)區(qū)與彈頭壁面相關(guān)聯(lián),隨壁面以相同角速度旋轉(zhuǎn);內(nèi)部旋轉(zhuǎn)區(qū)與外部固定區(qū)通過交界面?zhèn)鬟f數(shù)據(jù),交界面采用滑移邊界條件;外部固定區(qū)的外邊界采用壓力遠(yuǎn)場(chǎng)邊界條件.

    1.4 6自由度剛體外彈道方程

    由于12.7 mm彈頭的長(zhǎng)徑比較小(彈長(zhǎng)與彈徑之比約為5),其外彈道過程可以近似為6自由度的剛體外彈道過程,彈頭在外彈道飛行過程中受到的外力主要有重力、空氣阻力、升力、馬格努斯力等,外力矩主要有俯仰力矩、馬格努斯力矩、極阻尼力矩、赤道阻尼力矩等.彈頭在諸外力、外力矩的作用下運(yùn)動(dòng)過程滿足的方程為

    (4)

    (5)

    (6)

    式中:m為彈頭質(zhì)量;v為彈頭質(zhì)心速度矢量;x為彈頭質(zhì)心空間位置矢量;F為合外力矢量;G為彈頭角動(dòng)量矢量;M為合外力矩矢量.建立合適的輔助坐標(biāo)系并將各矢量分解即可建立彈頭的6自由度剛體外彈道運(yùn)動(dòng)方程,將彈頭氣動(dòng)特性數(shù)值模型計(jì)算得到的各氣動(dòng)參數(shù)代入式(4)~(6),結(jié)合彈頭的初始運(yùn)動(dòng)狀態(tài)即可基于4階龍格-庫塔法編程求解彈頭的外彈道運(yùn)動(dòng)過程.

    2 損傷槍管所發(fā)射彈頭表面形貌的獲得方法

    要研究槍管內(nèi)膛損傷對(duì)彈頭氣動(dòng)特性的影響,首先需要獲得各壽命階段槍管所發(fā)射彈頭出膛時(shí)的表面形貌,具體流程如下:

    ① 根據(jù)系統(tǒng)的槍管壽命試驗(yàn)中獲得的多根相同材料、工藝和結(jié)構(gòu)的槍管在相同射擊規(guī)范下完成不同射彈量的射擊試驗(yàn)過程后的內(nèi)膛尺寸數(shù)據(jù)、內(nèi)窺和解剖過程的視頻及照片等,總結(jié)槍管在各壽命階段其內(nèi)膛損傷的主要形式及損傷沿槍管軸向的分布規(guī)律;

    ② 根據(jù)槍管的內(nèi)膛損傷形式及其分布規(guī)律因射彈量增加和軸向位置不同而導(dǎo)致的差異性,在槍管內(nèi)膛預(yù)置諸如裂紋、燒蝕坑及鍍鉻層剝落等損傷,完成不同壽命階段損傷槍管有限元模型的建立.

    ③ 考慮內(nèi)膛損傷導(dǎo)致的彈后空間的改變,建立彈頭在損傷槍管內(nèi)完成擠進(jìn)和沿線膛運(yùn)動(dòng)直至出膛的熱力耦合模型,并將計(jì)算結(jié)果與試驗(yàn)結(jié)果進(jìn)行對(duì)比以驗(yàn)證模型的準(zhǔn)確性.

    ④ 基于該熱力耦合模型對(duì)彈頭沿射彈量分別為0,1 400,3 000和6 000發(fā)(分別對(duì)應(yīng)了槍管的無損傷狀態(tài)、壽命中前期狀態(tài)、壽命中期狀態(tài)和壽終狀態(tài))的槍管內(nèi)膛運(yùn)動(dòng)直至出膛的過程進(jìn)行仿真分析,獲得4個(gè)壽命階段槍管所發(fā)射彈頭出膛時(shí)的表面形貌狀態(tài)和運(yùn)動(dòng)姿態(tài)等.

    內(nèi)膛損傷數(shù)據(jù)的分析及損傷形式和分布規(guī)律的總結(jié)過程、有限元模型的建立過程及彈頭沿4個(gè)壽命階段槍管運(yùn)動(dòng)的內(nèi)彈道計(jì)算過程詳見參考文獻(xiàn)[11].最終獲得了上述4個(gè)壽命階段的槍管所發(fā)射彈頭的表面形貌狀態(tài)如圖4所示(右側(cè)為彈頭頭部方向)以及彈頭表面刻槽尺寸如表1所示.

    表1 彈頭出膛時(shí)表面刻槽尺寸

    將有限元分析結(jié)果中的彈頭網(wǎng)格導(dǎo)入Hypermesh或ICEM-CFD等軟件進(jìn)行流場(chǎng)計(jì)算域網(wǎng)格劃分后,即可利用Fluent對(duì)4種不同表面形貌的彈頭的氣動(dòng)特性進(jìn)行仿真分析.

    3 數(shù)值計(jì)算模型的實(shí)驗(yàn)驗(yàn)證

    為驗(yàn)證本文所建立的高速旋轉(zhuǎn)彈頭氣動(dòng)特性數(shù)值計(jì)算方法的有效性和準(zhǔn)確性,選取DeSpirito在文獻(xiàn)[4]及Plostins在文獻(xiàn)[12]風(fēng)洞實(shí)驗(yàn)中采用的彈頭外形,對(duì)其進(jìn)行氣動(dòng)特性的數(shù)值仿真,并將仿真結(jié)果與風(fēng)洞實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比.DeSpirito及Plostins實(shí)驗(yàn)中采用的彈頭外形結(jié)構(gòu)與本文中的12.7 mm彈頭相似,如圖5所示.

    數(shù)值計(jì)算條件為:來流馬赫數(shù)Ma在0.6~4.5之間,攻角分別為0°和3°,彈頭量綱一化轉(zhuǎn)速ω*=0.17,力矩參考點(diǎn)取在彈頭質(zhì)心,距彈頭頭部49.9 mm.選取數(shù)值計(jì)算結(jié)果中攻角為0°時(shí)的阻力系數(shù)以及攻角為3°時(shí)的俯仰力矩系數(shù)導(dǎo)數(shù)及法向力壓力中心位置等氣動(dòng)特性參數(shù),其數(shù)值模擬結(jié)果與實(shí)驗(yàn)結(jié)果的比較情況如圖6~圖8所示.從3個(gè)氣動(dòng)特性參數(shù)的對(duì)比曲線圖上可以看出,數(shù)值計(jì)算得到的氣動(dòng)參數(shù)的變化規(guī)律與風(fēng)洞試驗(yàn)得到的結(jié)果相一致,并且數(shù)值計(jì)算結(jié)果與風(fēng)洞試驗(yàn)所獲得的相應(yīng)氣動(dòng)參數(shù)值的最大誤差均在10%以內(nèi),這說明本文所建立的高速旋轉(zhuǎn)彈頭的氣動(dòng)特性數(shù)值計(jì)算模型是準(zhǔn)確可信的.

    4 氣動(dòng)特性計(jì)算結(jié)果及外彈道過程分析

    為了更深入地研究槍管內(nèi)膛損傷對(duì)彈頭外彈道過程的影響,本節(jié)基于上述數(shù)值計(jì)算方法獲得了4個(gè)壽命階段的槍管所發(fā)射的具有不同表面形貌的彈頭的氣動(dòng)特性參數(shù)之間的差異及變化規(guī)律,并結(jié)合6自由度剛體外彈道方程,分析了隨機(jī)擾動(dòng)作用下,12.7 mm機(jī)槍槍管處于4個(gè)壽命階段時(shí),其發(fā)射的彈頭在100 m靶處的散布圓半徑和橢圓彈孔率(彈孔的長(zhǎng)軸與短軸之比大于1.25即為橢圓彈孔).

    4.1 內(nèi)膛損傷對(duì)彈頭的氣動(dòng)特性的影響

    將4個(gè)不同表面形貌的彈頭依次編號(hào)為1,2,3和4號(hào),分別對(duì)應(yīng)于射彈量0,1 400,3 000和6 000發(fā)的槍管所發(fā)射的彈頭.

    圖9給出了4個(gè)彈頭在Ma=2.2,ω*=0.215的飛行狀態(tài)下,法向力壓力中心距離彈頭質(zhì)心的距離(以質(zhì)心為原點(diǎn),指向彈頭頭部方向?yàn)檎?隨攻角α的變化規(guī)律.從圖9中可以看出,4個(gè)不同表面形貌的彈頭的壓心均在質(zhì)心前方,且隨著攻角的增加,壓心總體均呈現(xiàn)后移的趨勢(shì).此外還可以看出,損傷程度越嚴(yán)重的槍管所發(fā)射的彈頭,其壓心位置就越靠近彈頭頭部,這一現(xiàn)象在身管壽命的中后期(對(duì)應(yīng)于3,4號(hào)彈頭)尤為明顯,根據(jù)文獻(xiàn)[13]中彈丸急螺穩(wěn)定系數(shù)計(jì)算公式(見式(7))可知,由內(nèi)膛損傷的加劇導(dǎo)致的彈頭壓心的前移會(huì)加大壓心質(zhì)心的距離,降低彈頭的飛行穩(wěn)定行,使散布直徑和橢圓彈孔率增加.

    (7)

    式中:σ為急螺穩(wěn)定系數(shù);h為壓心距質(zhì)心距離;D為彈頭直徑;H(y)為氣重函數(shù),取決于彈道高度;v為彈頭速度;kmz(Ma)為翻轉(zhuǎn)力矩的速度函數(shù);α為攻角;A、B分別為彈頭的極轉(zhuǎn)動(dòng)慣量和赤道轉(zhuǎn)動(dòng)慣量;ω為彈頭轉(zhuǎn)速;g為重力加速度.

    圖10~圖12分別對(duì)應(yīng)為來流馬赫數(shù)為2.2、量綱一轉(zhuǎn)速為0.215下的4個(gè)不同表面形貌彈頭的阻力系數(shù)、升力系數(shù)和馬格努斯力矩系數(shù)(規(guī)定馬格努斯力矩方向豎直向上為正)隨攻角變化規(guī)律.由圖10可以看出,隨著彈頭表面形貌較初始彈頭的改變量的增大,阻力系數(shù)也逐漸增大,且攻角越大,阻力系數(shù)的增大量也越大,這主要是因?yàn)殡S著槍管內(nèi)膛損傷的不斷發(fā)展,其發(fā)射的彈頭表面刻槽的深度在不斷降低[11],導(dǎo)致了彈頭阻力面積的增大,進(jìn)而引起了阻力系數(shù)的上升.由圖11、圖12可以看出,從絕對(duì)數(shù)值上來看,各壽命階段槍管所發(fā)射彈頭的升力系數(shù)和馬格努斯力矩系數(shù)在各攻角下均隨射彈量增加而增加,且攻角越大不同彈頭間氣動(dòng)參數(shù)差值越大,其中升力系數(shù)的增加原因與阻力系數(shù)相同,即刻槽變淺導(dǎo)致的受力面積的增加導(dǎo)致了氣動(dòng)參數(shù)的增加;馬格努斯力矩系數(shù)增加的主要原因?yàn)閮?nèi)膛損傷槍管所發(fā)射的彈頭,其表面除規(guī)則的陽線刻槽外,還會(huì)出現(xiàn)由內(nèi)膛損傷導(dǎo)致的不規(guī)則分布的局部刻痕,內(nèi)膛損傷還會(huì)劃破彈頭被甲材料,破壞其表面完整性,這會(huì)使得彈頭高速旋轉(zhuǎn)時(shí)邊界層厚度的非對(duì)稱畸變、徑向壓力梯度的非對(duì)稱畸變等增大,從而引起馬格努斯力矩系數(shù)的增大.

    結(jié)合圖10~圖12及壽終槍管和無損傷槍管所發(fā)射彈頭氣動(dòng)參數(shù)在α=3°時(shí)對(duì)比表(表2)可以看出,4號(hào)彈頭的各氣動(dòng)參數(shù)相對(duì)1號(hào)彈頭的改變量在α=3°達(dá)到最大,相對(duì)變化量均在10%以上,對(duì)彈頭外彈道過程的影響不可忽略.

    表2 攻角為3°時(shí)1號(hào)、4號(hào)彈頭氣動(dòng)參數(shù)對(duì)比

    4.2 外彈道過程仿真分析

    彈頭發(fā)射過程中,造成彈頭出膛時(shí)運(yùn)動(dòng)狀態(tài)產(chǎn)生隨機(jī)波動(dòng)的原因有很多,對(duì)于某一支特定的自動(dòng)武器來說,將其架在固定槍架上進(jìn)行射擊時(shí),槍械自身的隨機(jī)因素對(duì)彈頭的影響已降至較低,使彈頭出膛狀態(tài)產(chǎn)生隨機(jī)波動(dòng)進(jìn)而產(chǎn)生射彈散布的主要原因在于彈藥自身參數(shù)在一定范圍內(nèi)的波動(dòng),主要有彈頭外形尺寸偏差、彈頭質(zhì)量及質(zhì)量偏心偏差、彈頭殼材料力學(xué)性能偏差、發(fā)射藥主控參數(shù)偏差等.

    本文針對(duì)上述5個(gè)會(huì)對(duì)彈頭出膛擾動(dòng)隨機(jī)響應(yīng)產(chǎn)生較大影響的因素,依據(jù)設(shè)計(jì)要求和工程實(shí)踐經(jīng)驗(yàn),給出了各因素隨機(jī)波動(dòng)的范圍如表3所示.5個(gè)因素分別取11個(gè)水平值,若按照正交試驗(yàn)方法進(jìn)行因素不同水平值間的組合,則每個(gè)壽命階段的槍管均需要進(jìn)行112次理論及仿真計(jì)算.為避免參數(shù)組合情況太多帶來的時(shí)間成本的增加,同時(shí)又要充分考慮各因素的各水平值對(duì)彈頭初始擾動(dòng)的影響,采用均勻設(shè)計(jì)方法進(jìn)行各因素間的組合[14],均勻設(shè)計(jì)能在較好地反映試驗(yàn)體系主要特征的前提下大大地減少試驗(yàn)次數(shù),且相對(duì)于正交設(shè)計(jì)不會(huì)產(chǎn)生大的偏差.基于均勻設(shè)計(jì)方法獲得了11組彈頭主要隨機(jī)因素的組合狀態(tài),并將這11個(gè)不同初始狀態(tài)的彈頭分別代入第2節(jié)建立的4個(gè)壽命階段槍管內(nèi)彈道過程計(jì)算模型中,獲得了4個(gè)壽命階段的槍管所發(fā)射彈頭出膛擾動(dòng)的隨機(jī)響應(yīng)結(jié)果.

    表3 彈藥隨機(jī)因素取值范圍表

    將1~4號(hào)彈頭的各11組膛口隨機(jī)擾動(dòng)狀態(tài)作為初值代入式(4)~(6),結(jié)合氣動(dòng)特性參數(shù)對(duì)各階段彈頭的外彈道過程進(jìn)行求解,獲得了4個(gè)壽命階段槍管所發(fā)射的各11發(fā)彈頭在100 m處的散布密集度R50和橢圓彈孔率.判斷彈孔是否為橢圓彈孔的方法如圖13所示,圖13截面為彈頭軸線和質(zhì)心速度所在平面(橢圓彈孔長(zhǎng)軸所在平面),由6自由度外彈道方程可以得到100 m處彈頭的入射角值(質(zhì)心速度和靶面法線夾角)及攻角值(質(zhì)心速度和彈頭軸線夾角),根據(jù)幾何關(guān)系即可求得彈孔的長(zhǎng)軸長(zhǎng)度,短軸長(zhǎng)與彈頭圓柱部直徑相等,二者之比大于1.25即為橢圓彈孔.

    各階段槍管R50、橢圓彈孔率與槍管壽命試驗(yàn)結(jié)果的對(duì)比情況如表4所示(壽命試驗(yàn)的R50和橢圓彈孔率取多次試驗(yàn)的平均值,每次發(fā)射20發(fā)彈頭),以壽終槍管為例,其計(jì)算得到100 m靶處彈著點(diǎn)分布情況如圖14所示,可以看出橢圓彈孔往往也偏離散布中心較遠(yuǎn).槍管壽終的判斷標(biāo)準(zhǔn)為R50≥30 cm或者橢圓彈孔率≥50%.

    表4 各階段槍管散布密集度和橢圓彈孔率

    Tab.4 Dispersion intensity and ratio of elliptical bullet hole at different life period

    槍管射彈量/發(fā)R50/cm橢圓彈孔率計(jì)算值試驗(yàn)平均值計(jì)算值試驗(yàn)平均值012.814.600140013.013.400300015.115.600600021.722.854.5%60%

    從表4可以看出,散布密集度R50和橢圓彈孔率的計(jì)算值均與槍管壽命試驗(yàn)實(shí)測(cè)值較好地符合,這進(jìn)一步驗(yàn)證了損傷彈頭的獲得方法及氣動(dòng)特性參數(shù)的計(jì)算模型的準(zhǔn)確性.從表4中還可以看到,試驗(yàn)過程中無損傷槍管的R50值反而大于射彈量1 400發(fā)的槍管的相應(yīng)值,而數(shù)值仿真結(jié)果則恰恰相反,這也導(dǎo)致了無損傷槍管R50的數(shù)值計(jì)算結(jié)果與試驗(yàn)結(jié)果誤差較大,這是因?yàn)樵囼?yàn)時(shí)的無損傷槍管內(nèi)膛表面存在制造和鍍鉻時(shí)留下的細(xì)小毛刺等缺陷,導(dǎo)致了彈頭出膛擾動(dòng)增大,經(jīng)過一定次數(shù)的射擊后,細(xì)小毛刺被磨去,而且內(nèi)膛的損傷還未充分形成,因此在身管壽終的中前期,R50往往會(huì)出現(xiàn)不升反降的現(xiàn)象.計(jì)算值與試驗(yàn)情況均表明,該12.7 mm重機(jī)槍槍管壽終表現(xiàn)為橢圓彈孔率的超標(biāo),槍管在壽命中前期R50的增大量較小,且基本不會(huì)出現(xiàn)橢圓彈孔,在槍管壽命的后期,R50顯著增加以及橢圓彈孔率的急劇增加導(dǎo)致了槍管的壽終,這是因?yàn)闃尮軌勖笃趶楊^出膛時(shí)的初始攻角等膛口擾動(dòng)量會(huì)顯著增加[11],各氣動(dòng)特性參數(shù)在大攻角時(shí)的明顯改變(如壓心前移,馬格努斯力矩系數(shù)的增大等)也會(huì)顯著降低彈頭外彈道過程的飛行穩(wěn)定性,從而進(jìn)一步增大了彈頭的攻角,形成了惡性循環(huán),導(dǎo)致了槍管的迅速壽終.

    5 結(jié) 論

    基于槍管壽命試驗(yàn)獲得的內(nèi)膛損傷數(shù)據(jù),建立了損傷槍管的有限元模型,獲得了12.7 mm槍管一個(gè)壽命周期內(nèi)4個(gè)壽命階段下所發(fā)射的彈頭出膛時(shí)的表面形貌狀態(tài),研究獲得了各階段彈頭的氣動(dòng)特性參數(shù)隨攻角的變化規(guī)律,結(jié)合彈頭的初始擾動(dòng)狀態(tài)對(duì)其6自由度外彈道過程進(jìn)行了理論計(jì)算與分析,獲得的結(jié)論如下.

    ① 內(nèi)膛損傷導(dǎo)致的彈頭表面形貌的改變會(huì)使得越靠近壽命末期的槍管,其發(fā)射的彈頭壓心越靠前,阻力系數(shù)、升力系數(shù)及馬格努斯力矩系數(shù)等也越大,且攻角越大,不同表面形貌的彈頭間氣動(dòng)參數(shù)的差異也越明顯.

    ② 壽終槍管所發(fā)射彈頭的主要?dú)鈩?dòng)參數(shù)相對(duì)于無損傷槍管的改變量最大可達(dá)10%以上,內(nèi)膛損傷導(dǎo)致的彈頭表面形貌的改變對(duì)彈頭氣動(dòng)特性及外彈道過程的影響不可忽略.

    ③ 試驗(yàn)和計(jì)算結(jié)果表明,彈頭在100 m處散布圓半徑R50和橢圓彈孔率在槍管壽命中前期不會(huì)明顯增加,槍管壽終往往表現(xiàn)為壽命后期散布圓半徑R50和橢圓彈孔率的急劇增大.

    ④ 彈頭初始擾動(dòng)在槍管壽命末期的迅速增大,以及大攻角下彈頭所受氣動(dòng)力壓力中心的前移和馬格努斯力矩的增大,是造成壽命末期槍管所發(fā)射彈頭飛行穩(wěn)定性降低及槍管迅速壽終的主要原因.

    ⑤ 本文計(jì)算獲得的100 m處散布圓半徑R50和橢圓彈孔率與多次試驗(yàn)的平均值吻合度較高,表明本文建立的彈頭氣動(dòng)特性模型和內(nèi)、外彈道模型合理、有效.

    猜你喜歡
    內(nèi)膛彈孔彈頭
    果樹整形修剪技術(shù)
    熱身管下小口徑槍彈彈頭殼材料對(duì)其膛內(nèi)運(yùn)動(dòng)的影響規(guī)律
    克服獼猴桃架面內(nèi)膛光禿的技術(shù)措施
    最后的彈孔
    基于視頻圖像處理的彈孔檢測(cè)提取算法研究?
    導(dǎo)彈突防概率解析與仿真計(jì)算結(jié)果差異分析
    柑橘內(nèi)膛腹接增產(chǎn)技術(shù)
    怎樣防止桃樹內(nèi)膛過早光禿
    閑趣一刻
    讀者(2010年2期)2010-02-11 11:51:19
    智力大魔方
    久久影院123| 精品人妻在线不人妻| 国产av码专区亚洲av| 亚洲美女黄色视频免费看| 大片电影免费在线观看免费| 亚洲国产精品成人久久小说| 一级爰片在线观看| 亚洲成色77777| 亚洲精品,欧美精品| 亚洲内射少妇av| 国语对白做爰xxxⅹ性视频网站| 特大巨黑吊av在线直播| 精品酒店卫生间| 欧美日韩在线观看h| av一本久久久久| 中文乱码字字幕精品一区二区三区| 欧美成人午夜免费资源| 国产av精品麻豆| 久久久久精品久久久久真实原创| 热re99久久国产66热| 三上悠亚av全集在线观看| 亚洲无线观看免费| 精品国产一区二区三区久久久樱花| 99久久精品一区二区三区| 亚洲精品乱久久久久久| 最新的欧美精品一区二区| 亚洲精品日韩在线中文字幕| 国产精品偷伦视频观看了| 国产欧美日韩综合在线一区二区| 桃花免费在线播放| 黑人猛操日本美女一级片| 校园人妻丝袜中文字幕| 午夜老司机福利剧场| 午夜激情av网站| 美女视频免费永久观看网站| 老女人水多毛片| 中国美白少妇内射xxxbb| 久久久久久久久久人人人人人人| 伊人亚洲综合成人网| 国产精品国产三级国产av玫瑰| 美女内射精品一级片tv| 一个人免费看片子| 国产av码专区亚洲av| 街头女战士在线观看网站| 国产精品久久久久久久久免| 一二三四中文在线观看免费高清| 最近中文字幕高清免费大全6| 91精品伊人久久大香线蕉| 日本欧美国产在线视频| 国产乱来视频区| 亚洲av在线观看美女高潮| h视频一区二区三区| 日本黄色日本黄色录像| 中文字幕人妻丝袜制服| 18在线观看网站| 一区二区av电影网| 国国产精品蜜臀av免费| videossex国产| 日韩一区二区三区影片| 精品久久久噜噜| 熟女电影av网| 日韩av免费高清视频| 五月天丁香电影| 黄色欧美视频在线观看| 人妻夜夜爽99麻豆av| 国产男人的电影天堂91| 久久久久国产精品人妻一区二区| 日日摸夜夜添夜夜添av毛片| 五月伊人婷婷丁香| av国产久精品久网站免费入址| 亚洲精品乱码久久久v下载方式| 女人久久www免费人成看片| 26uuu在线亚洲综合色| 人人妻人人添人人爽欧美一区卜| 蜜桃久久精品国产亚洲av| 五月开心婷婷网| 久久久久精品久久久久真实原创| 亚洲av免费高清在线观看| 男人添女人高潮全过程视频| 最新的欧美精品一区二区| 黑人欧美特级aaaaaa片| 欧美成人午夜免费资源| 大话2 男鬼变身卡| 国产视频首页在线观看| 亚洲国产成人一精品久久久| 久久久久久久久久人人人人人人| 久久久久视频综合| 97精品久久久久久久久久精品| 午夜免费男女啪啪视频观看| 制服诱惑二区| 2018国产大陆天天弄谢| 在现免费观看毛片| 涩涩av久久男人的天堂| 欧美激情极品国产一区二区三区 | 亚洲av福利一区| 久久久久网色| 久久亚洲国产成人精品v| 如日韩欧美国产精品一区二区三区 | 9色porny在线观看| 91在线精品国自产拍蜜月| 国产精品人妻久久久久久| av.在线天堂| 母亲3免费完整高清在线观看 | 91在线精品国自产拍蜜月| 热re99久久国产66热| 熟女人妻精品中文字幕| 丁香六月天网| 乱人伦中国视频| 亚洲国产毛片av蜜桃av| kizo精华| 99视频精品全部免费 在线| 如日韩欧美国产精品一区二区三区 | 18禁裸乳无遮挡动漫免费视频| 久久久久久久精品精品| 亚洲国产精品成人久久小说| av线在线观看网站| 尾随美女入室| 中文天堂在线官网| 精品久久久久久电影网| 99九九线精品视频在线观看视频| 丝袜脚勾引网站| 国产亚洲精品第一综合不卡 | 免费观看av网站的网址| av一本久久久久| 999精品在线视频| 性高湖久久久久久久久免费观看| 午夜久久久在线观看| 日本黄色片子视频| 人体艺术视频欧美日本| 亚洲精品国产色婷婷电影| 在线观看www视频免费| 九草在线视频观看| 久久久国产一区二区| 三级国产精品片| 久久鲁丝午夜福利片| 亚洲国产精品专区欧美| 亚洲欧美中文字幕日韩二区| 国产淫语在线视频| 婷婷色综合大香蕉| 色婷婷久久久亚洲欧美| 日韩,欧美,国产一区二区三区| 色婷婷av一区二区三区视频| 在线免费观看不下载黄p国产| 午夜福利,免费看| 国产精品秋霞免费鲁丝片| 国产亚洲精品久久久com| 精品久久久噜噜| 成年人免费黄色播放视频| 97精品久久久久久久久久精品| 亚洲精品乱码久久久v下载方式| a级片在线免费高清观看视频| 亚洲成人av在线免费| 草草在线视频免费看| 一区二区三区乱码不卡18| 精品人妻一区二区三区麻豆| 国产精品久久久久久久电影| 精品久久蜜臀av无| 免费久久久久久久精品成人欧美视频 | 亚洲国产精品国产精品| 国产在线视频一区二区| 久久99精品国语久久久| 欧美另类一区| 国产片特级美女逼逼视频| 91成人精品电影| 伦理电影大哥的女人| 亚洲天堂av无毛| 丰满迷人的少妇在线观看| 国产精品99久久99久久久不卡 | 永久免费av网站大全| 免费不卡的大黄色大毛片视频在线观看| 夜夜骑夜夜射夜夜干| 日韩欧美一区视频在线观看| 国产成人91sexporn| 亚洲国产色片| 国产免费现黄频在线看| 丰满乱子伦码专区| av在线老鸭窝| 国产精品无大码| videossex国产| 久热这里只有精品99| 亚洲熟女精品中文字幕| 春色校园在线视频观看| 久久精品国产亚洲av涩爱| 各种免费的搞黄视频| 青春草国产在线视频| 免费高清在线观看日韩| 丰满迷人的少妇在线观看| 十八禁高潮呻吟视频| av网站免费在线观看视频| 一区二区三区四区激情视频| 人妻 亚洲 视频| 日本午夜av视频| 亚洲四区av| 午夜福利,免费看| 最近最新中文字幕免费大全7| 人妻系列 视频| 国产高清不卡午夜福利| 亚洲情色 制服丝袜| 精品亚洲成a人片在线观看| 久久久午夜欧美精品| av国产久精品久网站免费入址| 午夜视频国产福利| 日韩大片免费观看网站| 18禁观看日本| 久久久午夜欧美精品| 欧美日韩国产mv在线观看视频| 18+在线观看网站| 亚洲欧美成人综合另类久久久| 久久久久网色| 99视频精品全部免费 在线| 国产精品一区二区在线不卡| 在线观看一区二区三区激情| 亚洲国产av影院在线观看| 国产成人aa在线观看| 大香蕉久久网| 国产精品女同一区二区软件| 久久久国产欧美日韩av| 免费人妻精品一区二区三区视频| 国产精品嫩草影院av在线观看| 国产探花极品一区二区| 欧美另类一区| 日韩视频在线欧美| www.色视频.com| 人妻系列 视频| 夜夜看夜夜爽夜夜摸| 久久精品人人爽人人爽视色| 亚洲国产精品专区欧美| 三上悠亚av全集在线观看| 久久热精品热| 国内精品宾馆在线| 性高湖久久久久久久久免费观看| 久久综合国产亚洲精品| 欧美97在线视频| av福利片在线| 少妇熟女欧美另类| 高清黄色对白视频在线免费看| 日日爽夜夜爽网站| 中国美白少妇内射xxxbb| 18禁裸乳无遮挡动漫免费视频| 蜜桃国产av成人99| 伊人亚洲综合成人网| 王馨瑶露胸无遮挡在线观看| www.av在线官网国产| 夫妻午夜视频| 国产成人精品一,二区| 美女主播在线视频| 一边摸一边做爽爽视频免费| 美女大奶头黄色视频| 国产精品久久久久久av不卡| 亚洲国产日韩一区二区| 黄片播放在线免费| 欧美变态另类bdsm刘玥| 精品久久久久久久久av| 少妇精品久久久久久久| 搡女人真爽免费视频火全软件| 欧美日韩国产mv在线观看视频| 亚洲欧美色中文字幕在线| 国产成人精品福利久久| 久久久午夜欧美精品| 91久久精品国产一区二区三区| 亚洲精品第二区| 99国产精品免费福利视频| 国产熟女午夜一区二区三区 | 亚洲精品美女久久av网站| 亚洲av成人精品一二三区| 亚洲综合精品二区| 男女无遮挡免费网站观看| 日韩一本色道免费dvd| 免费大片18禁| 亚洲av不卡在线观看| 国产高清不卡午夜福利| 午夜久久久在线观看| 亚洲av综合色区一区| 波野结衣二区三区在线| 精品少妇久久久久久888优播| 日韩欧美一区视频在线观看| xxx大片免费视频| 蜜桃久久精品国产亚洲av| 大香蕉久久成人网| 在线免费观看不下载黄p国产| 亚洲第一区二区三区不卡| 少妇猛男粗大的猛烈进出视频| 80岁老熟妇乱子伦牲交| 国产成人免费观看mmmm| 97在线人人人人妻| 菩萨蛮人人尽说江南好唐韦庄| 性色av一级| 在线看a的网站| 能在线免费看毛片的网站| 亚洲在久久综合| 伦理电影免费视频| 午夜免费男女啪啪视频观看| 一级毛片电影观看| 久热这里只有精品99| 精品久久蜜臀av无| 一级爰片在线观看| av又黄又爽大尺度在线免费看| 免费黄网站久久成人精品| 亚洲美女黄色视频免费看| 久久久久国产精品人妻一区二区| 久久久亚洲精品成人影院| 国产精品熟女久久久久浪| 亚洲av二区三区四区| 国产黄片视频在线免费观看| 99re6热这里在线精品视频| 国产精品欧美亚洲77777| 免费黄频网站在线观看国产| 日日撸夜夜添| av在线app专区| 国产一区二区三区av在线| 少妇熟女欧美另类| 日韩av不卡免费在线播放| 青青草视频在线视频观看| 国产在线一区二区三区精| 精品久久久久久久久亚洲| 午夜老司机福利剧场| 亚洲精品久久成人aⅴ小说 | 极品少妇高潮喷水抽搐| 日韩制服骚丝袜av| 乱人伦中国视频| 99热国产这里只有精品6| 国产69精品久久久久777片| 美女大奶头黄色视频| 色婷婷久久久亚洲欧美| 亚洲精品aⅴ在线观看| 国语对白做爰xxxⅹ性视频网站| 亚洲av二区三区四区| 免费观看无遮挡的男女| 欧美激情极品国产一区二区三区 | 飞空精品影院首页| 精品国产乱码久久久久久小说| 丰满饥渴人妻一区二区三| 午夜日本视频在线| 亚洲精品第二区| 久久免费观看电影| 狂野欧美激情性bbbbbb| 久久久久久久久久久久大奶| 久久久久国产精品人妻一区二区| 18在线观看网站| 少妇人妻精品综合一区二区| 亚洲精品亚洲一区二区| 在线免费观看不下载黄p国产| 22中文网久久字幕| 一级黄片播放器| 22中文网久久字幕| 少妇熟女欧美另类| av播播在线观看一区| 少妇 在线观看| 中文天堂在线官网| 免费黄色在线免费观看| 亚洲精品,欧美精品| 人妻系列 视频| 亚洲av国产av综合av卡| 国产精品国产三级国产专区5o| 国产黄色免费在线视频| 日本黄色片子视频| 婷婷色综合大香蕉| 国产深夜福利视频在线观看| 七月丁香在线播放| tube8黄色片| 亚洲精品色激情综合| 日本91视频免费播放| 精品一区二区免费观看| 老女人水多毛片| 99久久中文字幕三级久久日本| 亚洲av在线观看美女高潮| 精品一区二区免费观看| 又黄又爽又刺激的免费视频.| 各种免费的搞黄视频| 久久久久久久亚洲中文字幕| 精品久久久久久久久av| 免费看不卡的av| 男女高潮啪啪啪动态图| 在线播放无遮挡| 丰满迷人的少妇在线观看| 色5月婷婷丁香| 欧美日韩视频高清一区二区三区二| 亚洲无线观看免费| 一区二区三区四区激情视频| 欧美最新免费一区二区三区| 国语对白做爰xxxⅹ性视频网站| 熟女电影av网| 曰老女人黄片| 男女边摸边吃奶| 99热全是精品| 69精品国产乱码久久久| 国产爽快片一区二区三区| 26uuu在线亚洲综合色| 九九爱精品视频在线观看| 精品国产一区二区久久| 精品亚洲成a人片在线观看| 国产精品国产三级国产av玫瑰| 国产精品麻豆人妻色哟哟久久| 日韩中文字幕视频在线看片| 男女无遮挡免费网站观看| 男女边吃奶边做爰视频| 亚洲无线观看免费| 91久久精品国产一区二区成人| 国产免费现黄频在线看| 免费观看的影片在线观看| 亚洲第一区二区三区不卡| 久久久久网色| 免费高清在线观看视频在线观看| 女的被弄到高潮叫床怎么办| 亚洲美女视频黄频| 丝袜喷水一区| 十八禁网站网址无遮挡| 26uuu在线亚洲综合色| 18禁观看日本| 97超视频在线观看视频| 久久久精品免费免费高清| 制服人妻中文乱码| 岛国毛片在线播放| 午夜老司机福利剧场| 永久网站在线| 日韩 亚洲 欧美在线| 午夜日本视频在线| 如何舔出高潮| www.色视频.com| 亚洲av成人精品一区久久| 亚洲综合色网址| 男女无遮挡免费网站观看| 久久久精品94久久精品| 精品人妻熟女毛片av久久网站| 啦啦啦中文免费视频观看日本| 亚洲经典国产精华液单| 99久久精品国产国产毛片| av网站免费在线观看视频| 一边摸一边做爽爽视频免费| 最黄视频免费看| 国产深夜福利视频在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 考比视频在线观看| 久久精品久久精品一区二区三区| 大陆偷拍与自拍| 中文字幕精品免费在线观看视频 | 精品久久国产蜜桃| 国产精品99久久久久久久久| 国产精品久久久久久av不卡| 高清黄色对白视频在线免费看| 在线免费观看不下载黄p国产| 国产 一区精品| 欧美成人午夜免费资源| 欧美精品国产亚洲| 国产毛片在线视频| 久久午夜福利片| 制服丝袜香蕉在线| 三级国产精品片| 99热国产这里只有精品6| 又黄又爽又刺激的免费视频.| 五月玫瑰六月丁香| 赤兔流量卡办理| 日本爱情动作片www.在线观看| 久久国产精品大桥未久av| 色婷婷久久久亚洲欧美| 在线观看免费日韩欧美大片 | 亚洲av福利一区| 国产不卡av网站在线观看| 久久99热这里只频精品6学生| videosex国产| 日韩亚洲欧美综合| 日产精品乱码卡一卡2卡三| 美女视频免费永久观看网站| 成人国产av品久久久| 啦啦啦在线观看免费高清www| 人妻 亚洲 视频| 黄色怎么调成土黄色| 国产免费现黄频在线看| 曰老女人黄片| av免费观看日本| 在线观看人妻少妇| 狂野欧美激情性xxxx在线观看| 看免费成人av毛片| av不卡在线播放| 男女边摸边吃奶| 欧美日韩亚洲高清精品| 国产一区亚洲一区在线观看| 国产精品人妻久久久久久| 日日摸夜夜添夜夜爱| 全区人妻精品视频| 久久久久国产精品人妻一区二区| 少妇丰满av| 精品久久久噜噜| 亚洲国产精品一区三区| 黑丝袜美女国产一区| 水蜜桃什么品种好| 欧美bdsm另类| 美女中出高潮动态图| 久久影院123| 亚洲av二区三区四区| a级毛片黄视频| 免费人妻精品一区二区三区视频| 日日啪夜夜爽| 黑人猛操日本美女一级片| 国产欧美亚洲国产| 男人爽女人下面视频在线观看| 亚洲美女黄色视频免费看| 国产日韩欧美亚洲二区| 成人影院久久| 亚洲成人av在线免费| 狂野欧美白嫩少妇大欣赏| 久久精品久久久久久噜噜老黄| 国产色婷婷99| 啦啦啦中文免费视频观看日本| av.在线天堂| 99热这里只有是精品在线观看| 人人妻人人澡人人看| 日韩精品免费视频一区二区三区 | 日本欧美国产在线视频| 亚洲色图 男人天堂 中文字幕 | a级片在线免费高清观看视频| 嫩草影院入口| 亚洲精品亚洲一区二区| 久久精品国产a三级三级三级| 亚洲欧洲日产国产| 女人久久www免费人成看片| av女优亚洲男人天堂| 欧美人与性动交α欧美精品济南到 | 精品一区在线观看国产| 国产精品.久久久| 五月天丁香电影| 国产熟女午夜一区二区三区 | 国产成人av激情在线播放 | 美女福利国产在线| 日韩av不卡免费在线播放| 国产爽快片一区二区三区| 日韩av免费高清视频| 岛国毛片在线播放| 丝袜脚勾引网站| 少妇丰满av| 精品一品国产午夜福利视频| av在线播放精品| 永久网站在线| 欧美 日韩 精品 国产| 又粗又硬又长又爽又黄的视频| 国产在线视频一区二区| a级毛片黄视频| 26uuu在线亚洲综合色| 在线观看免费视频网站a站| 国产高清有码在线观看视频| 中文天堂在线官网| 欧美成人午夜免费资源| 亚洲欧洲精品一区二区精品久久久 | 一区二区av电影网| 国产成人午夜福利电影在线观看| 午夜老司机福利剧场| 大香蕉久久成人网| 亚洲av成人精品一区久久| 精品少妇久久久久久888优播| 久久99一区二区三区| 国产一级毛片在线| 亚洲综合色网址| 精品酒店卫生间| 狂野欧美白嫩少妇大欣赏| 少妇人妻久久综合中文| 亚洲无线观看免费| 日本欧美国产在线视频| 久久精品国产自在天天线| 人妻 亚洲 视频| 午夜91福利影院| 亚洲精品久久久久久婷婷小说| 高清在线视频一区二区三区| 美女xxoo啪啪120秒动态图| 免费观看av网站的网址| 飞空精品影院首页| 男女边吃奶边做爰视频| 色网站视频免费| 啦啦啦啦在线视频资源| 国产成人精品婷婷| videosex国产| 99热这里只有是精品在线观看| 久久久久视频综合| 岛国毛片在线播放| 午夜激情福利司机影院| 十分钟在线观看高清视频www| 欧美激情国产日韩精品一区| 在现免费观看毛片| 日韩不卡一区二区三区视频在线| 欧美变态另类bdsm刘玥| 日韩免费高清中文字幕av| 亚洲天堂av无毛| 一区二区日韩欧美中文字幕 | 一边摸一边做爽爽视频免费| 亚州av有码| 大香蕉久久网| 麻豆乱淫一区二区| 精品午夜福利在线看| 九草在线视频观看| 秋霞在线观看毛片| 国产欧美另类精品又又久久亚洲欧美| 成年人免费黄色播放视频| 最近的中文字幕免费完整| 亚洲综合色惰| 高清视频免费观看一区二区| 老司机影院成人| 亚洲国产av新网站| 中国美白少妇内射xxxbb| 五月天丁香电影| 亚洲国产av新网站| 欧美激情国产日韩精品一区| 黄色视频在线播放观看不卡| 桃花免费在线播放| 黄色配什么色好看| 久久精品国产亚洲网站| 亚洲综合色网址| 国产高清不卡午夜福利| 亚洲精品456在线播放app| 精品久久久久久久久av| 美女脱内裤让男人舔精品视频| 超碰97精品在线观看| 国产成人精品婷婷| 欧美激情国产日韩精品一区| 性色avwww在线观看| 汤姆久久久久久久影院中文字幕| 日本欧美视频一区|