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

不同熱裝方式下鑄坯溫度分布研究

2017-09-28 03:17:54葛建華王明林趙文博
上海金屬 2017年4期

葛建華 王明林 趙文博 張 慧

(1.鋼鐵研究總院連鑄技術國家工程研究中心,北京 100081; 2.西安建筑科技大學冶金工程學院,陜西西安 710055)

不同熱裝方式下鑄坯溫度分布研究

葛建華1王明林1趙文博2張 慧1

(1.鋼鐵研究總院連鑄技術國家工程研究中心,北京 100081; 2.西安建筑科技大學冶金工程學院,陜西西安 710055)

采用有限元軟件對步進式加熱爐內鋼坯加熱過程進行分析計算,建立鑄坯加熱溫度場模型,計算不同熱裝條件下鑄坯中心與表面的溫度變化曲線。以鑄坯1/4斷面為研究對象,分析鑄坯斷面全流程溫度分布,得到連鑄坯表面溫度、角部溫度、心部溫度的變化曲線。模擬分析不同熱裝方式裝爐的鑄坯在加熱爐內的加熱升溫情況,提取鑄坯表面和中心位置的加熱曲線。模擬加熱爐內加熱過程,為降低鑄坯斷面溫差、加熱爐能耗以及提高生產效率提供依據。

熱裝方式 連鑄坯 加熱爐 溫度分布 能耗

冶金企業是能源消耗大戶,各工序幾乎都有熱量消耗,若能在上一道工序中保存余熱供下一道工序使用,可以明顯降低該工序的能源消耗[1]。隨著連鑄技術的發展成熟,充分利用連鑄坯自身熱能的熱送熱裝技術已成為鋼鐵生產中節能降耗的重要手段之一[2]。鑄坯從進入結晶器至出加熱爐過程中溫度變化、相變過程十分復雜,實際生產中無法時刻測量鑄坯表面溫度在各生產工序的變化情況。尤其對于Ti和Nb等微合金鋼而言,直裝、熱裝的應用會導致析出物的析出歷程改變,給鑄坯的組織、性能和表面質量等帶來不利影響。鑄坯的溫控是實現連鑄全流程智能化的重要因素之一;同時,合理的熱裝制度也有利于鋼鐵企業的節能降耗。近年來,越來越多的科技工作者開始采用數值模擬方法探究各類鑄坯溫度的變化過程,為鑄坯的直裝、熱裝工藝制定提供指導[3]。

1 鑄坯加熱過程數學模型

為了研究熱裝坯和冷裝坯在加熱爐內的升溫過程,需得到鑄坯進入加熱爐時的溫度分布情況。采用有限元數值模擬軟件結合現場生產流程,對237 mm×1 600 mm斷面的連鑄坯進行加熱過程數值模擬計算。建立鋼水經過結晶器、二冷區、空冷區的溫降過程模型,得到直裝坯的溫度場,可以提高數學模型的計算精度和準確度[4]。采用有限元數值模擬軟件建立鑄坯在加熱爐內的加熱模型,通過非穩態計算獲得鑄坯加熱溫度變化曲線。

1.1 連鑄坯溫度場數學模型

鋼水的主要化學成分見表1。

表1 鋼水的主要化學成分(質量分數)Table 1 Main chemical composition of the liquid steel (mass fraction) %

其他工藝參數如下:大包鋼水溫度為1 596 ℃;澆注溫度為1 556~1 561 ℃,本文澆注溫度取tP=1 559 ℃;連鑄機恒拉速為1.2 m/min;板坯斷尺寸為237 mm×1 600 mm。

鑄坯1/4斷面模型邊界條件如下:

(1)結晶器換熱采用綜合換熱系數,見式(1)[5]:

(1)

式中:u為拉坯速度,m/min;X為彎液面到結晶器底部的距離,cm;Fi為結晶器四面冷卻水的流量,L/min;Δt為結晶器冷卻水溫差,℃;W為結晶器寬度或厚度,mm;Ts為鋼的固相線溫度,℃;取k=0.029 274 kJ/(m·>s·>℃);c=0.752 4 kJ/(kg·℃);ρ=7 000 kg/m3;ρw=1 000 kg/m3;CW=4.182 kJ/(kg·>℃)。

鋼的液相線和固相線溫度計算經驗公式見式(2)、式(3)[6]:

