999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

TA7電火花加工實驗研究與工藝參數(shù)優(yōu)化

2017-08-08 03:01:25張立新史衛(wèi)朝
制造業(yè)自動化 2017年7期
關(guān)鍵詞:工藝優(yōu)化實驗

張立新,路 雄,史衛(wèi)朝

ZHANG Li-xin, LU Xiong, SHI Wei-chao

(西安理工大學(xué),西安 710082)

TA7電火花加工實驗研究與工藝參數(shù)優(yōu)化

Experimental research and the optimization of process parameters in the electrical discharge of the TA7

張立新,路 雄,史衛(wèi)朝

ZHANG Li-xin, LU Xiong, SHI Wei-chao

(西安理工大學(xué),西安 710082)

針對TA7鈦合金難切削以及加工表面質(zhì)量差的問題,在通用電火花機床上探究了空載電壓、峰值電流、脈沖寬度以及脈沖間隔對TA7鈦合金電火花加工中材料去除率與表面粗糙度的影響。設(shè)計中心復(fù)合實驗,在實驗數(shù)據(jù)的基礎(chǔ)上采用多線回歸技術(shù)建立材料去除率和表面粗糙度的二階預(yù)測模型。方差分析結(jié)果表明預(yù)測模型可以正確映射出TA7鈦合金電火花加工的工藝規(guī)律。以提高材料去除率降低表面粗糙度為目的建立工藝目標(biāo)優(yōu)化模型,根據(jù)實際加工條件對工藝參數(shù)的約束,設(shè)計粒子群算法求取優(yōu)化模型的解集,通過實際加工對優(yōu)化結(jié)果進(jìn)行驗證,結(jié)果表明該算法可以正確可靠的獲得多約束條件下的最優(yōu)加工參數(shù)。

電火花加工;TA7鈦合金;材料去除率;表面粗糙度;優(yōu)化

0 引言

TA7是一種中等強度的單向α型鈦合金,其成分為Ti-5Al-2.5Sn,由于TA7具有比強度高[1,2]、熱穩(wěn)定性好、抗腐蝕能力強、化學(xué)活性高等眾多優(yōu)點,因此被廣泛應(yīng)用與航空航天、兵器制造、醫(yī)療器械等領(lǐng)域。采用傳統(tǒng)的加工方式對TA7加工時,當(dāng)被加工件的硬度大于HB350時,刀具對工件的切割較為困難,而當(dāng)被加工工件的硬度小于HB300時,由于材料較軟會出現(xiàn)粘刀現(xiàn)象,因此TA7屬于機械難加工材料。而TA7的機械難加工性不僅由其硬度決定,同時還受到TA7本身物理、化學(xué)以及力學(xué)間的綜合性能的影響。電火花成形加工[3]是一種非傳統(tǒng)的特種加工技術(shù)[4],其加工對象必須具有導(dǎo)電能力,由于電火花加工無宏觀加工應(yīng)力,且具有加工能量密度高、加工不受工件材料硬度限制等優(yōu)點,因此適用于加工TA7或其它硬質(zhì)合金材料。對于金屬的電火花放電加工過程及加工工藝相對比較成熟[5,6]。加工效率和表面質(zhì)量是電火花加工中極為重要的工藝目標(biāo),而這兩項工藝目標(biāo)在很大程度上受到電源參數(shù)[7,8]和機械參數(shù)的影響。文獻(xiàn)[9]采用誤差反向傳遞神經(jīng)網(wǎng)絡(luò)模型研究加工的電流、脈沖寬度和脈沖間隔作對材料去除率和工具損耗的影響。R.Rajesh等人采用多元回歸和遺傳算法[10]建立了工作電流、工作電壓、油壓、脈沖寬度、脈沖間隔和材料去除率以及表面粗糙度的灰色關(guān)聯(lián)模型。M.R. Shabgard等人[11]采用模糊控方法分析了電流、脈沖寬度和刀具的超聲振動對材料的去除率和表面粗糙度的影響。然而對于合金材料尤其是鈦合金材料的電火花加工研究的相對較少,文獻(xiàn)[12]在電火花超聲輔助加工WC-Co合金過程中使用模糊方法建立了電極損耗預(yù)測模型,進(jìn)一步優(yōu)化了峰值電流、脈沖寬度等工藝參數(shù),有效的降低了電極損耗。P Sengottuvel等人[13]在正交實驗的基礎(chǔ)上采用模糊邏輯方法對電火花放電加工鉻鎳鐵合金的電極損耗進(jìn)行了實驗建模。

