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

平板玻璃大型熔窯碹頂結(jié)構(gòu)受力仿真分析

2024-03-27 08:02:20趙智均高博文宋清超張志勇劉世民1許世清1
玻璃 2024年3期
關(guān)鍵詞:變形分析模型

趙智均 高博文 宋清超 張志勇 劉世民1, 許世清1,

(1.沙河市安全實(shí)業(yè)有限公司創(chuàng)新研發(fā)中心 沙河 054100;2.望美實(shí)業(yè)集團(tuán)有限公司 沙河 054100;3.燕山大學(xué) 秦皇島 066004)

0 引言

隨著平板玻璃行業(yè)節(jié)能減排的要求越來(lái)越高,浮法玻璃生產(chǎn)技術(shù)水平也在不斷提升,很多企業(yè)在技術(shù)改造中,通過(guò)改建大型熔窯生產(chǎn)線來(lái)降低能耗的同時(shí)提高產(chǎn)量,以追求更高的經(jīng)濟(jì)效益。浮法玻璃熔窯一般為橫火焰池窯,熔化部分因熔化能力決定玻璃生產(chǎn)線的產(chǎn)能,玻璃熔窯的噸位受熔化面積和熔化率控制。熔化率反映了熔窯熔化能力的大小,玻璃窯爐設(shè)計(jì)過(guò)程中先選定好熔化率,可進(jìn)行窯池的結(jié)構(gòu)設(shè)計(jì)。在熔窯設(shè)計(jì)過(guò)程中,推導(dǎo)計(jì)算熔化部尺寸與熔化率的關(guān)系有一定的 “經(jīng)驗(yàn)公式”[1]可以使用,根據(jù)公式可計(jì)算出最高1000 t/d典型浮法玻璃熔窯熔化部的主要數(shù)據(jù),但更大噸位下的窯池寬數(shù)據(jù)是否可靠,需要進(jìn)一步進(jìn)行科學(xué)探討。大噸位玻璃熔窯所需的窯池寬度越大,相應(yīng)的大碹跨度也會(huì)增加,為保證熔化部溫度、氣氛等穩(wěn)定,大跨度大碹的中心角也偏大,其碹弧彎曲度會(huì)較大形成更高碹頂?shù)膯?wèn)題,增加大碹的跨度還會(huì)導(dǎo)致大碹的自重增加,加大碹磚之間的應(yīng)力不均問(wèn)題,而且對(duì)碹碴和胸墻的壓力也會(huì)增大,會(huì)影響到整個(gè)熔窯的窯齡。熔化部大碹的使用從理論計(jì)算和結(jié)構(gòu)設(shè)計(jì)上有很多可供參考的力學(xué)體系和工藝參數(shù),但在實(shí)際砌筑過(guò)程中還需要安裝施工人員保證預(yù)組裝、組裝碹胎的可靠性和穩(wěn)定性。本文就利用工程仿真技術(shù)對(duì)平板玻璃大型熔窯碹頂結(jié)構(gòu)的受力情況進(jìn)行仿真分析,評(píng)估碹頂在大跨度下應(yīng)力、變形方面的情況,為后續(xù)大噸位玻璃熔窯的設(shè)計(jì)提供參考。

1 模型建立與仿真模擬分析

1.1 熔化部碹頂模型

浮法玻璃熔窯大碹結(jié)構(gòu)由碹頂、保溫層以及力傳導(dǎo)的鋼結(jié)構(gòu)、拉條、胸墻等組成,熔化部碹頂模型見(jiàn)圖1。碹頂?shù)淖饔檬桥c胸墻、前臉墻組成火焰空間,還可以吸收燃料燃燒時(shí)釋放的熱量,再輻射到玻璃液面上,熔化部的碹頂是由楔形優(yōu)質(zhì)硅磚采用錯(cuò)縫磚砌的方式排列組合,磚縫按照砌筑砂漿的標(biāo)準(zhǔn)進(jìn)行確定。大碹的長(zhǎng)度方向需要進(jìn)行分節(jié),一般至少為三節(jié)以上[2]。砌筑時(shí)碹間預(yù)留100~120 mm的膨脹縫,前后山墻處的碹頂膨脹縫可以多預(yù)留一些。

圖1 熔化部碹頂模型示意圖

