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

        ?

        車削刀具幾何參數低碳優(yōu)化建模與求解

        2015-12-02 01:25:34劉志杰劉曉龍林成新孫德平
        計算機集成制造系統(tǒng) 2015年11期
        關鍵詞:切削液偏角切削力

        劉志杰,劉曉龍,林成新,孫德平

        (1.大連海事大學 交通運輸裝備與海洋工程學院,遼寧 大連 116026;2.大連海事大學 輪機工程學院,遼寧 大連 116026)

        0 引言

        當前,大力發(fā)展并實現低碳經濟,應對經濟衰退、環(huán)境惡化和氣候變化,已經成為世界各國的發(fā)展目標。國際能源署(International Energy Agency,IEA)的調查表明,制造業(yè)CO2排放量占全球排放量的36%[1]。麻省理工Gutowski教授研究表明:一臺主軸功率22kW 的數控機床運行一年消耗電能所產生的CO2和SO2排放量,分別相當于61 輛和248 輛SUV汽車的排放量[2]。目前,我國機床保有量已達800余萬臺,機床作為制造加工系統(tǒng)主體,其能耗大、能效低,低碳節(jié)能潛力巨大[3]。

        切削加工中,刀具幾何參數選擇影響機床加工質量、穩(wěn)定性、生產效率和能量消耗等。如Shinozuka等通過穩(wěn)態(tài)金屬切削機理的有限元分析方法模擬斷屑的形狀、溫度和屈服應力,獲得了最優(yōu)刀具切削角度[4]。近年來,低碳制造引起了企業(yè)和學者的關注[5]。例如,韓國Song等[6]、美國國家標準與技術研究所(National Institute of Standards and Technology,NIST)[7]提出將碳排放信息集成于物料清單(Bill of Material,BOM)中,分別建立了g-BOM(green house gas-BOM)和碳清單(Bill of Carbon,BoC)以評估產品制造碳排放;唐任仲等[2]研究了基于活動的產品車間制造過程碳排放量計算模型,用于產品加工工藝過程的優(yōu)化;孫良峰等[8]利用功能分解樹映射結構,研究構建了復雜裝備制造過程中的碳排放分層遞階模型;李先廣等[9]基于Petri網建立了機床制造過程的碳排放模型。在制造工藝優(yōu)化方面,李聰波等建立了以最少加工時間和最低碳排放為優(yōu)化目標的切削參數多目標優(yōu)化模型,并采用復合型法獲得最優(yōu)切削參數[10]。本文在分析車削加工過程中碳排放組成的基礎上,綜合考慮機床設備、切削溫度和加工質量等約束條件,建立了以碳排放為優(yōu)化目標的刀具幾何參數優(yōu)化模型,并利用自適應遺傳算法進行優(yōu)化求解。

        1 車削刀具幾何參數低碳優(yōu)化建模

        車削刀具幾何參數主要包括前角γ、后角α和主偏角κ、副偏角、刃傾角、刀尖圓弧半徑等。其中:刀具前角(Cutting Tool Rake Angle,CTRA)對切削力和切削溫度的影響較大,研究表明,刀具前角每改變1°,主切削力約改變1.5%;刀具主偏角(Cutting Tool Edge Angle,CTEA)影響工件表面粗糙度、切屑層形狀、切削刃的工作長度和負荷,主偏角的合理選取將改善刀具使用壽命和散熱條件;副偏角主要影響表面粗糙度和刀具強度;刀尖圓弧半徑和刃傾角對切削力的影響不大;后角主要影響刀具后面與工件磨損[11]。因此,本文選擇對切削力和切削溫度影響較大的刀具前角和主偏角為優(yōu)化變量,研究其優(yōu)化選擇對加工過程中碳排放的影響。

        1.1 車削加工過程碳排放建模

        機床車削加工過程碳排放是指加工過程因生產和消耗能源、物料等產生的直接和間接碳排放,包括能源碳(切削過程、切削液開啟等消耗能源產生的碳排放)、物料碳(刀具、切削液、工件材料等物料消耗產生的碳排放)等[5]。本文主要研究刀具幾何參數對碳排放的影響,輔助設備、工件材料、消耗材料后期處理階段產生的碳排放與刀具幾何參數的相關性較小,為簡化計算模型,本文不考慮。車削加工過程碳排放Cp主要包括車削消耗電能碳排放Ce、刀具碳排放Ct和切削液碳排放Cc,即

        (1)車削消耗電能碳排放Ce

        車削加工消耗電能引起的碳排放表示為

        式中:Fe為電能碳排放因子(kgCO2/kWh),Ee為車削加工過程消耗的電能。

        車削加工過程的電能消耗包括空載能耗和切削過程能耗。研究表明,機床在負載切削加工時會產生附加載荷損耗功率Pa[12]。由于同一臺機床在加工過程中以某一固定轉速穩(wěn)態(tài)運行且負載一定時,輸入功率、空載功率Pu、切削功率Pc、附加載荷損耗功率Pa的波動較小,車削加工過程能耗Ee可表示為[12]

        式中:Tp為加工工時;tm為切削加工時間;Pu,Pc,Pa和Tp分別計算如下:

        1)空載功率Pu

        機床空載功率與機床的自身情況、轉速n等相關,通常表示為[13]

        式中:Pu0為機床的最低空載功率;A1和A2為機床主軸轉速系數。

        2)切削功率Pc

        實際加工過程切深抗力Fp不消耗電能,進給力Ff所消耗的功率占總功率的1%~5%,可忽略不計,因此加工過程車削功率Pc可表示為[11]

        式中:vc為切削速度;切削力Fc=AcKs,Ac為切屑的截面積,Ks為單位切削力(N/mm2)[14]。

        式中:b為切屑截面寬;l為切屑截面長;h為切削深度;f為切削進給量。

        式中:ks11為額定單位切削力(N/mm2),z為曲線斜率;Kv為切削速度影響因子;Kγ為刀具影響因子,Kγ=(C-1.5γ)/100,對于鋼材料,C=109;Ka為刀具磨損影響因子;Kt為工具材料和金屬去除率影響因子。

        3)附加載荷損耗功率Pa

        附加載荷損耗功率Pa可表示為[12]

        式中bm為附加載荷損耗系數,在實際加工過程中常取0.15~0.25。

        4)加工過程工時Tp

        本文以車削加工單工步單工序的加工過程為例,建立工時Tp的數學模型,即[10]

        式中:tct為換刀一次所需時間;T為刀具耐用度;tot為使用同一把刀具時除換刀時間外的其他輔助時間;tm為車削加工時間。

        綜合切削要素v,f,a和其他因素對刀具壽命的影響規(guī)律,刀具壽命可表示為[11]

        式中:xT,yT為背吃刀量ap和進給量f對刀具壽命T的影響程度指數;CK為表示加工材料、加工形式、工件材料和進給量影響因素的修正系數;KT為與刀具前角和主偏角有關的修正系數。

        切削加工時間tm可表示為

        式中:lw為刀具一次走刀長度(mm);d為待加工工件直徑(mm);Δ為工件單邊加工余量(mm);C0為耐用度系數,與刀具、工件材料和切削條件有關,根據經驗通常取204。

        將式(13)代入式(3),可計算車削過程電能消耗Ee,進而由式(2)計算電能消耗產生的碳排放。

        (2)刀具碳排放Ct

        機械加工過程中,刀具碳排放包括刀具材料載能和刀具制備產生的碳排放,零件加工過程刀具碳排放Ct可表示為

        式中:Ft1和Ft2分別為刀具的材料碳排放因子和刀具制備碳排放因子;wt和Et分別為刀具質量和刀具制備過程消耗的能量;Tt為刀具壽命,Tt=(N+1)T,N為重磨次數。RAJEMI[15]研究表明,常用刀具單個刀片的平均質量為9.5g,結合全國電網電能平均碳排放因子0.067 47kgCO2/kWh[10],當計算只考慮刀具材料引起的碳排放而不考慮刀具制造過程能耗引起的碳排放時,刀具碳排放因子為75kg-CO2/kg;兩者都考慮時,刀具碳排放因子為104.6 kgCO2/kg。

        (3)切削液碳排放Cc

        切削液在使用過程中存在蒸發(fā)和泄漏等損耗現象。切削液的碳排放包括純礦物油制備和廢切削液處理引起的碳排放。生產中常用的水溶性(水基)切削液碳排放可表示為

        式中:F0為純礦物油排放因子;Fw為廢棄切削液處理碳排放因子;Cc0為初始切削油用量;Ac為附加切削油用量;δ為切削液濃度;Tc為切削液更換周期?;诘V物油的內含能值可得礦物油的碳排放因子為2.85kgCO2/L[10]。由于水基切削液濃度較低(3%~5%),其碳排放因子可取廢水處理的碳排放因子,即0.2kgCO2/L[16]。

        綜上所述,車削加工過程碳排放計算數學模型可表示為

        1.2 切削溫度建模

        車削加工切削溫度影響加工質量、加工效率、刀具使用壽命和切削液損耗量,進而影響碳排放。因此,切削溫度是選擇切削工藝參數時需考慮的重要因素之一。切削溫度主要包括兩部分[14,17]:①剪切面(即第一變形區(qū))的平均溫度,稱為剪切溫度,用θs表示;②刀—屑接觸面(即第二變形區(qū))的平均溫度,稱為摩擦溫度,用θf表示。切削溫度可表示為

        將切削面看成均勻熱源,剪切面單位時間內產生的熱量為Qs,其中一部分傳入工件,剩余熱量完全被切屑帶走。刀具導熱系數較小,傳入刀具的熱量忽略不計,則θs可表示為[14,17]

        式中:R1為傳入切屑的比例系數;Qs為剪切面單位時間內產生的熱量;ρ為工件材料的密度;c為工件的比熱;v為切削速度;ac為切削厚度;aw為切削寬度。

        對刀—屑之間的摩擦溫度采用Jaeger的動熱源方法進行計算,得[14,17]

        式中:ρ為工件材料的密度;c為工件的比熱;τ為剪切面上的剪切應力;As為剪切面上的截面積;φ為剪切角;τε為材料的剪切屈服極限;nθ為材料強化系數;K為導熱系數;Lf為刀具與切削區(qū)間的接觸長度;a0為切屑厚度;β為摩擦角。

        1.3 約束條件

        車削加工過程切削參數和刀具參數等工藝參數的選取受所選機床設備主軸轉速、進給量、最大切削力、最大切削功率和加工質量等條件的限制。對于刀具參數的選擇,本文考慮以下約束條件:

        (1)刀具前角約束 刀具前角影響加工質量、切削溫度和加工效率等,其取值范圍為-15°≤γ≤20°。

        (2)刀具主偏角約束 刀具主偏角影響刀具耐用度和切削溫度,其取值范圍為30°≤κr≤90°。

        (3)切削力約束 數控加工過程的最大進給抗力應小于機床進給機構的最大進給力,即vnfKFf≤Fmax。其中:Fmax為機構最大進給力;KFf,CFf,xf,yf和nf為與工件材料及加工條件相關的參數,可查閱機械工程手冊得到。

        (4)功率約束 機床加工功率不能大于規(guī)定的最大有效切削功率,即,其中η為機床功率系數。

        (5)加工質量約束 加工零件的表面粗糙度應滿足質量要求,即·Rmax。其中:rε為刀具刀尖圓弧半徑;Rmax為零件表面粗糙度要求的最大值;κ′為刀具副偏角。

        (6)切削溫度約束 加工過程切削溫度影響刀具使用壽命、工件加工質量、切削液蒸發(fā)速度等,切削溫度滿足θ≤θmax。其中:θ可由式(17)計算;θmax為加工過程允許的最大切削溫度,與刀具材料和被加工工件材料有關。切削溫度的限定影響因素復雜,本文僅從有利于工件加工質量和切削液消耗的角度給出切削溫度限定條件。

        1.4 車削刀具幾何參數低碳優(yōu)化模型

        綜合上述因素和模型,針對具體工程問題采用不同的性能指標,例如以最低碳排放或最低切削溫度或兩者同時滿足為優(yōu)化目標,建立刀具幾何參數低碳優(yōu)化模型。當以最低碳排放為優(yōu)化目標、以切削溫度等為約束條件時,建立的車削加工刀具幾何參數低碳優(yōu)化模型表示為:

        其中:式(22)為機床功率約束,式(23)為切削力約束,式(24)為表面加工質量約束,式(25)為刀具角度約束,式(26)為溫度約束。

        2 車削刀具幾何參數低碳模型優(yōu)化求解

        由于自適應遺傳算法能根據環(huán)境的變化過程發(fā)現環(huán)境規(guī)律和特性,具有較好的全局尋優(yōu)能力,本文采用該方法求解。針對本文優(yōu)化問題,對自適應遺傳算法交叉和變異概率進行改進,交叉概率Pc、變異概率Pm的表達式分別為:

        式中:f′為兩個交叉?zhèn)€體中的較大適應度值;f″為變異個體的適應度值;fmax為群體的最大適應度值;ˉf為適應度平均值。

        3 實例分析

        以如圖1所示的某數控機床(規(guī)格參數如表1)車削加工的軸類零件為例,分析本文模型和方法的應用。

        表1 數控機床規(guī)格參數

        3.1 模型參數確定

        工件和加工參數條件如下:工件為45 鋼棒,σb=0.637GPa,切削深度ap=1 mm,切削速度v=1.57m/s,進給量f=0.34mm/r,刀具為硬質合金車刀,刀具壽命60min,刃傾角λs=5°,刀具的刀尖圓弧半徑rε=0mm,加工表面粗糙度Ra要求不得超過6.4μm,加工時使用切削液。

        通過查閱相關資料,單位切削力和相關參數如表2所示,刀具壽命相關參數如表3所示,其他相關系數如表4和表5所示。值得說明的是,表3中KT是與刀具前角和主偏角有關的修正系數。根據相關資料,當刀具材料為硬質合金、高速鋼,工件材料為結構碳鋼及合金鋼、耐熱鋼、鑄鐵等,刀具耐用度為60min時,KT取值均為1.0[18],因此本算例中KT=1.0。

        表2 單位切削力及相關參數

        表3 刀具壽命及切削力系數

        表4 碳排放計算相關系數1

        表5 碳排放計算相關系數2

        通常,不同刀具在切削不同材料工件時,刀具允許切削溫度的范圍不同。考慮切削溫度對加工精度(熱變形)和切削液消耗等的影響,本文硬質合金刀具允許的最大切削溫度θmax=500 ℃。其他計算切削溫度的相關參數如表6和表7所示。

        表6 切削溫度相關系數1

        表7 切削溫度相關系數2

        3.2 優(yōu)化計算及結果分析

        采用交叉及變異概率如式(22)和式(23)所示的自適應遺傳算法,式中系數Pc1=0.6,Pc2=0.9,Pc3=0.2,Pm1=0.1,Pm2=0.001。迭代次數選取50,映射精度取10-3,利用MATLAB編程進行優(yōu)化計算。為對比分析考慮不同性能指標時的最優(yōu)刀具幾何參數,分別以低溫(以切削溫度最低為優(yōu)化目標)、低碳(以碳排放最小為優(yōu)化目標,不考慮切削溫度)和低碳限溫(以碳排放最小為優(yōu)化目標、以最高允許切削溫度θmax=500 ℃為約束條件)建立優(yōu)化模型,進行刀具參數的優(yōu)化選擇,其優(yōu)化結果如表8所示。

        表8 實例優(yōu)化結果

        如表8所示,以低溫為優(yōu)化目標時,最優(yōu)刀具前角為19.832 8°,主偏角為30.672 5°;以低碳為優(yōu)化目標時,最優(yōu)刀具前角為19.563 7°,主偏角為89.460 2°;低碳限溫時的最優(yōu)刀具前角為18.956 5°,主偏角為60.205 1°。這一結果表明:當刀具前角增大時,切削刃鋒利,切屑變形小,前刀面摩擦減小,產生的熱量減少,因此切削溫度隨前角的增大先降低;當前角過大時,由于刀具楔角變小,刀具散熱體積減少,切削溫度反而會升高;當主偏角增加時,在背吃刀量不變的條件下主切削刃的工作長度減小,散熱面積減小,散熱條件變差,因此切削溫度上升。在實際加工過程中,取低碳限溫時的最優(yōu)刀具前角20°、主偏角60°作為最優(yōu)刀具參數。優(yōu)化前后,刀具參數及其碳排放和車削溫度計算結果如表9所示。由表9可知切削溫度和碳排放分別降低了17%和12%,從而有助于實現低碳降溫加工。

        表9 優(yōu)化前后結果對比分析

        4 結束語

        本文以刀具前角和主偏角為優(yōu)化變量,考慮機床主軸、刀具前角、刀具主偏角、最大切削功率、最大切削力、切削溫度和加工質量等約束條件,建立了刀具幾何參數低碳優(yōu)化模型,并利用自適應遺傳算法進行了優(yōu)化計算。車削加工階梯軸實例優(yōu)化分析表明,優(yōu)化獲得的最優(yōu)刀具前角和主偏角,使車削過程碳排放和切削溫度分別比優(yōu)化前降低12%和17%。本文可為制造企業(yè)優(yōu)化選擇刀具幾何參數、降低碳排放提供理論方法支持。在建模中,多工序、多工步的加工工時模型,以及建立切削液、刀具消耗、加工質量、加工效率與切削溫度之間的關系模型是需進一步研究的問題。此外,除切削溫度外,工序成本等因素也是刀具幾何參數低碳優(yōu)化模型中需進一步考慮的因素。

        [1]International Energy Agency.Tracking industrial energy efficiency and CO2emissions[EB/OL].[2014-10-15].http://www.global-greenhouse-warming.com/industry-CO2-emissions.html.

        [2]TANG Renzhong,JIA Shun,LYU Jingxiang.Activity-based carbon-computing method for manufacturing process[J].Journal of Mechanical &Electrical Engineering,2011,28(6):639-644(in Chinese).[唐任仲,賈 順,呂景祥.基于活動的產品車間制造過程碳排量計算方法[J].機電工程,2011,28(6):639-644.]

        [3]LI Dongru.Prospect and development of green manufacturing technology in China[J].Mechanical Engineering Trends,2014(4):6-12(in Chinese).[李冬茹.中國綠色制造技術發(fā)展與展望[J].機械工程導報,2014(4):6-12.]

        [4]SHINOZUKA J,OBIKAWA T,SHIRAKASHI T.Chip breaking analysis from the viewpoint of the optimum cutting too geometry design[J].Journal of Materials Processing Technology,1996,62(4):345-351.

        [5]CAO Huajun,LI Hongcheng,DU Yanbin,et al.Current situation and development trend of low-carbon manufacturing[J].Aeronautical Manufacturing Technology,2012(9):26-31(in Chinese).[曹華軍,李洪丞,杜彥斌,等.低碳制造研究現狀、發(fā)展趨勢及挑戰(zhàn)[J].航空制造技術,2012(9):26-31.]

        [6]SONG J S,LEE K M.Development of a low-carbon product design system based on embedded GHG emissions[J].Resources,Conservation and Recycling,2010,54(9):547-556.

        [7]AMETA G,MANI M,RACHURI S.Carbon weight analysis for machining operation and allocation for redesign[J].International Journal of Sustainable Engineering,2009,2(4):241-251.

        [8]SUN Liangfeng,QIU Lemiao,ZHANG Shuyou,et al.Hierarchical model for carbon emissions of complex equipment for low-carbon design[J].Computer Integrated Manufacturing Systems,2012,18(11):2381-2390(in Chinese).[孫良峰,裘樂淼,張樹有,等.面向低碳設計的復雜裝備碳排放分層遞階模型[J].計算機集成制造系統(tǒng),2012,18(11):2381-2390.]

        [9]LI Xianguang,LI Congbo,LIU Fei,et al.Modeling and quantization methods for carbon emission in manufacturing processes of machine tools based on Petri net[J].Computer Integrated Manufacturing Systems,2012,18(12):2723-2735(in Chinese).[李先廣,李聰波,劉 飛,等.基于Petri網的機床制造過程碳排放建模與量化方法[J].計算機集成制造系統(tǒng),2012,18(12):2723-2735.]

        [10]LI Congbo,CUI Longguo,LIU Fei,et al.Multi-objective NC machining parameters optimization model for high efficiency and low carbon[J].Journal of Mechanical Engineering,2013,49(9):87-96(in Chinese).[李聰波,崔龍國,劉飛,等.面向高效低碳的數控加工參數多目標優(yōu)化模型[J].機械工程學報,2013,49(9):87-96.]

        [11]WU Wenge,XIN Zhijie.Principles of metal cutting and cutting tools[M].Beijing:National Defence Industry Press,2009(in Chinese).[武文革,辛志杰.金屬切削原理及刀具[M].北京:國防工業(yè)出版社,2009.]

        [12]XU Zongjun,LIU Fei,DAN Bin.Energy characteristics of maching systems and its application[M].Beijing:China Machine Press,1995(in Chinese).[徐宗俊,劉 飛,但 斌.機械加工系統(tǒng)能量特性及其應用[M].北京:機械工業(yè)出版社,1995.]

        [13]HU Shaohua,LIU Fei,HE Yan,et al.No-load energy parameter characteristics of computerized numerical control machine tool main transmission system[J].Computer Integrated Manufacturing Systems,2012,18(2):326-331(in Chinese).[胡韶華,劉 飛,何 彥,等.數控機床變頻主傳動系統(tǒng)的空載能量參數特性研究[J].計算機集成制造系統(tǒng),2012,18(2):326-331.]

        [14]SAGLAM H,UNSAEAR F,YALDIZ S.Investigation of the effect of rake angel and approaching angle on main cutting force and tool tip temperature[J].International Journal of Machine Tools &Manufacture,2006,46(2):1-10.

        [15]RAJEMI M F,MATIVENGA P T,ARAMCHROEN A.Sustainable machining:Selection of optimum turning conditions based on minimum energy considerations[J].Manufacturing Technology,2011,60(1):145-148.

        [16]China National Institute of Standardization.Enterprise greenhouse gas accounting and reporting[M].Beijing:China Quality Inspection Press,2011(in Chinese).[中國標準化研究院.企業(yè)溫室氣體核算與報告[M].北京:中國質檢出版社,2011.]

        [17]LI Huiyun.Investigation on cutting temperatures in metal cutting process[J].Journal of Hebei Institute of Mechano-electric Engineering,1996,3(2):35-39(in Chinese).[李惠云.關于切削溫度理論計算公式的探討[J].河北機電學院學報,1996,3(2):35-39.]

        [18]YANG Shuzi.Mechanical process engineer handbook[M].Beijing:China Machine Press,2011(in Chinese).[楊叔子.機械加工工藝師手冊[M].北京:機械工業(yè)出版社,2011.]

        猜你喜歡
        切削液偏角切削力
        軸承磨加工水基切削液的合理選擇
        哈爾濱軸承(2021年2期)2021-08-12 06:11:08
        基于PLC切削液自動勾兌機的設計
        翼吊長涵道發(fā)動機短艙內偏角優(yōu)化和機理研究
        2018全國Ⅱ卷選修3-4中偏角的解法探討
        30CrMnSiNi2A插銑加工的切削力研究
        圓刀片切削力計算方法
        歐姆表偶然誤差分析
        一種深孔鉆削切削力測量裝置
        集成對稱模糊數及有限元法的切削力預測
        山推公司切削液使用情況
        一本大道久久精品 东京热 | 丰满少妇高潮惨叫久久久| 久久久久亚洲精品无码网址蜜桃| 亚洲av无码一区二区三区网站| 国产精品久久久亚洲第一牛牛| 国产偷闻女邻居av在线观看| 日本伊人精品一区二区三区| 亚洲精品aa片在线观看国产| 久久精品国产亚洲5555| 五月激情在线观看视频| 欧美成人精品第一区| 99香蕉国产精品偷在线观看| 日韩精品国产自在欧美| 青青草好吊色在线视频| 精品无码av一区二区三区不卡| 大肉大捧一进一出好爽视色大师| 日本一区二区啪啪视频| 日本在线无乱码中文字幕| 97一期涩涩97片久久久久久久| 厕所极品偷拍一区二区三区视频| 国产交换精品一区二区三区| 国产精品无码一区二区在线观一 | 韩国免费一级a一片在线| 一区二区亚洲精品在线| 夜鲁很鲁在线视频| 好爽受不了了要高潮了av| av免费在线观看网站大全| 亚洲精品国产精品乱码在线观看| 中国a级毛片免费观看| 五月天无码| av成人资源在线播放| 亚洲高清在线天堂精品| 色一情一乱一伦一区二区三区日本| 国产亚洲高清不卡在线观看| 91国产熟女自拍视频| 丰满少妇a级毛片| 日韩A∨精品久久久久| 国产自产自现在线视频地址| 欧美亅性猛交内射| 精品福利视频一区二区三区| 日韩精品区欧美在线一区|