由于現(xiàn)有電火花加工設(shè)備中沒有針對鈦合金或其它合金材料的加工工藝參數(shù),而且鈦合金的電火花加工研究相對較少。在通用電火花成形機床上研究空載電壓、峰值電流、脈沖寬度以及脈沖間隔對TA7材料去除率和表面粗糙度的影響,采用響應(yīng)曲面法(Response surface methodology,RSM)在四因素五水平實驗數(shù)據(jù)的基礎(chǔ)上,通過多線回歸方法建立材料去除率、表面粗糙度與相關(guān)工藝參數(shù)的數(shù)學(xué)模型,進(jìn)一步采用粒子群優(yōu)化算法(particleswarm optimization,PSO)對工藝目標(biāo)進(jìn)行優(yōu)化,獲得加工TA7鈦合金電火花加工的最優(yōu)工藝參數(shù)。

1 所用設(shè)備及實驗加工對象

加工實驗平臺采用ACT SA20電火花成形機床,機床主要結(jié)構(gòu)參數(shù)和電源參數(shù)如表1所示。

表1 P機床規(guī)格及電源參數(shù)

電火花成形加工中工具電極的性能對工藝目標(biāo)有一定的影響,鑒于紫銅不但具備良好的導(dǎo)電性、導(dǎo)熱性,而且還同時兼顧了熔點高、磨損小、穩(wěn)定性好、易于加工等特點,因此本文選擇紫銅作為電極材料,紫銅電極直徑為Φ10、長度為150mm。

工件質(zhì)量選擇ESJ-B精密電子分析天平進(jìn)行測量,其最小稱重為0.05g,最大稱重為220g,測量精度為0.0001g,重復(fù)性誤差小于±0.0001g。

加工后的TA7鈦合金表面粗糙度采用Leica DCM 3D白光干涉儀對的進(jìn)行測量,其測量精度可達(dá)0.01μm。

實驗加工對象為長方體的TA7鈦合金,尺寸為80×60×50(mm),TA7鈦合金材料的物理屬性如表2所示。

2 試驗設(shè)計

響應(yīng)曲面法是研究多輸入多輸出問題的一種非線性回歸方法[14,15]。它在實驗的基礎(chǔ)上結(jié)合數(shù)學(xué)、統(tǒng)計學(xué)原理,對影響因素與響應(yīng)輸出之間的函數(shù)關(guān)系進(jìn)行描述[16,17]。本文研究TA7鈦合金電火花成形加工中材料去除率、表面粗糙度與工藝參數(shù)之間的量化關(guān)系,因此響應(yīng)輸出為材料去除率(MRR)和表面粗糙度(Ra)。MRR計算公式為:

通過測量加工前后工件的質(zhì)量差,再除以工件的密度得到工件蝕除的體積V(mm3),再除以加工時間T(min)得到MRR。Ra采用Leica DCM 3D白光干涉儀進(jìn)行測量。