Tl= 1 536-[90W(C)+6.2W(Si)+1.7W(Mn)+

28W(P)+40W(S)+2.6W(Cu)+

2.9W(Ni)+1.8W(Cr)+5.1W(Al)]

(2)

Ts=1 536-[415.3W(C)+12.3W(Si)+

6.8W(Mn)+124.5W(P)+183.9W(S)+

4.3W(Ni)+1.4W(Cr)+4.1W(Al)]

(3)

(2)二冷區寬面采用對流換熱邊界,見式(4)[4,7]:

h=A·Wn

(4)

式中:h為對流換熱系數,kW/(m·℃);W為水流密度,L/(m2·s);n為常數,取0.5~1.0;A為常數,取120~420。

(3)二冷區窄面采用輻射換熱邊界,見式(5)[8]:

q=εσ[(Tb+273)4-(Ta+273)4]

(5)

式中:ε為輻射系數,取0.8;σ為玻爾茲曼常數,取5.67×10-8W/(m·℃);Ta為噴水冷卻區環境溫度,40 ℃;Tb為鑄坯表面溫度。

(4)空冷區寬面和窄面采用輻射換熱,公式同式(5)。

建立鑄坯1/4斷面二維非穩態傳熱模型,倒角與鑄坯寬面呈30°,倒角斜邊55 mm[9]。

鋼液流經結晶器用時40 s,凝固形成坯殼。鑄坯在二冷區快速冷卻,經1 200 s到達火焰切割處切割,切割后在輥道上運輸至加熱爐前,用時約1 000 s。到達加熱爐前輥道上的鑄坯根據生產計劃進行調度,或下線堆垛緩冷,或直接進入加熱爐加熱。直接進入加熱爐加熱的鑄坯初始溫度場為經過結晶器、二冷區、空冷區后到達加熱爐前輥道的鑄坯的溫度場,如圖1所示。此時鑄坯中心溫度為1 010 ℃,角部溫度為650 ℃,表面中心溫度為850 ℃,等溫線呈橢圓環狀分布,溫度由內向外逐漸降低,角部溫度與鑄坯中心溫差約360 ℃,斷面溫差較大。

圖1 鑄坯經過空冷區后的溫度云圖Fig.1 Temperature contour of section in casting slab out of air cooling zone

下線的鑄坯經堆垛區堆垛緩冷,待溫度達到溫裝或冷裝標準時送入加熱爐加熱。下文分別模擬經堆垛后平均溫度為30、400、500、600 ℃的鑄坯在加熱爐中的加熱情況,以及不下線直接進加熱爐的鑄坯的加熱情況。

1.2 加熱爐加熱制度模型

鑄坯到達加熱爐前輥道后,一般采用三種生產制度進行生產。其一是直裝,從連鑄機出來的鑄坯不經過下線處理,直接送至加熱爐進行加熱,或者進行邊角補熱后直接熱軋,這種方式可以最大化地利用鑄坯的熱量;其二是溫裝,即從連鑄機生產的鑄坯經過下線堆垛,待鑄坯溫度冷卻到500 ℃左右時送至加熱爐進行加熱,待溫度達到設定值后進行熱軋;其三是冷裝,即鑄坯下線堆垛至常溫后送至加熱爐,經過加熱爐加熱到目標溫度后進行軋制,此方法所需能耗最大。鑄坯在加熱爐達到出爐溫度的標準是鑄坯最高溫度和最低溫度的差值≤20 ℃,且鑄坯整體溫度達到目標溫度[10]。

加熱爐爐膛分預熱段、一加熱段、二加熱段和均熱段。各段對應的長度和溫度見表2。

表2 加熱爐結構分布Table 2 Structure distribution of heating furnace

利用有限元數值模擬方法模擬鑄坯在加熱爐內的加熱過程,基本假設如下[11]:

(1)鑄坯內部沒有內熱源,同一時刻鑄坯各節點的熱物性參數相同;

(2)加熱爐內溫度分布穩定,爐膛各段介質溫度分布均勻,忽略爐長方向各段之間的輻射換熱;

(3)忽略鑄坯表面氧化鐵皮對傳熱的影響;

(4)加熱爐內相鄰兩塊鑄坯之間有一定的間隙,鑄坯上下受熱條件相同,鑄坯左右兩側受熱條件相同;

(5)忽略板坯加熱過程中相變熱,忽略板坯長度方向的熱傳導。