碹頂所使用的硅磚主要由鱗石英、方石英、少量殘余石英和玻璃相組成,優(yōu)質(zhì)硅磚的SiO2含量在96%以上,真密度為2.35~2.38 g/cm3,在高溫環(huán)境下也能保持較高的強(qiáng)度,在烤窯階段會(huì)有1.5%~2.2%的體積膨脹,可將硅磚之間的縫隙密合。硅磚的荷重軟化開(kāi)始溫度為1680 ℃,正常窯內(nèi)使用不會(huì)出現(xiàn)問(wèn)題[3]。硅磚屬于酸性材料,在加料過(guò)程中,堿性粉料或堿氣體及其揮發(fā)物在高溫狀態(tài)下會(huì)對(duì)酸性的硅磚進(jìn)行侵蝕。這個(gè)過(guò)程中粉料在耐火材料表面冷凝形成熔體,熔體在硅磚表面磚縫、氣孔處積聚,降低磚的表面熔點(diǎn),從而與磚內(nèi)組分發(fā)生反應(yīng)形成新的玻璃相。在溫度波動(dòng)過(guò)大時(shí)碹頂內(nèi)側(cè)玻璃相會(huì)剝落,導(dǎo)致碹頂變薄變輕結(jié)構(gòu)強(qiáng)度降低,發(fā)生破頂塌陷或者斷磚而影響生產(chǎn)。

1.2 熔化部碹頂受力特點(diǎn)

熔化部碹頂在理想工作狀態(tài)時(shí)的受力情況如圖2所示。

圖2 碹頂受力情況示意圖

圖2中碹頂受到自身重力G,兩側(cè)碹碴對(duì)稱,受到與碹半徑R垂直的斜向支撐力FA和FB,F(xiàn)A和FB分 別可分解為垂直向上的托力FA2、FB2和水平推力FA1、FB1, 為 保 證 平 衡,碹 重 力G與FA2和FB2相抵,方向相反;FA1=FB1,方向相對(duì),與拉條張力相關(guān),進(jìn)而平衡掉兩側(cè)碹碴的水平推力。單獨(dú)的每塊磚均受到三個(gè)外力的作用并相互平衡,將獨(dú)立硅磚的合力作用點(diǎn)相連,連線是作用在硅磚排放布置方向的多段折線,折線相連為一條近似圓弧線。碹頂上的每一列硅磚受的合力作用點(diǎn)都落在這條近似圓弧線上,這條近似圓弧線可以稱為碹頂?shù)摹案G碹推力作用線”[4]。

推力作用線合力作用在碹弧中心線的1/3碹厚范圍內(nèi)最為合理和安全,這樣才不會(huì)對(duì)碹頂產(chǎn)生過(guò)大的壓應(yīng)力,減少出現(xiàn)碹頂硅磚破碎等問(wèn)題。冷態(tài)時(shí)碹頂通過(guò)圖紙?jiān)O(shè)計(jì)和合理的施工布置,可以保證碹頂內(nèi)部推力作用線的穩(wěn)定;熱態(tài)時(shí)硅磚會(huì)出現(xiàn)受熱膨脹擠壓磚縫間砂漿,使得碹頂硅磚相互契合提高整體氣密性和結(jié)構(gòu)完整性。因?yàn)楣璐u的溫差膨脹原因碹弧會(huì)發(fā)生變化,出現(xiàn)整體推力作用線上移現(xiàn)象。這種狀態(tài)下需要靠拉條張緊裝置控制拉條的松緊度,拉條在受熱過(guò)程中出現(xiàn)一定的松弛,需要松開(kāi)一定量減少拉條的張力,平衡因硅磚受熱膨脹而出現(xiàn)的力。

1.3 仿真模擬分析

ANSYS Workbench數(shù)值模擬分析軟件可用來(lái)模擬復(fù)雜的多物理場(chǎng)環(huán)境的實(shí)際工程問(wèn)題,它引入項(xiàng)目流程圖通過(guò)各分析系統(tǒng)間的連接,將分析過(guò)程結(jié)合在一起。分析系統(tǒng)一般需要簡(jiǎn)化假定或真實(shí)的物理模型,將模型網(wǎng)格化分為有限元網(wǎng)格模型,再通過(guò)施加載荷和邊界條件后,運(yùn)行求解得到分析結(jié)果。軟件具有強(qiáng)大的結(jié)構(gòu)、流體、熱、電磁及其相互耦合等分析功能,可以實(shí)現(xiàn)結(jié)構(gòu)靜力分析、結(jié)構(gòu)動(dòng)力學(xué)分析、結(jié)構(gòu)熱分析和流體動(dòng)力學(xué)等分析技術(shù)[5]。