電火花成形加工的加工效率和工件表面質(zhì)量直接由極間放電能量的大小決定,而極間放電能量跟極間電壓、放電電流以及火花持續(xù)時間有著非常緊密的聯(lián)系,因此將極間空載電壓(U)、峰值電流(I)、脈沖寬度(Tw)以及脈沖間隔(Ti)作為影響MRR和Ra的主要因素。通過大量單因素實驗研究規(guī)律發(fā)現(xiàn)U低于80V放電通道不易被擊穿,U高于120V會因放電能量過大而使TA7鈦合金表面產(chǎn)生微裂紋;I在2A以下不能產(chǎn)生正常的火花放電,隨著I的增大火花放電處于穩(wěn)定狀態(tài),但當(dāng)I超過6A時TA7鈦合金表面有燒傷的痕跡,并出現(xiàn)翻邊現(xiàn)象,使工件表面質(zhì)量急劇惡化。正常的花火放電有一定的擊穿延時存在,擊穿延時時間一般在幾微秒到十幾微秒之間,因此Tw的最小值取20μs,而當(dāng)Tw大于420μs時,用于蝕除工件的能量有一部分通過熱交換的方式消耗掉了,造成了蝕除飽和現(xiàn)象,在這種情況下繼續(xù)增大Tw,工件的蝕除量也不會繼續(xù)增加。Ti較小時不能充分消電離,容易造成二次放電及拉弧現(xiàn)象。Ti較大時雖能有效的對極間進(jìn)行消電離,但是影響加工效率。針對所選4個影響因素,每個影響因素取5個水平,這樣細(xì)分因素可使實驗結(jié)果更加精確,中心復(fù)合實驗的因素及水平如表3所示。

實驗中采用電火花專用油作為極間絕緣工作液介質(zhì),實驗過程中其他加工條件都保持不變。

表2 TA7鈦合金參數(shù)

表3 因素水平

3 實驗建模與結(jié)果分析

3.1 實驗建模

根據(jù)表3中的影響因素與因素水平,設(shè)計四因素五水平中心組合實驗。實驗共計30組,按照中心組合實驗中參數(shù)在電火花成形機床上對TA7進(jìn)行加工,由于放電加工存在很大的隨機性,因此每組實驗重復(fù)進(jìn)行3次,分別求取MRR和Ra的平均值,中心復(fù)合實驗參數(shù)與實驗結(jié)果如表4所示。

通過對實驗數(shù)據(jù)的多元回歸擬合處理可以得到響應(yīng)與影響因素的函數(shù)關(guān)系模型,根據(jù)模型階次的不同,則有不同的表達(dá)形式。若采用1階模型很難對表4中的實驗數(shù)據(jù)進(jìn)行精確描述,雖然3階模型可以對實驗數(shù)據(jù)進(jìn)行精確描述,但是各影響因素的三次方對模型沒有任何貢獻(xiàn),因此選擇二階模型對實驗數(shù)據(jù)進(jìn)行描述,二階模型可描述為:

表4 四因素五水平實驗結(jié)果

式中:Y為響應(yīng)值;Xi為影響因素;β0為常數(shù)項;βi為一次項系數(shù);βii為二次項系數(shù);βij為交互項系數(shù);ε為誤差(其中包括實驗誤差和擬合不足引起的誤差)。

采用Design-Expert對表4中的實驗數(shù)據(jù)進(jìn)行回歸擬合,得到TA7鈦合金電火花成形加工MRR和Ra的二階數(shù)學(xué)模型:

由于所建立MRR和Ra的二階數(shù)學(xué)模型是在實驗的基礎(chǔ)上通過對實驗數(shù)據(jù)進(jìn)行回歸擬合得到的,因此模型中參數(shù)的取值應(yīng)在實驗所選參數(shù)范圍內(nèi),模型參數(shù)取值范圍如式(5)所示。

3.2 結(jié)果分析

