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

    某大口徑機(jī)槍內(nèi)膛損傷對彈頭擠進(jìn)過程的影響研究

    2019-01-02 08:36:12沈超周克棟陸野喬自平
    兵工學(xué)報 2018年12期
    關(guān)鍵詞:內(nèi)膛射彈彈頭

    沈超, 周克棟, 陸野, 喬自平

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

    0 引言

    機(jī)槍連發(fā)射擊過程中,槍管受到較高頻率的火藥燃?xì)馑矐B(tài)熱作用、化學(xué)作用、機(jī)械摩擦作用(包括彈頭和火藥燃?xì)庵械墓滔囝w粒)以及膛壓作用,槍管內(nèi)膛會出現(xiàn)諸如裂紋、燒蝕坑等損傷。大量試驗顯示,即使在槍管內(nèi)壁鍍上一層甚至多層耐高溫、耐腐蝕和耐磨損的鉻層以保護(hù)槍管基體材料,鉻層仍會在射擊壽命試驗的前期就出現(xiàn)裂紋,并隨著射擊過程的進(jìn)行,裂紋不斷加深加寬,并出現(xiàn)剝落,基體在失去鉻層保護(hù)之后便會快速燒蝕磨損,導(dǎo)致槍管壽終。壽命試驗表明,大口徑機(jī)槍的壽終常表現(xiàn)為橢圓彈孔率的超標(biāo),而坡膛的結(jié)構(gòu)損傷是造成槍管產(chǎn)生橢圓彈孔的根本原因[1],彈頭擠進(jìn)過程是一個具有高瞬態(tài)、高溫強(qiáng)沖擊、高速摩擦、大變形等特點的瞬態(tài)動力學(xué)過程[2],內(nèi)膛損傷導(dǎo)致的內(nèi)膛結(jié)構(gòu)及尺寸的細(xì)微變化都會對彈頭擠進(jìn)過程的受力及運(yùn)動姿態(tài)產(chǎn)生較大影響,進(jìn)而影響整個內(nèi)彈道過程。

    近年來,國內(nèi)外多位學(xué)者對彈- 槍耦合過程及彈槍耦合退化內(nèi)彈道進(jìn)行了研究,并取得了一些成果。文獻(xiàn)[3-5]對槍管有無磨損下的內(nèi)彈道過程進(jìn)行了理論研究,但均未涉及擠進(jìn)過程的大變形塑性力學(xué)及彈頭材料的損傷失效。文獻(xiàn)[6]對火炮身管失效機(jī)理進(jìn)行了研究并提出了基于火炮身管外壁面應(yīng)變變化量預(yù)測身管壽命的方法。文獻(xiàn)[7]建立了考慮彈頭塑性大變形及斷裂失效的彈帶擠進(jìn)身管坡膛的有限元模型,研究了擠進(jìn)過程中彈帶塑性變形流動及彈帶的斷裂失效模式。文獻(xiàn)[8]對某狙擊步槍的彈- 槍耦合作用過程進(jìn)行了研究,建立了彈- 槍相互作用數(shù)值計算模型,針對彈- 槍結(jié)構(gòu)參數(shù)變化對彈頭膛內(nèi)運(yùn)動過程的影響進(jìn)行了分析。文獻(xiàn)[9]基于試驗獲得了火炮不同壽命階段內(nèi)膛的燒蝕磨損情況,并建立理論模型研究了由于內(nèi)膛燒蝕磨損而引起的火炮彈道壽命變化情況。文獻(xiàn)[10-11]建立了磨損內(nèi)膛的彈- 炮耦合模型,采用參數(shù)化方法對身管內(nèi)膛進(jìn)行建模,分析了火炮內(nèi)彈道膛壓和彈丸初速隨內(nèi)膛損傷的退化過程。

    以上對彈- 炮、彈- 槍耦合作用的研究中,身管內(nèi)膛損傷的表現(xiàn)形式主要局限在內(nèi)膛的磨損量上,不論是以磨損最大處磨損量(文獻(xiàn)[9]中稱該處為磨損特征點)代替整個內(nèi)膛的磨損量,抑或是對身管分段建模,考慮身管不同截面處磨損量的差異,均只考慮了內(nèi)膛磨損對彈- 槍耦合過程的影響,并不能反映身管真實的內(nèi)膛損傷形式(如裂紋、燒蝕坑、鉻層剝落及磨損等),且未對內(nèi)膛損傷對彈頭內(nèi)彈道過程運(yùn)動姿態(tài)、彈頭表面形貌等的影響進(jìn)行深入研究,不能解釋壽終身管在彈頭初速未出現(xiàn)明顯下降情況下精度就已不滿足指標(biāo)要求的現(xiàn)象。

    為分析某大口徑機(jī)槍真實內(nèi)膛損傷情況下,不同壽命階段槍管內(nèi)膛損傷對彈頭擠進(jìn)過程的影響,本文基于槍管壽命試驗獲得的不同壽命階段槍管內(nèi)膛的損傷數(shù)據(jù),建立了各階段損傷槍管內(nèi)膛的三維模型并對其進(jìn)行了精確的網(wǎng)格劃分。采用有限元分析軟件Abaqus建立了擠進(jìn)過程的熱力耦合有限元分析模型,考慮內(nèi)膛損傷導(dǎo)致的彈后空間變大,并通過編寫彈底壓力推力子程序(VUAMP)實現(xiàn)了內(nèi)彈道方程組和彈頭擠進(jìn)過程的耦合求解。研究并得到了槍管在4個壽命階段、不同的內(nèi)膛損傷對彈頭擠進(jìn)過程中擠進(jìn)阻力、運(yùn)動姿態(tài)及被甲表面刻痕形成的影響規(guī)律。

    1 擠進(jìn)過程有限元模型的建立

    1.1 基本假設(shè)

    1) 彈頭各部分材料為各向同性材料,比熱容、導(dǎo)熱系數(shù)等熱物理性能參數(shù)隨溫度變化而變化;

    2) 忽略槍管變形,假設(shè)槍管為剛體;彈頭材料考慮塑性變形和損傷失效,其屈服強(qiáng)度服從Mises屈服準(zhǔn)則;

    3) 不考慮槍管后坐,忽略彈頭前端空氣對彈頭運(yùn)動的影響;

    4) 不考慮彈頭對流換熱和輻射放熱。

    1.2 有限元網(wǎng)格

    1.2.1 槍管與彈頭有限元網(wǎng)格

    在基于試驗獲得的內(nèi)膛損傷數(shù)據(jù)基礎(chǔ)上建立彈- 槍三維模型如圖1所示。內(nèi)膛部分由彈膛、坡膛及線膛組成,彈頭由銅被甲、鉛套、鋼芯及頭部填充物組成,初始狀態(tài)下彈頭和坡膛間有一定的間隙,彈頭走過一段自由行程后,弧形部與坡膛接觸開始擠進(jìn)過程。

    對槍管和彈頭進(jìn)行必要的簡化后利用Hyper-Mesh軟件對槍管和彈頭進(jìn)行六面體網(wǎng)格劃分,如圖2所示。槍管及彈頭均采用C3D8RT三維八節(jié)點熱力耦合單點積分單元,由于要準(zhǔn)確地還原內(nèi)膛損傷,必須對損傷處網(wǎng)格進(jìn)行加密,導(dǎo)致局部很小的網(wǎng)格尺寸,若采用顯式算法求解動力學(xué)方程,其穩(wěn)定時間步長也會很小,因此為了節(jié)約計算成本,對含內(nèi)膛損傷的槍管進(jìn)行網(wǎng)格劃分后,將槍管設(shè)定為剛體,并固定在空間一定位置處。對彈頭圓柱部網(wǎng)格進(jìn)行局部細(xì)化,多次計算結(jié)果表明圓柱部網(wǎng)格尺寸為0.1 mm時,可以在不降低計算精度的前提下縮短計算時間,最終彈頭總網(wǎng)格數(shù)為463 175個。

    1.2.2 材料模型

    被甲在擠進(jìn)線膛過程中出現(xiàn)高應(yīng)變率、材料溫度升高、熱軟化、應(yīng)變強(qiáng)化等一系列復(fù)雜的非線性力學(xué)形態(tài),故本文對有較大塑性變形的被甲和鉛套材料采用能描述材料熱黏性行為的Johnson-Cook本構(gòu)模型[12]:

    (1)

    材料有塑性變形時,以等效塑性臨界應(yīng)變εf作為損傷的判活依據(jù):

    (2)

    式中:D1~D5為材料斷裂失效參數(shù);σ*為應(yīng)力三軸度。

    采用線性損傷演化規(guī)律描述材料剛度下降,材料累積損傷參數(shù)為

    (3)

    式中:D為材料內(nèi)單元的損傷值,當(dāng)單元損傷值D=1時,單元失效失去承載能力并從仿真中刪除;Δε為等效塑性應(yīng)變變化量。

    1.2.3 接觸算法及邊界條件

    彈頭被甲與槍管內(nèi)膛間的接觸設(shè)定為基于罰函數(shù)法的通用接觸,基于文獻(xiàn)[13]中所述,對于高速、高壓擠進(jìn)過程,動摩擦系數(shù)應(yīng)不大于0.02,本文取動摩擦系數(shù)為0.01. 彈頭底部加以壓力載荷pd,載荷值來自于內(nèi)彈道方程的計算,并由推力子程序VUAMP施加于彈尾。邊界條件為槍管和彈殼的全自由度約束。

    1.3 熱力耦合模型

    文獻(xiàn)[11,14]研究發(fā)現(xiàn),由于溫度升高引起的彈頭被甲材料變軟導(dǎo)致的彈頭內(nèi)彈道性能(擠進(jìn)阻力、膛壓、彈頭初速等)的變化不可忽略。因此本文建立的有限元分析模型為考慮彈頭被甲材料塑性變形熱Qε及被甲與槍管內(nèi)膛間摩擦生熱Qf導(dǎo)致的被甲材料溫度改變的熱力耦合模型,二者引起的溫度改變量分別為Tε、Tf. 由摩擦熱和塑性變形熱引起的溫度改變?nèi)缦拢?/p>

    (4)

    式中:ρj和cj分別為被甲材料的密度和比熱容;σj和εj分別為被甲的應(yīng)力和塑性應(yīng)變;β為塑性熱轉(zhuǎn)化系數(shù),一般為0.85~0.95,本文取β=0.9.

    熱量由彈頭表面向內(nèi)傳導(dǎo),采用三維非穩(wěn)態(tài)、變物性熱傳導(dǎo)微分方程:

    (5)

    式中:T為彈帶材料溫度;kX、kY、kZ分別為3個方向的熱傳導(dǎo)系數(shù);Q為內(nèi)熱源,一般指塑性熱。

    2 含內(nèi)膛損傷槍管有限元模型的建模方法

    槍管內(nèi)膛損傷的主要形式有:裂紋、燒蝕坑、鍍層剝落及內(nèi)膛表面磨損。其中:內(nèi)膛表面磨損會使槍管對彈頭的導(dǎo)向及導(dǎo)轉(zhuǎn)作用減小,使得彈頭在內(nèi)彈道過程中擾動增大,導(dǎo)轉(zhuǎn)力減?。涣鸭y、燒蝕坑等一方面會對彈頭擠進(jìn)及沿膛運(yùn)動的穩(wěn)定性產(chǎn)生不利影響,另一方面也會影響彈頭嵌入膛線過程中被甲表面刻痕的形貌,嚴(yán)重時甚至?xí)潅患撞牧?,破壞其表面完整性,?dǎo)致彈頭氣動性能發(fā)生改變。在上述因素的綜合影響下,隨著射彈量的增加彈頭出膛口時的初始擾動及氣動性能均與從新槍管中射出的彈頭狀態(tài)有較大差異,最終導(dǎo)致了槍管壽終。

    針對某大口徑機(jī)槍進(jìn)行的大量壽命試驗表明:槍管內(nèi)膛損傷的形式及分布規(guī)律隨射彈數(shù)的增加表現(xiàn)出極大的規(guī)律性[1],這為使用數(shù)值方法模擬不同壽命階段槍管的擠進(jìn)乃至整個內(nèi)彈道過程提供了依據(jù)。

    通過對槍管內(nèi)膛直徑進(jìn)行測量,獲得了4支不同射彈量的槍管陽線直徑沿槍管軸向分布情況如圖3所示。射擊試驗中通過內(nèi)窺觀測得到的槍管坡膛處損傷主要有裂紋、燒蝕坑及鉻層剝落,4支槍管坡膛處損傷形式相同,區(qū)別在于損傷的嚴(yán)重程度及分布位置有所不同。以射彈量800發(fā)的槍管為例,其坡膛處損傷情況內(nèi)窺圖如圖4所示,對圖4中3種主要的內(nèi)膛損傷形式分別劃分網(wǎng)格如圖5所示(圖5同時也是射彈數(shù)800發(fā)槍管的網(wǎng)格劃分圖)。

    以射彈數(shù)800發(fā)(壽命前期)的槍管為例,分析其有限元模型的建模方法:首先根據(jù)圖3中試驗測得的槍管內(nèi)膛直徑數(shù)據(jù)建立其三維模型并劃分網(wǎng)格;然后根據(jù)由試驗獲得的坡膛部位內(nèi)窺照片及視頻所得到的坡膛段損傷分布規(guī)律及嚴(yán)重程度在內(nèi)膛預(yù)置不同數(shù)量及分布位置的損傷,即可建立該壽命階段含內(nèi)膛損傷槍管的有限元模型。由于不同射彈量的槍管坡膛處損傷形式相同,因此同理可以由圖3中其他壽命階段槍管內(nèi)膛陽線直徑沿軸向分布情況,結(jié)合已由試驗獲得對應(yīng)射彈量下的內(nèi)膛形貌特征,建立其他各壽命階段的槍管有限元模型。

    3 內(nèi)彈道過程耦合計算

    本文采用VUAMP子程序的方法,將內(nèi)彈道方程的求解和彈頭運(yùn)動過程耦合,能夠更好地獲得內(nèi)膛有損傷時彈后的膛壓數(shù)據(jù),提高計算結(jié)果的準(zhǔn)確度。實現(xiàn)過程如下:使用Fortran語言編寫內(nèi)彈道方程的VUAMP子程序,子程序根據(jù)t時刻的火藥已燃相對厚度zt等內(nèi)彈道參數(shù)計算得到該時刻彈底壓力pd,t,并將彈底壓力傳入有限元軟件中,彈頭在pd,t作用下沿膛線向前運(yùn)動,Δt時間后有限元軟件將t+Δt時刻的彈頭運(yùn)動相關(guān)參數(shù)傳回VUAMP子程序,由子程序計算得到此時的彈底壓力pd,t+Δt并再次傳入有限元軟件,對彈頭運(yùn)動參數(shù)進(jìn)行求解。如此往復(fù)循環(huán)直至彈頭出膛。

    3.1 基本假設(shè)

    內(nèi)彈道方程組滿足以下基本假設(shè):

    1)火藥燃燒服從幾何燃燒定律;

    2)火藥燃燒在平均壓力條件下進(jìn)行;

    3)火藥燃燒速度與壓力呈指數(shù)關(guān)系;

    4)火藥燃燒期間和燃燒結(jié)束后,火藥燃燒生成物始終保持不變,即把火藥力f和火藥氣體余容α當(dāng)作常量處理。

    3.2 內(nèi)膛磨損導(dǎo)致彈后空間增加的計算

    槍管內(nèi)膛有損傷后,彈后空間會比初始槍管有所增加,并隨著彈頭沿槍管軸線向前運(yùn)動,彈后空間的增加量也在不斷增加。在彈頭運(yùn)動dt時間內(nèi),彈頭向前運(yùn)動ds距離,t時刻時彈尾位置處槍管內(nèi)膛截面積增加量為ΔS,故彈后空間增加量可以表示為

    dV=ΔSds,

    積分得到:

    (6)

    式中:ΔV為彈后空間增加量;N、K、ΔR分別為膛線條數(shù)、陽線角度系數(shù)、陽線半徑增量[11]。

    由于陰線磨損量較小且不易測量,以及由裂紋、燒蝕坑等引起的彈后空間變化很小,本文僅考慮由于陽線磨損導(dǎo)致彈后空間的增大。使用復(fù)合辛普森公式求解該積分項以便于編程計算,并將結(jié)果加入內(nèi)彈道方程的彈后空間項中。

    3.3 內(nèi)彈道方程組

    對經(jīng)典內(nèi)彈道方程組進(jìn)行改寫,考慮擠進(jìn)過程及由于內(nèi)膛損傷導(dǎo)致彈后藥室容積的增加,并將次要功系數(shù)具體為各項能量的形式,同時為便于使用FORTRAN編程,將內(nèi)彈道方程寫為1個主要方程和4個輔助方程的形式,采用4階龍格- 庫塔法求解(7)式中的主要方程得到火藥已燃相對厚度z,然后由輔助方程(8)式計算得到膛內(nèi)平均壓力p,根據(jù)文獻(xiàn)[15]中彈底壓力與膛內(nèi)平均壓力的轉(zhuǎn)換公式,計算得到彈底壓力pd.

    (7)

    式中:u1、e1、n1分別為火藥燃燒速度系數(shù)、火藥弧厚及燃速指數(shù);ω為裝藥量;mb、S分別為彈質(zhì)量和無損傷線膛的截面積;ψ為火藥已燃百分比;θ為熱力指數(shù);v為彈頭速度;s、Lψ分別為彈頭軸向行程、藥室容積縮徑長;Eq為擠進(jìn)過程中的總次要功。

    輔助方程:

    (8)

    式中:χ、λ、μ為火藥形狀特征量;Er、Ef、Eε、Em分別為擠進(jìn)過程中的彈頭旋轉(zhuǎn)動能、摩擦耗散能、彈塑性應(yīng)變能及火藥氣體運(yùn)動功。

    綜上所述,擠進(jìn)過程內(nèi)彈道方程和有限元程序的耦合計算流程設(shè)計如圖6所示。

    4 有限元模型驗證及結(jié)果分析

    本文對某大口徑機(jī)槍射彈數(shù)分別為0發(fā)(無損傷)、1 400發(fā)、3 000發(fā)、6 000發(fā)(壽終)的槍管建立了上述彈- 槍耦合有限元模型,對彈頭擠進(jìn)含內(nèi)膛損傷槍管的過程進(jìn)行了分析。為便于敘述,將射彈量0發(fā)、1 400發(fā)、3 000發(fā)、6 000發(fā)的槍管分別編號為1號、2號、3號、4號槍管。

    4.1 模型驗證

    為了驗證所構(gòu)建有限元模型的準(zhǔn)確性,本文將4支槍管彈- 槍耦合模型的內(nèi)彈道計算結(jié)果和試驗結(jié)果進(jìn)行了對比。對于無損傷槍管,有限元模型計算得到的內(nèi)彈道過程最高膛壓為328 MPa、彈頭出膛口時速度為806 m/s、擠進(jìn)阻力最大值為10 954 N,與試驗測得的最高膛壓平均值320 MPa相比升高了2.5%、與平均實測初速800 m/s相比升高了0.75%,與靜態(tài)擠進(jìn)試驗測得最大擠進(jìn)阻力12 000 N相比降低了8.7%,其中擠進(jìn)阻力降低較多是因為高速高壓動態(tài)擠進(jìn)過程的摩擦系數(shù)比靜態(tài)擠進(jìn)試驗低[13]。此外,還將理論計算得到的4支槍管彈頭初速與試驗值進(jìn)行了對比,如表1所示。從表1可以看出,各槍管發(fā)射彈頭初速的計算值與試驗值誤差較小。綜上所述,本文建立的彈- 槍耦合有限元計算模型是合理與準(zhǔn)確的。

    表1 彈頭初速對比

    4.2 內(nèi)膛損傷對彈頭被甲表面形貌的影響

    彈頭表面形貌和運(yùn)動姿態(tài)綜合影響了彈頭外彈道性能,決定了武器射擊精度,對擠進(jìn)過程,尤其是槍管內(nèi)膛不同程度損傷形貌下的擠進(jìn)過程進(jìn)行分析,可以從根源上分析武器精度不達(dá)標(biāo)而壽終的機(jī)理。

    由裝配條件及彈- 槍結(jié)構(gòu)可知,整個擠進(jìn)過程對應(yīng)的彈頭軸向行程s為5~50 mm,圖7~圖9所示分別為彈頭被甲擠進(jìn)4支槍管過程中對應(yīng)彈頭軸向行程分別為10 mm、30 mm、50 mm時彈頭圓柱部的von Mises應(yīng)力云圖(右側(cè)為彈頭頭部),也同時可以間接反映被甲的變形狀態(tài)。圖7對應(yīng)彈頭擠進(jìn)的初期(被甲剛開始嵌入膛線),從中可以看出隨著槍管內(nèi)膛磨損量的增加,相同軸向位移時彈頭表面刻痕不斷變淺,即彈頭擠進(jìn)起始位置向槍管口部方向移動,導(dǎo)致了擠進(jìn)前彈頭自由運(yùn)動行程(彈頭初始位置到初始擠進(jìn)位置之間的軸向距離)增加。自由行程的增加會使彈頭在脫離彈殼后缺少軸向定位,彈頭運(yùn)動的不確定性增大,彈頭擾動增大,對彈頭擠進(jìn)過程中的運(yùn)動姿態(tài)及表面形貌的規(guī)則性及完整性均會產(chǎn)生不利影響。圖8、圖9分別為彈頭剛?cè)壳度肫绿盘啪€以及彈頭完全嵌入全深線膛的時刻,從圖8(a)、圖9(a)可以看出,彈頭擠進(jìn)無損傷槍管時,其表面刻痕的對稱性及表面的完整性保持得較好,表現(xiàn)為除規(guī)則由陽線形成的刻槽外,其他位置無不規(guī)則刻痕形成,且導(dǎo)轉(zhuǎn)側(cè)與非導(dǎo)轉(zhuǎn)側(cè)處刻痕寬度、深度也基本相同,以及彈頭所受應(yīng)力值及應(yīng)力的軸向梯度也較小。隨著內(nèi)膛損傷的不斷發(fā)展,陽線高度及寬度隨著磨損及鉻層的剝落不斷減小,使得被甲表面刻槽的寬度、深度及長度都在不斷減小,如表2所示,這會使得彈頭所受導(dǎo)轉(zhuǎn)力及轉(zhuǎn)速隨損傷的發(fā)展而不斷降低,不利于彈頭飛行穩(wěn)定。

    類別無損傷槍管射彈數(shù)1400發(fā)射彈數(shù)3000發(fā)射彈數(shù)6000發(fā)刻槽寬度/mm2.031.991.921.71刻槽深度/mm0.1690.1570.1150.052刻槽長度/mm28.0327.5727.0526.61

    此外,由于陽線導(dǎo)轉(zhuǎn)側(cè)損傷程度較非導(dǎo)轉(zhuǎn)側(cè)更為嚴(yán)重,使得彈頭在這兩個位置處刻槽的形狀隨射彈數(shù)增加表現(xiàn)出明顯的非對稱性,如圖9(b)、圖9(c)、圖9(d)所示。結(jié)合圖9(b)、圖9(c)可以看出:陽線表面不規(guī)則的鉻層剝落及燒蝕坑會使得被甲刻槽的中間出現(xiàn)不規(guī)則的刻痕,并隨著損傷的擴(kuò)展而擴(kuò)展;而在槍管壽命末期,由于陽線表面鉻層已經(jīng)基本全部剝落,陽線表面被磨損得比較規(guī)則,被甲的刻槽中間反而不會出現(xiàn)明顯的不規(guī)則刻痕,如圖9(d)所示。

    4.3 內(nèi)膛損傷對彈頭擠進(jìn)過程運(yùn)動姿態(tài)的影響

    除彈頭表面形貌外,內(nèi)膛損傷也會對彈頭擠進(jìn)過程的運(yùn)動姿態(tài)產(chǎn)生較大影響,如彈頭質(zhì)心在垂直槍管軸線平面內(nèi)的擺動位移、彈頭軸線與槍管軸線之間的夾角(擺動角)等。彈頭擠進(jìn)4根槍管過程質(zhì)心在垂直槍管軸線平面內(nèi)的豎直和水平偏移量如圖10所示。

    從圖10可以看出:擠進(jìn)初始階段,對應(yīng)圖10中Ⅰ區(qū),彈頭運(yùn)動完自由行程剛開始擠進(jìn),此階段由于彈頭運(yùn)動缺少約束且沖擊較大,彈頭擠進(jìn)4支槍管時質(zhì)心偏移量均較大,隨著彈頭在槍管約束下逐漸歸正,質(zhì)心偏移量均在逐漸減小,槍管內(nèi)膛損傷越嚴(yán)重,彈頭完全歸正所需的彈頭軸向行程越長;Ⅱ區(qū)為平穩(wěn)擠進(jìn)階段,此階段彈頭已大部分?jǐn)D進(jìn)膛線,彈頭運(yùn)動較平穩(wěn),對比4條曲線的質(zhì)心偏移量可以明顯發(fā)現(xiàn),隨著槍管內(nèi)膛損傷的擴(kuò)展,彈頭質(zhì)心在垂直槍管軸線平面內(nèi)的偏移量呈明顯的增加趨勢,彈頭完全擠進(jìn)膛線時,2號~4號槍管對應(yīng)的質(zhì)心總偏移量,即圖10中擠進(jìn)終了時刻(對應(yīng)曲線終點位置)彈頭質(zhì)心在豎直方向和水平方向偏移量的矢量和,分別為1號槍管對應(yīng)總偏移量的2.90倍、3.01倍、3.61倍,內(nèi)膛損傷對彈頭質(zhì)心運(yùn)動姿態(tài)的影響十分嚴(yán)重。

    槍管內(nèi)膛損傷也會對彈頭軸線與槍管軸線的夾角(即擺動角)產(chǎn)生較大影響,而彈頭出膛口后的擺動角和彈頭表面形貌綜合影響了彈頭外彈道過程的氣動參數(shù),是影響彈頭飛行穩(wěn)定性和自動武器精度的決定性因素。本文在彈頭質(zhì)心建立局部坐標(biāo)系,并在質(zhì)心局部坐標(biāo)系和地面坐標(biāo)系之間建立無自由度約束的連接器,獲得了彈頭擠進(jìn)4支槍管過程中彈頭軸線與槍管軸線的夾角隨彈頭擠進(jìn)過程的變化情況,并將其分解到過槍管軸線的豎直平面和水平平面內(nèi),分別對應(yīng)為鉛垂擺動角與側(cè)向擺動角,如圖11所示。從圖11中可以看出:在擠進(jìn)的初始階段,彈頭由于剛嵌入膛線,存在著動態(tài)沖擊并且缺少內(nèi)膛的有效約束,兩個方向的擺動角均較大;隨著彈頭不斷向前擠進(jìn),擺動角呈下降趨勢,且內(nèi)膛損傷越嚴(yán)重的槍管,彈頭擺動角下降得越慢。在擠進(jìn)的中后期,彈頭運(yùn)動趨于平穩(wěn),擺動角均較小,且隨著槍管射彈數(shù)的增加,兩個方向擺動角均呈增加的趨勢。與圖10不同的是,圖11中Ⅲ區(qū)對應(yīng)的擠進(jìn)末期內(nèi),4支槍管的擺動角均較之前開始增加,這是因為在擠進(jìn)末期,彈頭已經(jīng)絕大部分完成了擠進(jìn)并開始了沿線膛運(yùn)動過程,彈頭所受軸向阻力迅速下降,彈頭軸向速度及轉(zhuǎn)速開始快速上升,使得彈頭軸向運(yùn)動的擺動量增大,這與文獻(xiàn)[8]中對內(nèi)彈道時期彈頭擺動角變化規(guī)律的研究結(jié)果是一致的。

    4.4 內(nèi)膛損傷對彈頭擠進(jìn)阻力的影響

    (9)

    圖12為彈頭擠進(jìn)4支槍管過程的擠進(jìn)阻力隨彈頭軸向位移的曲線,彈頭擠進(jìn)每支槍管的最大擠進(jìn)阻力及此時對應(yīng)的彈頭軸向位移列于表3中。結(jié)合圖12及表3可以發(fā)現(xiàn):擠進(jìn)阻力的峰值出現(xiàn)在彈頭完全擠進(jìn)全深膛線之前,這是因為隨著彈頭被甲不斷擠進(jìn)膛線,被甲塑性變形量不斷增大,擠進(jìn)阻力迅速上升,在塑性變形量最大時擠進(jìn)阻力達(dá)到峰值,此后隨著塑性變形量的減小而降低,其變化規(guī)律與文獻(xiàn)[7]關(guān)于擠進(jìn)阻力的變化規(guī)律研究結(jié)果相符;彈頭擠進(jìn)過程擠進(jìn)阻力隨槍管射彈數(shù)增加而不斷下降,且擠進(jìn)阻力達(dá)到峰值時彈頭軸向位移在不斷增加,這是因為隨著射彈數(shù)的增加,內(nèi)膛磨損不斷加重及坡膛處鉻層逐漸剝落均導(dǎo)致了坡膛陰陽線直徑的擴(kuò)大,使得被甲材料完全嵌入膛線時的塑性變形量降低,達(dá)到相同塑性變形量時的彈頭軸向位移增加,這與文獻(xiàn)[16]中不同坡膛角下的擠進(jìn)阻力變化規(guī)律相似。

    5 結(jié)論

    本文針對大口徑槍械槍管壽命較短這一嚴(yán)重問題,在試驗基礎(chǔ)上建立了彈頭擠進(jìn)不同射彈數(shù)槍管的數(shù)值分析模型,研究了內(nèi)膛損傷對彈頭表面形貌及運(yùn)動姿態(tài)等的影響規(guī)律,得出以下結(jié)論:

    1)內(nèi)膛損傷會顯著影響彈頭表面形貌及運(yùn)動姿態(tài)等與外彈道性能緊密相關(guān)的參數(shù),彈頭內(nèi)彈道過程擾動增加,彈頭所受導(dǎo)轉(zhuǎn)力和轉(zhuǎn)速下降以及氣動力參數(shù)發(fā)生改變等是導(dǎo)致彈頭外彈道飛行穩(wěn)定性降低,橫彈孔率上升,進(jìn)而導(dǎo)致槍管壽終的主要原因。

    2)彈頭擠進(jìn)無損傷槍管時,其表面刻痕的對稱性及表面的完整性保持得較好,表現(xiàn)為除規(guī)則由陽線形成的刻槽外,其他位置無不規(guī)則刻痕形成,且導(dǎo)轉(zhuǎn)側(cè)與非導(dǎo)轉(zhuǎn)側(cè)處刻痕寬度、深度也基本相同,彈頭所受應(yīng)力值及應(yīng)力軸向梯度也較小。

    3)隨著槍管射彈數(shù)增加,彈頭擠進(jìn)后表面形貌的改變明顯,表現(xiàn)為被甲表面刻槽的寬度、深度及長度在不斷減?。粚?dǎo)轉(zhuǎn)側(cè)與非導(dǎo)轉(zhuǎn)側(cè)刻槽的形狀呈明顯的非對稱性;被甲表面除陽線形成的刻槽外還出現(xiàn)了不規(guī)則分布的刻痕。

    4)彈頭擠進(jìn)過程質(zhì)心的非軸向位移及彈頭擺動角隨著擠進(jìn)過程彈頭的歸正而逐漸減??;隨著射彈數(shù)增加,擠進(jìn)過程彈頭歸正所需要的彈頭軸向行程以及歸正后的非軸向位移量、擺動角等都在增加。

    5)隨著槍管射彈數(shù)增加,擠進(jìn)過程擠進(jìn)阻力不斷下降,且擠進(jìn)阻力達(dá)到峰值時彈頭軸向位移在不斷增加。

    猜你喜歡
    內(nèi)膛射彈彈頭
    果樹整形修剪技術(shù)
    熱身管下小口徑槍彈彈頭殼材料對其膛內(nèi)運(yùn)動的影響規(guī)律
    克服獼猴桃架面內(nèi)膛光禿的技術(shù)措施
    煙臺果樹(2021年2期)2021-07-21 07:18:48
    高速射彈并聯(lián)入水過程空泡演化特性試驗
    并列超空泡射彈彈道特性研究
    水下高速超空泡射彈串行運(yùn)動流體動力特性研究
    導(dǎo)彈突防概率解析與仿真計算結(jié)果差異分析
    柑橘內(nèi)膛腹接增產(chǎn)技術(shù)
    怎樣防止桃樹內(nèi)膛過早光禿
    質(zhì)心位置對超空泡射彈穩(wěn)定性影響的規(guī)律研究*
    色精品久久人妻99蜜桃| 97人妻天天添夜夜摸| 中文字幕人妻熟女乱码| 亚洲性夜色夜夜综合| 欧美在线黄色| 中文字幕制服av| 黑人巨大精品欧美一区二区mp4| 人妻丰满熟妇av一区二区三区 | 69av精品久久久久久| 免费看a级黄色片| 少妇 在线观看| 99久久99久久久精品蜜桃| 国产人伦9x9x在线观看| 热99re8久久精品国产| 国产成人精品久久二区二区91| 久久香蕉国产精品| tube8黄色片| 一级黄色大片毛片| 久久人妻av系列| 国产亚洲精品久久久久5区| 波多野结衣av一区二区av| 亚洲五月色婷婷综合| 18禁黄网站禁片午夜丰满| 亚洲性夜色夜夜综合| 黄色 视频免费看| 性少妇av在线| 国产三级黄色录像| 亚洲成人免费av在线播放| 可以免费在线观看a视频的电影网站| 99riav亚洲国产免费| 欧美日韩乱码在线| 久久精品国产清高在天天线| 99久久人妻综合| 一区二区三区激情视频| 欧美日韩瑟瑟在线播放| 啦啦啦在线免费观看视频4| av天堂久久9| 首页视频小说图片口味搜索| 如日韩欧美国产精品一区二区三区| 久久精品成人免费网站| 亚洲七黄色美女视频| 久久久精品免费免费高清| 久久久久久人人人人人| 午夜两性在线视频| 国产精品久久久av美女十八| 亚洲第一青青草原| 国产av精品麻豆| 亚洲精品国产精品久久久不卡| 天堂√8在线中文| 欧美日韩亚洲综合一区二区三区_| 亚洲人成77777在线视频| 天天操日日干夜夜撸| 一区在线观看完整版| 国产成人一区二区三区免费视频网站| 久久这里只有精品19| 久久九九热精品免费| 一边摸一边抽搐一进一小说 | 精品久久久久久久毛片微露脸| 日韩大码丰满熟妇| 国产精品亚洲av一区麻豆| 我的亚洲天堂| 国产成人免费无遮挡视频| 夜夜躁狠狠躁天天躁| 中文字幕av电影在线播放| aaaaa片日本免费| 制服诱惑二区| 狂野欧美激情性xxxx| 高清在线国产一区| 啦啦啦视频在线资源免费观看| 91精品国产国语对白视频| av国产精品久久久久影院| xxxhd国产人妻xxx| 国产精品99久久99久久久不卡| 国产精品久久久久成人av| 亚洲欧美激情在线| 国产精品免费一区二区三区在线 | 精品无人区乱码1区二区| 免费观看精品视频网站| 欧美乱色亚洲激情| 12—13女人毛片做爰片一| 国产一区在线观看成人免费| 女人久久www免费人成看片| 如日韩欧美国产精品一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 免费观看精品视频网站| 亚洲国产欧美网| 无限看片的www在线观看| ponron亚洲| 在线免费观看的www视频| 精品国产美女av久久久久小说| 午夜福利乱码中文字幕| 国产不卡av网站在线观看| 久久精品国产亚洲av香蕉五月 | 精品久久久久久久毛片微露脸| 黄片小视频在线播放| 国产又色又爽无遮挡免费看| 国产又爽黄色视频| 黄色女人牲交| 精品人妻在线不人妻| 欧美激情 高清一区二区三区| 国产精品二区激情视频| 国产一区二区三区在线臀色熟女 | 日本a在线网址| 如日韩欧美国产精品一区二区三区| tocl精华| 777米奇影视久久| www.熟女人妻精品国产| 淫妇啪啪啪对白视频| 中文字幕色久视频| 婷婷丁香在线五月| 午夜免费鲁丝| 亚洲色图综合在线观看| 99国产精品一区二区三区| 老汉色∧v一级毛片| 精品一品国产午夜福利视频| 亚洲精品美女久久久久99蜜臀| 欧美日韩亚洲综合一区二区三区_| 麻豆国产av国片精品| 校园春色视频在线观看| 国产99久久九九免费精品| 欧美人与性动交α欧美软件| 免费看a级黄色片| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲精品国产精品久久久不卡| 欧美日韩成人在线一区二区| 久久久久久久久久久久大奶| 超碰97精品在线观看| 每晚都被弄得嗷嗷叫到高潮| 国产亚洲精品第一综合不卡| 国产精品.久久久| 无人区码免费观看不卡| 欧美最黄视频在线播放免费 | 老司机福利观看| 久久久久久久久免费视频了| 一级黄色大片毛片| 免费在线观看完整版高清| 少妇裸体淫交视频免费看高清 | 一本大道久久a久久精品| 欧美不卡视频在线免费观看 | 无人区码免费观看不卡| 欧美丝袜亚洲另类 | 操美女的视频在线观看| 亚洲avbb在线观看| 老司机深夜福利视频在线观看| 自线自在国产av| 久热这里只有精品99| 亚洲成人国产一区在线观看| 久久99一区二区三区| 国产单亲对白刺激| 怎么达到女性高潮| 国产亚洲一区二区精品| 每晚都被弄得嗷嗷叫到高潮| 国产一区二区三区综合在线观看| 交换朋友夫妻互换小说| 老司机福利观看| 在线十欧美十亚洲十日本专区| 国产男靠女视频免费网站| 免费黄频网站在线观看国产| 露出奶头的视频| 成人黄色视频免费在线看| 视频区图区小说| 两个人看的免费小视频| 黑人猛操日本美女一级片| 欧美乱码精品一区二区三区| 亚洲国产欧美日韩在线播放| 高清视频免费观看一区二区| 香蕉国产在线看| 黑人欧美特级aaaaaa片| 亚洲va日本ⅴa欧美va伊人久久| 日韩成人在线观看一区二区三区| 久久人妻熟女aⅴ| 首页视频小说图片口味搜索| 精品国产国语对白av| 人人妻,人人澡人人爽秒播| 国产免费现黄频在线看| 国产欧美日韩精品亚洲av| 国内毛片毛片毛片毛片毛片| 久久精品国产亚洲av高清一级| 在线观看午夜福利视频| 别揉我奶头~嗯~啊~动态视频| 欧美日韩一级在线毛片| 亚洲人成伊人成综合网2020| 极品人妻少妇av视频| 婷婷丁香在线五月| 人妻丰满熟妇av一区二区三区 | av天堂久久9| 无限看片的www在线观看| 伊人久久大香线蕉亚洲五| 欧美黑人精品巨大| 老鸭窝网址在线观看| 女警被强在线播放| 免费久久久久久久精品成人欧美视频| 自拍欧美九色日韩亚洲蝌蚪91| 一本大道久久a久久精品| 亚洲成人手机| 精品久久久精品久久久| 亚洲人成电影免费在线| 久久精品成人免费网站| 少妇的丰满在线观看| 黄色丝袜av网址大全| 脱女人内裤的视频| 亚洲va日本ⅴa欧美va伊人久久| 黄网站色视频无遮挡免费观看| 亚洲国产欧美日韩在线播放| 欧美av亚洲av综合av国产av| 成人亚洲精品一区在线观看| cao死你这个sao货| 女警被强在线播放| 日本撒尿小便嘘嘘汇集6| 亚洲欧美激情综合另类| www.精华液| 国产精品久久电影中文字幕 | 91大片在线观看| 真人做人爱边吃奶动态| 91九色精品人成在线观看| av不卡在线播放| 国产精品98久久久久久宅男小说| 超碰97精品在线观看| 亚洲 国产 在线| 老熟妇仑乱视频hdxx| 国产精品香港三级国产av潘金莲| 免费观看人在逋| 国产区一区二久久| 啦啦啦视频在线资源免费观看| 三上悠亚av全集在线观看| 美女高潮喷水抽搐中文字幕| 亚洲av日韩在线播放| 狂野欧美激情性xxxx| 在线十欧美十亚洲十日本专区| 999久久久国产精品视频| 国产精品.久久久| 超色免费av| 亚洲七黄色美女视频| 视频区图区小说| 免费av中文字幕在线| 777久久人妻少妇嫩草av网站| 成年动漫av网址| 老鸭窝网址在线观看| 午夜精品久久久久久毛片777| 亚洲精品乱久久久久久| 亚洲av片天天在线观看| 国产精品乱码一区二三区的特点 | 国产欧美亚洲国产| 国产黄色免费在线视频| 人人妻人人澡人人看| 校园春色视频在线观看| 精品久久久久久久久久免费视频 | 国产精品秋霞免费鲁丝片| 法律面前人人平等表现在哪些方面| 亚洲一区二区三区不卡视频| 国产在线观看jvid| 黄频高清免费视频| 99久久精品国产亚洲精品| 精品高清国产在线一区| 伊人久久大香线蕉亚洲五| 午夜福利免费观看在线| 亚洲成a人片在线一区二区| 无限看片的www在线观看| av欧美777| 亚洲一区二区三区欧美精品| 欧美大码av| 国产成人欧美| 90打野战视频偷拍视频| 狂野欧美激情性xxxx| 成年版毛片免费区| 亚洲第一青青草原| 99香蕉大伊视频| 欧美日韩国产mv在线观看视频| 日本欧美视频一区| 久久亚洲真实| 人人妻人人添人人爽欧美一区卜| 好男人电影高清在线观看| 天堂中文最新版在线下载| 欧美中文综合在线视频| 亚洲三区欧美一区| 久久香蕉国产精品| 亚洲成a人片在线一区二区| 成年人午夜在线观看视频| 黄色女人牲交| 后天国语完整版免费观看| 91精品国产国语对白视频| 免费看a级黄色片| 啪啪无遮挡十八禁网站| 国产aⅴ精品一区二区三区波| 午夜视频精品福利| 亚洲av电影在线进入| 久久久国产一区二区| 久久国产精品影院| 欧美激情 高清一区二区三区| 久久久国产成人免费| 啦啦啦视频在线资源免费观看| 日韩制服丝袜自拍偷拍| 久久人妻熟女aⅴ| 多毛熟女@视频| 水蜜桃什么品种好| av欧美777| 中国美女看黄片| 亚洲成人手机| 国产精品香港三级国产av潘金莲| 亚洲七黄色美女视频| 制服人妻中文乱码| 午夜视频精品福利| 香蕉丝袜av| 在线播放国产精品三级| 精品国产超薄肉色丝袜足j| 亚洲人成77777在线视频| 欧洲精品卡2卡3卡4卡5卡区| 9热在线视频观看99| 欧美黑人欧美精品刺激| 色婷婷av一区二区三区视频| 亚洲欧美日韩另类电影网站| 亚洲性夜色夜夜综合| 国产高清国产精品国产三级| 亚洲 国产 在线| 91av网站免费观看| 亚洲精品国产一区二区精华液| 亚洲熟女毛片儿| 国产精品一区二区在线观看99| 国产精品免费大片| 午夜视频精品福利| 欧美在线一区亚洲| 在线播放国产精品三级| 免费在线观看影片大全网站| 国产精品久久视频播放| 国产成人精品久久二区二区免费| 精品人妻1区二区| 嫁个100分男人电影在线观看| 亚洲自偷自拍图片 自拍| 男女午夜视频在线观看| 91av网站免费观看| 国产精品亚洲av一区麻豆| 777久久人妻少妇嫩草av网站| 亚洲久久久国产精品| 午夜福利在线观看吧| 亚洲色图综合在线观看| 动漫黄色视频在线观看| 久久精品国产清高在天天线| 国产成人av教育| 一二三四在线观看免费中文在| 成人亚洲精品一区在线观看| 男女高潮啪啪啪动态图| 18禁观看日本| 久久久久久久午夜电影 | 日韩一卡2卡3卡4卡2021年| 色老头精品视频在线观看| 日本精品一区二区三区蜜桃| 欧美乱妇无乱码| 99精品在免费线老司机午夜| 精品亚洲成a人片在线观看| 人成视频在线观看免费观看| 最近最新中文字幕大全电影3 | 久久久久国产精品人妻aⅴ院 | 日韩欧美在线二视频 | 无人区码免费观看不卡| 日韩免费高清中文字幕av| 人妻一区二区av| 亚洲免费av在线视频| 亚洲色图av天堂| 日本黄色日本黄色录像| 欧洲精品卡2卡3卡4卡5卡区| 嫩草影视91久久| 中亚洲国语对白在线视频| 在线观看免费高清a一片| 精品国产一区二区三区四区第35| 夜夜爽天天搞| 国产成人精品久久二区二区免费| 黄色女人牲交| 久久久国产成人精品二区 | 亚洲av成人不卡在线观看播放网| 午夜亚洲福利在线播放| 国产激情欧美一区二区| 天天影视国产精品| 1024视频免费在线观看| 一级毛片精品| 免费黄频网站在线观看国产| 在线观看免费日韩欧美大片| 两性午夜刺激爽爽歪歪视频在线观看 | 国产高清国产精品国产三级| 国产亚洲精品第一综合不卡| 久久精品亚洲av国产电影网| 欧美乱色亚洲激情| 亚洲精品美女久久久久99蜜臀| 成年女人毛片免费观看观看9 | 成人亚洲精品一区在线观看| 欧美大码av| 欧美精品高潮呻吟av久久| 国产精华一区二区三区| 高清欧美精品videossex| 18禁国产床啪视频网站| 日本wwww免费看| 亚洲一区中文字幕在线| 欧美黑人欧美精品刺激| 最新在线观看一区二区三区| 下体分泌物呈黄色| 女人久久www免费人成看片| 大码成人一级视频| 黄网站色视频无遮挡免费观看| 国产亚洲精品久久久久5区| 日日爽夜夜爽网站| 亚洲视频免费观看视频| 校园春色视频在线观看| 99国产极品粉嫩在线观看| 久久99一区二区三区| 一级毛片精品| 热99久久久久精品小说推荐| 久久久久久亚洲精品国产蜜桃av| 亚洲欧美激情在线| av网站免费在线观看视频| 咕卡用的链子| 国产欧美亚洲国产| 久久热在线av| 丁香欧美五月| 国产麻豆69| 欧美亚洲日本最大视频资源| 大陆偷拍与自拍| 国产成人av激情在线播放| 80岁老熟妇乱子伦牲交| 国产单亲对白刺激| 成人国产一区最新在线观看| 成人手机av| 国产精品98久久久久久宅男小说| 久久人妻熟女aⅴ| 91在线观看av| 天堂动漫精品| 国产精品一区二区在线不卡| 国产精品99久久99久久久不卡| 精品国产乱码久久久久久男人| 精品国产乱子伦一区二区三区| 两个人免费观看高清视频| 欧美一级毛片孕妇| 国产野战对白在线观看| 国产av又大| a级毛片黄视频| 亚洲av电影在线进入| 91大片在线观看| 80岁老熟妇乱子伦牲交| 中文字幕色久视频| 老熟妇仑乱视频hdxx| 国产黄色免费在线视频| 搡老熟女国产l中国老女人| 亚洲精品美女久久av网站| 久久ye,这里只有精品| 巨乳人妻的诱惑在线观看| 亚洲人成电影免费在线| 国产精品98久久久久久宅男小说| 老汉色av国产亚洲站长工具| 在线av久久热| 国产成人欧美| 桃红色精品国产亚洲av| 成年版毛片免费区| 亚洲久久久国产精品| 色综合欧美亚洲国产小说| 热re99久久精品国产66热6| 天堂中文最新版在线下载| 久久ye,这里只有精品| 男女午夜视频在线观看| 嫁个100分男人电影在线观看| 亚洲av片天天在线观看| 一区福利在线观看| 久久久精品国产亚洲av高清涩受| av线在线观看网站| 亚洲午夜精品一区,二区,三区| 啦啦啦 在线观看视频| 久久久国产一区二区| 精品亚洲成国产av| 热99re8久久精品国产| 999久久久精品免费观看国产| 色精品久久人妻99蜜桃| 91精品国产国语对白视频| 免费在线观看完整版高清| √禁漫天堂资源中文www| 亚洲少妇的诱惑av| 无人区码免费观看不卡| 日韩免费av在线播放| 如日韩欧美国产精品一区二区三区| 亚洲专区国产一区二区| 欧美精品人与动牲交sv欧美| 亚洲精品中文字幕在线视频| 91成年电影在线观看| 欧美另类亚洲清纯唯美| 亚洲精品在线观看二区| 免费在线观看黄色视频的| 国产亚洲一区二区精品| 十八禁人妻一区二区| 侵犯人妻中文字幕一二三四区| 人成视频在线观看免费观看| 亚洲午夜精品一区,二区,三区| 男女免费视频国产| 欧美 亚洲 国产 日韩一| 黑丝袜美女国产一区| 国产1区2区3区精品| 欧美日韩中文字幕国产精品一区二区三区 | 91大片在线观看| 91麻豆精品激情在线观看国产 | 色综合婷婷激情| 国产精品一区二区在线观看99| 中文字幕高清在线视频| 久久精品人人爽人人爽视色| 9191精品国产免费久久| 亚洲人成77777在线视频| 精品高清国产在线一区| 黄色视频不卡| 中文字幕av电影在线播放| 日本黄色日本黄色录像| 国产aⅴ精品一区二区三区波| 国产精品免费一区二区三区在线 | 老司机靠b影院| 波多野结衣一区麻豆| 天堂俺去俺来也www色官网| 视频在线观看一区二区三区| 亚洲五月色婷婷综合| 十八禁人妻一区二区| 激情视频va一区二区三区| 久久午夜亚洲精品久久| 亚洲精品在线美女| 两性夫妻黄色片| 久久久久久久国产电影| 无人区码免费观看不卡| 涩涩av久久男人的天堂| 操美女的视频在线观看| 丝瓜视频免费看黄片| 久9热在线精品视频| 国产欧美日韩一区二区三区在线| 女警被强在线播放| 老熟女久久久| x7x7x7水蜜桃| 亚洲国产毛片av蜜桃av| 亚洲精品乱久久久久久| 91成人精品电影| 欧美丝袜亚洲另类 | 国产日韩一区二区三区精品不卡| 亚洲成人免费av在线播放| 久久久久久亚洲精品国产蜜桃av| 亚洲成人免费av在线播放| 少妇的丰满在线观看| av天堂在线播放| 99热网站在线观看| 亚洲av成人不卡在线观看播放网| 如日韩欧美国产精品一区二区三区| 美女视频免费永久观看网站| 久久九九热精品免费| 国产欧美日韩一区二区三区在线| 国产精品美女特级片免费视频播放器 | 97人妻天天添夜夜摸| 亚洲人成77777在线视频| 女人爽到高潮嗷嗷叫在线视频| 国产精品一区二区在线观看99| 男男h啪啪无遮挡| 久久香蕉国产精品| 黄色a级毛片大全视频| 一本一本久久a久久精品综合妖精| 久久人人爽av亚洲精品天堂| 宅男免费午夜| 9色porny在线观看| 在线天堂中文资源库| 国产成人免费无遮挡视频| 日韩视频一区二区在线观看| 国产在视频线精品| 18禁美女被吸乳视频| 亚洲国产精品合色在线| 夜夜爽天天搞| 国产亚洲av高清不卡| 老鸭窝网址在线观看| 色老头精品视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 窝窝影院91人妻| 亚洲熟妇熟女久久| 看黄色毛片网站| 亚洲人成77777在线视频| 黑人操中国人逼视频| 国产男女超爽视频在线观看| 亚洲av片天天在线观看| 久久久国产精品麻豆| 一级毛片精品| 五月开心婷婷网| 精品卡一卡二卡四卡免费| 91av网站免费观看| 精品福利观看| 99riav亚洲国产免费| ponron亚洲| 国产97色在线日韩免费| 一区二区日韩欧美中文字幕| 伦理电影免费视频| 日本黄色日本黄色录像| 亚洲av片天天在线观看| 精品免费久久久久久久清纯 | 老司机影院毛片| 亚洲va日本ⅴa欧美va伊人久久| 这个男人来自地球电影免费观看| 精品国产超薄肉色丝袜足j| 精品国产一区二区久久| 啦啦啦免费观看视频1| 最新美女视频免费是黄的| 亚洲精品国产精品久久久不卡| 午夜福利影视在线免费观看| 电影成人av| 老司机靠b影院| e午夜精品久久久久久久| 久热这里只有精品99| 精品国内亚洲2022精品成人 | 两人在一起打扑克的视频| 一a级毛片在线观看| 国产精品一区二区在线不卡| 亚洲精品美女久久av网站| 精品高清国产在线一区| 久久精品国产亚洲av高清一级| 亚洲国产精品合色在线| 黄片大片在线免费观看| ponron亚洲| 欧美丝袜亚洲另类 | 自拍欧美九色日韩亚洲蝌蚪91| 国产亚洲精品一区二区www | 欧美人与性动交α欧美软件| 欧美在线一区亚洲| 女同久久另类99精品国产91|