鑄坯斷面在不停的運動中,隨時間的變化其溫度分布不斷變化。根據單元體的能量平衡原理,在直角坐標系下建立的控制微分方程為[12- 13]:

(6)

式中:CP為鑄坯的質量熱容,J/(kg·℃);λ為鑄坯的導熱系數,W/(m·℃);T為鑄坯的溫度,℃;ρ為鑄坯的密度,kg/m3;x、y為鋼坯寬度和厚度方向的坐標。

鑄坯加熱時,鑄坯入爐溫度可認為是初始溫度,即鑄坯斷面初始溫度為:

T(x,y,t)>|t=0=T入爐(x,y)

(7)

鑄坯在加熱爐內的邊界條件主要分為兩種,對流、輻射邊界和絕熱邊界。鑄坯外表面與爐氣接觸,是對流及輻射邊界;鑄坯對稱面是絕熱邊界。由于加熱爐中溫度高達1 200 ℃左右,輻射傳熱量占總傳熱量的90%以上,而加熱爐內爐氣的對流換熱系數很難確定,邊界條件常通過輻射換熱量乘以系數加以補償處理[14- 16]。邊界條件如式(8)、式(9)所示:

(8)

(9)

式中:k1、k2為系數;σgws為導來輻射系數,W/m2·K4,取σgws為=4.85。

設定加熱爐爐膛溫度為1 250 ℃,鑄坯加熱目標溫度為1 200 ℃。計算鑄坯采用三種不同裝爐方式加熱到目標溫度耗用時長以及鑄坯表面和心部溫度的變化曲線。

2 不同熱裝方式鑄坯加熱過程分析

圖2為鑄坯1/4斷面關鍵點位置分布,“1”為鑄坯1/4斷面中心,“7”為寬面中心,“6”為寬面1/4,“5”為角部上端點,“4”為角部中心,“3”為角部下端點,“2”為窄面中心。分析7個關鍵點在加熱爐內的溫度變化曲線,以得到鑄坯加熱溫度變化規律,為加熱爐加熱制度的制定和完善提供依據。冷裝鑄坯要求鑄坯入爐溫度為室溫~400 ℃,此處模擬30 ℃鑄坯在加熱爐內升溫過程。熱裝鑄坯溫度為400~700 ℃,本文模擬400、500、600 ℃鑄坯入爐加熱過程。直裝鑄坯溫度為700~900 ℃。

圖2 鑄坯1/4斷面7個關鍵點示意圖Fig.2 Seven key points at a quarter billet cross section

2.1 鑄坯加熱過程斷面溫度云圖

圖3(a)、3(b)~3(d)和3(e)分別為30 ℃鑄坯、溫裝鑄坯和直裝鑄坯加熱3 600 s的斷面溫度云圖。

圖3 不同熱裝方式鑄坯加熱3 600 s的 斷面溫度云圖Fig.3 Section temperature contour of casting slab under different heat loading conditions after being heated for 3 600 s

采取相同的溫度標尺,經過相同的加熱時間,冷裝坯升溫最為劇烈,中心溫度上升620 ℃,表面中心溫度上升740 ℃,表面1/4位置上升770 ℃,角部溫度上升920 ℃,斷面溫差為300 ℃;400 ℃溫裝鑄坯加熱3 600 s溫度迅速上升,中心溫度上升430 ℃,表面中心溫度上升490 ℃,表面1/4位置上升550 ℃,角部溫度上升640 ℃,斷面溫差為210 ℃;500 ℃溫裝鑄坯加熱3 600 s溫度迅速上升,中心溫度上升390 ℃,表面中心溫度上升450 ℃,表面1/4位置上升480 ℃,角部溫度上升570 ℃,斷面溫差為180 ℃;600 ℃溫裝鑄坯加熱3 600 s,中心溫度上升320 ℃,表面中心溫度上升380 ℃,表面1/4位置上升410 ℃,角部溫度上升500 ℃,斷面溫差為180 ℃。可見冷裝坯和溫裝坯在加熱開始階段溫度上升較為明顯,而冷裝坯升溫速率明顯大于溫裝坯;鑄坯斷面溫差則是冷裝坯最大,最小為600 ℃溫裝坯。從整體溫度分布看,冷裝坯和溫裝坯不同溫度區域呈橢圓環狀分布,從鑄坯角部到鑄坯中心溫度上升。若將30 ℃冷裝坯繼續加熱,將會得到和溫裝坯大致相同的溫度場,即若要達到相同的溫度狀態,冷裝坯將需加熱更長時間。