方差分析是用來驗證模型擬合實驗數(shù)據(jù)的可靠性及其顯著性。其本質(zhì)是將實驗因素(變量和水平)引起的響應(yīng)值與試驗隨機誤差造成響應(yīng)值的變化進(jìn)行比較,做出因素對響應(yīng)值的影響是否顯著的判斷。實驗數(shù)據(jù)方差分析結(jié)果如表5所示,MRR和Ra模型的P值都小于0.0001,失擬項均大于0.05且F值均大于F0.01(9,5)=10.16,則表示所建立MRR和Ra的模型整體極其顯著。各影響因素所對應(yīng)的P值若小于0.0001,則說明該影響因素對模型會產(chǎn)生極其顯著的影響。由于表5中各因素所對應(yīng)的P值均小于0.0001,意味這四個因素對模型都產(chǎn)生了極其顯著的影響,從而進(jìn)一步驗證了實驗設(shè)計中影響因素選擇的正確性。

表5 方差分析結(jié)果

通過影響因素F的值大小可進(jìn)一步判斷各因素對模型影響顯著的次序。分析表5中各因素F值的大小可知,對MRR和Ra影響顯著的順序均為:I、U、Tw、Ti。由MRR和Ra模型的回歸系數(shù)可以看出影響顯著的次序也為:I、U、Tw、Ti。

利用復(fù)相關(guān)系數(shù)R2的值對模型做進(jìn)一步分析,MRR模型的復(fù)相關(guān)系數(shù)R2=0.9854,Ra模型的復(fù)相關(guān)系數(shù)R2=0.9852,說明模型值和實驗值有著高度的相關(guān)性,而且模型具有良好的擬合程度。MRR模型的校正系數(shù)AdjR2=0.9718,說明該模僅有總變異的2.82%不能用此模型來解釋。Ra模型的校正系數(shù)AdjR2=0.9715,說明該模型僅有總變異的2.85%不能用此模型來解釋。模型的殘差正態(tài)分布圖可以檢驗?zāi)P偷膬?yōu)劣,圖1為MRR模型殘差正態(tài)積累分布圖,圖2為Ra模型殘差正態(tài)積累分布圖,從圖中可以看出,殘差都在正態(tài)分布直線附近,沒有異常點,所以可以判斷殘差呈正態(tài)分布,模型具有良好的適應(yīng)性。MRR模型的信噪比Adeq precision=31.578,Ra模型的信噪比Adeq precision=32.814,同樣也表明模型有良好的適應(yīng)性。

圖1 材料去除率模型殘差正態(tài)分布圖

圖2 表面粗糙度模型殘差正態(tài)分布圖

4 粒子群(PSO)優(yōu)化算法

TA7電火花成形加工中MRR和Ra是兩個非常重要但又相互矛盾的工藝目標(biāo)。由于目前電火花成形機床中沒有針對TA7鈦合金的加工工藝參數(shù),多目標(biāo)優(yōu)化就是在同時兼顧這兩個工藝目標(biāo)的前提下尋求合適的加工參數(shù),為后續(xù)加工提供參數(shù)可選依據(jù)。

粒子群算法[18,19](particleswarm optimization,PSO)是被Kennedy和Eberhart提出的,該算法基于Hepper的模擬鳥群(魚群)模型進(jìn)行修正,從而得到了粒子群優(yōu)化算法。針對MRR和Ra這兩個工藝目標(biāo)所建立的多目標(biāo)優(yōu)化模型為:

4.1 PSO算法流程

對于多目標(biāo)優(yōu)化問題,通常存在一個解集。就目標(biāo)函數(shù)而言,這些解之間是無法比較優(yōu)劣的,其特點是:無法在改進(jìn)任何目標(biāo)函數(shù)的同時不削弱至少一個其他目標(biāo)函數(shù),這些解稱作Pareto最優(yōu)解或非支配解。求解多目標(biāo)優(yōu)化問題的主要任務(wù)是:毫無偏好地找到盡可能多的具有代表性的符合要求的Pareto最優(yōu)解,在計算得到均勻分布的Pareto最優(yōu)解之后,根據(jù)設(shè)計要求和工程實際經(jīng)驗,從中客觀地選擇最滿意的優(yōu)化結(jié)果。本文采用基于Pareto解的PSO算法進(jìn)行求解多目標(biāo)優(yōu)化問題,同遺傳算法類似,是一種基于群體的迭代,通過粒子在解空間追隨最優(yōu)的粒子進(jìn)行求解。該算法的流程圖如圖3所示。