碹頂模擬計(jì)算利用ANSYS Workbench軟件包括3個(gè)過(guò)程:前處理,創(chuàng)建幾何模型并劃分網(wǎng)格;加載和求解,施加材料載荷、邊界條件并求解計(jì)算;后處理,基本數(shù)據(jù)和導(dǎo)出數(shù)據(jù)[6]。

本文在模型建立過(guò)程中進(jìn)行了一定的簡(jiǎn)化假定,整體式3D幾何模型如圖3所示,假定砂漿彌散于整個(gè)硅磚單元中,將碹頂簡(jiǎn)化為各項(xiàng)同性的均勻整體,整體式模型的各種應(yīng)力狀態(tài)體現(xiàn)在結(jié)構(gòu)主軸上,圖4為幾何模型自動(dòng)網(wǎng)格劃分示意圖。為比較大型熔窯碹頂結(jié)構(gòu)差異化,建立了中心角60° 、碹厚500 mm、跨度14、15和16 m的三個(gè)尺寸的局部碹頂模型并分別進(jìn)行計(jì)算。

圖3 幾何模型圖

圖4 幾何模型自動(dòng)網(wǎng)格劃分示意圖

加載和求解部分先確定碹頂硅磚參數(shù),如彈性模量、泊松比、密度、膨脹系數(shù)和比熱等材料特性再進(jìn)行計(jì)算。對(duì)碹頂力的加載按照1.2中提到的碹頂斜向支撐力FA和FB的分力方向進(jìn)行施加。溫度通過(guò)以施工時(shí)的冷態(tài)溫度和工作熱態(tài)溫度兩種狀態(tài)進(jìn)行控制。后處理部分對(duì)不同尺寸的碹頂結(jié)構(gòu)進(jìn)行總變形和等效應(yīng)力的數(shù)據(jù)采集和導(dǎo)出,進(jìn)而討論分析。

2 結(jié)果與討論

2.1 碹頂總變形分析

2.1.1 冷態(tài)碹頂總變形情況分析

圖5中(a)~(c)是冷態(tài)22 ℃時(shí)跨度分別為14、15和16 m的碹頂總變形變化情況云圖。

圖5 冷態(tài)時(shí)不同跨度碹頂總變形圖

云圖顯示變形量從中部向兩側(cè)對(duì)稱性的逐漸減少,最大的變形情況出現(xiàn)在碹頂?shù)闹虚g位置,最大變形量數(shù)值分別為1.8721e-5m、3.7677e-5m、5.925e-5m,跨度增加最大變形量數(shù)值也逐漸增大,說(shuō)明跨度的增加導(dǎo)致碹頂更容易出現(xiàn)因自重加大而產(chǎn)生變形的情況。最小變形量出現(xiàn)在碹碴支撐的位置,主要是平衡水平推力和自重較為穩(wěn)定,所以變形量小。

按照最大變形量到1.65e-5m變形量區(qū)間劃定變形量區(qū)域,可以找出不同跨度此區(qū)域的節(jié)點(diǎn)數(shù)量分別為21041、53327和67321。隨跨度的增大變形量區(qū)域逐漸增加,變形量區(qū)域占整體結(jié)構(gòu)區(qū)域比例分別為23.2%、55.1%和65.4%,最大變形量隨跨度的增大會(huì)增加約2.02e-5m,增大碹頂結(jié)構(gòu)變形的風(fēng)險(xiǎn)。

2.1.2 熱態(tài)碹頂總變形情況分析

圖6中(a)~(c)是熱態(tài)1600 ℃時(shí)跨度分別為14、15和16 m的碹頂總變形變化情況云圖。

圖6 熱態(tài)時(shí)不同跨度碹頂總變形圖

