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

煤巖破落過(guò)程中螺旋滾筒的可靠性研究*

2017-04-16 01:37:51趙麗娟樊志海周文潮
關(guān)鍵詞:采煤機(jī)

趙麗娟,樊志海,周文潮

(遼寧工程技術(shù)大學(xué) 機(jī)械工程學(xué)院,遼寧 阜新 123000)

數(shù)字出版日期: 2017-07-14

0 引言

螺旋滾筒是采煤機(jī)的工作機(jī)構(gòu),其承擔(dān)著破煤、裝煤及除塵等任務(wù),采煤機(jī)裝機(jī)功率絕大部分消耗在螺旋滾筒截割煤巖過(guò)程當(dāng)中[1]。因此,螺旋滾筒的設(shè)計(jì)是否合理將直接影響截齒受力及其波動(dòng)、截割比能耗、塊煤率、裝煤率及生產(chǎn)效率等性能指標(biāo)[2-7]。通過(guò)離散元方法[8-10]對(duì)滾筒的裝煤性能進(jìn)行分析,研究采煤機(jī)運(yùn)動(dòng)學(xué)參數(shù)對(duì)滾筒裝煤性能以及對(duì)螺旋滾筒可靠性的影響,不僅大大的降低了企業(yè)的研發(fā)成本,縮短其產(chǎn)品設(shè)計(jì)時(shí)間,提高滾筒的截煤和裝煤效率,也可以更好的提升煤礦機(jī)械生產(chǎn)的安全性,對(duì)提升企業(yè)的經(jīng)濟(jì)效益具有重要意義[11-15]。

以某新型采煤機(jī)截割部為研究對(duì)象,利用EDEM對(duì)螺旋滾筒的裝煤性能進(jìn)行了數(shù)值模擬。分析了不同牽引速度下裝煤率的變化規(guī)律,得到了不同工況下螺旋滾筒受到的非線(xiàn)性沖擊載荷曲線(xiàn),并對(duì)沖擊載荷作用下滾筒的動(dòng)應(yīng)力分布及其工作的可靠性進(jìn)行了研究,分析結(jié)果為薄煤層采煤機(jī)螺旋滾筒的設(shè)計(jì)及采煤機(jī)的定型生產(chǎn)提供了參考。

1 煤巖性質(zhì)的測(cè)定及離散元仿真模型的建立

1.1 煤巖性質(zhì)的測(cè)定

煤巖參數(shù)準(zhǔn)確與否直接影響著仿真結(jié)果的可靠性,因此對(duì)某煤礦5-1和5-2煤巖試樣進(jìn)行標(biāo)準(zhǔn)化測(cè)試得出其煤巖的物理力學(xué)性質(zhì),實(shí)驗(yàn)所用煤巖樣本及部分測(cè)試系統(tǒng)分別如圖1所示,利用DQ-1型巖石切割機(jī)將煤巖切割長(zhǎng)方體,而后分別進(jìn)行了單軸抗壓強(qiáng)度試驗(yàn)、三軸試驗(yàn)以及硬度塑性系數(shù)試驗(yàn),利用型微機(jī)控制電子試驗(yàn)機(jī)檢測(cè)切割試樣的物理力學(xué)性質(zhì)參數(shù)并將其傳到計(jì)算機(jī)上,得到的煤巖各項(xiàng)參數(shù)如表1所示。

圖1 煤樣性質(zhì)測(cè)試Fig.1 The property test of coal sample

參數(shù)名稱(chēng)5-15-2真密度/(kg·m-3)13321319孔隙率/%9.5310.28彈性模量/MPa43884112泊松比μ0.230.24抗壓強(qiáng)度/MPa17.7115.78抗拉強(qiáng)度/MPa1.080.83黏聚力/MPa1.851.45內(nèi)摩擦角/(°)5958堅(jiān)固性系數(shù)2.01.9

1.2 螺旋滾筒裝煤模型的建立