直裝坯加熱3 600 s斷面溫度較入爐時有較大變化,最低溫度區域位于X正方向鑄坯1/4位置處,圍繞最低溫度區域呈橢圓環狀分布,中心溫度上升90 ℃,表面中心溫度上升260 ℃,表面1/4位置上升250 ℃,角部溫度上升500 ℃,斷面溫差為60 ℃。直裝坯的加熱升溫過程與冷裝和溫裝坯的規律不同,與其加熱初始溫度分布密切相關。直裝坯入爐時中心溫度高,內部熱量向外傳導,中心溫度降低;而爐膛溫度高于鑄坯溫度,因此鑄坯表面溫度迅速上升并向鑄坯內部傳遞熱量,鑄坯溫度由表面向中心逐漸上升。鑄坯中心1/4位置溫度上升的速率低于鑄坯中心和表面的升溫速率,因此成為溫度最低區域。隨著加熱時間的延長,直裝坯的斷面溫度分布趨勢與冷裝、溫裝坯的溫度分布趨勢一致。

2.2 鑄坯加熱過程關鍵點升溫曲線

圖4(a)、4(b)~4(d)和4(e)分別為冷裝坯、溫裝坯和直裝坯在加熱爐內關鍵節點的升溫曲線,圖中平行X軸的直線為Y=1 200,為鑄坯出爐的目標溫度。

整體來看,冷裝與溫裝鑄坯的溫度上升過程趨勢相似,幾個關鍵節點溫度的高低關系基本相同,即角部溫度>窄面溫度>中心溫度。升溫速率:角部位置>窄面位置>中心位置。繼續加熱,幾個節點處溫差逐漸減小,到達加熱終點時溫度基本相同。不同裝爐方式鑄坯各關鍵節點的初始溫度相同,但升溫速率不同。冷裝和溫裝方式下,窄面中心、角部下端點、角部中點、角部上端點的升溫速率明顯高于寬面中心點和寬面1/4點的升溫速率, 斷面中心升溫速率最慢。鑄坯表面與中心的溫差短時間內迅速增大,表面熱量迅速向中心傳遞,鑄坯中心溫度上升。由于鑄坯的熱傳導系數低于加熱爐爐膛與鑄坯表面的換熱系數,因此中心溫度的上升速率小于鑄坯表面溫度的上升速率。

圖4 不同熱裝方式鑄坯加熱關鍵點的升溫曲線Fig.4 Heating curves of the key points of casting slab under different heat loading conditions

圖4(e)所示為直裝坯加熱過程的升溫趨勢。可見加熱初期,斷面中心溫度下降,其他關鍵點溫度迅速升高,窄面中心、角部位置升溫速率高于寬面中心、寬面1/4位置升溫速率。將圖4(e)局部放大如圖4(f)所示,加熱192 s后,斷面中心溫度與寬面中心溫度同為974 ℃。加熱400 s時斷面中心溫度降到最低為968 ℃。加熱1 008 s時斷面中心溫度與角部下端點溫度相同為996 ℃,寬面中心溫度約1 030 ℃,其余位置溫度較低。加熱1 160 s時,斷面中心溫度與角部中心溫度同為1 004 ℃,寬面中心溫度為1 050 ℃,窄面中心溫度986 ℃。加熱1 200 s時,斷面中心溫度與角部上端點溫度同為1 006 ℃,窄面中心991 ℃,其余關鍵點位置溫度較低。加熱1 472 s時,斷面中心溫度與窄面中心溫度同為1 020 ℃,低于其他位置處溫度。斷面中心升溫速率較小,繼續加熱斷面中心溫度會低于其他關鍵位置溫度,而窄面中心和角部位置溫度會大于寬面溫度。繼續加熱至1 520 s時,角部下端點與寬面1/4位置溫度同為1 047 ℃;加熱至2 040 s時,角部下端點溫度與寬面中心溫度相等。加熱至1 704 s時,角部中心溫度與寬面1/4溫度同為1 054 ℃;加熱到2 384 s時,角部中心溫度與寬面中心溫度同為1 095 ℃。加熱1 768 s時,角部上端點溫度與寬面1/4溫度相同為1 056 ℃;加熱到2 408 s時,角部上端點溫度與寬面中心溫度相同為1 096 ℃。加熱到2 024 s時,窄面中心溫度與寬面1/4位置溫度同為1 065 ℃;加熱至2 672 s時,窄面中心溫度與寬面中心溫度同為1 104 ℃。