云圖顯示變形量與冷態(tài)碹頂變形量分布趨勢(shì)相同,均是從中部向兩側(cè)對(duì)稱性的逐漸減少,最大的變形情況出現(xiàn)在碹頂?shù)闹虚g位置,隨著跨度的增加最大變形量數(shù)值逐漸增大,最小變形量還是出現(xiàn)在碹碴支撐的位置。但與冷態(tài)對(duì)比,在模擬玻璃熔窯內(nèi)高溫條件下,最大變形量較冷態(tài)提高了5個(gè)數(shù)量級(jí),碹頂受熱膨脹從而形成了更大程度的變形,這與實(shí)際溫度情況時(shí)碹頂中部的鎖磚位置有較大抬升的情況相符。按照最大變形量到0.46 m變形量區(qū)間劃定變形量區(qū)域,可以找出不同跨度此區(qū)域節(jié)點(diǎn)數(shù)量分別為21001、28563和35219,此區(qū)域占整體結(jié)構(gòu)區(qū)域比例分別為23.2%、29.5%和34.2%,最大變形量隨跨度的增大會(huì)增加約0.037 m,碹頂會(huì)出現(xiàn)更大面積的內(nèi)部硅磚相互擠壓?jiǎn)栴},硅磚更易出現(xiàn)裂紋或破碎,從而危害到整體碹頂?shù)膲勖?/p>

2.2 碹頂?shù)刃?yīng)力分析

2.2.1 冷態(tài)碹頂?shù)刃?yīng)力情況分析

圖7中(a)~(c)是冷態(tài)22 ℃時(shí)跨度分別為14、15和16 m的碹頂?shù)刃?yīng)力變化情況云圖。

圖7 冷態(tài)時(shí)不同跨度碹頂?shù)刃?yīng)力圖

云圖顯示整個(gè)碹頂結(jié)構(gòu)積聚的應(yīng)力主要為正值的壓應(yīng)力,較小的壓應(yīng)力主要集中在碹厚的中部下表面及其兩側(cè)上表面;較大的壓應(yīng)力主要集中在碹頂中部碹厚的2/3和碹角位置。隨著跨度的增加,最大應(yīng)力值和最小應(yīng)力值都逐漸增大,說(shuō)明自重增加時(shí),拉條的水平推力平衡分力作用減弱,使得壓應(yīng)力集中位置下移至碹角位置。

圖8顯示的是冷態(tài)碹頂?shù)撞康刃?yīng)力局部分布情況云圖。

圖8 冷態(tài)碹頂?shù)撞康刃?yīng)力局部放大圖

圖中顯示出碹頂?shù)膽?yīng)力集中在位于碹頂兩側(cè)的底部以及下表面,碹角位置的最大應(yīng)力值隨跨度增加分別為1.3827e5Pa、2.3837e5Pa和3.3191e5Pa,本研究模型建立碹角磚使用鋯英石磚,為中性耐火磚,常溫耐壓強(qiáng)度可達(dá)到100 MPa。碹頂頂部的應(yīng)力值分布隨跨度增加分別為71332 Pa、1.2893e5Pa和1.8638e5Pa,此位置的優(yōu)質(zhì)硅磚常溫耐壓強(qiáng)度為34 ~ 40 MPa,綜合數(shù)據(jù)可以看出碹頂用磚可以滿足條件。下表面壓應(yīng)力的集中會(huì)導(dǎo)致硅磚之間擠壓程度變大,在升溫過(guò)程會(huì)出現(xiàn)倒“V”字形裂縫。通過(guò)對(duì)比不同碹頂?shù)膽?yīng)力變化情況,可以發(fā)現(xiàn)隨跨度的增加,下表面壓應(yīng)力集中區(qū)域向碹頂兩側(cè)的底部靠近,出現(xiàn)中部支撐力減少,底部應(yīng)力過(guò)于集中現(xiàn)象。在這種現(xiàn)象下出現(xiàn)裂縫的位置也會(huì)逐漸集中在碹頂兩側(cè)的底部,如果施工過(guò)程控制不當(dāng),以及碹頂跨度增加自重變大等問(wèn)題,碹頂?shù)姆€(wěn)定性就會(huì)受到影響,從而出現(xiàn)垮塌的危險(xiǎn)。

2.2.2 熱態(tài)碹頂?shù)刃?yīng)力情況分析

圖9中(a)~(c)是熱態(tài)1600 ℃時(shí)跨度分別為14、15和16 m的碹頂?shù)刃?yīng)力變化情況云圖。