首先定義模型之間粘結(jié)模型,包括設(shè)置煤巖顆粒與煤巖顆粒粘結(jié)模型為Hertz-Mindlin with bonding、煤巖顆粒與采煤機(jī)滾筒的粘結(jié)模型為Hertz-Mindlin(no slip)、煤巖顆粒體積力的類(lèi)型為Electrostatics。對(duì)材料屬性進(jìn)行賦值,煤巖顆粒的泊松比及密度分別為0.235及1.325 5 g/cm3,剪切模量為1 720 MPa。滾筒材料的泊松比、密度、剪切模量分別為0.31,7.85 g/cm3和70 000 MPa。定義煤巖顆粒之間恢復(fù)系數(shù)、靜摩擦系數(shù)、動(dòng)摩擦系數(shù)分別為0.5,0.8及0.58;煤巖顆粒與滾筒之間恢復(fù)系數(shù)、靜摩擦系數(shù)、動(dòng)摩擦系數(shù)分別為0.5,0.9及0.58。

對(duì)顆粒之間的黏結(jié)平行鍵進(jìn)行設(shè)置,設(shè)置的主要參數(shù)涉及法向剛度kn,切向剛度ks,法向應(yīng)力σ、切向應(yīng)力τ以及黏結(jié)半徑。其中,kn,ks可通過(guò)Hertz接觸理論及Mindlin的研究成果按公式(1),(2)進(jìn)行計(jì)算:

(1)

ks=λkn

(2)

式中:kn為法向剛度,N/m;ks為切向剛度,N/m;ν1,ν2為2個(gè)顆粒的泊松比;E1,E2為2個(gè)顆粒的彈性模量,MPa;r1,r2為2個(gè)顆粒的半徑,mm;λ為剛度系數(shù)折算值,一般取1/2~2/3。σ及τ可通過(guò)莫爾-庫(kù)倫理論計(jì)算得到,當(dāng)應(yīng)力值超過(guò)σ或τ,則發(fā)生相應(yīng)的拉壓破壞或剪切破壞[9],如公式(3)所示。

(3)

式中:σ為破壞面上的法向應(yīng)力,MPa;τ為破壞面的剪切應(yīng)力,MPa;σ1為最大主應(yīng)力,MPa;σ3為最小主應(yīng)力,MPa;α為剪切破壞角,(°);φ為內(nèi)摩擦角,(°);C為煤巖的黏聚力,MPa。

而其中的σ1,σ3可通過(guò)McClintock和Walsh修正后的格里菲斯公式(4)進(jìn)行計(jì)算。

(4)

式中:σt為材料的抗拉強(qiáng)度,MPa;σc為材料的抗壓強(qiáng)度,MPa;μ是裂隙間的摩擦系數(shù)。

將表1中相關(guān)的煤巖力學(xué)物理性能參數(shù)代入公式(1)~(4)中,計(jì)算出kn,ks,σ,τ分別為9.5×109N/m,5.5×109N/m,0.955 MPa以及3.21 MPa。根據(jù)仿真模型需要,設(shè)置煤巖顆粒的半徑設(shè)置為10 mm。在Geometry界面中經(jīng)過(guò)空間壓縮等系列的操作最后生成的煤壁長(zhǎng)寬高為12 400 mm×880 mm×1 150 mm,顆粒的總質(zhì)量為7 834 kg。形成煤巖顆粒煤壁。之后將Pro/E中建立好的采煤機(jī)截割部三維實(shí)體模型存成iges格式后導(dǎo)入EDEM中,通過(guò)調(diào)整滾筒的旋轉(zhuǎn)約束,可實(shí)現(xiàn)對(duì)煤壁截割裝煤的模擬。建立好的EDEM仿真模型如圖2所示,建立2個(gè)區(qū)域,其中統(tǒng)計(jì)區(qū)Ⅰ為未落到刮板輸送機(jī)上的煤巖顆粒數(shù),統(tǒng)計(jì)區(qū)Ⅱ?yàn)榻芈涞焦伟遢斔蜋C(jī)上的煤巖顆粒數(shù),則裝煤率可通過(guò)Ⅱ區(qū)顆粒累計(jì)質(zhì)量與Ⅰ、Ⅱ區(qū)之和的比進(jìn)行計(jì)算。