直裝坯各關鍵點的升溫曲線不同于冷裝和溫裝坯,直裝坯加熱初期各區域溫度高低不一,升溫速率各不相同才會產生各區域溫度在加熱過程中產生交點的情況;對于冷裝和溫裝坯,加熱初期各區域溫度相同,但鑄坯各位置升溫速率不同,各位置溫度的高低趨勢在加熱過程中不改變。

3 結論

(1)冷裝和溫裝鑄坯在加熱爐加熱過程中溫度分布趨勢相似,即以斷面中心低溫區域為中心,向角部高溫區域溫度呈橢圓環狀分布。

(2)直裝坯加熱過程溫度場變化趨勢與冷裝、溫裝坯不同,在加熱初期,直裝坯中心區域為高溫區域,角部為低溫區域,溫度場呈橢圓環狀分布。

(3)直裝坯加熱后中心溫度上升速率小于角部溫度上升速率,一度呈現以X正方向鑄坯1/4位置為低溫區域,呈橢圓環狀溫度向外逐漸升高。繼續加熱,低溫區域轉移至斷面中心,與冷裝和溫裝坯溫度分布趨勢相似。

(4)冷裝和溫裝坯關鍵位置處升溫趨勢相似,而在加熱初期直裝坯各關鍵位置溫度上升趨勢較為復雜,加熱一段時間后,冷裝、溫裝坯溫度分布相似。

(5)提高加熱爐預熱段溫度及鑄坯入爐溫度,縮短鑄坯在高溫區內的停留時間,可以縮小鑄坯斷面溫差,提高加熱效率,從而降低能耗。

[1] 魏軍廣. 帶鋼步進梁式加熱爐節能技術的應用[J]. 山東冶金, 2000, 22(3): 6- 7.

[2] 李永中. 連鑄坯熱裝熱送技術[J].山西冶金, 2007(3): 10- 12.

[3] 于洋,王暢,王林,等.有限元模擬在板坯直裝及熱裝中的應用[J].中國冶金, 2014, 24(1): 41- 46.

[4] 王衛華, 王文軍, 劉洋, 等. 250~320 mm連鑄板坯堆垛緩冷過程仿真計算[J].特殊鋼, 2013, 34(1): 13- 15.

[5] 張炯明,王立峰,王新華,等.板坯連鑄結晶器傳熱系數[J].金屬學報, 2003, 39(12): 1281- 1284.

[6] 鄒達基,鄒宗樹.關于連鑄凝固傳熱數值模擬中鋼液有效導熱系數的探討[J].連鑄, 2009(6): 5- 8.

[7] 鐘婧. 低合金高強度鋼連鑄板坯熱送過程中溫度控制的模擬研究[D]. 重慶:重慶大學,2014.

[8] 甄圣明. 0Cr13不銹鋼板坯連鑄溫度場及應力場的數值模擬[D].包頭:內蒙古科技大學,2008.

[9] 王偉寧, 胡鵬, 任飛飛, 等. 倒角坯熱軋過程角部變形和溫度分布研究[J]. 上海金屬, 2015, 37(2): 57- 62.

[10] CHEN W H, CHUNG Y C, LIU J L. Analysis on energy consumption and performance of reheating furnaces in a hot strip mill[J]. International Communications in Heat & Mass Transfer, 2005, 32(5): 695- 706.

[11] 徐東, 祭程, 朱苗勇, 等. 步進式加熱爐內鋼坯加熱過程的模擬研究[J]. 東北大學學報 (自然科學版), 2013, 34(2): 244- 247.

[12] 青格勒,程素森,楊天鈞,等. 步進梁式加熱爐內的板坯溫度場數值模擬[J].北京科技大學學報, 2004, 26(2): 164- 168.

[13] 李俊洪,林霞,劉勇,等.熱軋板坯加熱溫度場數值模擬及應用[J].鋼鐵, 2009, 44(1): 43- 46.

[14] 朱正海,仇圣桃,廖鵬,等.熱送熱裝過程大板坯傳熱數值模擬研究[J].煉鋼, 2009(5): 43- 47.