圖9 熱態(tài)時(shí)不同跨度碹頂?shù)刃?yīng)力圖

云圖顯示整個(gè)碹頂結(jié)構(gòu)積聚的應(yīng)力主要為正值的壓應(yīng)力,較小的壓應(yīng)力主要集中在碹厚的中間位置;較大的壓應(yīng)力主要集中在碹角兩側(cè)底部。在工作溫度下碹頂硅磚受熱膨脹,使硅磚之間的砌筑更加緊密,整體應(yīng)力值分布也變得均勻。碹頂?shù)膽?yīng)力集中在位于碹頂兩側(cè)的底部,對(duì)碹角磚及立柱的推力方向呈傾斜向下,此時(shí)碹角頂鐵與碹角磚的接觸面是斜面向上形狀。集中在碹厚中部的應(yīng)力值隨跨度增加分別為4.5999e5Pa、8.7099e5Pa和1.2737e6Pa,仍然在硅磚的應(yīng)力范圍(34 MPa)內(nèi)。

圖10顯示的是熱態(tài)碹頂?shù)撞康刃?yīng)力局部分布情況云圖。

圖10 熱態(tài)碹頂?shù)撞康刃?yīng)力局部放大圖

為研究極端情況,溫度設(shè)置接近硅磚的荷重軟化溫度1680 ℃,發(fā)現(xiàn)碹角底部的應(yīng)力值在5.30e9Pa以上,超過(guò)底部碹角的鋯英石磚耐壓強(qiáng)度,就會(huì)出現(xiàn)碹角磚的碎裂失效,導(dǎo)致碹頂坍塌,但浮法玻璃熔窯內(nèi)部溫度通常不會(huì)達(dá)到這個(gè)溫度,但極端實(shí)驗(yàn)也得出要注意碹頂兩側(cè)底部的耐火材料選用,而且烤窯過(guò)程中的升溫也需要對(duì)拉條進(jìn)行放松操作,減少碹角頂鐵和碹角磚契合位置發(fā)生較大移動(dòng)導(dǎo)致應(yīng)力過(guò)于集中的問(wèn)題。隨著跨度的增加中部的下表面應(yīng)力集中面積逐漸增大,而且中部位置受到窯內(nèi)溫度和堿性粉料的侵蝕可能性也更大,更易出現(xiàn)損壞。

2.3 碹頂極端情況分析

為研究碹頂因?yàn)槭芰Φ牟黄胶舛霈F(xiàn)坍塌的問(wèn)題,設(shè)置跨度14 m和16 m碹頂?shù)撞坷瓧l力減少,碹頂中部因自重向下坍塌的模型,并進(jìn)行對(duì)比分析討論。

圖11和圖12是冷態(tài)22 ℃時(shí)跨度分別為14 m和16 m的碹頂總變形和等效應(yīng)力變化情況對(duì)比云圖,拉條的水平推力減少80%時(shí),碹頂出現(xiàn)坍塌的位置分別從總變形2.21e-6m和3.60e-6m位置開(kāi)始,此位置的等效應(yīng)力值分別為16872 Pa和22749 Pa,出現(xiàn)坍塌的范圍為整個(gè)碹頂結(jié)構(gòu)的90.3%和95.5%。在拉條松開(kāi)水平推力同時(shí)減少的情況下,跨度增加2 m時(shí),坍塌程度增加了5%左右,說(shuō)明在碹頂砌筑過(guò)程中拉條的張緊控制要隨跨度的增加而增大,而且拉條力不足時(shí)會(huì)出現(xiàn)碹頂坍塌問(wèn)題。

圖11 冷態(tài)時(shí)跨度14 m與跨度16 m碹頂總變形圖

圖12 冷態(tài)時(shí)跨度14 m與跨度16 m碹頂?shù)刃?yīng)力圖

3 結(jié)語(yǔ)

(1)總變形結(jié)果顯示,碹頂變形量從中部向兩側(cè)對(duì)稱性的逐漸減少。對(duì)比冷態(tài)22 ℃與熱態(tài)1600 ℃的總變形量,熱態(tài)碹頂受熱膨脹形成了更大程度的變形,冷態(tài)時(shí)跨度每增加1 m最大變形量增加約2.02e-5m,熱態(tài)時(shí)增加約0.37 m。