圖2 截割部的EDEM仿真模型Fig.2 Cutting of EDEM simulation model

2 螺旋滾筒裝煤過(guò)程仿真

2.1 螺旋滾筒裝煤率分析

采煤機(jī)落煤量與其牽引速度有著直接的影響關(guān)系,牽引速度對(duì)采煤機(jī)滾筒的截割性能和滾筒葉片空間內(nèi)煤流量都有影響,為了得出牽引速度與滾筒裝煤效果之間的關(guān)系,選擇滾筒運(yùn)動(dòng)學(xué)參數(shù)如下:滾筒的截深為800 mm,滾筒轉(zhuǎn)速為58 m/min,仿真時(shí)間為2 s,牽引速度依次設(shè)置為6,7,8,9,10 m/min。牽引速度為8 m/min時(shí)的螺旋滾筒的裝煤情況如圖3所示。通過(guò)對(duì)不同牽引速度下滾筒拋射的裝煤率進(jìn)行仿真,最后得到的各個(gè)工況下的裝煤率如表2所示。

圖3 牽引速度8 m/min下顆粒速度Fig.3 Under the traction speed is 8 m/min speed particle cloud

牽引速度/(m·min-1)678910裝煤率/%58.760.861.359.657.5

牽引速度逐漸增大的過(guò)程中,煤巖體在滾筒中分布逐漸增多,煤巖體沿螺旋葉片向輸送帶端運(yùn)移的能力增強(qiáng),裝煤率逐漸上升;當(dāng)牽引速度繼續(xù)增大,落煤量超過(guò)了滾筒容積,部分煤巖掉落到了統(tǒng)計(jì)區(qū)Ⅰ致使裝煤率下降。

2.2 螺旋滾筒載荷的提取與分析

通過(guò)菜單欄中的Test data按鈕將截割煤壁的過(guò)程受到的載荷提取出來(lái),牽引速度為8 m/min的三向力如圖4(a)可知,截割阻力(Z向)最大,其次是牽引阻力(X向),側(cè)向力(Y方向)最小。滾筒Z向切削合力方向與滾筒轉(zhuǎn)動(dòng)方向相反,滾筒X向截齒合力與牽引方向相反,滾筒的Y向合力在零線(xiàn)上下波動(dòng),但其平均值不為零。牽引速度8 m/min螺旋葉片受到裝煤反力曲線(xiàn)如圖4(b)所示, 0~0.65 s期間,裝煤反力逐漸增大,這是由于螺旋葉片起到裝煤作用,隨著滾筒轉(zhuǎn)動(dòng)與進(jìn)給,參與截割煤巖的截齒數(shù)量將增大,截掉的煤巖質(zhì)量也將增大。當(dāng)工作到0.65 s之后,滾筒處于穩(wěn)定截割,裝煤反力將不再持續(xù)增大。

圖4 螺旋滾筒受到的載荷Fig.4 The load of spiral drum

圖5 滾筒載荷理論值與仿真值對(duì)比Fig.5 Drum load comparison of theoretical value and simulation value

利用數(shù)值模擬技術(shù),得到了不同牽引速度下螺旋滾筒受到的載荷,并基于破煤理論分別將相同工況下計(jì)算所得數(shù)據(jù)與本次仿真所得數(shù)據(jù)進(jìn)行比較(如圖5),以驗(yàn)證本次模擬的正確性,X向載荷最大誤差為8.968%,Y向載荷最大誤差為9.793%,Z載荷最大誤差為9.482%,所有誤差均小于10%在誤差允許范圍內(nèi),故仿真結(jié)果可靠。

3 不同牽引速度對(duì)裝煤過(guò)程中螺旋滾筒的可靠性分析

3.1 剛?cè)狁詈夏P偷慕?/h3>

