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

    82B小方坯連鑄提高拉速的工藝優(yōu)化實(shí)踐

    2024-06-29 09:38:51赫英利羅洪彥裘文李天坤陽(yáng)劍
    冶金信息導(dǎo)刊 2024年3期

    赫英利 羅洪彥 裘文 李天坤 陽(yáng)劍

    摘 要:為提高82B生產(chǎn)效率,在北營(yíng)煉鋼廠2號(hào)小方坯鑄機(jī)上開(kāi)展提高拉速的生產(chǎn)試驗(yàn),特別針對(duì)鑄坯芯部質(zhì)量進(jìn)行多工藝的綜合優(yōu)化,以實(shí)現(xiàn)高效高質(zhì)生產(chǎn)。首先,建立該鑄機(jī)的凝固傳熱數(shù)學(xué)模型,并結(jié)合比色測(cè)溫方法對(duì)模型中待定參數(shù)進(jìn)行修正;針對(duì)中心偏析和縮孔質(zhì)量問(wèn)題,結(jié)合修正后凝固傳熱模型和相應(yīng)冶金準(zhǔn)則對(duì)拉速、二次冷卻、電磁攪拌等參數(shù)進(jìn)行了優(yōu)化計(jì)算;最后通過(guò)生產(chǎn)試驗(yàn)比較,確定工藝優(yōu)化的方向和效果。最終優(yōu)化確定拉速2.25 m/min,比水量為1.7 L/Kg,結(jié)晶器電磁攪拌參數(shù)為260 A/3Hz,末端電磁攪拌參數(shù)320 A/8Hz。工藝實(shí)踐及分析表明,優(yōu)化工藝后鑄坯低倍中心碳偏析、縮孔均有較明顯的改善。

    關(guān)鍵詞:小方坯;芯部質(zhì)量;凝固傳熱模型;二冷;電磁攪拌

    TOPTIMIZATION AND PRACTICE OF BILLET CONTINUOUS CASTING PROCESS FOR IMPROVING CASTING SPEED OF 82B STEEL

    He Yingli1? ? Luo Hongyan1? ? Qiu Wen1? ? Li Tiankun1? ? Yang Jian2

    (1.Bensteel Group Coperation Limited Benxi 117017, China;2. School of Information Science and Engineering, Northeastern University Shenyang 110819,China.)

    Abstract:In order to improve the production efficiency of 82B, production tests were carried out on the No.2 billet continuous casting machine at Beiying Steelmaking Plant to increase the casting speed, with a special focus on the comprehensive optimization of the quality of central part, to achieve both high-efficiency and high-quality. Firstly, a heat transfer and solidification model has been established for the continuous casting billets, and then the undetermined parameters in the model were calibrated using a colorimetric pyrometer. Secondly, in response to the issues of central segregation and shrinkage quality, optimization calculations were carried out on parameters such as casting speed, secondary cooling, and electromagnetic stirring by combining both the calibrated solidification heat transfer model and corresponding metallurgical criteria. Finally, by comparing production tests, the direction and effectiveness of process optimization are determined. The final optimization results indicates that it is proper with casting speed of 2.25 m/min, a specific water flow of 1.7 L/Kg, electromagnetic stirring parameters of 260 A/3Hz for the mold electromagnetic stirring, and 320 A/8Hz for the final electromagnetic stirring. Process practice and analysis both indicate that the optimization of the process has significantly improved the center carbon macro-segregation and shrinkage of the billet products.

    Key words:billets; macro-segregation; solidification and heat transfer model; secondary cooling; electromagnetic stirring

    0? ? 前? ? 言

    鑄坯芯部質(zhì)量問(wèn)題(包括中心偏析、縮孔及疏松)是生產(chǎn)高碳硬線鋼所面臨的普遍而突出的難題[1],容易造成后續(xù)拉拔過(guò)程中的杯錐狀斷裂[2],并引起相應(yīng)的質(zhì)量異議和生產(chǎn)加工損失。同時(shí),該問(wèn)題又難以完全克服。實(shí)際上,由于選分結(jié)晶的存在和對(duì)流、擴(kuò)散對(duì)溶質(zhì)的傳輸作用,連鑄凝固過(guò)程中宏觀偏析必然形成,且在后序軋制過(guò)程中也難以消除。特別是對(duì)于過(guò)共析高碳鋼硬線品種如82B,嚴(yán)重的中心碳偏析在軋制過(guò)程中還導(dǎo)致晶界處網(wǎng)狀滲碳體的析出,造成后序拉拔加工過(guò)程中產(chǎn)生斷裂而影響加工過(guò)程或最終產(chǎn)品的質(zhì)量。由于高碳鋼兩相區(qū)較寬,凝固收縮較大,也較容易形成中心疏松和縮孔缺陷。因此,生產(chǎn)82B等硬線品種鋼時(shí),有效改善其芯部質(zhì)量對(duì)于滿足供貨質(zhì)量要求至關(guān)重要。

    在設(shè)備穩(wěn)定的條件下,連鑄工藝的優(yōu)化控制是改善其鑄坯質(zhì)量的主要手段和途徑[3]。其中,對(duì)于小方坯連鑄生產(chǎn),拉速、二冷和電磁攪拌是影響連鑄坯芯部質(zhì)量的核心要素。其中,拉速既決定生產(chǎn)效率又影響鑄坯質(zhì)量,合適的恒穩(wěn)拉速是改善鑄坯質(zhì)量的首要條件。二冷則通過(guò)影響溫度分布進(jìn)而影響凝固過(guò)程和凝固組織的形成,對(duì)中心偏析的形成有顯著的影響,對(duì)于小斷面,強(qiáng)冷條件下有利于坯殼與心部的同步收縮,同時(shí)較大的冷卻速率能減小枝晶間距,有利于縮孔和偏析的形成。電磁攪拌則電磁感應(yīng)產(chǎn)生電磁力驅(qū)動(dòng)鋼水做旋轉(zhuǎn)運(yùn)動(dòng)改變鋼水流動(dòng)狀態(tài),改變傳熱、傳質(zhì)過(guò)程,從而影響凝固過(guò)程,細(xì)化和提高等軸晶率,改善芯部質(zhì)量。合理優(yōu)化拉速、二冷和電磁攪拌是提高鑄坯質(zhì)量的關(guān)鍵。

    本文針對(duì)82B生產(chǎn),基于凝固傳熱模型的開(kāi)發(fā),通過(guò)有效結(jié)合工藝優(yōu)化計(jì)算和生產(chǎn)實(shí)踐,有效地改善了連鑄小方坯的芯部質(zhì)量。

    1? ? 凝固傳熱數(shù)學(xué)模型的建立與修正

    1.1? ? 模型的建立

    凝固傳熱模型是進(jìn)行連鑄工藝優(yōu)化控制的核心和基礎(chǔ)[4-5]。而由于連鑄過(guò)程涉及對(duì)流、凝固相變等復(fù)雜現(xiàn)象,模型建立和求解均較為困難,通常進(jìn)行必要的簡(jiǎn)化或轉(zhuǎn)化。

    1)將連鑄過(guò)程中的對(duì)流換熱等效為導(dǎo)熱。

    2)將凝固過(guò)程中潛熱釋放過(guò)程等效為比熱的增加。采用跟隨鑄坯移動(dòng)的隨動(dòng)坐標(biāo)系,建立連鑄坯三維凝固傳熱模型:

    ?T? ? ? ??? ? ? ? ? ?T? ? ? ? ?? ? ? ? ? ?T? ? ? ? ?? ? ? ? ? ?T

    ρceff ——=——(keff——)+——(keff——)+——(keff——)

    ?t? ? ? ?x? ? ? ? ??x? ? ? ??y? ? ? ? ??y? ? ? ??z? ? ? ? ??z

    (1)

    式中:ρ為鋼液密度,kg/m3; ceff為等價(jià)比熱容,J/(kg·K),ceff=dH/dT,H為熱焓;keff為等效導(dǎo)熱系數(shù),W/(m·K) ,keff=fsks+m(1-fs)kl[6],ks、kl分別為固液相導(dǎo)熱系數(shù),fs為固相分率,m表征對(duì)流強(qiáng)度;T為鑄坯溫度,℃。其中,物性參數(shù)ρ、H、k是溫度和鋼種成分的函數(shù),本文通過(guò)偽二元相圖方法進(jìn)行計(jì)算[7]。

    模型的初始條件和邊界條件如下.。

    1)初始條件,忽略中間包到彎月面溫降

    T(x,y,z,0)=Tc? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?(2)

    式中:Tc為中間包鋼水溫度。

    2)結(jié)晶器邊界條件

    ?T

    -k——=A-B√t? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? (3)

    ?n

    式中:t為鑄坯自彎月面開(kāi)始在結(jié)晶中停留的時(shí)間,A、B為常數(shù),通過(guò)測(cè)量結(jié)晶器進(jìn)出水溫差計(jì)算得到。

    3)二冷區(qū)邊界條件

    ?T

    -k——=h(T-Twater) +εσ[(T+273)4- (Tair+273)4]? ? (4)

    ?n

    式中:hi為二冷區(qū)第i冷卻段換熱系數(shù),W/(m2·K);Twater,Tair分別冷卻水溫度和環(huán)境空氣溫度,℃;ε為鑄坯表面綜合輻射系數(shù),σ為斯蒂芬-波爾茲曼常數(shù),W/(m2·K4)。

    4)空冷區(qū)邊界條件

    ?T

    -k——=εσ[(T+273)4- (Tair+273)4]? ? ? ? ? ? ? ? ? ? ? ?(5)

    ?n

    式中:數(shù)值模型采用有限容積法進(jìn)行離散化并采用交差隱式算法進(jìn)行求解,考慮方坯的對(duì)稱性,只計(jì)算1/4斷面。同時(shí),為了在保證模型精度的條件下提高模型的計(jì)算速度,考慮網(wǎng)格劃分與溫度梯度分布相匹配,通過(guò)優(yōu)化網(wǎng)格劃分和時(shí)間步長(zhǎng),采用基于非均勻變網(wǎng)格變步長(zhǎng)的實(shí)時(shí)算法對(duì)模型進(jìn)行求解[4,8]。

    1.2? ? 模型中的待辨識(shí)參數(shù)及修正

    依據(jù)傳熱機(jī)理所建立的模型,由于物性參數(shù)和邊界條件的不確定性引入一些待定參數(shù),它們是影響模型準(zhǔn)確性進(jìn)而可靠應(yīng)用的關(guān)鍵。根據(jù)噴淋方式的差異,二冷區(qū)各段換熱系數(shù)與水流密度的關(guān)系采用不同的經(jīng)驗(yàn)公式確定:

    1)噴水冷卻[5]

    1 570w2 i(1-0.007 5Tw)

    hi=—————————— ? ? ? ? ? ? ? ?(6)

    αi

    2)氣霧冷卻[4]

    hi=1 000wi/ αi? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?(7)

    式中:wi為二冷第i段水流密度,L/(m2·s),αi為相應(yīng)冷卻段對(duì)應(yīng)的待辨識(shí)參數(shù).

    共有(n+1)個(gè)參數(shù)需要在模型應(yīng)用前進(jìn)行辨識(shí),其中包括n個(gè)冷卻段參數(shù)αi及有效導(dǎo)熱系數(shù)中參數(shù)m,共同記為α, α=[α1, α2, …, αn, m]。

    凝固傳熱模型參數(shù)辨識(shí)的基本原理是通過(guò)不斷尋優(yōu)待辨識(shí)的參數(shù)值,使得模型計(jì)算值與相應(yīng)檢測(cè)值偏差不斷減小,以滿足應(yīng)用需求。本文采用比色測(cè)溫儀(HWSG-2H)多點(diǎn)測(cè)溫進(jìn)行模型參數(shù)的辨識(shí)修正。

    Tcal,j-Tmea,j? ?2

    minJ(α)=∑〔————〕 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?(8)

    Tmea,j

    式中:Tmea,j是鑄坯表面溫度,Tcal,j則是根據(jù)第j組測(cè)量時(shí)工藝條件及傳熱模型得到的相應(yīng)計(jì)算值。優(yōu)化目標(biāo)函數(shù)的最優(yōu)解α*作為模型參數(shù)辨識(shí)的結(jié)果。在優(yōu)化方法上,考慮凝固傳熱模型參數(shù)辨識(shí)問(wèn)題具有明顯的非線性特征,采用混沌粒子群算法[9]對(duì)目標(biāo)函數(shù)(7)進(jìn)行優(yōu)化,其中,α是作為粒子群算法中的“粒子”。其中,所使用的測(cè)溫?cái)?shù)據(jù)和相應(yīng)鑄機(jī)參數(shù)分別見(jiàn)表1和表2,修正后模型參數(shù)見(jiàn)表3。圖1給出模型修正前后模型計(jì)算表面溫度與實(shí)際測(cè)溫?cái)?shù)據(jù)的比較。

    2? ? 基于模型的工藝優(yōu)化計(jì)算

    2.1? ? 拉速的優(yōu)化

    拉速是影響鑄坯產(chǎn)量和質(zhì)量的重要因素之一。拉速工藝的制定首先是在保證爐機(jī)匹配、生產(chǎn)順行的條件下盡可能提高拉坯速度,以保證較高的生產(chǎn)效率;其次,合理的拉坯速度應(yīng)保證鑄坯具有較好的質(zhì)量。即拉速是保證高效高質(zhì)生產(chǎn)的關(guān)鍵。特別對(duì)于小方坯連鑄機(jī),一般拉速上限要保證不帶液芯矯直以避免產(chǎn)生中心裂紋或中間裂紋;同時(shí),為了充分發(fā)揮凝固末端電磁攪拌的作用,其對(duì)應(yīng)位置的中心凝固分率推薦為0.3 ~ 0.4[10]。根據(jù)凝固傳熱模型計(jì)算得到表4,則對(duì)應(yīng)拉速上限和合適的拉速范圍分別為2.4 m/min、2.2 ~2.4 m/min。優(yōu)化后拉速較原來(lái)拉速2.0~2.2 m/min提高10%~20%。

    2.2? ? 二次冷卻的優(yōu)化

    連鑄過(guò)程中,二冷對(duì)鑄坯溫度場(chǎng)具有決定性影響,而鑄坯表面溫度是鑄坯溫度場(chǎng)和凝固狀態(tài)的重要表征[9]。鑄坯溫度的變化與裂紋的產(chǎn)生擴(kuò)展有關(guān),優(yōu)化的二冷制度應(yīng)該保證鑄坯表面溫度處于較高的塑性區(qū)間。本文基于目標(biāo)溫度的方法[11]對(duì)二冷進(jìn)行優(yōu)化。對(duì)于小斷面生產(chǎn)高碳鋼,一般采用強(qiáng)冷模式,通過(guò)強(qiáng)冷可以減小二次枝晶間距,減輕鑄坯偏析程度。對(duì)應(yīng)于強(qiáng)冷模式,二冷區(qū)鑄坯表面溫度設(shè)置于高溫塑性區(qū)(900 ~ 1 100 ℃)低溫區(qū)段,通常接近900 ℃,本文中二冷分為四個(gè)區(qū),各區(qū)目標(biāo)溫度依次設(shè)定為900 ℃、900 ℃、900 ℃、900 ℃。以82B為例,過(guò)熱度(30 ℃)、拉速條件下,優(yōu)化調(diào)整二冷各區(qū)水量,使二冷各段末溫度達(dá)到目標(biāo)溫度,見(jiàn)圖2。不同拉速下可得到不同的二冷各段水量,見(jiàn)表5。優(yōu)化后的比水量由1.28 ~ 1.49 L/kg增強(qiáng)到1.52 ~ 1.81 L/kg。

    2.3? ? 電磁攪拌的優(yōu)化

    電磁攪拌是目前針對(duì)改善小方坯生產(chǎn)內(nèi)部質(zhì)量的重要措施,其作用主要是提高中心等軸晶率、改善中心偏析等[12]。電磁攪拌包含結(jié)晶器電磁攪拌和凝固末端電磁攪拌。要充分發(fā)揮電磁攪拌的作用,主要是兩點(diǎn):一是在合適的位置進(jìn)行攪拌,關(guān)鍵是凝固末端電磁攪拌的位置與拉速匹配,其中,對(duì)于新安裝電磁攪拌的鑄機(jī),凝固末端電磁攪拌安裝的位置要進(jìn)行合理設(shè)計(jì),而對(duì)于電磁攪拌安裝位置已固定的鑄機(jī),則通過(guò)調(diào)節(jié)拉速使中心凝固分率在安裝位置處取得合適的值;二是設(shè)定合適的攪拌強(qiáng)度(凝固前沿?cái)嚢杷俣龋┖头绞?,其中影響攪拌?qiáng)度的參數(shù)包括攪拌電流和頻率,攪拌方式則包括連續(xù)攪拌和交替攪拌方式。本文中結(jié)晶器電磁攪拌和凝固末端電磁攪拌均采用連續(xù)攪拌。連續(xù)攪拌下,凝固前沿?cái)嚢杷俣扰c磁感應(yīng)強(qiáng)度、頻率、液芯大小以及鋼水粘度和電導(dǎo)率等有關(guān),具體如式(9)[10]。

    uθmax=0.1KB√fr(√σ/ρμ-0.44)≈0.064 4B√fr? ? ? ? ? ? ? (9)

    式中:K是與攪拌器有關(guān)的常數(shù),B是鑄坯中心磁感應(yīng)強(qiáng)度(mT),f為頻率(Hz),r為液芯半徑(指兩相區(qū)+液相區(qū)部分,m),括號(hào)內(nèi)為與鋼種相關(guān)的項(xiàng),σ是電導(dǎo)率(S/cm),ρ為密度(g/cm3),μ為鋼種動(dòng)力粘度(Pa·s)。根據(jù)式(9),實(shí)踐中可通過(guò)磁感應(yīng)強(qiáng)度推算凝固前沿?cái)嚢杷俣龋ǔUJ(rèn)為凝固前沿合適的攪拌速度為0.3 ~ 0.6 m/s,不超過(guò)1 m/s。本文中磁感應(yīng)強(qiáng)度通過(guò)特斯拉計(jì)(HT201)測(cè)量獲得。

    3? ? 工藝優(yōu)化實(shí)踐與分析

    根據(jù)模型優(yōu)化測(cè)算結(jié)果,結(jié)合現(xiàn)場(chǎng)生產(chǎn)經(jīng)驗(yàn),設(shè)計(jì)相關(guān)試驗(yàn)以確定工藝調(diào)整方向。主要涉及拉速、二冷、結(jié)晶器電磁攪拌、凝固末端電磁攪拌4個(gè)工藝變量。選擇相同流進(jìn)行不同工藝的比較研究,見(jiàn)表6。

    根據(jù)試驗(yàn)得到的結(jié)果,f工藝條件下中心縮孔和中心碳偏析指數(shù)均表現(xiàn)良好,后續(xù)擬以該工藝為基礎(chǔ)進(jìn)一步開(kāi)展批量生產(chǎn)試驗(yàn)。針對(duì)所取得的試驗(yàn)效果,有以下分析。

    1)提高拉速對(duì)縮孔和偏析均有所改善,其主要原因?yàn)椋哼m當(dāng)提高拉速有利于更充分發(fā)揮凝固末端電磁攪拌的作用,有利于打破“小鋼錠”凝固,增強(qiáng)中心區(qū)域的補(bǔ)縮。

    2)適當(dāng)增強(qiáng)MEMS也有利于同時(shí)改善偏析和縮孔,其主要原因?yàn)椋涸鰪?qiáng)結(jié)晶器電磁攪拌:①促進(jìn)中心形核的增加,增加并細(xì)化中心等軸晶區(qū),有利于減輕中心疏松和縮孔;②由于等軸晶區(qū)增加,中心偏析有所改善。

    3)適當(dāng)減弱凝固末端電磁攪拌有利于改善偏析,原因可能是過(guò)強(qiáng)的攪拌產(chǎn)生的溶質(zhì)沖刷作用將溶質(zhì)帶到中心會(huì)引起中心碳偏析的加重。

    4)強(qiáng)冷具有多方面的作用:①促進(jìn)柱狀晶生長(zhǎng)更容易搭橋,加重中心偏析和縮孔;② 鑄坯表面冷卻加強(qiáng)起到一定熱壓下作用,鑄坯芯部致密度增加,中心疏松有所減輕,對(duì)縮孔影響不大;③冷卻速率增加,二次枝晶間距減小,液相流動(dòng)減弱微觀偏析減輕,中心偏析減輕;④液芯變短,不利于發(fā)揮凝固末端電磁攪拌作用增強(qiáng)補(bǔ)縮,縮孔和V型偏析加重。從結(jié)果分析來(lái)看,②③④的作用更為主導(dǎo)。

    4? ? 結(jié)? ? 論

    本文針對(duì)小方坯82B連鑄提高拉速條件下質(zhì)量?jī)?yōu)化問(wèn)題,開(kāi)展建模、優(yōu)化計(jì)算和試驗(yàn),得到以下結(jié)論:

    1)為取得良好質(zhì)量,針對(duì)高碳鋼合適的拉速應(yīng)與凝固末端電磁攪拌位置相匹配;

    2)適當(dāng)強(qiáng)冷有利于提高鑄坯芯部致密度,改善疏松和縮孔和偏析狀態(tài);

    3)采用較強(qiáng)的結(jié)晶器電磁攪拌和稍弱的凝固末端電磁攪拌有利于改善鑄坯芯部質(zhì)量。

    參考文獻(xiàn)

    [1]? ? 王長(zhǎng)順, 郭福建, 李廣龍,等. 中心偏析控制對(duì)型鋼低溫韌性的影響[J]. 鋼鐵, 2019, 54(8):202.

    [2]? ? 李桂英, 李麥麥. 鋼簾線用鋼C82D2盤(pán)條拉拔斷絲成因分析和工藝改進(jìn)[J]. 特殊鋼, 2021, 42(4):30.

    [3]? ? 賀洪印, 劉海濤, 陳偉慶,等. 高碳鋼小方坯連鑄參數(shù)優(yōu)化及鑄坯中心偏析控制[J]. 鋼鐵釩鈦, 2011, 32(1):84.

    [4]? ? Yang J, Xie Z, Ji Z P, et al. Real-time heat transfer model based on variable non-uniform grid for dynamic control of continuous casting billets[J]. ISIJ International, 2014, 54(2): 328.

    [5]? ? Ma J C, Xie Z, Jia G L. Applying of Real-time heat transfer and solidification model on the dynamic control system of billet continuous casting[J]. ISIJ International, 2008, 48(12): 1722.

    [6]? ? Louhenkilpi S, Laitinen E, Nieminen R. Real-time simulation of heat transfer in continuous casting[J]. Metallurgical Transactions B, 1993, 24(4): 685.

    [7]? ? Xie Z, Yang J. Calculation of Solidification-Related Thermophysical properties of steels based on Fe-C pseudobinary phase diagram[J]. Steel Research International, 2015, 86(7): 766.

    [8]? ? Jian Yang, Zhi Xie, Hongji Meng, et al. Multiple time steps optimization for real-time heat transfer model of continuous casting billets[J], International Journal of Heat and Mass Transfer, 2014, 76:492.

    [9]? ? 孟紅記, 鄭鵬, 梅國(guó)暉,等. 基于混沌序列的粒子群優(yōu)化算法[J]. 控制與決策, 2006, 21(3):263.

    [10]? ? Yang J, Ji Z P, Liu W H, et al. Digital-twin-based coordinated optimal control for steel continuous casting process[J]. Metals, 2023, 13(4): 816.

    [11]? ? 紀(jì)振平,高志強(qiáng). 基于SVM的連鑄二冷目標(biāo)控制模型研究[J]. 沈陽(yáng)理工大學(xué)學(xué)報(bào), 2018, 37(1):47.

    [12]? ? 王彪, 謝植, 賈光霖,等. 凝固末端電磁攪拌參數(shù)確定及其對(duì)中心偏析的影響[J]. 鋼鐵, 2007, 42(3):18.

    白带黄色成豆腐渣| www.av在线官网国产| 亚洲精品影视一区二区三区av| 高清av免费在线| 国产精品国产三级专区第一集| 国产精品一区二区三区四区免费观看| 欧美极品一区二区三区四区| 久久精品夜色国产| 97精品久久久久久久久久精品| 一二三四中文在线观看免费高清| 国产成人免费观看mmmm| 一区二区av电影网| 欧美日韩精品成人综合77777| 人妻夜夜爽99麻豆av| 伊人久久精品亚洲午夜| 日韩一区二区视频免费看| 欧美高清成人免费视频www| 国产成人精品一,二区| 美女被艹到高潮喷水动态| 人妻 亚洲 视频| 99精国产麻豆久久婷婷| 国产成人免费观看mmmm| 午夜福利视频1000在线观看| 精品99又大又爽又粗少妇毛片| 久久精品国产亚洲av天美| 欧美老熟妇乱子伦牲交| 亚洲av中文av极速乱| 午夜激情久久久久久久| 简卡轻食公司| 国产午夜精品一二区理论片| 午夜福利网站1000一区二区三区| 日本wwww免费看| 视频中文字幕在线观看| 嘟嘟电影网在线观看| 欧美变态另类bdsm刘玥| 亚洲精品视频女| 国产国拍精品亚洲av在线观看| 国产男女内射视频| 国产精品三级大全| 最近的中文字幕免费完整| 国产成人精品一,二区| 人体艺术视频欧美日本| 日日撸夜夜添| 日本一二三区视频观看| 欧美 日韩 精品 国产| 中国美白少妇内射xxxbb| 22中文网久久字幕| 麻豆国产97在线/欧美| 欧美激情在线99| 麻豆乱淫一区二区| 亚洲三级黄色毛片| 国内少妇人妻偷人精品xxx网站| 亚洲熟女精品中文字幕| 国产免费视频播放在线视频| 91精品国产九色| 亚洲,欧美,日韩| 午夜日本视频在线| 青青草视频在线视频观看| 天堂俺去俺来也www色官网| 亚洲三级黄色毛片| 久久久午夜欧美精品| 99久久精品热视频| 亚洲色图综合在线观看| 日本三级黄在线观看| 网址你懂的国产日韩在线| 婷婷色麻豆天堂久久| 啦啦啦在线观看免费高清www| 久久久久久久精品精品| 欧美成人精品欧美一级黄| 亚洲,欧美,日韩| 国产一区二区三区av在线| 国产精品精品国产色婷婷| 国产亚洲精品久久久com| 亚洲国产成人一精品久久久| 国产伦理片在线播放av一区| 精品国产乱码久久久久久小说| 舔av片在线| 精品久久久久久久末码| 国产在线一区二区三区精| 免费大片黄手机在线观看| 日本午夜av视频| 久久久久久久精品精品| 亚洲av在线观看美女高潮| 三级经典国产精品| 精品久久久久久久久亚洲| 啦啦啦中文免费视频观看日本| 午夜福利视频1000在线观看| 蜜桃久久精品国产亚洲av| 日日摸夜夜添夜夜添av毛片| 亚洲一区二区三区欧美精品 | 激情 狠狠 欧美| 丝袜脚勾引网站| 午夜精品一区二区三区免费看| 国产黄片美女视频| 久久99热这里只有精品18| 国产av不卡久久| 欧美激情在线99| 韩国高清视频一区二区三区| 免费看不卡的av| 久久久亚洲精品成人影院| 一边亲一边摸免费视频| 亚洲人与动物交配视频| 国产69精品久久久久777片| tube8黄色片| 丝袜喷水一区| 嫩草影院入口| 一级片'在线观看视频| 99九九线精品视频在线观看视频| 成人国产av品久久久| 日本爱情动作片www.在线观看| 最后的刺客免费高清国语| av女优亚洲男人天堂| 三级国产精品片| 日本三级黄在线观看| 国产又色又爽无遮挡免| 久久女婷五月综合色啪小说 | 日韩一本色道免费dvd| 亚洲怡红院男人天堂| 国产视频首页在线观看| 能在线免费看毛片的网站| 欧美成人精品欧美一级黄| 18+在线观看网站| 在线看a的网站| 日韩精品有码人妻一区| 国产精品嫩草影院av在线观看| 水蜜桃什么品种好| 国产午夜精品一二区理论片| 精品国产三级普通话版| 免费大片18禁| 精品人妻熟女av久视频| 一个人看视频在线观看www免费| 亚洲av国产av综合av卡| 亚洲欧美日韩卡通动漫| 直男gayav资源| 精品国产露脸久久av麻豆| 亚洲国产欧美在线一区| 欧美丝袜亚洲另类| 伦理电影大哥的女人| 亚洲国产精品专区欧美| 国产亚洲精品久久久com| 国产综合懂色| 成人毛片a级毛片在线播放| 成人特级av手机在线观看| 国产精品99久久99久久久不卡 | 欧美另类一区| 久久久久久久大尺度免费视频| 久久久久久久久久成人| 最近最新中文字幕免费大全7| 亚洲婷婷狠狠爱综合网| 另类亚洲欧美激情| 亚洲精品视频女| 嫩草影院新地址| 丝瓜视频免费看黄片| av.在线天堂| 国产精品三级大全| 亚洲久久久久久中文字幕| 成人亚洲欧美一区二区av| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 午夜精品一区二区三区免费看| 亚洲电影在线观看av| xxx大片免费视频| 久久久久九九精品影院| 日本三级黄在线观看| 国产免费福利视频在线观看| 免费电影在线观看免费观看| 午夜福利视频精品| 小蜜桃在线观看免费完整版高清| 黄片无遮挡物在线观看| 亚洲av免费高清在线观看| 日韩视频在线欧美| 免费观看a级毛片全部| 精品久久久久久久末码| 免费少妇av软件| 蜜臀久久99精品久久宅男| 亚洲,一卡二卡三卡| 成年av动漫网址| 国产免费福利视频在线观看| 久久影院123| 国产精品久久久久久久久免| av天堂中文字幕网| 777米奇影视久久| 97精品久久久久久久久久精品| 日韩国内少妇激情av| 99久久中文字幕三级久久日本| 久久国内精品自在自线图片| 亚洲,欧美,日韩| 少妇人妻 视频| 五月天丁香电影| 91狼人影院| 草草在线视频免费看| 欧美成人一区二区免费高清观看| 另类亚洲欧美激情| 中国三级夫妇交换| 亚洲成人一二三区av| 久久久久久久大尺度免费视频| 一个人看的www免费观看视频| 最近最新中文字幕免费大全7| 在现免费观看毛片| 成人高潮视频无遮挡免费网站| 亚洲欧美成人综合另类久久久| 男插女下体视频免费在线播放| 国产日韩欧美在线精品| 少妇被粗大猛烈的视频| 2021少妇久久久久久久久久久| 国产免费又黄又爽又色| 日韩人妻高清精品专区| 一边亲一边摸免费视频| 久久久色成人| 国产人妻一区二区三区在| 99热全是精品| 中文在线观看免费www的网站| 亚洲性久久影院| av在线蜜桃| 99热国产这里只有精品6| av一本久久久久| 国产亚洲av嫩草精品影院| 超碰97精品在线观看| 久久久精品免费免费高清| 国产免费福利视频在线观看| 久久99热这里只有精品18| 欧美少妇被猛烈插入视频| 男人爽女人下面视频在线观看| 美女高潮的动态| 日韩电影二区| 午夜福利网站1000一区二区三区| 亚州av有码| 欧美激情在线99| 又爽又黄无遮挡网站| 美女被艹到高潮喷水动态| 久久久久久久久久人人人人人人| 夫妻午夜视频| av在线app专区| freevideosex欧美| 日韩欧美精品免费久久| 尾随美女入室| 欧美另类一区| 久久久久久久亚洲中文字幕| 国产一区二区在线观看日韩| 国产视频内射| .国产精品久久| 国国产精品蜜臀av免费| 别揉我奶头 嗯啊视频| 亚洲精品日韩在线中文字幕| 高清欧美精品videossex| 在线播放无遮挡| 亚洲av.av天堂| 国产69精品久久久久777片| 亚洲电影在线观看av| 男女边吃奶边做爰视频| 乱码一卡2卡4卡精品| freevideosex欧美| 男人爽女人下面视频在线观看| kizo精华| 九色成人免费人妻av| 日本黄色片子视频| 国产欧美亚洲国产| 我要看日韩黄色一级片| 中文在线观看免费www的网站| 亚洲电影在线观看av| 少妇的逼好多水| 嘟嘟电影网在线观看| 日韩欧美 国产精品| 亚洲自偷自拍三级| 亚洲国产精品国产精品| 黄片无遮挡物在线观看| 韩国av在线不卡| 亚洲第一区二区三区不卡| 精品少妇久久久久久888优播| 亚洲精品乱码久久久v下载方式| 寂寞人妻少妇视频99o| 69人妻影院| 一级毛片aaaaaa免费看小| 日韩不卡一区二区三区视频在线| 男人和女人高潮做爰伦理| 国产成人午夜福利电影在线观看| 日本wwww免费看| 日本黄色片子视频| 人人妻人人澡人人爽人人夜夜| 国产成人福利小说| 日本三级黄在线观看| 午夜福利网站1000一区二区三区| 久久6这里有精品| 波多野结衣巨乳人妻| 亚洲成人精品中文字幕电影| 高清欧美精品videossex| 亚洲国产日韩一区二区| av在线天堂中文字幕| 亚洲国产成人一精品久久久| 亚洲自拍偷在线| 国产高清三级在线| av天堂中文字幕网| 久久久久久久午夜电影| 美女高潮的动态| 欧美激情久久久久久爽电影| 久久人人爽人人爽人人片va| 亚洲人成网站在线播| 亚洲精品自拍成人| 久久国内精品自在自线图片| 精品国产露脸久久av麻豆| 国产免费视频播放在线视频| 久久久久国产网址| 亚洲精品久久午夜乱码| 熟女av电影| 亚洲内射少妇av| av在线app专区| 欧美bdsm另类| 99热6这里只有精品| 三级男女做爰猛烈吃奶摸视频| 国产黄a三级三级三级人| 亚洲国产精品国产精品| 日韩伦理黄色片| 免费在线观看成人毛片| 国产成人精品久久久久久| 亚洲高清免费不卡视频| 黄色一级大片看看| 一级二级三级毛片免费看| 亚洲三级黄色毛片| 色吧在线观看| 美女高潮的动态| 亚洲在久久综合| 国产精品爽爽va在线观看网站| 欧美 日韩 精品 国产| 亚洲真实伦在线观看| 少妇猛男粗大的猛烈进出视频 | 国产精品三级大全| 精品久久久噜噜| 精品一区二区免费观看| 亚洲人成网站在线播| 国产真实伦视频高清在线观看| 亚洲国产最新在线播放| 国产黄片美女视频| 久久ye,这里只有精品| 欧美97在线视频| 国产成人91sexporn| 中文字幕亚洲精品专区| 免费观看a级毛片全部| 国产精品一二三区在线看| av网站免费在线观看视频| 国产亚洲av嫩草精品影院| 国产一区二区三区综合在线观看 | videossex国产| 黄色配什么色好看| 亚洲欧美日韩卡通动漫| 亚洲丝袜综合中文字幕| 一级片'在线观看视频| 国产亚洲av片在线观看秒播厂| 精品人妻偷拍中文字幕| 男人和女人高潮做爰伦理| 大香蕉久久网| 啦啦啦中文免费视频观看日本| 久久这里有精品视频免费| 18禁裸乳无遮挡动漫免费视频 | 国产成人91sexporn| 少妇人妻精品综合一区二区| 国产女主播在线喷水免费视频网站| 校园人妻丝袜中文字幕| 亚洲精华国产精华液的使用体验| 一级毛片 在线播放| 99热这里只有精品一区| 亚洲精品aⅴ在线观看| 国产免费视频播放在线视频| 亚洲在线观看片| av福利片在线观看| 又大又黄又爽视频免费| 日本av手机在线免费观看| av在线app专区| 街头女战士在线观看网站| 亚洲欧美清纯卡通| 久久久久久久久久成人| 一级毛片 在线播放| 男人狂女人下面高潮的视频| 寂寞人妻少妇视频99o| 日本av手机在线免费观看| 国产欧美日韩精品一区二区| 亚洲精品,欧美精品| 美女cb高潮喷水在线观看| av在线亚洲专区| 日韩三级伦理在线观看| 美女视频免费永久观看网站| 天堂中文最新版在线下载 | 国产成人aa在线观看| 色5月婷婷丁香| xxx大片免费视频| 蜜桃久久精品国产亚洲av| 中文字幕av成人在线电影| 春色校园在线视频观看| 亚洲av二区三区四区| 日韩成人av中文字幕在线观看| 亚洲三级黄色毛片| 日韩制服骚丝袜av| 亚洲人成网站高清观看| 97超视频在线观看视频| 最近最新中文字幕大全电影3| 国产精品久久久久久精品电影| 亚洲国产成人一精品久久久| 国产精品久久久久久精品电影| 亚洲不卡免费看| 亚洲综合色惰| 国产色婷婷99| 91久久精品国产一区二区三区| 午夜免费观看性视频| av一本久久久久| 中国美白少妇内射xxxbb| 亚洲av一区综合| 亚洲最大成人中文| 91精品伊人久久大香线蕉| 久久久久久久久久久免费av| 一区二区三区乱码不卡18| 亚洲电影在线观看av| 日韩免费高清中文字幕av| 亚洲欧美精品自产自拍| 黄色视频在线播放观看不卡| 自拍偷自拍亚洲精品老妇| 超碰97精品在线观看| 国产一区二区在线观看日韩| 涩涩av久久男人的天堂| 日本猛色少妇xxxxx猛交久久| 亚洲丝袜综合中文字幕| 欧美一区二区亚洲| 黄色一级大片看看| 日本一二三区视频观看| 国产免费一级a男人的天堂| 欧美xxxx性猛交bbbb| 免费看a级黄色片| 蜜桃亚洲精品一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 尤物成人国产欧美一区二区三区| 亚洲av免费在线观看| videos熟女内射| 亚洲国产欧美在线一区| 亚洲精品一二三| 最近中文字幕高清免费大全6| 一个人看视频在线观看www免费| 69人妻影院| 免费av不卡在线播放| 能在线免费看毛片的网站| 免费av不卡在线播放| 18禁在线无遮挡免费观看视频| 免费av不卡在线播放| av卡一久久| 五月开心婷婷网| 一二三四中文在线观看免费高清| 欧美最新免费一区二区三区| 亚洲欧洲国产日韩| 亚洲精品中文字幕在线视频 | 国产一区亚洲一区在线观看| 亚洲精品国产av成人精品| 黄色欧美视频在线观看| 久久久久久久午夜电影| 亚洲性久久影院| 亚洲aⅴ乱码一区二区在线播放| 成人亚洲精品一区在线观看 | 国产亚洲精品久久久com| av在线蜜桃| 最近2019中文字幕mv第一页| 国产一区二区三区av在线| 亚洲美女搞黄在线观看| av又黄又爽大尺度在线免费看| 大又大粗又爽又黄少妇毛片口| 久久久国产一区二区| 国产淫语在线视频| 天堂网av新在线| 欧美日韩一区二区视频在线观看视频在线 | 国产一区二区三区av在线| 免费av毛片视频| 久久久a久久爽久久v久久| 国产成人福利小说| 少妇熟女欧美另类| 久久精品国产自在天天线| 免费黄频网站在线观看国产| 涩涩av久久男人的天堂| 又粗又硬又长又爽又黄的视频| 精品人妻一区二区三区麻豆| 日产精品乱码卡一卡2卡三| 蜜桃亚洲精品一区二区三区| 一本色道久久久久久精品综合| 日韩亚洲欧美综合| 国产乱人视频| 亚洲欧美日韩另类电影网站 | av在线亚洲专区| 麻豆国产97在线/欧美| 欧美变态另类bdsm刘玥| 亚洲,一卡二卡三卡| 国产精品国产av在线观看| 久久97久久精品| 国产免费又黄又爽又色| 久久久久久伊人网av| 久久国内精品自在自线图片| 日本一本二区三区精品| 午夜福利高清视频| 六月丁香七月| 国产日韩欧美亚洲二区| eeuss影院久久| 国产美女午夜福利| 欧美高清性xxxxhd video| 亚洲精品成人久久久久久| 日韩成人av中文字幕在线观看| 一级av片app| 国产乱人视频| eeuss影院久久| 极品教师在线视频| 国产精品99久久久久久久久| 精品久久久久久电影网| 成年免费大片在线观看| 亚洲av欧美aⅴ国产| 日韩三级伦理在线观看| 中文字幕人妻熟人妻熟丝袜美| 亚洲欧美一区二区三区黑人 | 亚洲av不卡在线观看| av福利片在线观看| 日韩亚洲欧美综合| 91午夜精品亚洲一区二区三区| 男男h啪啪无遮挡| 天堂中文最新版在线下载 | 黄色欧美视频在线观看| 视频区图区小说| 国产日韩欧美在线精品| 国产亚洲91精品色在线| av在线亚洲专区| 国产高清国产精品国产三级 | 国产探花极品一区二区| 国产成人精品婷婷| 黄色配什么色好看| 亚洲国产精品成人综合色| 成人高潮视频无遮挡免费网站| 日韩强制内射视频| 成年版毛片免费区| 一二三四中文在线观看免费高清| 人人妻人人澡人人爽人人夜夜| 舔av片在线| 亚洲欧美成人精品一区二区| 五月玫瑰六月丁香| 亚洲欧美日韩另类电影网站 | 国产日韩欧美在线精品| 男女国产视频网站| 国产亚洲av片在线观看秒播厂| 五月玫瑰六月丁香| av线在线观看网站| 国产女主播在线喷水免费视频网站| 干丝袜人妻中文字幕| 亚洲精品日韩在线中文字幕| 国产日韩欧美亚洲二区| 日本午夜av视频| 国产成人精品一,二区| 亚洲精品色激情综合| 国产欧美日韩精品一区二区| 亚洲综合色惰| 免费黄网站久久成人精品| 两个人的视频大全免费| 一本久久精品| 亚洲熟女精品中文字幕| 免费黄网站久久成人精品| 免费av毛片视频| 午夜精品国产一区二区电影 | av女优亚洲男人天堂| 99热国产这里只有精品6| 人妻 亚洲 视频| 午夜福利视频精品| 久久久午夜欧美精品| 老女人水多毛片| 精品一区在线观看国产| av天堂中文字幕网| 亚洲欧美日韩另类电影网站 | 日韩精品有码人妻一区| 国产午夜精品一二区理论片| 午夜福利视频精品| 天堂网av新在线| 一级毛片电影观看| 夫妻性生交免费视频一级片| 韩国高清视频一区二区三区| 欧美精品国产亚洲| 日韩国内少妇激情av| 国产成人免费无遮挡视频| 人人妻人人澡人人爽人人夜夜| 国产精品99久久99久久久不卡 | 高清日韩中文字幕在线| 男女无遮挡免费网站观看| 亚洲综合色惰| 久久精品国产亚洲av涩爱| 国产精品熟女久久久久浪| 亚洲欧美一区二区三区黑人 | 小蜜桃在线观看免费完整版高清| 欧美人与善性xxx| 亚洲欧美中文字幕日韩二区| av播播在线观看一区| 大又大粗又爽又黄少妇毛片口| 在线观看一区二区三区| 国产久久久一区二区三区| 黑人高潮一二区| 天天一区二区日本电影三级| 一级毛片久久久久久久久女| 波多野结衣巨乳人妻| av在线播放精品| 午夜日本视频在线| 中文字幕免费在线视频6| 亚洲精品视频女| 久久久久久久亚洲中文字幕| 69av精品久久久久久| 精品一区二区三卡| 一区二区三区免费毛片| 日产精品乱码卡一卡2卡三| 午夜精品一区二区三区免费看| 各种免费的搞黄视频| 搞女人的毛片| 亚洲国产精品成人久久小说| 色综合色国产| 天堂俺去俺来也www色官网| 免费不卡的大黄色大毛片视频在线观看| 午夜视频国产福利| 久久久成人免费电影| 国产精品国产av在线观看| 大码成人一级视频| 白带黄色成豆腐渣| 亚洲欧美中文字幕日韩二区| 肉色欧美久久久久久久蜜桃 |