亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于鉻層與基體結合部疲勞損傷的轉管機槍身管壽命預測研究

        2016-11-09 10:08:59徐寧吳永海王永娟徐誠
        兵工學報 2016年10期
        關鍵詞:身管機槍壽命

        徐寧, 吳永海, 王永娟, 徐誠

        (1.南京理工大學 機械工程學院, 江蘇 南京 210094; 2.南京理工大學 格萊特納米科技研究所, 江蘇 南京 210094;3.淮陰工學院 機電工程實驗教學中心, 江蘇 淮陰 223003)

        ?

        基于鉻層與基體結合部疲勞損傷的轉管機槍身管壽命預測研究

        徐寧1,2, 吳永海1,3, 王永娟1, 徐誠1

        (1.南京理工大學 機械工程學院, 江蘇 南京 210094; 2.南京理工大學 格萊特納米科技研究所, 江蘇 南京 210094;3.淮陰工學院 機電工程實驗教學中心, 江蘇 淮陰 223003)

        為預測身管壽命,采用熱- 固耦合場理論、沖擊疲勞理論和疲勞累積損傷理論,提出基于身管鉻層與基體結合部疲勞損傷的身管壽命預測方法。以3管轉管機槍的身管壽命為例,建立有限元模型,計算鉻層與基體結合面的耦合應力,預測身管壽命,理論計算與試驗結果進行了對比,數據基本吻合;分析身管鍍鉻層厚度對身管壽命的影響。研究結果表明:鉻層與基體結合部的疲勞損傷壽命基本決定了身管的壽命,證明所提壽命預測方法是可行的;適當增加鍍鉻層厚度,可以提高身管壽命。

        兵器科學與技術; 自動武器; 身管; 耦合應力; 壽命; 疲勞; 預測模型

        0 引言

        自動武器發(fā)射時,身管內膛承受復雜、劇烈的熱流脈沖和膛壓載荷,載荷具有持續(xù)時間短、幅度大等特點。身管內膛薄層,特別是鉻層的溫度梯度、應力遠高于其他部位的應力。發(fā)射時在內彈道時期,內膛表面至內層0.92 mm左右的厚度范圍內,溫度差值和應力差值分別達到360 K和600 MPa,應力的幅值變化達到600~800 MPa,內部應力最高值達到820 MPa;在后效期和間隔期,身管內膛薄層的溫度和應力都迅速下降[1-4]。連續(xù)射擊時,身管內壁薄層承受如此劇烈的持續(xù)熱和膛壓載荷作用,內膛薄層特別是鉻層必然發(fā)生疲勞破壞,疲勞裂紋萌生和擴展將首先出現在接合力薄弱的兩種金屬接觸面——鉻鋼結合面部位,疲勞破壞最終將導致身管鉻層的脫落,身管失效、壽終。身管是自動武器的核心部件之一,同時也是自動武器中壽命最短的零件之一。在提高武器壽命的研究中,有關身管壽命的研究一直占有重要地位,其中大口徑轉管機槍的身管壽命問題尤為突出。目前,國內外對身管壽命預測作了大量研究,例如:基于隨機有限元法和身管熔化燒蝕理論模型[5]、基于灰色系統理論[6]等,但是國內外關于鉻層疲勞破壞對身管壽命影響尚未見報道。

        本文運用傳熱與變形(簡稱熱- 固或熱- 彈)耦合場理論、沖擊疲勞理論和疲勞累積損傷理論,研究轉管機槍的身管在高頻循環(huán)動態(tài)熱載荷和瞬態(tài)膛壓載荷耦合作用下基于鉻層與基體結合部疲勞損傷的身管壽命預估方法,分析身管的鉻層厚度對壽命的影響,為身管結構設計提供理論依據。

        1 轉管機槍身管數學模型

        根據熱力學第一定律[7-8],轉管機槍身管單元體熱量的平衡關系為

        (1)

        (2)

        式中:λ、G為Lame常數;ur為徑向變形;uz為軸向變形;r為半徑。

        2 轉管機槍鍍鉻身管有限元建模及熱- 固耦合應力計算

        2.1模型簡化及網格劃分

        身管是轉管機槍的核心零件之一,在發(fā)射過程中受到高溫、高壓及高速流動的火藥燃氣及彈帶的熱的、化學的及機械的作用,工作狀況很惡劣,若要對其進行完全真實的模擬,過于復雜。因此,在實際問題的研究中,對問題進行了簡化處理:1)忽略彈頭在運動過程中對身管的擠進力;2)在身管的線膛部分,熱- 固耦合計算不考慮膛線的影響,在幾何尺寸上,以陽線尺寸為準;3)忽略身管重力的影響;4)身管尾部受軸向約束;5)忽略身管定位面、孔、螺紋等局部細節(jié)的影響。

        身管由鍍鉻層和鋼本體構成,根據其結構和傳熱特點,建模時將其作為雙層管處理。為分析鉻鍍層對身管溫度和應力的影響,同時考慮到耦合分析中身管內壁薄層存在較大的溫度梯度,對內膛表面采用加密網格劃分,有限元模型采用HyperMesh工具生成,采用六面體單元劃分網格,如圖1為3管轉管武器的三維有限元模型,共1 082 032個節(jié)點, 991 038個單元。材料模型考慮了溫度對材料熱物理性能和力學性能的影響。為全面評價身管的壽命情況,綜合反映身管的熱和應力狀況,沿身管軸線取3個典型截面進行數據分析,其中A截面為膛線起始處截面,距身管尾端面127.5 mm,B截面為最大膛壓截面,距身管尾端面252 mm,C截面在膛口區(qū)域,距身管尾端面1 325.5 mm.

        圖1 3管有限元模型Fig.1    Finite element model of rotating barrels machine gun with 3 barrels

        2.2邊界條件

        2.2.1內壁邊界條件

        槍彈發(fā)射過程中,身管主要受到火藥燃氣瞬態(tài)壓力和瞬態(tài)熱沖擊作用。身管壁與火藥氣體熱交換的主要形式為強迫對流換熱,輻射換熱通過對火藥氣體強迫對流系數的修正來進行補償[9-12]。因此,內膛邊界條件為

        r=r1,

        (3)

        式中:T(z,r1)為膛內流場溫度函數,f(z,t)為內膛壓力函數,這兩個函數可由內彈道確定;T0為環(huán)境溫度;hi為火藥氣體強迫對流系數;g(z,t)為內膛流場溫度函數。

        2.2.2外壁邊界條件

        當武器有多根身管時,需要考慮外壁的邊界條件。武器發(fā)射時,空氣橫掠通過身管外壁,空氣在身管外壁上形成流動邊界層,沿與身管軸線垂直的方向流動[11-12]。身管外壁以強迫對流和輻射方式放熱,其外壁邊界條件為

        r=r2,

        (4)

        式中:ho為外壁氣流對流系數;ε為輻射率(黑度);σ為斯蒂芬- 波爾茲曼常數;A為輻射面積。

        2.3動態(tài)熱- 固耦合應力計算

        轉管機槍發(fā)射時,將熱載荷和膛壓載荷同時施加給槍膛內表面,加載見圖2所示,使用動態(tài)耦合應力計算法,可以得到身管應力情況。

        圖2 單發(fā)單管邊界條件Fig.2    The boundary conditions of single barrel during single shot

        圖3為連續(xù)射擊時,鉻鋼結合面上的Von Mises應力變化規(guī)律。從第2發(fā)開始,由熱載荷產生的應力峰值已經超過膛壓載荷產生的應力峰值,熱載荷居于主導地位。從第4發(fā)開始,由于殘余熱應力的存在,膛壓載荷的主要作用抑制熱載荷產生的應力,繼續(xù)連續(xù)發(fā)射,會出現類似現象。

        圖3 鉻鋼結合面的Von Mises應力(A截面)Fig.3Von Mises stress on chromium steel surface (Section A)

        3 身管壽命預測模型

        3.1身管鉻鋼結合面的疲勞方程

        身管鉻鋼結合面的疲勞可以看作熱和膛壓載荷共同作用下的低周疲勞,其壽命可以用Manson-Coffin疲勞方程[13-14]和Morrow方程等[15],可得當平均應力σm≠0和平均應變εm≠0時,身管疲勞方程為

        (5)

        材料在低周循環(huán)加載下,產生循環(huán)硬化或軟化,其循環(huán)應力- 應變關系與單調加載不同。低周循環(huán)加載時,剛開始,應力- 應變關系隨循環(huán)數而改變,但到達一定循環(huán)次數后,材料對變形的抗力趨于穩(wěn)定。由于循環(huán)穩(wěn)定階段占疲勞壽命的大部分,所以通常以穩(wěn)定的循環(huán)應力- 應變曲線來代表材料的循環(huán)應力、應變性質,其表達式為

        (6)

        式中:Δσ為應力幅;n′為循環(huán)應變硬化指數,n′=0.1~0.2.

        3.2結構局部應力、應變計算

        局部應力、應變近似計算的典型方法是Neuber法,為提高疲勞壽命預測精度,利用修正的Neuber公式[16-17]計算構件局部應力、應變。在循環(huán)加載過程中,修正的Neuber公式可寫為

        (7)

        3.3疲勞累積損傷理論

        Miner線性疲勞累積損傷理論具有很好可驗性和可行性,在工程上得到廣泛應用。設加載歷史由δ1,δ2,…,δl等l個不同的應力水平構成,各應力水平下的疲勞壽命依次為N1,N2,…,Nl,各應力水平下的循環(huán)次數依次為n1,n2,…,nl,則零件的疲勞壽命[18-19]為

        (8)

        4 3管轉管機槍的身管壽命預測

        基于鍍鉻3管轉管機槍的身管有限元模型,首先計算身管鉻層與基體結合面上每發(fā)彈的應力峰值、谷值及平均值,在此基礎上,采用本文所建的身管鉻層與基體結合部疲勞壽命的壽命模型推斷身管壽命,并分析了身管鍍層厚度對壽命的影響。所采用的射擊規(guī)范為:每組發(fā)射102發(fā):連射18發(fā),停5 s,連射18發(fā),停5 s,連射33發(fā),停5 s,連射33發(fā),停4 min,循環(huán)射擊。

        4.13管轉管機槍的身管壽命預測結果及試驗結果比較

        4.1.1計算結果

        圖4 鉻鋼結合面(A截面)計算結果Fig.4    Calculated results of Section A of chromium steel interface

        利用有限元模型,可以分別計算出身管3個典型截面A、B、C的每發(fā)Von Mises應力峰值、谷值和平均值,典型計算結構如圖4為鉻鋼結合面(A截面)的應力情況。本文研究的槍管材料是高強度合金鋼,經過膛線冷精鍛加工、熱處理、表面鍍鉻和表面激光強化后,材料屈服強度達1 360 MPa,抗拉強度達1 620 MPa. 由圖4(a)和圖4(b)可以看出,連續(xù)射擊時,鉻鋼結合面上的每發(fā)Von Mises應力峰值、谷值和平均值均大幅上升。3管轉管機槍連射18發(fā),對應于每個身管6發(fā)連續(xù)發(fā)射,每發(fā)應力幅值持續(xù)上升;3管連射33發(fā),對應于每個身管11連續(xù)發(fā)射,應力幅值呈現先單邊上升,在第7發(fā)時達到頂點,隨后每發(fā)應力幅值緩慢下降的規(guī)律。與應力規(guī)律對應的是,連發(fā)時,隨著平均應力、應力幅的大幅上升,身管等幅載荷對應的疲勞壽命迅速衰減,單個身管6發(fā)連續(xù)發(fā)射時,對應第6發(fā)彈的身管等幅載荷的疲勞壽命較第1發(fā)彈下降53.8%,單個身管11發(fā)連續(xù)發(fā)射時,對應第11發(fā)彈的身管等幅載荷下的疲勞壽命較第1發(fā)彈下降65%,見圖4(d)。由圖4(c)可以看出,射擊間隔對每發(fā)身管等幅載荷下的疲勞壽命影響很大。經過5 s停射散熱后,每發(fā)身管等幅載荷下的疲勞壽命迅速回升,對應第7發(fā)彈的身管等幅載荷下的疲勞壽命為6 700發(fā),僅與第1發(fā)彈相差1.3%. 經4 min停射散熱后,第2組(至68發(fā))、第3組(至102發(fā))與第1組(至34發(fā))在每發(fā)槍管等幅載荷下疲勞壽命均相差無幾。

        采用疲勞累計損傷理論及壽命預測規(guī)范[20],預測身管3個典型截面A、B、C的疲勞壽命計算結果,如表1所示。從表1中可以看出,A截面的壽命相對較小,這與A截面所在的位置(膛線起始區(qū)域)有關,由于A截面的壁厚較大,在強迫對流散熱的工況下其熱量存儲在壁內,不易散失。另外,該截面的熱載荷較大,連續(xù)射擊時其平均應力和應力幅值均較大。

        4.1.2身管壽命試驗驗證

        試驗結果見表2,從中可以看到,當身管發(fā)射累計1 700發(fā)開始,彈頭出槍口的初速呈明顯的直線下降趨勢,到3 000發(fā)時,初速的下降率已達到5.27%~6.09%,2支槍管的平均壽命為3 535發(fā)。由于槍管壁面承受多次劇烈的火藥氣體脈沖熱和膛壓載荷作用,“鉻層與鋼體結合部位” 疲勞破壞,導致槍管鉻層脫落,膛線破壞,橫彈產生,故試驗的壽終判據“橫彈率”與“鉻層與鋼體結合部位” 疲勞壽命直接相關。采用本文所提出的壽命預測放大的計算結果與試驗結果相比,誤差小于5%,證實了本文所建的模型和方法是正確、可行的。

        表1 槍管截面壽命

        表2 3管轉管機槍壽命試驗結果

        4.2鉻層厚度的影響分析

        為研究鉻層厚度對槍管壽命的影響,設定3種方案,分別取鉻層厚度為0.15 mm、0.25 mm、0.41 mm進行對比分析,計算結果見圖5和表3. 從圖5可以看出,鉻層厚度對于身管鉻鋼結合面上的每1發(fā)彈Von Mises最大應力、Von Mises平均應力和壽命有著顯著影響。鉻層厚度越大,每1發(fā)彈Von Mises最大應力、Von Mises平均應力越小,每1發(fā)彈的壽命越高。射擊中,第1個6連發(fā)中的第6發(fā)彈時,鉻層厚度為0.15 mm的方案每1發(fā)彈壽命比鉻層厚度為0.41 mm的方案低39.7%;第34發(fā)彈時,鉻層厚度為0.15 mm的方案每1發(fā)彈壽命比鉻層厚度為0.41 mm的方案低32.1%;第102發(fā)彈時,鉻層厚度為0.15 mm的方案每1發(fā)彈壽命比鉻層厚度為0.41 mm的方案低30.5%. 由表3可以看出,鉻層厚度為0.41 mm的方案身管總壽命比鉻層

        表3 鉻層厚度不同時的槍管總壽命(A截面)

        圖5 A截面的每發(fā)彈最大應力、平均應力及壽命Fig.5    The comparison of max stress, average stress and lifetime of Section A during each shot

        厚度為0.25 mm的方案高30.2%. 綜上所述,鉻層厚度的增加對于提高身管總壽命十分有利,在工藝條件允許的情況下,應合理增加鉻層厚度。

        5 結論

        本文建立了身管的溫度和應力耦合計算模型,提出了一種基于身管鉻層與基體結合部疲勞壽命的身管壽命預測模型,研究了轉管機槍身管在熱載荷和膛壓載荷耦合作用下的身管壽命預測問題,理論計算與試驗結果對比說明,鍍鉻層與基體結合部的疲勞壽命基本決定了身管的壽命,驗證了本文模型的正確性,為轉管機槍槍管壽命推斷提供了理論與方法。本文也分析了鉻層厚度對3管轉管機槍的身管壽命影響,分析結果表明,在工藝條件許可的條件下,適當增加鉻層厚度,可以提高身管壽命。

        References)

        [1]Bundy M L, Gerber N, Bradley J W. Gun barrel erosion and bore temperature modeling[C]∥Proceedings of the 1993 Summer Computer Simulation Conference. Boston,US: Society for Computer Simulation, 1993:397-402.

        [2]Johnston I A.Understanding and predicting gun barrel erosion, DSTO-TR-1757[R]. Australia: Weapons Systems Division,Defence Science and Technology Organisation,2005.

        [3]Kerwien S C, Rocha R. Metallurgical failure analysis of the M230 machine gun barrel-MIF-A-001-1999-I, ADA366133[R]. NJ, US:Army Armament Research Development and Engineering Center, 1999.

        [4]Tzeng J T. Dynamic fracture of composite gun tubes, ARL-TR-1869[R].MD, US:US Army Research Laboratory, 1999.

        [5]方峻,吳華清.融合理論退化模擬與試驗數據的身管壽命預測[J].機械科學及技術,2014,33(10): 1468-1472.

        FANG Jun, WU Hua-qing. Prediction of the barrel life based on the theoretical degradation simulation and the experimental data[J]. Mechanical Science and Technology for Aerospace Engineering, 2014, 33(10):1468-1472. (in Chinese)

        [6]孟祥飛,王昌明,何博俠, 等. 基于灰色線性回歸組合模型的火炮身管壽命預測[J]. 南京理工大學學報,2012,36(4):636-638.

        MENG Xiang-fei, WANG Chang-ming, HE Bo-xia, et al. Prediction ofgun barrel life based on grey-linear regression combined model[J]. Journal of Nanjing University of Science and Technology, 2012, 36(4):635-638. (in Chinese)

        [7]王洪綱. 熱彈性力學概論[M]. 北京:清華大學出版社, 1989.

        WANG Hong-gang. Introduction to thermoelasticity [M].Beijing:Tsinghua University Press, 1989.(in Chinese)

        [8]李維特, 黃保海, 畢仲波. 熱應力理論分析及應用[M]. 北京:中國電力出版社, 2004.

        LI Wei-te, HUANG Bao-hai, BI Zhong-bo.Theoretical analysis and application of thermal stress [M].Beijing: China Electric Power Press, 2004.(in Chinese)

        [9]吳永海, 徐誠, 李峰,等. 承受熱沖擊的大口徑機槍槍管的熱效應分析[J]. 南京理工大學學報, 2007, 31(1):1-5.

        WU Yong hai,XU Cheng, LI Feng, et al. Thermal effect of large-calibre machine-gun barrel subjected to thermal shock[J]. Journal of Nanjing University of Science and Technology, 2007, 31(1):1-5.(in Chinese)

        [10]吳永海, 徐誠, 張海兵. 某大口徑機槍槍管的瞬態(tài)熱彈耦合動力響應分析[J]. 彈道學報, 2007, 18(4):16-20.

        WU Yong-hai, XU Cheng, ZHANG Hai-bing. Dynamic response analysis of transient thermoelasticity couple of a large-calibre machine gun barrel[J]. Journal of Ballistics, 2007, 18(4):16-20.(in Chinese)

        [11]吳永海. 速射武器身管現代設計理論研究 [D].南京:南京理工大學,2010.

        WU Yong-hai. Modern design theory of rapid-fire weapon barrel [D].Nanjing: Nanjing University of Technology and Engineering, 2010.(in Chinese)

        [12]Cheng X, Wu Y, Guan X. The transient stress and life analysis of a gun tube[C]∥ASME 2014 Pressure Vessels and Piping Conference. Anaheim, CA, US:American Society of Mechanical Engineers, 2014.

        [13]Manson S S. Fatigue: a complex subject—some simple approximations[J]. Experimental Mechanics,1965,5(4): 193-226.

        [14]Wetzel R M, Bussa S. Fatigue under complex loading: analyses and experiments[M]. Warrendale, PA, US: Society of Automotive Engineers, 1997.

        [15]Morrow J D. Laboratory simulation of the low cycle fatigue behavior of the hook region of a steam turbine blade subjected to start stop cycles[C]∥Fourth National Congress on Pressure Vessel and Piping Technology.Portland, OR, US: Society of Automotive Engineers, 1983.

        [16]馬康民.某型飛機主起落架半輪叉疲勞壽命分析[J].強度與環(huán)境,1999(1):37-42.

        MA Kang-min. An analysis of fatigue life for half-fork of main landing gear of aircraft[J]. Structure & Environment Engineering, 1999(1):37-42.(in Chinese)

        [17]舒陶, 任宏光, 郭克平. 局部應力應變Neuber法與有限元求法的比較[J]. 彈箭與制導學報, 2009, 29(1):267-269.

        SHU Tao, REN Hong-guang, GUO Ke-ping. Comparison of Neuber’s rule and FEM for local stress-strain[J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2009, 29(1):267-269.(in Chinese)

        [18]劉維信.機械可靠性設計[M].北京:清華大學出版社, 1996.

        LIU Wei-xin.Mechanical reliability design[M]. Beijing:Tsinghua University Press, 1996.(in Chinese)

        [19]趙少汴. 常用累積損傷理論疲勞壽命估算精度的試驗研究[J]. 機械強度, 2000, 22(3):206-209.

        ZHAO Shao-bian. Study on the accuracy of fatigue life prediction by the generally used damage accumulation theory[J].Journal of Mechanical Strength, 2000, 22(3):206-209.(in Chinese)

        [20]呂文軍,王惠源,陳輝. 美軍轉管武器試驗規(guī)范分析[J]. 四川兵工學報,2015,36(3):60-63,75.

        LYU Wen-jun, WANG Hui-yuan, CHEN Hui. Analysis of test specification of US Gatling weapons[J]. Journal of Sichuan Ordnance, 2015, 36(3):60-63,75.(in Chinese)

        Barrel Life Prediction of Rotating Barrels Machine Gun Based on Fatigue Damage of Chromium-steel Interface

        XU Ning1,2, WU Yong-hai1,3, WANG Yong-juan1, XU Cheng1

        (1.School of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing 210094, Jiangsu, China; 2.Herbert Gleiter Institute of Nanoscience, Nanjing University of Science and Technology, Nanjing 210094, Jiangsu, China; 3.Experimental-teaching Center of Mechanical and Electrical Engineering, Huaiyin Institute of Technology, Huaiyin 223003, Jiangsu,China)

        To predict the life of a gun barrel, the theories of thermal-solid coupled stress, the impact fatigue and fatigue cumulative damage are used to derive the life prediction of a gun barrel based on fatigue damage of chromium-steel interface. A 3D finite element model for a 3-barrels rotating machine gun is built to calculate the coupling stress of chromium-steel interface. The barrel life of the rotating barrels machine gun is estimated. The estimated result is basically consistent with the experimental result. The analysis results show that the fatigue damage life at chromium-steel interface almost determines the life of the gun barrel. The effect of the thickness of chromium coating on the life of gun barrel is studied. The research result shows that the life of gun barrel can be improved by increasing the thickness of chromium coating properly.

        ordnance science and technology; automatic weapon; gun barrel; coupling stress; life; fatigue; prediction model

        2016-01-07

        國家自然科學基金項目(51575279)

        徐寧(1982—),男,博士研究生。E-mail:njustxn@163.com;

        王永娟(1972—),女,教授,碩士生導師。E-mail:13951643935@139.com

        TJ202

        A

        1000-1093(2016)10-1926-08

        10.3969/j.issn.1000-1093.2016.10.020

        猜你喜歡
        身管機槍壽命
        火炮身管壽命評估方法及其標準研究
        M240通用機槍
        人類壽命極限應在120~150歲之間
        中老年保健(2021年8期)2021-12-02 23:55:49
        M134 機槍
        玩轉百變機槍
        倉鼠的壽命知多少
        基于垂向固有振動仿真的身管口徑倍數設計研究
        那些被“黑”的機槍們
        馬烈光養(yǎng)生之悟 自靜其心延壽命
        華人時刊(2018年17期)2018-12-07 01:02:20
        基于數值仿真的四種類型身管強度比較
        機械工程師(2018年6期)2018-06-14 06:12:28
        亚洲国产日韩a在线乱码| 国产精品色内内在线播放| 国产精品国产午夜免费看福利| 一本久道久久综合婷婷五月| 久久久亚洲欧洲日产国码是AV| 亚洲国产精品色一区二区| 日本国产精品久久一线| 国产情侣一区二区| 不卡一卡二卡三乱码免费网站| 国产精品免费久久久久影院仙踪林 | 中文字幕乱码熟女人妻水蜜桃| 狠狠狠色丁香婷婷综合激情| 不卡无毒免费毛片视频观看| 一区二区三区国产天堂| 久久亚洲中文字幕乱码| 国产女主播白浆在线观看| 久久精品国产99国产精2020丨| 亚洲每天色在线观看视频| 日韩精品视频在线观看无| 一本大道熟女人妻中文字幕在线| 99久久er这里只有精品18| 高清无码一区二区在线观看吞精| 国产区高清在线一区二区三区| 国产精品亚洲二区在线看| 国产又色又爽又刺激在线播放| 日韩免费小视频| 久久久人妻丰满熟妇av蜜臀| 亚洲2022国产成人精品无码区| 亚洲欧美另类激情综合区| 亚洲香蕉毛片久久网站老妇人| 一区二区三区观看在线视频| 亚洲一区二区三区99| 亚洲精品无码久久久影院相关影片| 波多野结衣视频网址| 日韩精品有码中文字幕| 免费观看全黄做爰大片| 色婷婷综合中文久久一本| 久久无码中文字幕东京热| 国产草逼视频免费观看| 48沈阳熟女高潮嗷嗷叫| 国产av综合一区二区三区最新|