在軟件PRO/E中建立采煤機(jī)主要零部件并進(jìn)行虛擬裝配與干涉檢查,將采煤機(jī)三維實(shí)體模型通過(guò)接口軟件Mechanism/Pro導(dǎo)入到ADAMS中,設(shè)置零件密度、材料,添加約束等,同時(shí)在ANSYS中建立中性(mnf)文件,在ADAMS中進(jìn)行柔性件替換,然后驗(yàn)證模型是否有冗余約束。將EDEM仿真得到的動(dòng)態(tài)載荷保存成TXT文本數(shù)據(jù)樣條曲線(xiàn)導(dǎo)入到ADAMS中進(jìn)行加載,設(shè)置參數(shù),進(jìn)行仿真計(jì)算。采煤機(jī)截割部剛?cè)狁詈夏P?,如圖6所示。

圖6 采煤機(jī)剛?cè)狁詈咸摂M樣機(jī)模型Fig.6 Rigid-flex coupled virtual prototyping model of shearer

3.2 螺旋滾筒的可靠性分析

合金頭與齒體截割煤巖時(shí)應(yīng)力分布如圖7和圖8所示。由圖7可以看出,合金頭的最大應(yīng)力主要集中在齒尖的局部接觸區(qū)域,在前刀面及其兩側(cè)呈現(xiàn)非對(duì)稱(chēng)性分布,這說(shuō)明合金頭主要以磨損失效為主,且呈現(xiàn)非對(duì)稱(chēng)性磨損。截割過(guò)程中合金頭長(zhǎng)期處于高應(yīng)力狀態(tài)下,使得合金頭迅速磨鈍,磨鈍后的截齒會(huì)大大增加其工作機(jī)構(gòu)的截割阻力,進(jìn)而影響到采煤機(jī)整機(jī)的工作性能。位于齒尖處所受應(yīng)力明顯比前刀面兩側(cè)所受應(yīng)力大。由圖8可以看出,齒體高應(yīng)力基本集中在其頂端及前刃面、齒柄頭部的軸肩處,且應(yīng)力值遠(yuǎn)遠(yuǎn)小于合金頭的應(yīng)力。位于齒柄頭部的軸肩處,受到應(yīng)力較大,齒體有可能由于應(yīng)力過(guò)大而斷裂;位于齒體錐段前刃面的應(yīng)力波動(dòng)比較平穩(wěn)。由表3可知,牽引速度為10 m/min時(shí)中合金頭與齒體受到的應(yīng)力最大,受到的載荷變大,這是由于隨著牽引速度的增加,單位時(shí)間內(nèi)參與截割截齒的切削厚度也會(huì)增加,使截齒受到的載荷,從而導(dǎo)致截齒的應(yīng)力增大,即參與截割截齒的合金頭與齒體的應(yīng)力均有不同程度的增加。

圖7 合金頭的應(yīng)力Fig.7 The stress of the alloy head

圖8 齒體的應(yīng)力Fig.8 The stress maximum unit

牽引速度/(m·min-1)合金頭應(yīng)力值/MPa合金頭安全系數(shù)齒體應(yīng)力值/MPa齒體安全系數(shù)6970.431.799669.051.86571061.121.560707.581.73881119.571.877765.482.54791209.561.723812.572.323101321.831.598872.302.098

截割煤巖時(shí)螺旋滾筒的應(yīng)力分布如圖9所示。螺旋滾筒上的應(yīng)力主要集中于工作截齒的齒座根部,并且端盤(pán)上齒座受到的應(yīng)力明顯大于葉片上齒座受到的應(yīng)力,并且螺旋滾筒上工作截齒的齒座根部的應(yīng)力遠(yuǎn)大于其他部位的應(yīng)力。齒座受到的應(yīng)力變化規(guī)律與截齒受到的應(yīng)力變化類(lèi)似,這是由于運(yùn)動(dòng)參數(shù)的不同匹配改變了單位時(shí)間內(nèi)參與截割截齒的切削厚度,從而改變了截齒受到的載荷,進(jìn)而影響了齒座的受力。由截齒齒座的結(jié)構(gòu)決定其薄弱位置處于根部,故齒座根部易發(fā)生損傷、破壞斷裂等失效形式。