圖3 PSO算法流程圖

其算法實現(xiàn)步驟為:

步驟1:初始化迭代次數(shù),設(shè)置交叉概率、變異概率及其種群的規(guī)模,隨機初始化種群P;

步驟2:執(zhí)行種群進(jìn)化操作,對種群P進(jìn)行選擇、交叉和變異,得到新的進(jìn)化種群R;

步驟3:構(gòu)造種群PUR的非支配NDSet。當(dāng)NDSet的大于或小于N時,需要對NDSet按照策略進(jìn)行調(diào)整;

步驟4:判別終止條件,如果滿足終止條件,算法結(jié)束,否則將NDSet中的個體復(fù)制到種群P中,作為新的種群P,重復(fù)步驟2。

4.2 PSO算法優(yōu)化結(jié)果及實驗驗證

針對所建立TA7鈦合金電火花成形加工的多目標(biāo)優(yōu)化模型,本文采用粒子群算法(PSO)對優(yōu)化模型進(jìn)行求解。粒子群規(guī)模為500,精英集規(guī)模為200,最大迭代次數(shù)為200,加速度常量c1,c2為1.05,慣性權(quán)重ω為0.44。對優(yōu)化模型進(jìn)行求解,圖4為第10代Pareto前沿,圖5為第100代Pareto前沿,其中橫坐標(biāo)為MRR,縱坐標(biāo)為Ra。表6為20組由粒子群算法(PSO)優(yōu)化得到的電火花成形加工多目標(biāo)優(yōu)化問題的Pareto最優(yōu)解。其中每組解的兩個目標(biāo)函數(shù)值均對應(yīng)圖5中Pareto前沿上的一個點。

為了進(jìn)一步證明優(yōu)化結(jié)果的正確性以及可靠性,采用優(yōu)化所得參數(shù)在電火花成形機床上對優(yōu)化結(jié)果進(jìn)行驗證。對表6中其中5組優(yōu)化結(jié)果進(jìn)行實驗驗證,為了盡可能消除加工過程中隨機性引起的誤差,每組實驗重復(fù)進(jìn)行3次求取MRR和Ra的平均值。驗證結(jié)果如表7所示。

圖4 第10代的Pareto前沿

圖5 第100代的Pareto前沿

表6 粒子群算法優(yōu)化結(jié)果

表7 實例驗證結(jié)果

5 結(jié)論

以TA7鈦合金為實驗研究對象,在通用電火花成形機床上采用四因素五水平實驗研究空載電壓、峰值電流、脈沖寬度以及脈沖間隔對TA7電火花加工MRR和Ra的影響,利用響應(yīng)曲面法建立MRR和Ra的數(shù)學(xué)模型,進(jìn)一步設(shè)計粒子群算法對所建立的工藝目標(biāo)優(yōu)化模型進(jìn)行優(yōu)化求解,并進(jìn)行了實驗驗證,給出以下結(jié)論:

1)方差分析中峰值電流所對應(yīng)的F值比其他影響因素所對應(yīng)的F都大,因此各影響因素中峰值電流對MRR和Ra產(chǎn)生的影響最大。

2)MRR和Ra模型的P值、F值及失擬項均滿足方差分析結(jié)果的要求,因此所建立的二階模型整體極其顯著。二階模型的復(fù)相關(guān)系數(shù)、校正系數(shù)、信噪比以及殘差分布均在正常范圍內(nèi),表明模型具有較好擬合程度和適應(yīng)性。

3)建立了工藝目標(biāo)優(yōu)化模型,結(jié)合實際加工條件采用粒子群算法對所建立的模型進(jìn)行了優(yōu)化求解,為保證優(yōu)化結(jié)果的正確性以及可靠性,選擇優(yōu)化解中的加工參數(shù)進(jìn)行了實驗驗證,驗證結(jié)果表明優(yōu)化結(jié)果具有很好的可重現(xiàn)性。

