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

        ?

        掃描電子束能量峰值對表面溫度場的影響規(guī)律

        2022-07-27 02:59:12李新凱王榮胡露瑤任旭隆王喜社
        表面技術(shù) 2022年7期
        關(guān)鍵詞:電子束熱源峰值

        李新凱,王榮,胡露瑤,任旭隆,王喜社

        掃描電子束能量峰值對表面溫度場的影響規(guī)律

        李新凱1,王榮1,胡露瑤2,任旭隆1,王喜社1

        (1.桂林電子科技大學,廣西 桂林 541004;2.桂林旅游學院,廣西 桂林 541004)

        掃描電子束;能量;溫度場;熱源;表面改性

        近年來,電子束以其非接觸、可控性強、真空無污染、能量利用率高等優(yōu)點而廣泛應(yīng)用于材料表面改性領(lǐng)域中[1-2]。常見的電子束表面處理技術(shù)包括電子束表面合金化、表面淬火、表面非晶態(tài)、表面熔凝等,以上電子束技術(shù)均是對金屬材料表層進行處理,需實現(xiàn)電子束大面域、均勻、穩(wěn)定下束[3-4]。

        目前,國內(nèi)外學者已對強流脈沖電子束與大面域輻照兩種電子束表面改性加工方法進行了深入研究,主要考察電子束類型[5]及工藝參數(shù)[6]、電子束能量分布規(guī)律[7]以及表面改性優(yōu)化方法[8]等方面。研究發(fā)現(xiàn)脈沖電子束能量密度較大,適用于處理高熔點、大深寬比金屬表面改性,但處理過程中表層金屬會發(fā)生濺射產(chǎn)生“熔坑”缺陷,同時高能作用下表層經(jīng)歷驟熱急冷過程,易造成內(nèi)應(yīng)力與組織應(yīng)力集中,產(chǎn)生結(jié)構(gòu)裂紋缺陷[9-11]。大面積電子束輻照是通過電子束散焦的方式實現(xiàn)最大面積為60 mm范圍的改性,因下束面積大所以更適用于表面微熔處理,然而該方式對電子槍功率要求較高,較難實現(xiàn)電子束能量的均勻分布[12]。本團隊在此基礎(chǔ)上,針對電子束表面微熔處理能量密度均勻且穩(wěn)定分布的需求,開發(fā)了一種新型連續(xù)掃描電子束技術(shù)(Continuous Scanning Electron Beam Technique Process,CSEBP),通過聚焦線圈與偏轉(zhuǎn)線圈共同作用,實現(xiàn)電子束聚焦的同時以高頻率旋轉(zhuǎn)實現(xiàn)環(huán)狀下束效果[13-15]。研究發(fā)現(xiàn)環(huán)狀掃描下束下電子槍無需高功率即可實現(xiàn)較高能量密度的均勻分布。另外,CSEBP能量密度分布除了受到電子束束流以及加速電壓等電子槍參數(shù)影響外,還受到能量分布特征參數(shù)影。因此,有必要對電子束下束過程中能量分布規(guī)律進行探索,以指導CSEBP在表面改性領(lǐng)域的應(yīng)用。

        本文以45鋼為電子束表面改性為研究對象,通過數(shù)值求解、仿真計算與試驗驗證相結(jié)合的方式探究能量峰值系數(shù)對電子束能量分布的影響規(guī)律,詳細探討了多種峰值系數(shù)下45鋼表層溫度的變化規(guī)律,并通過電子束微熔試驗進行驗證,力圖為電子束大面域掃描提供新的方法和理論依據(jù)。

        1 掃描電子束微熔處理數(shù)學物理模型的建立

        1.1 CSEBP下束方式原理

        本課題組自主研發(fā)的掃描電子束下束方式如圖1所示,其中,為工件移動速度,為掃描帶半徑,為掃描長度,掃描帶縱向5個均分點為后續(xù)熱循環(huán)曲線仿真取樣點。其原理是將編輯好的電子束掃描軌跡及運動方式的控制程序輸入信號發(fā)生器,用方程x+y=2來描述圓形的電子束掃描軌跡,其中方向的分量分別為cos,sin(0≤≤π),通過所產(chǎn)生的偏轉(zhuǎn)磁場實現(xiàn)電子槍內(nèi)部束流沿固定角度傾斜與高頻旋轉(zhuǎn),從而實現(xiàn)環(huán)狀掃描電子束軌跡。

        圖1 連續(xù)掃描電子束(CSEBP)示意圖

        1.2 環(huán)狀下束電子束能量分布數(shù)學模型

        環(huán)狀電子束熱源與常用的高斯熱源、雙橢球熱源有較大不同,其作用形式為高斯熱源在環(huán)形區(qū)域的疊加。為此以高斯熱源為基礎(chǔ)搭建環(huán)狀電子束能量密度函數(shù),如圖2和圖3所示。

        圖2 電子束能量密度分布圖

        設(shè)表面中心處的熱流密度()為:

        式中:為任意一點的熱流密度(W);m為最大熱流密度(W/m2);為任意一點到加熱中心的距離(m);為熱源集中系數(shù)。

        掃描電子束微熔處理時在掃描面的功率為:

        由式(2)得:

        將式(3)代入(1)得:

        通常情況下掃描電子束處理時,一般取95%的有效能量范圍[18]。

        由此可得:

        因電子束作用過程中,束斑內(nèi)能量峰值可調(diào),為此引入能量峰值位置參數(shù),綜合以上計算可將掃描電子束能量分布數(shù)學模型表征為:

        式中:Rx為電子束下束環(huán)外徑(mm);Rp為電子束能量峰值位置距離z軸的距離(mm);rx為圓環(huán)內(nèi)徑(mm);為束流偏轉(zhuǎn)角(°);為能量峰位置系數(shù),;z為離焦量(mm),R0為束斑寬度(mm)。

        1.3 能量峰值對電子束熱源模型的影響

        圖4 不同入射角度和能量峰值系數(shù)下的電子束能量分布模擬圖

        上述所建立的電子束熱源數(shù)學模型中引入能量峰位置系數(shù),可對環(huán)形區(qū)域內(nèi)的電子束能量分布狀態(tài)作更加詳細的定義,該系數(shù)對平衡環(huán)形束斑內(nèi)外側(cè)溫度補償與中心熱傳導有重要意義。

        2 建立掃描電子束有限元模型

        2.1 模型假設(shè)

        電子束表面處理過程是一個驟熱急冷的非穩(wěn)態(tài)過程,表層的金屬的熔融與凝固均在極短時間內(nèi)完成,為簡化模型、減少運算量,對模型作出以下假設(shè):45鋼的熱物性參數(shù)為溫度的函數(shù);電子束掃描過程中熱傳導處于穩(wěn)態(tài);樣為各向同性的均勻介質(zhì);試樣被處理前的溫度和所處的工作室溫度均為300 K;不考慮熱對流;忽略組織相變引起的塑性變形[19-20]。

        2.2 幾何建模與網(wǎng)格劃分

        采用COMSOL軟件對45鋼掃描電子束微熔過程進行模擬,模型尺寸與實際試樣尺寸一致,為50 mm× 50 mm×50 mm,在工作面選取50 mm×8 mm區(qū)域作為電子束加工區(qū)域,并進行網(wǎng)格加密處理,電子束掃描區(qū)域采用六面體的單元類型進行較細的網(wǎng)格劃分,其他區(qū)域采用四面體的單元類型進行智能網(wǎng)格劃分,網(wǎng)格劃分模型如圖5所示[21-22]。模型內(nèi)部為均勻介質(zhì),所以將試樣中某處定義為微元控制體積ddd。

        電子束掃描試樣表面時,熱傳導會通過控制體積的各個面發(fā)生。利用泰勒公式展開控制表面的導熱速度,直角坐標系中的熱擴散方程為:

        式中:為密度(kg/m3);為定壓比熱容(J·kg/K);k為導熱系數(shù)(W/(m·K))。

        2.3 電子束熱源模型及邊界條件

        掃描電子束是以高能電子束轟擊金屬表面產(chǎn)生的高溫為熱源,并以指定速度平移實現(xiàn)大面域掃描,由上述搭建能量分布數(shù)學模型,其移動熱源與時間的函數(shù)表達式為:

        掃描電子束微熔處理是在真空環(huán)境中進行,真空度為10?2Pa,故可忽略空氣熱對流造成的熱量損失,而熱輻射的傳遞不需借助任何介質(zhì)[23]。因此可認為基體的導熱與表面的熱輻射是工件主要散熱方式,電子束熱源作用區(qū)域的傳熱方程為:

        熱輻射滿足第三類邊界條件,可由斯蒂芬–波爾茲曼方程來計算:

        3 掃描電子束微熔處理溫度場仿真分析

        根據(jù)前期研究結(jié)果[24-25],溫度場仿真工藝參數(shù)為:=60 kV,=5 mA,=3 mm/s,=230 mm,R+0/2=2 mm,=5°,=400 Hz,整個電子束加熱時間為16.7 s。

        3.1 電子束能量峰值對熱循環(huán)過程的影響

        圖6 掃描帶上不同點熱循環(huán)曲線

        3.2 電子束能量峰值對溫度場的影響

        圖7 掃描區(qū)溫度云圖

        4 掃描電子束微熔處理試驗

        試驗選用45鋼作為實驗用原材料,使用銑床將原材料加工成50 mm×50 mm×50 mm的立方體,加工過程中通過控制進刀量、銑削速度、主軸轉(zhuǎn)速等參數(shù)恒定,將試樣表面粗糙度控制在1.9~2.0 μm內(nèi),電子束加工前使用酒精擦拭表面并風干,去除表面油污。使用HDZ–6F型高壓數(shù)控真空電子束機進行表面處理試驗,實驗參數(shù)與數(shù)值模擬參數(shù)一致。采用光學顯微鏡對熔融層進行觀測,采用OLS4100激光顯微鏡測試處理后表面粗糙度。每個待測面均勻測量5次粗糙度,取其均值作為該面粗糙度值。

        圖8 不同峰值系數(shù)下的表面形貌

        5 結(jié)論

        [1] LEE T, BIAN Hua-kang, AOYAGI K, et al. Fabricating 9-12 Cr Ferritic/Martensitic Steels Using Selective Electron Beam Melting[J]. Materials Letters, 2020, 271: 127747.

        [2] GUO Shun, ZHOU Qi, KONG Jian, et al. Effect of Beam Offset on the Characteristics of Copper/304stainless Steel Electron Beam Welding[J]. Vacuum, 2016, 128: 205-212.

        [3] LI Xin-kai, WANG Rong, XIN Zhe, et al. Changes in Surface Roughness and Microstructure of 45 Steel after Irradiation by Electron Beam[J]. Materials Letters, 2021, 296: 129934.

        [4] IVANOV Y F, ZAGULIAEV D V, GLEZER A M, et al. Changes in Surface Structure and Mechanical Chara-cteristics of Al-5 wt%Si Alloy after Irradiation by Ele-ctron Beam[J]. Materials Letters, 2020, 275: 128105.

        [5] LU Jian, SUI Xin-meng, YANG Bo-hang, et al. Ultrafast In-Situ Transformation of Graphite into Graphene Nano-sheets by High Current Pulsed Electron Beam Direct Irradiation[J]. Applied Surface Science, 2022, 572: 151498.

        [6] LEI Shuang, LI Xian-feng, DENG Ya-qi, et al. Micro-structure and Mechanical Properties of Electron Beam Freeform Fabricated TiB2/Al-Cu Composite[J]. Materials Letters, 2020, 277: 128273.

        [7] LI Xin-kai, WANG Rong, WANG Qi-chao, et al. Scanning Electron Beam Polishing and Defect Analysis of 65 Steel[J]. Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions With Materials and Atoms, 2021, 490: 34-38.

        [8] 陳軍, 李偉, 賀冬云, 等. 強流脈沖電子束表面改性FeCrAl涂層的顯微組織及耐高溫腐蝕性能研究[J]. 表面技術(shù), 2020, 49(5): 200-206.

        CHEN Jun, LI Wei, HE Dong-yun, et al. Surface Micro-structure and High-Temperature Erosion Resistance of FeCrAl Coating after High Current Pulsed Electron Beam Treatment[J]. Surface Technology, 2020, 49(5): 200-206.

        [9] LU Jian, WEI De-qiang, WANG Rong, et al. Surface Polishing and Modification of 3Cr2Mo Mold Steel by Electron Beam Irradiation[J]. Vacuum, 2017, 143: 283- 287.

        [10] LV Peng, SUN Xiao, CAI Jie, et al. Microstructure and High Temperature Oxidation Resistance of Nickel Based Alloy GH4169 Irradiated by High Current Pulsed Electron Beam[J]. Surface and Coatings Technology, 2017, 309: 401-409.

        [11] HAO Sheng-zhi, WANG Hui-hui, ZHAO Li-min. Surface Modification of 40CrNiMo7Steel with High Current Pulsed Electron Beam Treatment[J]. Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions With Materials and Atoms, 2016, 368: 81-85.

        [12] FU Yu-lei, HU Jing, SHEN Xian-feng, et al. Surface Hardening of 30CrMnSiA Steel Using Continuous Electron Beam[J]. Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions With Materials and Atoms, 2017, 410: 207-214.

        [13] WEI De-qiang, WANG Xiao-bing, WANG Rong, et al. Surface Modification of 5CrMnMo Steel with Continuous Scanning Electron Beam Process[J]. Vacuum, 2018, 149: 118-123.

        [14] CUI Hong-yang, WANG Rong, WEI De-qiang, et al. Surface Modification of the Carbon Tool Steel by Continuous Scanning Electron Beam Process[J]. Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions With Materials and Atoms, 2019, 440: 156-162.

        [15] LI Xin-kai, WANG Rong, DONG Yu-jian, et al. Surface Morphology and Grain Size of 45 Steel after Scanning by Electron Beam[J]. Materials Letters, 2021, 297: 129884.

        [16] 李廣琪, 王麗芳, 朱剛賢, 等. 掃描方式對中空環(huán)形激光熔覆層殘余應(yīng)力及基板變形的影響研究[J]. 表面技術(shù), 2021, 50(3): 158-170.

        LI Guang-qi, WANG Li-fang, ZHU Gang-xian, et al. Influence of Scanning Patterns on Residual Stress of Cladding Layer and Substrate Deformation Produced by Hollow-Ring Laser Cladding[J]. Surface Technology, 2021, 50(3): 158-170.

        [17] 胡興, 彭昭成, 馮廣杰, 等. SUS310S不銹鋼局部真空電子束焊接接頭殘余應(yīng)力及變形研究[J]. 機械工程學報, 2020, 56(21): 38-47.

        HU Xing, PENG Zhao-cheng, FENG Guang-jie, et al. Numerical Simulation of Residual Stress and Deformation of SUS310S Stainless Steel Local Vacuum Electron Beam Welded Joint[J]. Journal of Mechanical Engineering, 2020, 56(21): 38-47.

        [18] 魏德強, 李新凱, 王曉冰. 電子束拋光技術(shù)的研究進展[J]. 表面技術(shù), 2016, 45(4): 175-182.

        WEI De-qiang, LI Xin-kai, WANG Xiao-bing. Research Progress of Electron Beam Polishing Technology[J]. Surface Technology, 2016, 45(4): 175-182.

        [19] 郭順, 羅添元, 彭勇, 等. Ti/Cu異種金屬電子束焊接界面行為[J]. 焊接學報, 2019, 40(8): 26-32, 162.

        GUO Shun, LUO Tian-yuan, PENG Yong, et al. Interface Behavior and Mechanical Properties of Ti/Cu Dissimilar Metals Welding by Electron Beam[J]. Transactions of the China Welding Institution, 2019, 40(8): 26-32, 162.

        [20] 房玉超, 楊子酉, 何景山. 電子束點焊熔池的液態(tài)金屬沖刷效應(yīng)作用規(guī)律[J]. 焊接學報, 2019, 40(6): 137-142, 167.

        FANG Yu-chao, YANG Zi-you, HE Jing-shan. Study on Liquid Metal Flushing Effect during Electron Beam Spot Welding[J]. Transactions of the China Welding Institu-tion, 2019, 40(6): 137-142, 167.

        [21] 趙桐, 唐振云, 劉巧沐, 等. GH4065A合金電子束焊接工藝及接頭組織性能[J]. 材料導報, 2020, 34(22): 22105-22110.

        ZHAO Tong, TANG Zhen-yun, LIU Qiao-mu, et al. Electron Beam Welding Process and Microstructure and Properties of Joint of GH4065A Alloy[J]. Materials Reports, 2020, 34(22): 22105-22110.

        [22] 郭超, 林峰, 葛文君. 電子束選區(qū)熔化成形316L不銹鋼的工藝研究[J]. 機械工程學報, 2014, 50(21): 152-158.

        GUO Chao, LIN Feng, GE Wen-jun. Study on the Fab-rication Process of 316L Stainless Steel via Electron Beam Selective Melting[J]. Journal of Mechanical Engi-neering, 2014, 50(21): 152-158.

        [23] REISGEN U, OLSCHOK S, KRICHEL T, et al. Deter-mination of the Influence of Welding Parameters on the Efficiency of Electron Beam Welding by Measurement of Backscattered Electrons[J]. Vacuum, 2019, 159: 182-185.

        [24] 李新凱, 王榮, 王啟超, 等. 掃描電子束微熔拋光臨界功率密度規(guī)律及實驗研究[J]. 表面技術(shù), 2021, 50(7): 386-393.

        LI Xin-kai, WANG Rong, WANG Qi-chao, et al. Research on Critical Power Density and Experiment of Scanning Electron Beam Micro-Melting Polishing[J]. Surface Tech-nology, 2021, 50(7): 386-393.

        [25] Lu Jian, Wei De-qiang, Wang Rong, et al. Surface Polishing and Modification of 3Cr2Mo Mold Steel by Electron Beam Irradiation[j]. Vacuum, 2017, 143: 283-287.

        [26] NIU Shao-qiang, YOU Qi-fan, YOU Xiao-gang, et al. Mechanism of Impurities Reduction and Evaporation of Alloying Elements for a Multi-Elements Ni-Based Superalloy during Electron Beam Remelting[J]. Vacuum, 2018, 156: 345-350.

        [27] WANG Rong, YU Jie, WEI De-qiang, et al. Surface Microstructures and Improved Mechanical Property of 40CrMn Steel Induced by Continuous Scanning Electron Beam Process[J]. Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions With Materials and Atoms, 2019, 459: 130-136.

        The Influence of the Peak Energy of Scanning Electron Beam on the Surface Temperature Field

        1,1,2,1,1

        (1. Guilin University of Electronic Technology, Guangxi Guilin 541004, China; 2. Guilin Tourism University, Guangxi Guilin 541004, China)

        The mathematical model of energy distribution in the circular downward beam mode of scanning electron beams is clearly defined. Obtain the influence law of energy crest factor on the surface temperature field of 45 steel. Based on the Gaussian heat source model, the energy peak position parameter is introduced to calculate the mathematical model of the energy distribution in the scanning electron beam downward beam mode. The COMSOL software was used to simulate the thermal cycle curve and temperature field of the scanning zone. Revised the electron beam heat source model in the ring-shaped downward beam mode. The results show that the electron beam energy distribution was symmetrically distributed along the center line. The surface energy distribution was related to the deflection angle and the energy peak parameter. When the energy peak parameter was within 0 to 1, the value becomes higher and higher, the larger the first and second energy peaks at point, the larger the difference between the two. When the parameter was 0, the maximum temperature difference at the sampling point was 1 065 K. The smaller the temperature difference between the longitudinal points of the scanning belt, the smaller the distance between the thermal cycle curves. At the same time, it can be seen from the heat source model that the energy peak has a greater impact on the beam diameter of the ring electron beam, and the maximum ring diameter can be up to 8 mm under the selected basic parameters. When the parameter was 1, the temperature curves of the sampling points are the closest, which indicates that the surface heat distribution under this parameter was uniform. It can be seen from the temperature field simulation diagram that the beam spot temperature varies greatly during the down and converging phases of the electron beam, while the temperature in the middle of the scan was relatively stable, and the temperature difference was basically stable within 20 K. The larger the energy peak parameter, the larger the radius of the high temperature area on the surface of 45 steel, and the maximum temperature will increase accordingly. After 45 steel was subjected to different energy peak coefficients, the width of the scanning zone and the sub-high temperature zone were different. Finally, based on the simulation parameters, the scanning electron beam micro-melting polishing experiment was carried out. It was found that the surface roughness of 45 steel was reduced under this scanning mode, and the surface showed a bright white scanning area relative to the substrate. The scanning area width increased with the increase of the energy peak parameter. This was in full agreement with the simulation results. The surface roughness after scanning electron beam treatment was as low as 0.36 μm relative to the substrate. In the end, the following conclusion can be drawn that the energy peak parameter has a great influence on the energy distribution of the ring electron beam. When= 1, the energy gradient of each position on the surface of the scanning area was the smallest, which was beneficial to the uniform energy distribution under the surface modification of the large area electron beam.

        scanning electron beams; energy; temperature field; heat source; surface modification

        V261.6

        A

        1001-3660(2022)07-0306-08

        10.16490/j.cnki.issn.1001-3660.2022.07.030

        2021–04–02;

        2021–11–20

        2021-04-02;

        2021-11-20

        廣西自然科學基金項目(2020GXNSFBA297079,2022GXNSFAA035585);國家自然科學基金資助項目(52165057,51665009);桂林市重點研發(fā)計劃(20211B032068)

        Guangxi Natural Science Foundation Project (2020GXNSFBA297079, 2022GXNSFAA035585); National Natural Science Foundation of China (52165057, 51665009); Guilin City Key Research and Development Plan (20211B032068)

        李新凱(1993—),男,博士研究生,主要研究方向為電子束表面改性。

        LI Xin-kai (1993-), Male, Ph. D. candidate, Research focus: electron beam surface modification.

        王喜社(1966—),女,高級實驗師,主要研究方向為數(shù)控加工與工藝優(yōu)化。

        WANG Xi-she (1966-), Female, Senior experimenter, Research focus: CNC machining and process optimization.

        李新凱, 王榮, 胡露瑤, 等. 掃描電子束能量峰值對表面溫度場的影響規(guī)律[J]. 表面技術(shù), 2022, 51(7): 306-313.

        LI Xin-kai, WANG Rong, HU Lu-yao, et al. The Influence of the Peak Energy of Scanning Electron Beam on the Surface Temperature Field[J]. Surface Technology, 2022, 51(7): 306-313.

        責任編輯:萬長清

        猜你喜歡
        電子束熱源峰值
        “四單”聯(lián)動打造適齡兒童隊前教育峰值體驗
        少先隊活動(2022年9期)2022-11-23 06:55:52
        橫流熱源塔換熱性能研究
        煤氣與熱力(2021年3期)2021-06-09 06:16:20
        基于PCI-1721電子束磁掃描焊接的軟件設(shè)計
        電子測試(2017年12期)2017-12-18 06:35:21
        基于啟發(fā)式動態(tài)規(guī)劃的冷熱源優(yōu)化控制
        電子制作(2017年19期)2017-02-02 07:08:31
        寬占空比峰值電流型準PWM/PFM混合控制
        基于峰值反饋的電流型PFM控制方法
        Mo-3Nb單晶管定向電子束焊接的組織及性能
        焊接(2016年7期)2016-02-27 13:05:04
        5A90鋁鋰合金電子束焊接接頭顯微分析
        焊接(2015年7期)2015-07-18 10:59:18
        聚變堆用CLF-1鋼電子束焊接缺陷分析及控制
        焊接(2015年7期)2015-07-18 10:59:18
        中部槽激光-MAG復合熱源打底焊焊接工藝研究
        焊接(2015年8期)2015-07-18 10:59:13
        狠狠的干性视频| 久久亚洲精品一区二区| 国产一区二区三区涩涩涩| 日本综合视频一区二区| 日本三级片在线观看| 极品少妇被猛的白浆直喷白浆| 97免费人妻在线视频| 欧美成人网视频| 亚洲一区在线二区三区| 国产自拍高清在线观看| 无码中文字幕免费一区二区三区| 欧美性videos高清精品| 亚洲欧美成人在线免费| av有码在线一区二区三区| 福利视频一区二区三区| 欧美变态另类刺激| 国产中老年妇女精品 | 91精品啪在线观看国产色| 亚州终合人妖一区二区三区| 日本精品少妇一区二区三区| 中文无码一区二区不卡αv| 99久久国产亚洲综合精品| 丰满少妇又爽又紧又丰满动态视频 | 久久综合给合综合久久| 玩弄放荡人妻一区二区三区| 91久久精品国产性色tv | 日本一道dvd在线中文字幕| 白白色日韩免费在线观看| 欧美a级毛欧美1级a大片| 伊在人天堂亚洲香蕉精品区| 国产AV无码一区精品天堂| 青青草免费在线视频导航| 9久久婷婷国产综合精品性色| 亚洲人成人无码www影院| 国产免费又色又爽又黄软件 | 一区二区三区日本大片| 日本视频一区二区这里只有精品| 国产成人精品一区二区20p| 亚洲精品国产av天美传媒| 亚洲精品理论电影在线观看| 东京热东京道日韩av|