圖9 齒座與葉片的應(yīng)力云圖分布Fig.9 Drum stress cloud distribution

葉片的應(yīng)力變化規(guī)律與齒座的應(yīng)力變化不同,由圖9及表4可得,工況為牽引速度10 m/min時(shí)葉片的應(yīng)力值最大,最大值為131.52 MPa,小于葉片的屈服強(qiáng)度,說(shuō)明葉片的強(qiáng)度可以滿(mǎn)足使用要求。當(dāng)滾筒牽引速度增大時(shí),葉片的應(yīng)力隨運(yùn)動(dòng)參數(shù)的增大而變大。這是由于運(yùn)動(dòng)參數(shù)的不同匹配改變了單位時(shí)間內(nèi)參與截割截齒的切削厚度,從而改變了單位時(shí)間截掉煤巖的質(zhì)量,牽引速度越大單位時(shí)間內(nèi)截掉的煤巖質(zhì)量越大,截掉的煤巖由螺旋葉片裝到運(yùn)輸機(jī)上,轉(zhuǎn)速越高對(duì)裝煤效果越明顯,但對(duì)葉片的沖擊越大。由圖10可知當(dāng)滾筒其他影響因素不變的情況下,牽引速度越小,截掉煤巖的質(zhì)量也較少,滾筒葉片內(nèi)的部分煤巖將向葉片外緣運(yùn)動(dòng),不能及時(shí)裝到運(yùn)輸機(jī)上,不利于滾筒裝煤,同時(shí)對(duì)葉片沖擊較小。牽引速度增大到某一值后,單位時(shí)間內(nèi)采煤機(jī)截掉的煤巖量大于滾筒葉片內(nèi)的存煤空間大小,煤巖流動(dòng)將發(fā)生阻塞,不能及時(shí)排出,大大減小了滾筒的裝煤率,并且葉片易造成破壞。結(jié)合煤層的賦存條件和采煤機(jī)工作的可靠性,該采煤機(jī)的牽引速度應(yīng)控制在8 m/min左右。因此,運(yùn)動(dòng)參數(shù)的合理匹配對(duì)滾筒可靠性至關(guān)重要。

表4 齒座、葉片的應(yīng)力值及安全系數(shù)

圖10 牽引速度、葉片應(yīng)力、裝煤率三者的關(guān)系Fig.10 The three relationship between traction speed, >blade stress,loading rate

4 結(jié)論

1)利用EDEM軟件建立了采煤機(jī)截割部破落煤巖的耦合模型,得到了裝煤率隨牽引速度增大過(guò)程中,裝煤率出現(xiàn)先增大后減小的規(guī)律。

2)提取了不同牽引速度下截煤過(guò)程中螺旋滾筒受到得三向力與葉片受到的裝煤反力,并將其對(duì)剛?cè)狁詈夏P瓦M(jìn)行加載,通過(guò)仿真得到了螺旋滾筒受到的應(yīng)力隨牽引速度的增大而增大,在保證采煤機(jī)可靠性的前提下,牽引速度控制在8 m/min左右有較高的裝煤率。

3)通過(guò)協(xié)同仿真解決了EDEM離散元軟件運(yùn)用到采煤機(jī)滾筒的裝煤效率與煤巖破落過(guò)程中螺旋滾筒的可靠性問(wèn)題,為研究采煤機(jī)裝煤性能方面提供了一種新的仿真方法并可以大幅度的降低企業(yè)研發(fā)成本,更好的增加企業(yè)的經(jīng)濟(jì)效益。

[1] 劉春生,于信偉,任昌玉.滾筒式采煤機(jī)工作機(jī)構(gòu)[M].哈爾濱:哈爾濱工程大學(xué)出版社,2010:4-21.

[2] 劉送永,杜長(zhǎng)龍,崔新霞,等.采煤機(jī)滾筒運(yùn)動(dòng)參數(shù)的分析研究[J].煤炭科學(xué)技術(shù),2008,36(8):62-64.