[1] 龐洪,張海龍,王希哲,等.包覆疊軋TA7鈦合金薄板的組織和力學(xué)性能[J].中國有色金屬學(xué)報,2010,20(z1):66-69.

[2] 東赟鵬,于秋穎,方爽,等.TA7鈦合金高溫流變行為研究[J].航空材料學(xué)報,2015,35(1):13-19.

[3] 陳日,郭鐘寧,劉江文,等.電火花單脈沖放電加工的電熱模型與實驗研究[J].制造業(yè)自動化,2014,36(7):74-76.

[4] V Garbajs,Janez Peklenik. Statistical Model for an Adaptive Control of EDM-Process[J].CIRP Annals-Manufacturing Technology,1985,34(1):499-502.

[5] Y. S. Liao,M. P.Cheng,K. W. Liao.An on-line pulse tRains analysis system of the wire-EDM process[J].Journal of Materials Processing Technology,2009,209(9):4417-4422.

[6] Jean Bosco Byiringiro,Bernard W Ikua,George N Nyakoe.Fuzzy logic based controller for micro-electro discharge machining servo systems[A].IEEE[C].2009:1-6.

[7] Yih-fong Tzeng,F(xiàn)u-chen Chen. Multi-objective optimiRation of high-speed electrical discharge machining process using a Taguchi fuzzy-based approach[J].Materials & design,2007,28(4):1159-1168.

[8] GrzegorzSkrabala,Maria ZybuRa-SkRabalak,Adam Ruszaj.Building of rules base for fuzzy-logic control of the ECDM process[J].Journal of materials processing technology,2004,149(1):530-535.

[9] Debabrata Mandal,Surjya K. Pal,Partha Raha. Modeling of electrical discharge machining process using back propagation neuRal network and multi-objective optimization using nondominating sorting genetic algorithm-II[J].Journal of Materials Processing Technology,2007,186(1-3):154-162.

[10] R. Rajesh,M. Dev Anand. The Optimization of the Electro-Discharge Machining Process using Response Surface Methodology and Genetic Algorithms[J].Procedia Engineering,2012,38:3941-3950.

[11] M. R. Shabgard,M.A. Badamchizadeh,G. Ranjbary,etc. Fuzzy approach to select machining paRameters in electrical discharge machining(EDM) and ultRasonic-assisted EDM processes[J].Journal of Manufacturing Systems,2013,32(1):32-39.

[12] M. R. Shabgard,M.A.Badamchizadeh,G. Ranjbary,etc. Fuzzy approach to select machining paRameters in electrical discharge machining(EDM)and ultRasonic-assisted EDM processes[J].Journal of Manufacturing Systems,2013,32(1):32-39.

[13] P Sengottuvel,S Ratishkumar,D DinakaRan. Optimization of Multiple Characteristics of EDM Parameters Based on Desirability Approach and Fuzzy Modeling[J].Procedia Engineering,2013,64(4):1069-1078.

[14] Sameh S Habib. Study of the parameters in electrical discharge machining through response surface methodology approach[J].Applied Mathematical Modelling,2009,33(12):4397-4407.

[15] Sameh S Habib. Study of the parameters in electrical discharge machining through response surface methodology approach[J].Applied Mathematical Modelling,2009,33(12):4397-4407.

[16] 李淑娟,萬波,胡超.基于RSM的SiC單晶片切割過程多目標(biāo)優(yōu)化[J].機械科學(xué)與技術(shù),2014,33(3):452-459.

[17] 韓興國,祖曉琳.切削加工表面粗糙度的預(yù)測模型與參數(shù)優(yōu)化[J].制造業(yè)自動化,2014,36(10):39-42.

