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

    采煤機整機力學模型的預條件擬極小剩余算法

    2016-12-09 07:27:33劉春生李孝宇
    黑龍江科技大學學報 2016年5期
    關鍵詞:滑靴擺線采煤機

    劉春生, 李孝宇

    (1.黑龍江科技大學, 哈爾濱 150022; 2.黑龍江科技大學 機械工程學院, 哈爾濱 150022)

    ?

    采煤機整機力學模型的預條件擬極小剩余算法

    劉春生1, 李孝宇2

    (1.黑龍江科技大學, 哈爾濱 150022; 2.黑龍江科技大學 機械工程學院, 哈爾濱 150022)

    為了得到不同工況下采煤機各個滑靴的受力狀態(tài)及其運動特性,在考慮仰俯角和煤層傾角的條件下,建立了采煤機整機力學模型。采用Krylov子空間方法的擬極小剩余法,解算出不同工況下各滑靴受力的最佳逼近值,分析隨著煤層傾角和仰俯角的變化采煤機各滑靴受力的變化規(guī)律。結果表明:行走輪前一輪齒退出嚙合、后一輪齒進入嚙合時與銷齒嚙合區(qū)由圓弧段過渡到斜線段時,加速度值會出現(xiàn)兩次突變;后導向滑靴受正壓力較大,煤層傾角對各滑靴的受力狀態(tài)影響相對較小。該研究為采煤機整機力學特性的分析提供了參考。

    采煤機; 整機力學模型; 擬極小剩余法; 慣性力

    0 引 言

    采煤機是煤炭生產(chǎn)中主要的機械設備之一。隨著大功率、大采高采煤機的發(fā)展,整機裝機功率的增大,導向滑靴承受的負載力也越來越大,導致采煤機牽引部的一些關鍵部件因過載而損壞、失效如行走輪非正常斷齒、導向滑靴非正常損壞,以及行走輪和銷排銷齒的異常磨損等問題頻繁出現(xiàn)。諸多學者對采煤機行走輪動力學特性,陳偉等[1]利用Pro/E進行擺線行走輪的實體建模,通過ANSYS進行受力分析,證明了MG750/1815-WD采煤機擺線行走輪滿足強度。王振乾等[2]給出了行走輪五種齒形的強度和適應性。由于采煤機整機受力的復雜性,很難確定采煤機關鍵承載部件的參數(shù)。劉春生等[3]采用最小二乘解算方法對采煤機的力學模型進行解算,求解出采煤機實際工作過程中各個滑靴的力情況。采煤機在行走過程中,由于行走輪齒輪輪廓與銷排銷齒齒廓不能形成共軛嚙合,行走傳動比時刻發(fā)生變化,引起采煤機的速度和加速度變化,產(chǎn)生采煤機行走附加慣性力。筆者考慮附加慣性力的影響,建立采煤機整機力學模型,利用QMR算法求解非對稱大型稀疏線性系統(tǒng)的計算和存儲量均比較小的優(yōu)勢對采煤機整機力學模型進行求解,以期為分析采煤機的力學特性提供參考。

    1 整機力學模型

    在具有一定的仰俯角和煤層傾角的條件下,建立三維直角坐標系,如圖1所示,分析采煤機的力學狀態(tài)。

    圖1 采煤機受力分析

    由圖1可見,沿x、y、z三個坐標軸方向分別列出力平衡方程,ΣFx=0,ΣFy=0,ΣFz=0;以采煤機中心O點列三個力矩平衡方程,ΣFxoy=0,ΣFxoz=

    0,ΣFyoz=0,即

    (1)

    式中:a1=-F1-F2-Gcosαsinβ;

    a2=(2T+2FG-Gsinα-Fz1-Fz2)/μ;

    a3=Gcosαsinβ-2F-FJ1+FJ2;

    a4=-2Th0-2FG-FJ1(lcosθ1+l1+l0/2)-fz1(lsinθ1+h1-h0)-FJ2(lcosθ2+

    l1+l0/2)+Fz2(lsinθ2+h0-h1);

    a5=2Fb2+FJ2(b0+b1)-FJ1(b0+b1)-F1(lsinθ1+h1-h0)+F2(lsinθ2+h0-h1);

    a6=-2Tb2-2FGb2-(Fz1+Fz2)(b0+b1)+

    F1(lcosθ1+l1+l0/2+e)-F2(lcosθ2+

    l1+l0/2+e);

    T——單個牽引機構牽引力,kN;

    F——單個牽引機構的垂直分力,kN;

    μ——導向滑靴與導軌的摩擦系數(shù),μ=0.18;

    γ——行走輪與銷齒的嚙合角,(°);

    G——采煤機的整機重,kN;

    θ1、θ2——前后截割部的幾何擺角,(°);

    Ni——前、后支撐和導向滑靴所受支撐力,(i=1,z,…,6),kN;

    FJ1、FJ2——前、后螺旋滾筒截割阻力,kN;

    Fz1、Fz2——前、后螺旋滾筒推進阻力,kN;

    F1、F2——前、后螺旋滾筒軸向力,kN;

    α、β——分別為煤層傾角和仰俯角,(°);

    li、bi、hi——各受力點的主要位置尺寸,mm;

    e——軸向力作用位置距螺旋滾筒軸線的偏距,mm;

    FG——采煤機行走附加慣性力,kN;

    令C41=l0/2-μ(h0+h2),C42=l0/2-μ(h0+h2),C43=-μh0±l0/2,C44=-l0/2-μh0,C45=-μh0,C61=-μb1,C63=-μb2,C65=-μ(b3+b2)-l0/2,C66=-μ(b2-b3)-l0/2。

    將式(1)寫成矩陣的形式:

    (2)

    UX=A。

    (3)

    對式(3)首先判斷其解的存在情況,即獨立方程的個數(shù)。實際生產(chǎn)實踐中,N3的方向具有不確定性,當N3向上時取上行符號,當N3向下時取下行符號。分別對兩種情況的系數(shù)矩陣U作初等變換,當N3的方向向上是,U的秩是5;當N3的方向下時U的秩是6。N3方向的不確定性直接影響了系數(shù)矩陣U的秩。

    2 外部載荷的確定

    2.1 螺旋滾筒截割的確定

    根據(jù)文獻[4]前螺旋滾筒(截割阻力FJ1,推進阻力Fz1和軸向力F1)可按式(4)進行計算,后螺旋滾筒的按照前螺旋滾筒乘以工作條件系數(shù)k1,取k1=0.8。

    (4)

    式中:PH——單截割電機額定功率,kW;

    DC——螺旋滾筒直接,mm;

    ηj——截割部的傳動效率,ηj=0.85;

    J——螺旋滾筒有效截深,m;

    n——螺旋滾筒轉速,r/min;

    Lk——螺旋滾筒端盤部分截齒的截割寬度,m,取Lk=0.2 m;

    kq——截割磨損程度系數(shù),取kq=0.6~0.8;

    k2——螺旋滾筒端盤部分接近半封閉截割條件的系數(shù),取k2=2。

    2.2 附加慣性力的確定

    2.2.1 行走輪與銷齒輪廓線方程

    采煤機在行走過程中,由于行走輪與銷齒在嚙合時不共軛,會產(chǎn)生慣性力。文中研究的采煤機行走輪齒廓為擺線輪。擺線齒廓的行走輪,一般由齒頂處的外擺線和齒根處的內擺線兩部分組成,見2所示。

    a 外擺線

    b 內擺線

    (1)外擺線齒廓線方程

    由圖2可見,一個動圓外切與一個定圓作無滑動的純滾動,動圓圓周上一個定點的軌跡稱為普通外擺線,定圓稱為基圓,如圖2所示,外擺線的矢量-矩陣方程[5]為

    OP=OO1+O1P,

    即:

    r1(λ)=(R+ra)Kτ1j-raKτ2j。

    (5)

    直角坐標系中的方程為

    式中:R——采煤機行走輪節(jié)圓半徑,mm;

    ra——外擺線外滾圓半徑,mm;

    τ1——旋轉角,τ1=λ-ε,(°);

    λ——滾圓圓心和節(jié)圓圓心連線與起始線的夾角,(°);

    ε——行走輪節(jié)圓弦齒夾角的一半,(°)。

    (2) 內擺線齒廓線方程

    動圓內切與一個定圓作無滑動的純滾動時,動圓上一個定點形成的軌跡成為普通內擺線。如圖2所示。普通內擺線的矢量-矩陣方程為

    OP=OO1+O1P,

    r1(λ)=(R-rf)Kτ3j+rfKτ4j。

    (6)

    直角坐標系中的方程為

    式中:rf——內擺線內滾圓半徑,mm;

    τ3——旋轉角,τ3=-λ-ε,(°);

    (3) 銷齒齒廓線方程

    文中與行走輪配套的銷齒是截距為147mm的大截距銷齒。其齒廓線上面部分采用半徑為80mm的圓弧曲線,下面部分采用圓角曲線,中間用相切的直線過渡。上面圓弧段和過渡斜線段是主要的嚙合區(qū)。代入實際的相關尺寸后,嚙合區(qū)的矢量—矩陣方程圓弧段為

    r2(μ)=Kμri,

    (7)

    斜線段為

    (8)

    式中:μ——嚙合過程中銷齒滾動過的相對角度;

    τ5——旋轉角,τ1=40°;

    τ6——旋轉角,τ1=15°。

    2.2.2 行走輪與銷齒的嚙合傳動比

    采煤機實際行走過程中,銷排固定在地面上,行走輪一邊轉動一邊前進。在不改變此實際相對運動關系的條件下,假設行走輪作定軸轉動,銷排做平移運動。行走輪轉過ψ角度,原直角坐標oxy變?yōu)閛x1y1,銷排所在坐標系o2x2y2向前平移了S(ψ),如圖3所示。根據(jù)齒輪嚙合原理的接觸條件和相切條件,行走輪輪齒與銷齒的嚙合方程為:

    (9)

    式中:ψ——行走輪轉過的角度,(°);

    B——行走輪與銷齒嚙合中心距,mm。

    圖3 行走輪與銷齒嚙合矢量

    Fig. 3 Vector of meshing between traveling wheel and pin

    對行走輪與銷齒的嚙合方程兩邊對ψ求導,整理化簡后,得到行走輪與銷齒圓弧段嚙合的傳動比:

    (10)

    行走輪與銷齒的斜線段嚙合的傳動比為

    (11)

    將式(5)、(7)帶入式(10),整理可得到內、外擺線與圓形齒廓嚙合傳動比,內擺線為

    外擺線為

    將式(5)和(8)帶入式(11),簡化整理,可得外擺線與斜線段嚙合傳動比:

    斜線段

    根據(jù)傳動比即可求解采煤機行走速度以及加速度,進而求解采煤機行走附加慣性力。

    v=iω,

    FG=ma,

    式中:ω——行走輪角速度,rad/s;

    v——采煤機行走速度,mm/s;

    a——采煤機行走加速度,mm/s2;

    m——采煤機質量,kg。

    2.3 牽引力的確定

    根據(jù)文獻[3],牽引力可按式(17)進行計算:

    T=[μ(Gcosα-FJ1+FJ2+F1+F2)+FZ1+

    Fz2+Gsinα-FG]k3。

    (17)

    式中:k3——安全系數(shù),k3=1.20~1.25。

    3 解算方法

    3.1 擬極小剩余法

    采煤機受力力學模型的系數(shù)矩陣為奇異陣,因此方程無精確的解[6-10],但是采用Krylov子空間方法中擬極小剩余法(QMR)可以在求解某預條件下的近似解X。

    QMR算法是Krylov子空間迭代法。令x0為初始向量,r0為初始殘差,令Kk(A,r0)為k次krylov子空間,即,Kk(A,r0)=span{r0,Ar0,…,Ak-1r0},構造迭代序列{xk},使xk∈x0+Kk(A,r0),令β=‖r0‖2,v1=r0/β,取w1=v1,由Lanczos算法可得Lanczos向量v1,v2,…,vk,則

    xk=x0+[v1,v2,…,vk]zk=x0+Vkzk,zk∈k,

    對于線性方程組

    AX=B,

    (13)

    當矩陣A是奇異矩陣式,A-1是不存在的,故選用一個奇異矩陣M作為奇異線性系統(tǒng)的預條件陣可能更合理。設M是奇異的,求解線性系統(tǒng)AX=B可以通過求解下面的預條件方程:

    M+Ax=M+B,

    3.2 模型的計算

    采煤機力學模型的計算需要對式(3)增廣矩陣進行初等變換,首先考慮到N3的受力方向,當N3受力向下時,系數(shù)矩陣U的秩是6,通過正常求解方式即可求解。當N3方向向上是,系數(shù)矩陣U的秩為5,根據(jù)上述理論對奇異線性方程組UX=A進行求解。

    首先求解兩個初等變換陣P和Q,它們分別交換矩陣U的第ik和k行以及第jk和k列,k=1,2,…,r。即

    W11=W(i1,i2,…,ir;j1,j2,…,jr)是秩為5的W的子式,然后根據(jù)下面式子求解U+

    4 算例與結果分析

    以某厚煤層采煤機為例,計算不同仰俯角和煤層傾角的工況下,支撐于導向滑靴的受力情況?;拘阅軈?shù)為:機重約1 460kN,最大牽引力1 391kN,單截割電機功率1 000kW,螺旋滾筒轉速25.66r/min,螺旋滾筒直徑3.2mm,截深865mm,采煤機主要技術參數(shù)和位置尺寸,見表1。

    表1 整機技術參數(shù)

    根據(jù)整機技術參數(shù),通過式(4)計算采煤機前螺旋滾筒截割阻力FJ1為197.72 kN,推進阻力Fz1為158.17 kN和軸向力F1為91.43 kN,后螺旋滾筒的截割阻力FJ2為158.17 kN,推進阻力Fz1為126.53 kN和軸向力F1為73.14 kN。通過式(17)計算采煤機牽引力為572.14 kN。根據(jù)式(14)和(15)計算采煤機行走速度和加速度,如圖4所示。

    a 速度

    b 加速度

    擺線行走輪與嚙齒嚙合過程中,行走輪內擺線部分基本不參與嚙合,行走輪的主要嚙合區(qū)為外擺線部分。從圖4中可以看出,采煤機行走速度在74 mm/s處上下波動,速度的變化主要分為兩個階段:外擺線與銷齒圓弧段嚙合和外擺線與銷齒斜線段嚙合。這兩個階段速度的變化都呈遞減趨勢,外擺線與圓弧段嚙合到外擺線與斜線段嚙合的過渡區(qū)域速度發(fā)生突變。加速度變化曲線中出現(xiàn)兩個峰值,分別為114.07和-53.25 mm/s2,其中加速度最大值114.07 mm/s2出現(xiàn)在行走輪的前一個齒退出嚙合,后一個齒進入嚙合的時刻;加速度第二峰值-53.25 mm/s2出現(xiàn)在銷齒嚙合區(qū)由圓弧段到斜線段過渡的時刻。

    各滑靴所受變化如圖5所示,從圖5中可以看出,仰俯角β從-10°~10°變化,N3、N4、N6均有明顯增加其余滑靴受力均呈遞減趨勢, 當α=0°時,大約在β=6°處,N5的受力方向發(fā)生改變。當煤層傾角α由0°~14°變化,N3有略微下降,N4的值基本保持不變,其余滑靴受力均呈遞增趨勢。當β=6°時,大約α=10° 處,N5的受力方向發(fā)生改變。

    a 實線α=14°,虛線α=0°

    b 實線β=10°,虛線β=10°

    5 結 論

    (1)對采煤機整機力學模型進行分析時,考慮了采煤機行走過程中,行走輪與銷齒嚙合時不共軛產(chǎn)生的附加慣性力,采用Krylov子空間方法擬極小剩余法(QMR)對采煤機整機受力模型進行解算,通過迭代方法減小殘差,解出各滑靴受力的最佳逼近值。

    (2)行走輪與銷齒嚙合過程中,速度與加速度有兩個突變點,第一個突變點發(fā)生在前一個輪齒脫出嚙合下一個輪齒進入嚙合的時刻;第二個突變發(fā)生在銷齒的嚙合區(qū)由圓弧段過渡到斜線段的時刻。

    (3)隨著仰俯角和煤層傾角由小到大,各個支撐反力基本呈線性變化,變化趨勢受仰俯角的影響比較大。其中,N1和N2隨著仰俯角的增大而減??;N3、N4和N6隨著仰俯角的增大而增大;N5隨著仰俯角的增大先減小后增大。隨著煤層傾角的增大,N3呈遞減趨勢,其余滑靴的受力均呈遞增狀態(tài)。

    (4)N4受煤層傾角影響較小,其數(shù)值基本保持不變;N2數(shù)值較大,說明后支撐滑靴受力正壓力較大,磨損嚴重;隨著仰俯角的增大,N5先減小后增大,大約在β=6°時,受力方向發(fā)生改變,由原來滑靴與導軌的正壓力變?yōu)榛ハ鹿刺幱趯к壍睦?,這種受力方式對于滑靴非常不利,嚴重影響滑靴的使用壽命。

    (5)各滑靴的載荷受采煤機仰俯角影響較大。

    [1] 陳 偉, 杜長龍, 夏麗建, 等. 采煤機行走機構擺線輪設計研究[J]. 煤礦機械, 2010(1): 1-2.

    [2] 王振乾, 周常飛. 采煤機用行走輪的幾種齒形曲線初探[J]. 機械設計與研究, 2011(7): 39-42

    [3] 劉春生, 田 操. 采煤機整機力學模型的最小二乘解算方法[J]. 遼寧工程技術大學學報: 自然科學版, 2015, 34(4): 505-509.[4] 劉春生, 戴淑芝. 雙滾筒式采煤機整機力學模型與解算方法[J]. 黑龍江科技學院學報, 2012, 22(1): 33-38.

    [5] 王振乾. 滾筒式采煤機行走機構運動學分級及強度研究[D]. 上海: 煤炭科學研究總院上海分院, 2007.

    [6] 張世洪. 我國綜采采煤機技術的創(chuàng)新研究[J]. 煤炭學報, 2010, 11(5): 1898-1902.

    [7] 吳佳梁, 尹 力, 李 勇, 等. 智能化聯(lián)合采煤系統(tǒng)及其關鍵技術[J]. 遼寧工程技術大學學報: 自然科學版, 2014, 33(2): 226-231.

    [8] 劉春生, 李德根. 不同截割狀態(tài)下鎬型截齒側向力的實驗與理論模型[J]. 煤炭學報, 2016, 41(9): 2359-2366.

    [9] 李慶茹, 申 磊, 李玉標. 新型采煤機行走輪的研究與探討[J]. 煤礦機械, 2011, 32(12): 158-159.

    [10] 張勇濤, 王興文, 李 英, 等. 采煤機行走部導向滑靴分析[J]. 煤礦機械, 2012, 33(11): 210-211

    [11] 郎國軍, 苑雪濤, 庹文敏, 等. 采煤機行走輪斷齒現(xiàn)象分析[J]. 煤礦機電, 2010(3): 34-35.

    (編校 李德根)

    Shearer mechanics model preconditioned quasi-minimal residual algorithm

    LiuChunsheng1,LiXiaoyu2

    (1.Heilongjiang University Science & Technology, Harbin 150022, China;2.School of Mechanical Engineering, Heilongjiang University of Science & Technology, Harbin 150022, China)

    This paper is concerned with a better insight into the stress state of shearer and the motion characteristics of shearer under different working conditions; and the development of the mechanics model of shearer, along with the consideration of the conditions of the elevation angle and the inclination of coal seam. The research involves using Krylov subspace method of quasi minimal residual method, calculating the optimal approximation of slippers under force values due to different conditions, and delving into the mechanism by which the slippers are subjected to varying forces, depending on the variation of coal seam angle and pitch angle of the machine. Results show that there occur two sudden changes in the acceleration value when the engagement area of the pin tooth moves from the arc section to the oblique line segment, the front of the running gear is engaged with the gear wheel, and the rear gear tooth enters the engagement time. Rear guide sliding boots are subjected to a greater positive pressure and coal seam dip angle has a relatively small effect on the stress state of each shoe. This study might provide a reference for the analysis of mechanical properties of coal mining machines.

    coal winning machine; force of the whole machine; quasi-minimal residual method; inertial force

    2016-08-26

    國家自然科學基金項目(51274091);黑龍江省教育廳科學技術研究項目(12531006)

    劉春生(1961-),男,山東省牟平人,教授,研究方向:機械設計和液壓傳動與控制,E-mail:liu_chunsheng@163.com。

    10.3969/j.issn.2095-7262.2016.05.017

    TD421.61

    2095-7262(2016)05-0552-06

    A

    猜你喜歡
    滑靴擺線采煤機
    擺線的歷史和現(xiàn)狀
    煤礦用隔爆型采煤機電機FMECA分析
    防爆電機(2022年1期)2022-02-16 01:14:06
    采煤機支撐滑靴耳板斷裂問題分析與改進
    擺線銑削在硬零件槽銑加工中的應用技術
    防爆電機(2021年1期)2021-03-29 03:02:50
    采煤機滑靴失效分析
    擺線齒輪和傳動特點
    導向滑靴鑄鋼件的質量優(yōu)化
    基于平衡系數(shù)的滑靴優(yōu)化模型
    MGTY300/700-1.1D采煤機制動系統(tǒng)分析改進
    河南科技(2014年18期)2014-02-27 14:14:58
    AM500采煤機小支臂改造
    河南科技(2014年4期)2014-02-27 14:07:18
    久久99热6这里只有精品| 大陆偷拍与自拍| 国产午夜精品久久久久久一区二区三区| 色哟哟·www| 国产亚洲av嫩草精品影院| 丝瓜视频免费看黄片| 人妻一区二区av| 精品熟女少妇av免费看| 亚洲自偷自拍三级| 成人无遮挡网站| 国产欧美日韩一区二区三区在线 | 久久精品国产自在天天线| 色哟哟·www| 好男人在线观看高清免费视频| 免费黄网站久久成人精品| 亚洲在久久综合| 2021天堂中文幕一二区在线观| 九色成人免费人妻av| 春色校园在线视频观看| 一级爰片在线观看| 午夜福利网站1000一区二区三区| 五月玫瑰六月丁香| 中文乱码字字幕精品一区二区三区| av黄色大香蕉| 最近中文字幕2019免费版| 99热国产这里只有精品6| 黄色怎么调成土黄色| 99热这里只有是精品在线观看| av专区在线播放| 大话2 男鬼变身卡| 精品久久国产蜜桃| 秋霞在线观看毛片| 人人妻人人爽人人添夜夜欢视频 | 亚洲精品一二三| 久久热精品热| 搡老乐熟女国产| 国产成人aa在线观看| 久久精品熟女亚洲av麻豆精品| 午夜福利在线观看免费完整高清在| www.色视频.com| 欧美成人a在线观看| 狠狠精品人妻久久久久久综合| 日本熟妇午夜| 亚洲欧美成人综合另类久久久| 久久久久久久精品精品| 女的被弄到高潮叫床怎么办| 久久女婷五月综合色啪小说 | 色视频www国产| 中文资源天堂在线| 欧美日韩精品成人综合77777| 国产综合精华液| 麻豆乱淫一区二区| 日日摸夜夜添夜夜添av毛片| 哪个播放器可以免费观看大片| 天堂网av新在线| 成人高潮视频无遮挡免费网站| 久久99蜜桃精品久久| 久久影院123| 麻豆乱淫一区二区| 亚洲最大成人中文| 视频区图区小说| 69av精品久久久久久| 久热久热在线精品观看| 黄色配什么色好看| 一区二区av电影网| 亚洲高清免费不卡视频| 老女人水多毛片| 美女脱内裤让男人舔精品视频| 亚洲伊人久久精品综合| 国产淫片久久久久久久久| 联通29元200g的流量卡| 久久热精品热| 嫩草影院新地址| 欧美成人精品欧美一级黄| 日本熟妇午夜| 国产成人91sexporn| 王馨瑶露胸无遮挡在线观看| 国产在线一区二区三区精| 色网站视频免费| 七月丁香在线播放| 亚洲成人久久爱视频| 亚洲人成网站在线观看播放| 国产精品无大码| 观看美女的网站| 久久久久久久大尺度免费视频| 岛国毛片在线播放| 亚洲av免费在线观看| 成人无遮挡网站| 观看免费一级毛片| 中文字幕av成人在线电影| 丰满乱子伦码专区| 中文字幕亚洲精品专区| 插逼视频在线观看| 一级片'在线观看视频| 中文乱码字字幕精品一区二区三区| 国产欧美日韩精品一区二区| 国产精品无大码| 欧美激情在线99| 欧美高清性xxxxhd video| av福利片在线观看| xxx大片免费视频| 男人爽女人下面视频在线观看| 国产伦精品一区二区三区四那| 日韩av在线免费看完整版不卡| 激情 狠狠 欧美| 色视频www国产| www.av在线官网国产| 日本wwww免费看| 插阴视频在线观看视频| 看黄色毛片网站| 亚洲自拍偷在线| 69av精品久久久久久| 久久久久久久久久久丰满| 熟女电影av网| 我要看日韩黄色一级片| 99久久人妻综合| 精品熟女少妇av免费看| 18禁在线播放成人免费| 久久精品国产亚洲av天美| 你懂的网址亚洲精品在线观看| 亚洲av免费在线观看| 亚洲成人久久爱视频| av在线亚洲专区| 欧美日韩精品成人综合77777| 嫩草影院入口| 成人一区二区视频在线观看| 亚洲天堂av无毛| 好男人视频免费观看在线| av在线老鸭窝| 亚洲av欧美aⅴ国产| 国产亚洲最大av| 久久精品久久久久久久性| 免费大片黄手机在线观看| 亚洲怡红院男人天堂| 免费观看的影片在线观看| 国产淫语在线视频| 欧美日韩精品成人综合77777| 18禁在线无遮挡免费观看视频| 亚洲人成网站高清观看| 中文字幕av成人在线电影| 欧美日韩视频高清一区二区三区二| 久久久久久久大尺度免费视频| 最近中文字幕2019免费版| 日本与韩国留学比较| 亚洲精品国产成人久久av| 好男人在线观看高清免费视频| 欧美成人午夜免费资源| 男女啪啪激烈高潮av片| 国产成人免费无遮挡视频| 中文在线观看免费www的网站| 欧美精品一区二区大全| 久久97久久精品| 国产 一区 欧美 日韩| av网站免费在线观看视频| 青春草亚洲视频在线观看| 狠狠精品人妻久久久久久综合| av.在线天堂| 精品一区在线观看国产| 尤物成人国产欧美一区二区三区| 久久久久久九九精品二区国产| 在线观看国产h片| 国产91av在线免费观看| 日韩av不卡免费在线播放| 综合色丁香网| 欧美97在线视频| 亚洲精品久久久久久婷婷小说| 免费观看无遮挡的男女| 另类亚洲欧美激情| 我的女老师完整版在线观看| 成人鲁丝片一二三区免费| 网址你懂的国产日韩在线| 精品午夜福利在线看| 国产视频内射| 成年女人看的毛片在线观看| 一区二区三区四区激情视频| 日韩欧美一区视频在线观看 | 久久热精品热| 嫩草影院新地址| 观看免费一级毛片| 少妇人妻一区二区三区视频| 国产淫片久久久久久久久| 99re6热这里在线精品视频| 久久99热这里只有精品18| 亚洲国产色片| 精品酒店卫生间| 欧美精品国产亚洲| 99热这里只有精品一区| 亚洲一区二区三区欧美精品 | 91久久精品电影网| 日本av手机在线免费观看| 国产大屁股一区二区在线视频| 亚洲精品视频女| 亚洲国产精品成人综合色| 欧美97在线视频| 精品人妻视频免费看| 婷婷色综合www| 亚洲精华国产精华液的使用体验| 久久久久久伊人网av| 一级黄片播放器| 久久精品国产亚洲网站| 色综合色国产| av又黄又爽大尺度在线免费看| 国产精品一二三区在线看| 大陆偷拍与自拍| 嫩草影院精品99| 九九爱精品视频在线观看| 纵有疾风起免费观看全集完整版| 女人十人毛片免费观看3o分钟| 日韩电影二区| 免费av毛片视频| 欧美xxxx性猛交bbbb| 国产精品99久久久久久久久| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产一级毛片在线| av网站免费在线观看视频| 在线观看免费高清a一片| 精品视频人人做人人爽| 我的女老师完整版在线观看| 亚洲精品国产av蜜桃| 欧美精品一区二区大全| 亚洲精品国产av成人精品| 亚洲精品一二三| 免费电影在线观看免费观看| 女的被弄到高潮叫床怎么办| 免费高清在线观看视频在线观看| 插阴视频在线观看视频| 五月开心婷婷网| 久久久久久国产a免费观看| 午夜精品国产一区二区电影 | 一本久久精品| 亚洲色图av天堂| a级毛片免费高清观看在线播放| 99久久人妻综合| 国产精品国产三级专区第一集| 欧美成人一区二区免费高清观看| 国产探花极品一区二区| 亚洲精品国产色婷婷电影| 日日啪夜夜爽| 欧美日韩精品成人综合77777| 亚洲欧美中文字幕日韩二区| 波野结衣二区三区在线| 午夜免费鲁丝| a级一级毛片免费在线观看| 国产精品秋霞免费鲁丝片| 三级国产精品片| 欧美日韩一区二区视频在线观看视频在线 | 国产成人aa在线观看| .国产精品久久| 国产视频首页在线观看| 一级毛片aaaaaa免费看小| 国产伦精品一区二区三区视频9| 亚洲欧美日韩另类电影网站 | 国产黄色免费在线视频| 国产亚洲av嫩草精品影院| 人妻制服诱惑在线中文字幕| 美女xxoo啪啪120秒动态图| 国产日韩欧美在线精品| av免费在线看不卡| 欧美成人午夜免费资源| 99视频精品全部免费 在线| 国产伦精品一区二区三区视频9| 国产黄色免费在线视频| 国产一区有黄有色的免费视频| 成人免费观看视频高清| 国产欧美另类精品又又久久亚洲欧美| 久久久久九九精品影院| 嘟嘟电影网在线观看| 一级爰片在线观看| 22中文网久久字幕| 久久综合国产亚洲精品| 精品99又大又爽又粗少妇毛片| 在线看a的网站| 听说在线观看完整版免费高清| av黄色大香蕉| 国产精品福利在线免费观看| 在线观看免费高清a一片| 久久久a久久爽久久v久久| 在线看a的网站| 我要看日韩黄色一级片| 精品99又大又爽又粗少妇毛片| 全区人妻精品视频| 一边亲一边摸免费视频| a级毛色黄片| 大香蕉久久网| 97超视频在线观看视频| 亚洲欧美日韩另类电影网站 | 免费在线观看成人毛片| www.av在线官网国产| 亚洲av男天堂| 黄色视频在线播放观看不卡| av在线亚洲专区| 亚洲国产av新网站| 亚洲最大成人手机在线| 国产午夜精品一二区理论片| 男女下面进入的视频免费午夜| 亚洲四区av| av在线亚洲专区| 久久久成人免费电影| 好男人在线观看高清免费视频| 在线播放无遮挡| 最新中文字幕久久久久| 国产亚洲精品久久久com| 日韩成人伦理影院| 搡女人真爽免费视频火全软件| 久久久久国产精品人妻一区二区| 亚洲精品aⅴ在线观看| 亚洲国产最新在线播放| 免费av不卡在线播放| 天美传媒精品一区二区| 搞女人的毛片| 少妇裸体淫交视频免费看高清| 久久久精品免费免费高清| 五月玫瑰六月丁香| 国产高清三级在线| 国产精品人妻久久久影院| 久久精品国产亚洲av天美| 亚洲精品成人久久久久久| 国产亚洲一区二区精品| 2021天堂中文幕一二区在线观| 欧美xxxx性猛交bbbb| 成人毛片60女人毛片免费| 午夜老司机福利剧场| 国产欧美亚洲国产| 国精品久久久久久国模美| 尾随美女入室| 免费大片18禁| 高清视频免费观看一区二区| 久久综合国产亚洲精品| 精品少妇黑人巨大在线播放| 日本av手机在线免费观看| 日韩成人伦理影院| 欧美日韩视频高清一区二区三区二| 97超碰精品成人国产| 久久久久久国产a免费观看| 黄色欧美视频在线观看| 亚洲内射少妇av| 国产成人a区在线观看| 色吧在线观看| 亚洲天堂国产精品一区在线| 七月丁香在线播放| 又黄又爽又刺激的免费视频.| 免费av毛片视频| 别揉我奶头 嗯啊视频| 又大又黄又爽视频免费| 亚洲图色成人| 美女国产视频在线观看| 涩涩av久久男人的天堂| 成人国产av品久久久| 最新中文字幕久久久久| 高清毛片免费看| 又大又黄又爽视频免费| 午夜精品国产一区二区电影 | 亚洲精品成人av观看孕妇| 欧美精品一区二区大全| 少妇 在线观看| 国产乱人视频| 亚洲欧美日韩另类电影网站 | 免费播放大片免费观看视频在线观看| 亚洲不卡免费看| 久久精品久久久久久噜噜老黄| 欧美成人精品欧美一级黄| 亚洲精品自拍成人| 青春草国产在线视频| 1000部很黄的大片| 狂野欧美白嫩少妇大欣赏| 欧美精品一区二区大全| 一级黄片播放器| 神马国产精品三级电影在线观看| 国产视频内射| av国产免费在线观看| 免费看日本二区| av国产免费在线观看| 国产精品蜜桃在线观看| 亚洲欧洲国产日韩| 熟妇人妻不卡中文字幕| 搡女人真爽免费视频火全软件| 亚洲美女搞黄在线观看| 久久久久久伊人网av| 中文在线观看免费www的网站| av免费在线看不卡| 国产色婷婷99| 久久久久久久午夜电影| 亚洲无线观看免费| 久久精品久久精品一区二区三区| 不卡视频在线观看欧美| 天天躁日日操中文字幕| 亚洲国产高清在线一区二区三| 国产亚洲精品久久久com| 日本黄大片高清| 五月天丁香电影| 亚洲av日韩在线播放| 亚洲精品日本国产第一区| 两个人的视频大全免费| 欧美xxxx性猛交bbbb| 亚洲av中文av极速乱| 国产视频内射| 26uuu在线亚洲综合色| 99热这里只有精品一区| 99久久中文字幕三级久久日本| av在线老鸭窝| 男女边吃奶边做爰视频| 亚洲性久久影院| 国产黄片视频在线免费观看| 网址你懂的国产日韩在线| 欧美日韩国产mv在线观看视频 | 亚洲精品乱久久久久久| 欧美潮喷喷水| 舔av片在线| 搞女人的毛片| 亚洲aⅴ乱码一区二区在线播放| 网址你懂的国产日韩在线| 欧美日韩国产mv在线观看视频 | 中文精品一卡2卡3卡4更新| 国产淫语在线视频| av在线天堂中文字幕| 91狼人影院| 久久精品久久精品一区二区三区| 久久久久久久久久人人人人人人| 久久久色成人| 欧美亚洲 丝袜 人妻 在线| 国产精品伦人一区二区| 亚洲伊人久久精品综合| 青春草视频在线免费观看| 久久久久久国产a免费观看| 国产精品不卡视频一区二区| 欧美zozozo另类| 免费av不卡在线播放| 毛片女人毛片| 亚洲精品乱码久久久久久按摩| 人妻系列 视频| 少妇人妻一区二区三区视频| 高清在线视频一区二区三区| 2018国产大陆天天弄谢| 尾随美女入室| 嫩草影院入口| 两个人的视频大全免费| 性色av一级| eeuss影院久久| 国产在线一区二区三区精| 亚洲自偷自拍三级| 91精品一卡2卡3卡4卡| 一区二区三区四区激情视频| 熟女电影av网| 久久久久国产精品人妻一区二区| 日本午夜av视频| 久久精品国产a三级三级三级| 亚洲自偷自拍三级| 熟女av电影| 一二三四中文在线观看免费高清| 免费观看在线日韩| 国产中年淑女户外野战色| 亚洲精品日韩av片在线观看| 亚洲精品一区蜜桃| 国产男人的电影天堂91| 亚洲无线观看免费| 国产免费一区二区三区四区乱码| 日韩欧美 国产精品| 韩国av在线不卡| 亚洲天堂av无毛| 尤物成人国产欧美一区二区三区| 青春草国产在线视频| 国产精品福利在线免费观看| 直男gayav资源| 亚洲国产欧美在线一区| 免费av观看视频| 国精品久久久久久国模美| 国产精品久久久久久久电影| 少妇熟女欧美另类| 久久精品久久久久久噜噜老黄| 国产色爽女视频免费观看| 精品久久久久久久久av| 国产成人freesex在线| 王馨瑶露胸无遮挡在线观看| 国内少妇人妻偷人精品xxx网站| 欧美3d第一页| 日韩不卡一区二区三区视频在线| 中文字幕亚洲精品专区| 久久精品国产自在天天线| 老女人水多毛片| 麻豆国产97在线/欧美| 夜夜爽夜夜爽视频| 国产av不卡久久| 人人妻人人看人人澡| 日日摸夜夜添夜夜爱| 你懂的网址亚洲精品在线观看| 最近的中文字幕免费完整| 国产 一区精品| 国产综合懂色| 亚洲最大成人av| 三级经典国产精品| 午夜爱爱视频在线播放| 少妇被粗大猛烈的视频| 午夜精品一区二区三区免费看| 搡老乐熟女国产| 黄色日韩在线| 狂野欧美激情性bbbbbb| 成人国产麻豆网| 国产淫语在线视频| 久久97久久精品| 国内精品美女久久久久久| 国产精品一二三区在线看| 久久99热这里只有精品18| 日韩一区二区视频免费看| 啦啦啦中文免费视频观看日本| 免费观看在线日韩| 欧美丝袜亚洲另类| 高清日韩中文字幕在线| 欧美日本视频| 丰满少妇做爰视频| 蜜臀久久99精品久久宅男| 国产在线一区二区三区精| 97在线视频观看| 成人免费观看视频高清| 国产成人aa在线观看| 在线观看一区二区三区| 美女被艹到高潮喷水动态| 黄色怎么调成土黄色| 免费观看在线日韩| 亚洲成人一二三区av| 亚洲av.av天堂| 97超碰精品成人国产| 大片免费播放器 马上看| 天天躁夜夜躁狠狠久久av| 精品99又大又爽又粗少妇毛片| 少妇的逼好多水| 久久鲁丝午夜福利片| 在现免费观看毛片| 哪个播放器可以免费观看大片| 国产成人精品一,二区| 免费大片18禁| 永久网站在线| 特大巨黑吊av在线直播| 日本猛色少妇xxxxx猛交久久| 免费看光身美女| 大香蕉97超碰在线| 亚洲天堂av无毛| 成人综合一区亚洲| 亚洲色图综合在线观看| 日韩欧美 国产精品| 中文字幕免费在线视频6| 精品久久久久久久久av| 亚洲经典国产精华液单| 午夜激情福利司机影院| 精华霜和精华液先用哪个| 青春草视频在线免费观看| 久久97久久精品| 久久久久久九九精品二区国产| 免费黄网站久久成人精品| 欧美少妇被猛烈插入视频| 夫妻性生交免费视频一级片| 麻豆精品久久久久久蜜桃| 日韩大片免费观看网站| 国产伦在线观看视频一区| 大香蕉97超碰在线| eeuss影院久久| av专区在线播放| 18禁在线无遮挡免费观看视频| 2022亚洲国产成人精品| 久久久a久久爽久久v久久| 少妇裸体淫交视频免费看高清| 大话2 男鬼变身卡| 成年av动漫网址| 夫妻午夜视频| 欧美日韩综合久久久久久| 久久精品久久久久久久性| 777米奇影视久久| 99久久九九国产精品国产免费| 国产精品99久久99久久久不卡 | 日韩大片免费观看网站| 搞女人的毛片| 简卡轻食公司| 精品一区二区免费观看| 国产成人一区二区在线| 久久久久网色| 亚洲怡红院男人天堂| 熟女av电影| 一本一本综合久久| 99热国产这里只有精品6| 午夜免费男女啪啪视频观看| 亚洲欧美清纯卡通| 国产 一区精品| 秋霞伦理黄片| 精品久久久精品久久久| 国产精品无大码| 日韩电影二区| 王馨瑶露胸无遮挡在线观看| 亚洲成人精品中文字幕电影| 免费少妇av软件| 久久韩国三级中文字幕| 免费高清在线观看视频在线观看| 亚洲在线观看片| 丝袜美腿在线中文| 日产精品乱码卡一卡2卡三| 亚洲av成人精品一区久久| 91在线精品国自产拍蜜月| 91午夜精品亚洲一区二区三区| 涩涩av久久男人的天堂| 69人妻影院| 亚洲激情五月婷婷啪啪| 国产色爽女视频免费观看| 久久6这里有精品| 亚洲激情五月婷婷啪啪| 国产高清三级在线| 一二三四中文在线观看免费高清| 神马国产精品三级电影在线观看| 国产一区有黄有色的免费视频| 丝袜美腿在线中文| 免费看a级黄色片| 天堂网av新在线| 18禁在线无遮挡免费观看视频| 亚洲欧美日韩无卡精品| 午夜精品一区二区三区免费看| 亚洲精品影视一区二区三区av| 日韩欧美精品v在线| 国产欧美日韩精品一区二区|