LIU Songyong, DU Changlong, CUI Xinxia, et al. Analysis and research on cutting drum movement parameters of coal drum shearer[J]. Coal Science and Technology, 2008,36(8):62-64.

[3] 李曉豁,高健.采煤機(jī)運(yùn)動(dòng)參數(shù)對(duì)塊煤產(chǎn)量的影響[J].煤炭學(xué)報(bào),2009,34(9):1268-1270.

LI Xiaohuo, GAO Jian. Effect of kinematics parameters of shearer on lump coal production[J]. Journal of China Coal Society, 2009,34(9):1268-1270.

[4] Dolipski M, Jaszczuk M .Computer-aided determination of dynamic loads in a longwall shearer's cutting system[J]. Computer Applications in the Minerals Industries, 2001, (4):413-416.

[5] 趙麗娟,馬強(qiáng),劉旭南,等.新型薄煤層采煤機(jī)裝煤量與可靠性研究[J].機(jī)械設(shè)計(jì),2014(7):97-101.

ZHAO Lijuan, MA Qiang, LIU Xunan, et al. Study on the loading capacity and reliability of a new type of thin seam coal mining machine [J]. Mechanical Design, 2014(7): 97-101.

[6] 高魁東.薄煤層滾筒采煤機(jī)裝煤性能研究[D].徐州:中國(guó)礦業(yè)大學(xué),2014:1-15.

[7] 秦大同,王鎮(zhèn),胡明輝,等. 基于多目標(biāo)優(yōu)化的采煤機(jī)滾筒最優(yōu)運(yùn)動(dòng)參數(shù)的動(dòng)態(tài)匹配[J].煤炭學(xué)報(bào),2015(S2):532-539.

QIN Datong,WANG Zhen,HU Minghui, et al. Dynamic matching of optimal drum movement parameters of shearer based on multi-objective optimization[J]. Journal of China Coal Society,2015(S2):532-539.

[8] 雷強(qiáng).基于離散元的物料破碎機(jī)理研究[D].贛州:江西理工大學(xué),2012.

[9] Ferellec J.F, Mc Dowell G.R.A method to model realistic particle shape and inertia in DEM[J].Granular Matter,2010(12):459-467.

[10] 王國(guó)強(qiáng),郝萬(wàn)軍,王繼新. 離散單元法及其在EDEM上的實(shí)踐[M].西安:西北工業(yè)大學(xué)出版社,2010.

[11] 趙麗娟,田震,劉旭南. 薄煤層采煤機(jī)滾筒載荷特性仿真分析[J].系統(tǒng)仿真學(xué)報(bào),2015(12):3102-3108.

ZHAO Lijuan, TIAN Zhen, LIU Xunan. Simulation and analysis of the load characteristics of the thin coal seam shearer drum[J].Journal of System Simulation,2015(12):3102-3108.

[12] 馬聯(lián)偉. 薄煤層采煤機(jī)可靠性與疲勞壽命研究[D].葫蘆島:遼寧工程技術(shù)大學(xué),2013:30-45.

[13] 趙麗娟,劉旭南,劉鵬. 截齒排列方式對(duì)采煤機(jī)工作性能的影響[J]. 機(jī)械科學(xué)與技術(shù),2014(12):1838-1844.

ZHAO Lijuan, LIU Xunan, LIU Peng. Picks arrangement mode’s influence on shearer’s working performance[J]. Mechanical Science and Technology for Aerospace Engineering,2014(12):1838-1844.

[14] 趙麗娟,李明昊,范佳藝. 基于VP-MFC的采煤機(jī)搖臂系統(tǒng)齒輪性能退化可靠性分析[J]. 中國(guó)安全生產(chǎn)科學(xué)技術(shù),2016,12(9):92-97.

ZHAO Lijuan, LI Minghao, FAN Jiayi. Reliability analysis of shearer rocker gear system performance based on VP-MFC[J]. Chinese Degradation of Safety Science and Technology, 2016,12(9): 92-97.