[15] 蔣揚虎.武鋼連鑄坯熱送系統傳熱數學模型的研究[D]. 北京:北京科技大學,2000.

[16] 劉新忠,韓靜濤,余萬華,等.步進梁式加熱爐內鋼坯溫度場數值模擬[J].塑性工程學報, 2008, 15(3): 199- 203.

收修改稿日期:2016- 09- 18

StudyonTemperatureDistributioninContinuousCastingSlabunderDifferentHotChargingBehaviors

Ge Jianhua1Wang Minglin1Zhao Wenbo2Zhang Hui1

(1. National Engineering Research Center of Continuous Casting Technology, Central Iron & Steel Research Institute, Beijing 100081,China; 2. College of Metallurgy, Xi’an University of Architecture and Technology, Xi’an Shanxi 710055, China)

Finite element software was used to analyze the heating process of the billet in the walking beam reheating furnace.The model of the slab heating temperature field was established,and the temperature variation curves of the center and the surface of the slab under different heat loading conditions were calculated.Taking the 1/4 section slab as the research object,the temperature distribution in the whole slab section was analyzed,and the variation curves of surface temperature,corner temperature and core temperature were obtained.The heating condition of casting billet with different hot charging modes in heating furnace was simulated and analyzed,and the heating curve of the surface and center of the slab was extracted. The simulation of furnace heating process provided a basis for reducing the temperature difference between the strand section,furnace energy consumption and improving productivity.

hot charging mode,continuous casting slab,heating furnace,temperature distribution,energy consumption

葛建華,男,主要研究方向為連鑄坯熱送熱裝方式及鑄坯熱量損失,電話:15201461577,Email:gthustb@163.com

主站蜘蛛池模板: 国产乱肥老妇精品视频| 免费无遮挡AV| 99精品福利视频| 国产精品刺激对白在线| 国产浮力第一页永久地址| 无码AV高清毛片中国一级毛片| 精品一区二区三区波多野结衣 | 中文字幕av无码不卡免费| 欧美激情视频二区| 日a本亚洲中文在线观看| 久久综合色天堂av| av性天堂网| 国产原创自拍不卡第一页| 亚洲第一黄片大全| 狠狠色香婷婷久久亚洲精品| 精品国产一区91在线| 国产成人亚洲精品色欲AV| 538国产视频| 国产本道久久一区二区三区| 欧美a在线看| 国产精品短篇二区| 亚洲高清国产拍精品26u| 亚洲色偷偷偷鲁综合| 久久国产亚洲偷自| 欧美亚洲欧美区| 日韩精品一区二区三区中文无码| 成人在线视频一区| 国产成人久久777777| 国产精品成人观看视频国产 | 东京热av无码电影一区二区| 高清亚洲欧美在线看| 无遮挡国产高潮视频免费观看| 亚洲欧美自拍中文| 88av在线| 色天堂无毒不卡| 国产欧美又粗又猛又爽老| 亚洲综合九九| 全午夜免费一级毛片| 精品人妻一区无码视频| 久久综合伊人77777| 亚洲中久无码永久在线观看软件| 中文字幕无码中文字幕有码在线| aa级毛片毛片免费观看久| 国产日韩欧美在线播放| 91久久性奴调教国产免费| 国产av无码日韩av无码网站| 亚洲无码高清视频在线观看| 六月婷婷激情综合| 麻豆精品在线播放| 色综合国产| 丁香五月婷婷激情基地| 国产黄色免费看| 亚洲综合第一页| 成年女人18毛片毛片免费| 亚洲男人天堂2020| 日韩国产欧美精品在线| 狠狠色婷婷丁香综合久久韩国| 亚洲系列无码专区偷窥无码| 一本大道东京热无码av| 国产福利2021最新在线观看| 欧美特黄一级大黄录像| 国产精品林美惠子在线观看| 老司国产精品视频91| 中文字幕 欧美日韩| 99成人在线观看| 亚洲国产精品无码AV| 久精品色妇丰满人妻| 一级毛片免费播放视频| 亚洲男人天堂2018| 在线观看热码亚洲av每日更新| 国产视频一区二区在线观看| 亚洲精品片911| 国内黄色精品| 精品无码视频在线观看| 国内黄色精品| 国产网站在线看| 97国产在线视频| 久久永久精品免费视频| 精品黑人一区二区三区| 三区在线视频| 亚洲综合18p| 婷婷色婷婷|