[18] Xie X, Zhang W,Yang L. Particle swarm optimization[J].Control and Decision,2003,18:129-134.

[19] Dorigo M,de Oca M A M, Engelbrecht A.Particle swarm optimization[J].Scholarpedia,2008,3(11):1486-1491.

TG661

:A

1009-0134(2017)07-0054-07

2017-02-22

張立新(1968 -),男,陜西人,本科,研究方向為特種加工。

猜你喜歡
工藝優(yōu)化實驗
記一次有趣的實驗
超限高層建筑結(jié)構(gòu)設(shè)計與優(yōu)化思考
民用建筑防煙排煙設(shè)計優(yōu)化探討
關(guān)于優(yōu)化消防安全告知承諾的一些思考
一道優(yōu)化題的幾何解法
轉(zhuǎn)爐高效復(fù)合吹煉工藝的開發(fā)與應(yīng)用
山東冶金(2019年6期)2020-01-06 07:45:54
做個怪怪長實驗
5-氯-1-茚酮合成工藝改進(jìn)
NO與NO2相互轉(zhuǎn)化實驗的改進(jìn)
實踐十號上的19項實驗
太空探索(2016年5期)2016-07-12 15:17:55
主站蜘蛛池模板: 国产亚洲欧美在线专区| 亚洲成人黄色在线| 国产另类乱子伦精品免费女| 欧美不卡视频在线观看| 国产一区二区三区精品久久呦| 国产精品吹潮在线观看中文| 国产在线精品美女观看| 91精品亚洲| 91美女在线| 在线观看国产黄色| 久久亚洲高清国产| 制服丝袜一区二区三区在线| 欧美另类视频一区二区三区| 99久久国产精品无码| 中文字幕 91| 在线视频亚洲色图| 2021最新国产精品网站| 18禁高潮出水呻吟娇喘蜜芽| 国产成人无码Av在线播放无广告| 国产成人无码AV在线播放动漫 | 国产精品尤物铁牛tv| 欧美在线伊人| 国产91丝袜在线播放动漫| 日韩欧美中文| 亚洲制服丝袜第一页| 91精品网站| 四虎精品黑人视频| 亚洲日韩每日更新| 欧美一区中文字幕| 欧美一级高清片欧美国产欧美| 日韩欧美国产三级| 毛片免费视频| 亚洲女同欧美在线| 亚洲成a人片| 日本三级黄在线观看| 1769国产精品免费视频| 国产喷水视频| 免费国产福利| 国产成人综合亚洲欧美在| 沈阳少妇高潮在线| 欧美日韩专区| 久久精品国产电影| 亚洲成人精品久久| 九九九精品成人免费视频7| 国产亚洲欧美日韩在线一区| 国产美女在线观看| 丰满人妻中出白浆| 色婷婷国产精品视频| 久久人人爽人人爽人人片aV东京热| 热久久国产| 永久毛片在线播| 91美女视频在线| 99久久免费精品特色大片| 美女免费精品高清毛片在线视| 亚洲成A人V欧美综合天堂| 91精品国产自产在线老师啪l| 国产免费好大好硬视频| 亚洲中文字幕无码mv| 97人人模人人爽人人喊小说| 欧美精品啪啪| 欧美区国产区| 久久五月天国产自| 一级毛片不卡片免费观看| 成人免费午间影院在线观看| 四虎亚洲国产成人久久精品| 欧美色99| av尤物免费在线观看| 日本午夜影院| 国产在线97| 自拍偷拍一区| 就去吻亚洲精品国产欧美| 男女男精品视频| 2018日日摸夜夜添狠狠躁| 亚洲日韩精品综合在线一区二区| 欧美日韩激情在线| 韩日午夜在线资源一区二区| 日韩精品无码免费一区二区三区| 日韩欧美中文亚洲高清在线| 999国产精品永久免费视频精品久久| 久久久久久久久久国产精品| 情侣午夜国产在线一区无码| 国产精品网址你懂的|