[15] 張強(qiáng),王海艦,李立瑩,等. 采煤機(jī)多參數(shù)耦合匹配優(yōu)化智能降塵控制研究[J]. 中國(guó)安全科學(xué)學(xué)報(bào),2016(11):139-144.

ZHANG Qiang, WANG Haijian, LI Liying,et al. Research on intelligent multi parameter coupling matching optimization intelligent dust control [J]. Proceedings of the Chinese Safety Science, 2016(11):139-144.

猜你喜歡
采煤機(jī)
煤礦用隔爆型采煤機(jī)電機(jī)FMECA分析
采煤機(jī)故障分析及處理方法
煤礦井下采煤機(jī)截割驅(qū)動(dòng)控制系統(tǒng)的應(yīng)用研究
MG400/920-WD型采煤機(jī)防滑制動(dòng)力的分析
煤礦采煤機(jī)的常見(jiàn)故障與維修措施
綜采工作面采煤機(jī)安裝工藝優(yōu)化探析
采煤機(jī)截割部關(guān)鍵零件有限元優(yōu)化設(shè)計(jì)
固態(tài)繼電器在采煤機(jī)四象限變頻改造中的應(yīng)用
MGTY300/700-1.1D采煤機(jī)制動(dòng)系統(tǒng)分析改進(jìn)
河南科技(2014年18期)2014-02-27 14:14:58
AM500采煤機(jī)小支臂改造
河南科技(2014年4期)2014-02-27 14:07:18
主站蜘蛛池模板: 欧美日韩中文字幕在线| 亚洲精品老司机| 欧洲熟妇精品视频| 亚洲最黄视频| 99精品热视频这里只有精品7| 亚洲国产成人无码AV在线影院L| 国产二级毛片| 在线观看视频99| 婷婷在线网站| 免费一级无码在线网站| 精品久久久久成人码免费动漫| 精品五夜婷香蕉国产线看观看| 国产欧美高清| 毛片免费在线视频| 亚洲第一成年网| 五月天福利视频| 天天综合色网| 国产精品午夜福利麻豆| 波多野结衣一区二区三视频| 国产91小视频在线观看| 91小视频在线观看免费版高清| 浮力影院国产第一页| 久久不卡国产精品无码| 黄色网站不卡无码| 农村乱人伦一区二区| 免费无码又爽又刺激高| 在线免费不卡视频| 亚洲人成日本在线观看| 亚洲一区二区三区在线视频| 久青草网站| 为你提供最新久久精品久久综合| 这里只有精品在线播放| 欧美日韩在线成人| 国产精品自在自线免费观看| 国产成人综合欧美精品久久| 草逼视频国产| 成人精品区| 午夜激情婷婷| 一本久道久综合久久鬼色| 国产精品久久精品| 欧美日韩理论| 黄色成年视频| 中文字幕 91| 亚洲女人在线| 亚洲成a人片在线观看88| 中文字幕精品一区二区三区视频 | 97精品久久久大香线焦| 久久鸭综合久久国产| 成年免费在线观看| 国产成人1024精品| 久夜色精品国产噜噜| 日本久久网站| 性欧美在线| 国产视频a| 国产91蝌蚪窝| 国产乱肥老妇精品视频| 99热这里只有免费国产精品 | 老司机久久精品视频| 精品国产免费观看| 91网红精品在线观看| 国产激情无码一区二区免费| 亚洲精品自拍区在线观看| 国产精品视频导航| 在线观看国产精品第一区免费| 国产激情影院| 国产成人欧美| 一级成人a毛片免费播放| 亚洲激情区| 在线观看av永久| 亚洲综合在线最大成人| 99视频在线看| 日韩精品亚洲一区中文字幕| 天天爽免费视频| 亚洲天堂成人| 凹凸国产分类在线观看| 国产Av无码精品色午夜| 久久特级毛片| 91成人在线观看视频| 中文字幕调教一区二区视频| 日韩视频免费| 久久99国产乱子伦精品免| 2021天堂在线亚洲精品专区|