(2)碹頂內(nèi)部應(yīng)力結(jié)果顯示,最大應(yīng)力值位置主要集中在碹頂中部碹厚的2/3和兩側(cè)底部位置,隨著跨度的增加,最大應(yīng)力值增加。在熔窯設(shè)計(jì)過(guò)程要注意碹頂兩側(cè)底部耐火材料的選用,而且烤窯過(guò)程中的升溫也需要對(duì)拉條進(jìn)行放松操作,減少碹角頂鐵和碹角磚契合位置發(fā)生較大移動(dòng)應(yīng)力過(guò)于集中的問(wèn)題。

(3)在雙碳目標(biāo)的政策下,大跨度的浮法玻璃熔窯生產(chǎn)線要考慮熔化面積和成本的投入,實(shí)際結(jié)構(gòu)設(shè)計(jì)時(shí)碹厚和碹頂中心角都應(yīng)該與跨度相匹配,在合理的結(jié)構(gòu)設(shè)計(jì)下對(duì)施工、硅磚質(zhì)量以及熔化的燃燒制度等提出更高要求。

猜你喜歡
變形分析模型
一半模型
隱蔽失效適航要求符合性驗(yàn)證分析
重要模型『一線三等角』
談詩(shī)的變形
重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
電力系統(tǒng)不平衡分析
電子制作(2018年18期)2018-11-14 01:48:24
“我”的變形計(jì)
例談拼圖與整式變形
會(huì)變形的餅
電力系統(tǒng)及其自動(dòng)化發(fā)展趨勢(shì)分析
主站蜘蛛池模板: 成人夜夜嗨| 中国精品自拍| 亚洲第一成年人网站| 国产精品短篇二区| 黄色网在线| 国产制服丝袜91在线| 亚欧乱色视频网站大全| 女人一级毛片| 久久综合国产乱子免费| 成人福利在线视频| 黄片一区二区三区| 亚洲无码电影| 亚洲一区二区三区香蕉| 国产青青操| 亚洲自偷自拍另类小说| 黄片一区二区三区| 国产原创演绎剧情有字幕的| 中日无码在线观看| 成年av福利永久免费观看| 国产精品主播| 国产成人精品亚洲日本对白优播| 久久国产精品麻豆系列| 亚洲欧美成aⅴ人在线观看| 国产色偷丝袜婷婷无码麻豆制服| 欧美成人手机在线观看网址| 日韩精品无码不卡无码| 一边摸一边做爽的视频17国产| 谁有在线观看日韩亚洲最新视频| 18禁影院亚洲专区| 亚洲一区二区无码视频| 国产精品乱偷免费视频| 热这里只有精品国产热门精品| 日本国产在线| 成人年鲁鲁在线观看视频| 在线观看av永久| www欧美在线观看| 97se亚洲综合在线| 午夜精品区| 日韩无码视频网站| 国产精品入口麻豆| 视频一本大道香蕉久在线播放| 日韩一二三区视频精品| 蜜桃臀无码内射一区二区三区| 2018日日摸夜夜添狠狠躁| 国产97区一区二区三区无码| 国产成人久久777777| 国产高颜值露脸在线观看| 亚洲无限乱码一二三四区| 色国产视频| 波多野结衣一区二区三视频| 最新亚洲人成无码网站欣赏网| 欧美成人免费一区在线播放| 国产欧美日韩专区发布| 欧美日韩精品一区二区在线线| 91小视频在线观看| 国产成人无码播放| 国产精品一区在线麻豆| 激情综合网激情综合| 国产女人爽到高潮的免费视频| 在线播放真实国产乱子伦| 国产91av在线| 久久综合婷婷| 国产九九精品视频| 色香蕉网站| 亚洲永久色| 国产91丝袜| 亚洲热线99精品视频| 伊人精品视频免费在线| 四虎永久在线| 亚洲男人的天堂久久香蕉 | 美女无遮挡拍拍拍免费视频| 亚洲精品视频网| 色135综合网| 亚洲三级视频在线观看| 亚洲无码视频图片| 制服丝袜在线视频香蕉| 午夜性刺激在线观看免费| 欧美精品亚洲精品日韩专区va| 亚洲日韩第九十九页| 久久久久久久久久国产精品| 精品91视频| 国产无码在线调教|