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

    純電動拖拉機(jī)機(jī)電耦合系統(tǒng)非線性振動分析

    2018-10-18 07:15:16嚴(yán)國軍顧建華姚月琴
    關(guān)鍵詞:控制參數(shù)傳動系統(tǒng)機(jī)電

    嚴(yán)國軍, 楊 彥, 顧建華, 姚月琴,祁 淼,4

    (1. 鹽城工業(yè)職業(yè)技術(shù)學(xué)院 機(jī)電工程學(xué)院,江蘇 鹽城 224000; 2. 東南大學(xué) 機(jī)械工程學(xué)院,江蘇 南京 210000; 3. 南京理工大學(xué) 控制工程學(xué)院,江蘇 南京 210000;4. 江蘇大學(xué) 汽車工程學(xué)院,江蘇 鎮(zhèn)江 212000)

    0 引 言

    在我國,電動農(nóng)業(yè)機(jī)械還處于研究階段。無論從市場需求,還是能源發(fā)展方向等方面考慮,電動農(nóng)業(yè)機(jī)械都是我國新農(nóng)村建設(shè)和現(xiàn)代農(nóng)業(yè)發(fā)展的需要。我國“十二五”規(guī)劃曾明確提出:“提高設(shè)施裝備水平,提升產(chǎn)業(yè)發(fā)展質(zhì)量,促進(jìn)設(shè)施農(nóng)業(yè)可持續(xù)發(fā)展”。隨著新能源車輛的發(fā)展,作為農(nóng)業(yè)大國,純電動拖拉機(jī)正在慢慢替代傳統(tǒng)燃油拖拉機(jī),并且應(yīng)用在農(nóng)業(yè)的各個方面。

    純電動車輛動力傳動系統(tǒng)中普遍存在機(jī)電耦合現(xiàn)象,研究機(jī)電系統(tǒng)耦合振動,可減小系統(tǒng)非線性特性的影響。純電動車輛動力傳動系統(tǒng)的扭轉(zhuǎn)振動問題在電磁激勵和外部激勵的作用下變得復(fù)雜與突出。在機(jī)電耦合方面,林利紅等人對純電動車輛動力傳動系統(tǒng)進(jìn)行機(jī)電耦合動力學(xué)分析和實驗研究,建立了系統(tǒng)的機(jī)電耦合振動仿真模型,仿真分析了由于電流調(diào)節(jié)器參數(shù)、阻尼、諧波擾動、間隙以及負(fù)載擾動等因素引起的機(jī)電耦合振動動態(tài)過程,并將物理實驗研究與理論分析結(jié)果進(jìn)行對比和驗證[1]。在車用電機(jī)非線性扭轉(zhuǎn)振動方面,陳星等人采用多尺度法研究了純電動車輛動力傳動系統(tǒng)的扭轉(zhuǎn)振動系統(tǒng)特性,仿真結(jié)果表明合理的設(shè)計和控制電機(jī)的電磁參數(shù)和機(jī)械參數(shù),可以解決電機(jī)扭振系統(tǒng)出現(xiàn)的跳躍分岔等非線性特性現(xiàn)象[2-3]。

    針對低速重載起步工況下純電動拖拉機(jī)機(jī)電耦合傳動系統(tǒng)受到非周期沖擊干擾力作用,分析扭轉(zhuǎn)角與電磁轉(zhuǎn)矩的耦合影響。為此筆者通過建立純電動拖拉機(jī)機(jī)電耦合傳動系統(tǒng)的非線性扭轉(zhuǎn)振動模型,求解無擾動Hamilton系統(tǒng)的平衡點;利用Melnikov法求解分岔閾值曲線,分析電機(jī)控制參數(shù)對純電動拖拉機(jī)機(jī)電耦合傳動系統(tǒng)非線性動力學(xué)行為的影響;基于Melnikov理論分析研究傳動系統(tǒng)的分岔和混沌特性,從而求解系統(tǒng)穩(wěn)定性運行邊界條件。

    1 機(jī)電耦合扭轉(zhuǎn)振動模型

    1.1 電機(jī)模型的建立

    不同情況下,永磁同步電機(jī)的穩(wěn)定運行時的相量如圖1。

    圖1 永磁同步電機(jī)相量Fig. 1 Phasor diagram of permanent magnet synchronous motor

    由圖1得到如下關(guān)系式:

    (1)

    通過求解式(1)可得到電機(jī)定子電流的直、交軸分量:

    (2)

    定子相電流:

    (3)

    而電機(jī)的輸入功率:

    P1=mUI1cosφ

    (4)

    當(dāng)忽略定子電阻,由式(4)可得到電機(jī)的電磁功率:

    (5)

    則電機(jī)的電磁轉(zhuǎn)矩:

    (6)

    式中:m為電機(jī)相數(shù);Ω為機(jī)械角速度,且ω=pΩ,其中,ω為電機(jī)的電角速度,p為電機(jī)的極對數(shù);E0為空載反電動勢,可由式(7)計算得到。

    (7)

    式中:f為電源頻率;Kdp為基波繞組系數(shù);N為每相串聯(lián)匝數(shù);Φ10為永磁體基波磁通;Am為永磁體面積;Br為永磁體剩磁;σ0為空載漏磁系數(shù);bm0為永磁體空載工作點假設(shè)值;KΦ為氣隙磁通波形系數(shù);αi為計極弧系數(shù);αp為極弧系數(shù),δ為氣隙長度,D為定子內(nèi)徑。

    Te=b1sinθ+b2sin 2θ

    (8)

    由圖1可知,θ=φ+ψ,因此可得到sinθ和sin 2θ三角函數(shù)展開式。

    利用sinx≈x-1/6x3,cosx≈1-1/2x2,泰勒展開的近似式將含φ的三角函數(shù)展開并化簡,即可得到作用于轉(zhuǎn)子的電磁轉(zhuǎn)矩:

    Te=k0+k1φ-k2φ2-k3φ3

    (9)

    式中:k0=b1sinψ+b2cos 2ψ;k1=b1cosψ+2b2sinψ;k2=1/2b1sinψ+2b2sin 2ψ;k3=1/6b1cosψ+4/3b2cos 2ψ。

    1.2 機(jī)電耦合扭轉(zhuǎn)振動模型

    通過分析車載永磁同步電機(jī)機(jī)電耦合扭轉(zhuǎn)振動特性,將轉(zhuǎn)子系統(tǒng)和傳動系統(tǒng)進(jìn)行簡化,可以等效為二質(zhì)體機(jī)電耦合轉(zhuǎn)子模型[4-5],如圖2。圖中永磁同步電機(jī)轉(zhuǎn)子與機(jī)械轉(zhuǎn)子之間彈性聯(lián)軸器聯(lián)接,電機(jī)機(jī)座可視為剛性。圖中Te是電磁轉(zhuǎn)矩,TL是負(fù)載轉(zhuǎn)矩。

    圖2 二質(zhì)體機(jī)電耦合轉(zhuǎn)子系統(tǒng)Fig. 2 Two plastid electromechanical coupled rotor system

    根據(jù)牛頓定律,二質(zhì)體轉(zhuǎn)子系統(tǒng)扭轉(zhuǎn)振動動力學(xué)方程為:

    (10)

    (11)

    當(dāng)轉(zhuǎn)子系統(tǒng)處于穩(wěn)態(tài)運行時,式(10)通過方程變換并消除得到公式中的靜態(tài)項:

    (12)

    則機(jī)電耦合扭轉(zhuǎn)振動動力學(xué)方程可轉(zhuǎn)化為無量綱形式:

    (13)

    式中:v=J1/J2;r=1/(1+v);μ=(1/J1+1/J2)C;β=k2r2/J1;γ=k3r3/J1;η2=ω0-k1r/J1,ω0為不考慮機(jī)電耦合時轉(zhuǎn)子系統(tǒng)的機(jī)械固有頻率,可由式(14)進(jìn)行計算:

    (14)

    從式(14)可看出,動力學(xué)方程具有平方項和立方項,方程為典型的非線性動力學(xué)方程。因此,考慮到機(jī)電耦合效應(yīng),傳動系統(tǒng)轉(zhuǎn)子扭轉(zhuǎn)振動方程將有明顯的非線性特征。

    2 系統(tǒng)平衡點分析

    選取的拖拉機(jī)電機(jī)額定功率為7.5 kW,額定電壓為96 V,其他相關(guān)電機(jī)參數(shù)如表1。

    表1 電機(jī)參數(shù)表Table 1 Motor parameter list

    式中:ε為小參數(shù),且0<ε?1。

    無擾動系統(tǒng)的勢函數(shù)在研究非線性動力學(xué)行為中具有十分重要的作用,因此需要討論無擾動系統(tǒng)的函數(shù)。當(dāng)ε=0時,排除阻尼項和擾動項,函數(shù)方程變?yōu)闊o擾動系統(tǒng):

    (16)

    在排除阻尼項和擾動項后,轉(zhuǎn)子動力學(xué)系統(tǒng)為Hamilton系統(tǒng),其Hamilton方程為:

    (17)

    式(17)對應(yīng)的勢函數(shù)為:

    (18)

    Hamilton系統(tǒng)平衡點穩(wěn)定性分為如下幾種情況:

    1)情況A

    當(dāng)β2≤4γ(1-κ),Hamilton系統(tǒng)只有一個平衡點P0(0,0):

    ①情況A1:當(dāng)κ<1時,P0是中心點。

    ②情況A2:當(dāng)κ>1時,P0是鞍點。

    2)情況B

    當(dāng)β2>4γ(1-κ)和4γ(1-κ)>0時,有3個平衡點,分別為P0(0,0),P1(x1,0),P2(x2,0):

    ①情況B1:當(dāng)γ>0時,P1(x1,0)是鞍點,P0(0,0)和P2(x2,0)是中心點。當(dāng)β=0.8,κ=0.35,γ=0.2時,可以得到勢函數(shù)和相軌跡如圖3(a)和圖3(b),從圖3(a)和圖3(b)中可看出存在橫截同宿軌道。

    ②情況B2:當(dāng)γ<0時,P1(x1,0)是中心點,P0(0,0)和P2(x2,0)是鞍點。當(dāng)β=0.8,κ=1.65,γ=-0.2時,可得到勢函數(shù)和相軌跡如圖3(c)和3圖(d),從圖3(c)和3圖(d)中可以看出存在橫截異宿軌道。

    圖3 B1和B2情況下Hamilton系統(tǒng)勢函數(shù)和相軌跡Fig. 3 Potential functions and phase trajectories of Hamilton systems in the case of B1 and B2

    3)情況C

    當(dāng)β2>4γ(1-κ)和4γ(1-κ)<0時,有3個平衡0點,分別為P0(0,0),P1(x1,0),P2(x2,0):

    ①情況C1:當(dāng)γ>0時,P0(0,0)是鞍點,P1(x1,0)和P2(x2,0)是中心點。當(dāng)β=0,κ=1.5,γ=0.2時,可得到勢函數(shù)和相軌跡如圖4(a)和圖4(b),從圖4(a)和圖4(b)中可以看出存在橫截同宿軌道。

    ②情況C2:當(dāng)γ<0時,P0(0,0)是中心點,P1(x1,0)和P2(x2,0)是鞍點。當(dāng)β=0.8,κ=1.65,γ=-0.2時,可得到勢函數(shù)和相軌跡如圖4(c)和圖4(d),從圖中可以看出存在橫截異宿軌道。

    圖4 C1和C2情況下Hamilton系統(tǒng)勢函數(shù)和相軌跡Fig. 4 Potential functions and phase trajectories of Hamilton systems in the case of C1 and C2

    綜合以上分析可知:情況A中只存在一個平衡點,將不做討論;情況B和C均具有3個平衡點,而在情況B1和C1中,對較小的擾動幅值,Hamilton系統(tǒng)具有封閉軌道,此時將產(chǎn)生橫截同宿軌道,可能導(dǎo)致同宿混沌;在情況B2和C2中,增加激勵幅值,將會產(chǎn)生橫截異宿軌道,可能導(dǎo)致異宿混沌。

    從圖3和圖4中可以看出,考慮機(jī)電耦合作用的Hamilton系統(tǒng)是雙勢對稱或不對稱波動系統(tǒng)原型,同時由于機(jī)電耦合動力學(xué)方程系數(shù)的復(fù)雜多變,所以機(jī)電系統(tǒng)扭轉(zhuǎn)振動將會表現(xiàn)出復(fù)雜的非線性動力學(xué)行為。車用永磁同步電機(jī)傳動系統(tǒng)為獲得更高的傳動范圍和高功率密度,通常設(shè)計和控制于弱磁效應(yīng),因此,筆者只討論ψ>0的情況。當(dāng)內(nèi)功率因數(shù)角ψ在[0,π]內(nèi)變化時,Hamilton系統(tǒng)的情況如圖5。分析圖可知,系統(tǒng)只存在B1和C1兩種情況,筆者將以C1為例進(jìn)行進(jìn)一步討論(B1可采用相似的過程進(jìn)行計算討論)。

    圖5 機(jī)電傳動系統(tǒng)的平衡點的穩(wěn)定性分類Fig. 5 Stability classification of the equilibrium points ofelectromechanical transmission system

    已知例C1中,P0(0,0)是鞍點,P1(x1,0)和P2(x2,0)是中心點。系統(tǒng)Hamilton量為:

    (19)

    當(dāng)H(0,0)=0時,存在兩條連接鞍點(0,0)的同宿軌道{q0(t)}。同宿軌道求解過程如式(20):

    (20)

    由式(20)可得:

    (21)

    積分得:

    (22)

    當(dāng)t=0,x0=±1時,c=0。整理后得兩條同宿軌的參數(shù)方程為:

    (23)

    3 混沌閾值

    當(dāng)擾動項增加到Hamilton系統(tǒng)時,封閉的同宿軌將被破壞,此時可能會導(dǎo)致分岔,全局同宿分岔預(yù)示著混沌運動。Melnikov方法可以估算存在同宿分岔的參數(shù)空間,并可以轉(zhuǎn)化為混沌,從而確定機(jī)電耦合扭轉(zhuǎn)振動混沌運動的閾值[6-7]。因此,在分析同宿分岔時,將系統(tǒng)(15)寫成狀態(tài)方程的形式為:

    (24)

    Melnikov方法實際上是通過計算雙曲不動點的穩(wěn)定和不穩(wěn)定流形之間的距離來確定系統(tǒng)是否存在橫截同宿點。根據(jù)Smale-birkhoff同宿理論,當(dāng)μ≠0,f≠0,且足夠小時,系統(tǒng)(15)可能存在Smale馬蹄存在意義下的混沌。根據(jù)Melinkov理論,Melnikov方程為:

    (25)

    式中:q0(τ)=(x±,y±),代表Hamilton系統(tǒng)的同宿軌道;F∧G=f1g2-f2g1。

    將式(24)代入式(25)中,初始條件為t=0,x0=±1,y0=0時,y±(τ)是時間τ的奇函數(shù)。則可得到:

    (26)

    式中,I1和I2分別為部分積分:

    所以,對于一定頻率,當(dāng)滿足式(27)時,系統(tǒng)可能產(chǎn)生Smale馬蹄變換意義下的混沌。

    (27)

    式(26)第一項積分I1可通過數(shù)學(xué)軟件計算直接求得I1=0.372 7;第二項積分I2在求解過程相當(dāng)復(fù)雜,所以采用數(shù)值積分法,通過采用數(shù)值積分法畫出混沌閾值曲線,然后根據(jù)閾值曲線同樣能夠?qū)ο到y(tǒng)做定性和定量分析。Melnikov函數(shù)數(shù)值積分法求解過程如下:

    (28)

    對式(28)分離變量,并對兩邊同時積分,得到:

    (29)

    將式(28)和式(29)代入I2中得到:

    (30)

    式(30)可以計算得到對應(yīng)ω值的I2值,代入式(31)即可得出隨ω變化的Melnikov閾值曲線。以f作為控制參數(shù)時,Melnikov閾值曲線如圖6。在一定條件下,在f-ω平面上存在混沌區(qū)域。當(dāng)f在該區(qū)域取值時,系統(tǒng)存在Smale馬蹄變換意義下的混沌。

    圖6 同宿分岔的閾值曲線(μ=2.5)Fig. 6 Threshold curve of homoclinic bifurcation (μ=2.5)

    4 數(shù)值計算分析

    利用Runge-Kutta法對擾動系統(tǒng)(15)進(jìn)行數(shù)值仿真,來進(jìn)行驗證上述用Melnikov方法做出的判斷是否正確,從而對得到對應(yīng)的時域波形、相軌跡、龐加萊截面圖和幅值譜進(jìn)行分析。圖7為μ=2.5,ω=1.0,ε=0.1,時,F(xiàn)=εf,以f為控制參數(shù)時的系統(tǒng)分岔圖,初始條件為x=0,y=0。

    圖7 f為控制參數(shù)時的系統(tǒng)分叉Fig. 7 Bifurcation diagram of the system when f is thecontrol parameter

    從圖7中可以看出,當(dāng)f為控制參數(shù)時,機(jī)電耦合非線性扭振系統(tǒng)表現(xiàn)出復(fù)雜的非線性振動現(xiàn)象[8]:當(dāng)激勵值f較小時(f<3.6),系統(tǒng)為1周期振蕩;當(dāng)激勵值f增大到3.6~3.8區(qū)間,2周期振蕩取代1周期振蕩;當(dāng)f>3.8時,控制參數(shù)f進(jìn)一步增加,系統(tǒng)在倍周期分岔后進(jìn)入混沌運動;當(dāng)f>3.8時,機(jī)電耦合系統(tǒng)表現(xiàn)為周期運動與混沌運動交替存在。

    便于更好的理解,選取3組數(shù)值對周期運動變?yōu)榛煦邕\動的過程進(jìn)行具體分析系統(tǒng),分別選取控制參數(shù)f=2,f=3.2,f=4,如圖6(b)。當(dāng)ω=1.0時,控制參數(shù)f由0逐漸增加,橫截同宿點橫截交叉發(fā)生于f≥3.6處。圖8~圖10表示了3種情況下時域波形、相軌跡、龐加萊截面和幅值譜。f=2對應(yīng)的點處于閾值曲線之下,沒有達(dá)到混沌運動的前提條件,如圖6。

    圖8 f=2時1周期運動Fig. 8 1 period motion graph when f=2

    圖9 f=3.2時2周期運動Fig. 9 2 period motion graph when f=3.2

    圖10 f=4時混沌態(tài)Fig. 10 Chaotic state graph when f=4

    由圖8可以看出,f=2時,機(jī)電耦合扭轉(zhuǎn)系統(tǒng)為1周期運動;當(dāng)f開始增大并超過臨界閾值,系統(tǒng)開始失去穩(wěn)定性,并存在發(fā)生混沌運動可能性。從圖9可以看出,f=3.2時,此時運行工況點處于閾值曲線之上,系統(tǒng)運動為2周期運動,系統(tǒng)發(fā)生倍周期分岔,但并未發(fā)生混沌運動。從圖10可以看出:f=4時,運行工況點處于閾值曲線之上,進(jìn)入如圖7中的混沌區(qū)域,其計算結(jié)果如圖10;時域波形并無周期性,相軌跡為多個交叉的環(huán)形圈,龐加萊截面圖具有雜亂無章的多個點,頻譜分析結(jié)果具有能量集中區(qū)域,這種現(xiàn)象表明該情況下系統(tǒng)發(fā)生混沌運動[9-10],這與我們之前Melnikov方法解析得到的結(jié)果是一致的。

    5 結(jié) 論

    以永磁同步電機(jī)的機(jī)電復(fù)合傳動系統(tǒng)為研究對象,考慮機(jī)電耦合關(guān)系,建立永磁同步電機(jī)的電磁轉(zhuǎn)矩模型,并通過非線性動力學(xué)理論得到電磁參數(shù)作用下的機(jī)電耦合參數(shù)振動界限以及規(guī)律。分析得到如下結(jié)論:

    1)產(chǎn)生非線性振動的主要原因是由于車用永磁同步電機(jī)扭轉(zhuǎn)角引起的電磁激勵。

    2)考慮機(jī)電耦合作用,扭轉(zhuǎn)振動系統(tǒng)的固有特性與電機(jī)的運行狀態(tài)以及結(jié)構(gòu)參數(shù)相關(guān),且在某些情況下,固有頻率有較大的降幅。

    3)考慮到車用電機(jī)采用弱磁控制,并且固有頻率受到內(nèi)功率因數(shù)角的影響,因此車用永磁同步電機(jī)需要合理控制內(nèi)功率因數(shù)角的范圍。

    4)純電動拖拉機(jī)機(jī)電耦合傳動系統(tǒng)的控制參數(shù)激勵值f逐漸增大時,系統(tǒng)具有周期、擬周期和混沌等復(fù)雜的動力學(xué)行為。當(dāng)出現(xiàn)非線性振動時,機(jī)電耦合傳動系統(tǒng)受到非周期沖擊干擾力作用,產(chǎn)生較大強(qiáng)度動載荷并引起非線性扭振,嚴(yán)重影響純電動拖拉機(jī)機(jī)電耦合傳動系統(tǒng)可靠性和耐久性。因此,應(yīng)避免純電動拖拉機(jī)在低速重載工況下長時間運行。

    猜你喜歡
    控制參數(shù)傳動系統(tǒng)機(jī)電
    佛蘭德傳動系統(tǒng)有限公司
    中國造紙(2022年8期)2022-11-24 09:43:38
    高超聲速飛行器滑??刂茀?shù)整定方法設(shè)計*
    飛控與探測(2022年6期)2022-03-20 02:16:14
    拖拉機(jī)內(nèi)燃電傳動系統(tǒng)
    Birkhoff系統(tǒng)穩(wěn)定性的動力學(xué)控制1)
    簡述機(jī)電自動化在工程機(jī)械制造中的應(yīng)用
    高速公路站級機(jī)電維護(hù)管理模式創(chuàng)新探討
    機(jī)械制造過程中機(jī)電自動化的應(yīng)用
    電子制作(2018年10期)2018-08-04 03:24:44
    ACS6000中壓傳動系統(tǒng)在鋼管軋制中的應(yīng)用
    電子制作(2018年10期)2018-08-04 03:24:26
    基于PI與準(zhǔn)PR調(diào)節(jié)的并網(wǎng)逆變器控制參數(shù)設(shè)計
    黑龍江電力(2017年1期)2017-05-17 04:25:08
    一體化污泥干化焚燒裝置工藝控制參數(shù)的優(yōu)化及運行
    av女优亚洲男人天堂| 又黄又爽又刺激的免费视频.| 久久久久久久久中文| 成人特级黄色片久久久久久久| 久久人人精品亚洲av| 香蕉av资源在线| 国产高清有码在线观看视频| 亚洲成a人片在线一区二区| 老司机福利观看| 亚洲欧美日韩高清在线视频| 国产极品精品免费视频能看的| 亚洲乱码一区二区免费版| 国产男人的电影天堂91| 亚洲四区av| 如何舔出高潮| 韩国av一区二区三区四区| 国产女主播在线喷水免费视频网站 | 欧美最黄视频在线播放免费| 午夜福利在线观看免费完整高清在 | 国产精品99久久久久久久久| 午夜福利18| 亚洲国产欧美人成| 欧美xxxx性猛交bbbb| 国产大屁股一区二区在线视频| 老熟妇乱子伦视频在线观看| 午夜a级毛片| 免费观看人在逋| 国内精品久久久久久久电影| 中文资源天堂在线| 最好的美女福利视频网| 国内毛片毛片毛片毛片毛片| 女人被狂操c到高潮| 午夜精品一区二区三区免费看| 三级男女做爰猛烈吃奶摸视频| 美女免费视频网站| 国产视频内射| 免费黄网站久久成人精品| 两性午夜刺激爽爽歪歪视频在线观看| 日本色播在线视频| 女人十人毛片免费观看3o分钟| 亚洲成a人片在线一区二区| 少妇的逼水好多| 1000部很黄的大片| 成人特级av手机在线观看| 在现免费观看毛片| 99热精品在线国产| 午夜福利高清视频| 国产精品爽爽va在线观看网站| 麻豆国产av国片精品| 又黄又爽又刺激的免费视频.| 日本成人三级电影网站| 丰满人妻一区二区三区视频av| 国模一区二区三区四区视频| 国产高清三级在线| 国产精品一区二区性色av| 久久久成人免费电影| 亚洲欧美日韩卡通动漫| 人妻久久中文字幕网| 特级一级黄色大片| avwww免费| 最近中文字幕高清免费大全6 | 午夜福利在线观看免费完整高清在 | 国产精品98久久久久久宅男小说| 免费看美女性在线毛片视频| 精品午夜福利视频在线观看一区| 白带黄色成豆腐渣| 精华霜和精华液先用哪个| 免费观看人在逋| 性欧美人与动物交配| 中文字幕高清在线视频| 免费看a级黄色片| 中文字幕熟女人妻在线| 成年女人毛片免费观看观看9| 日韩欧美在线乱码| 国产精品无大码| 国产精品女同一区二区软件 | 一区福利在线观看| 午夜老司机福利剧场| 亚洲第一电影网av| 亚洲国产高清在线一区二区三| 日日撸夜夜添| 亚洲国产色片| 欧美日本视频| 少妇被粗大猛烈的视频| 日韩欧美精品免费久久| 黄色欧美视频在线观看| 国产三级中文精品| 一夜夜www| 男人狂女人下面高潮的视频| .国产精品久久| 久久久成人免费电影| www.www免费av| 国产在视频线在精品| 在线观看免费视频日本深夜| 少妇的逼好多水| 简卡轻食公司| 久久精品国产亚洲av涩爱 | 欧美一区二区精品小视频在线| 哪里可以看免费的av片| xxxwww97欧美| 国产白丝娇喘喷水9色精品| 欧美bdsm另类| 男人的好看免费观看在线视频| 九色国产91popny在线| 午夜福利在线观看免费完整高清在 | 夜夜看夜夜爽夜夜摸| 国产成年人精品一区二区| 亚洲四区av| 99热精品在线国产| 成人三级黄色视频| www.www免费av| 午夜精品一区二区三区免费看| 国产高潮美女av| 亚洲av日韩精品久久久久久密| 久久久久性生活片| 伊人久久精品亚洲午夜| 高清日韩中文字幕在线| 18禁黄网站禁片午夜丰满| 99久国产av精品| 色5月婷婷丁香| 欧美三级亚洲精品| 亚洲av.av天堂| 亚洲四区av| 99久久精品一区二区三区| 99国产极品粉嫩在线观看| 乱人视频在线观看| 欧美另类亚洲清纯唯美| 国产一级毛片七仙女欲春2| 久久久久久国产a免费观看| 色视频www国产| 又黄又爽又免费观看的视频| 亚洲av免费在线观看| 乱系列少妇在线播放| 国产欧美日韩精品一区二区| 午夜免费成人在线视频| 久久天躁狠狠躁夜夜2o2o| 精品久久久久久久久久免费视频| 欧美+亚洲+日韩+国产| 观看美女的网站| 久久精品国产亚洲av天美| 欧美三级亚洲精品| 成人特级av手机在线观看| 久久精品国产99精品国产亚洲性色| 国产精品国产高清国产av| 美女黄网站色视频| 欧美xxxx性猛交bbbb| 国产探花极品一区二区| 亚洲av免费在线观看| 少妇人妻精品综合一区二区 | 亚洲欧美精品综合久久99| 变态另类成人亚洲欧美熟女| 亚洲最大成人av| 精品日产1卡2卡| 小说图片视频综合网站| 舔av片在线| 亚洲人成网站在线播| 动漫黄色视频在线观看| 69人妻影院| 久久久久久久久久久丰满 | 性插视频无遮挡在线免费观看| 久久人妻av系列| 大型黄色视频在线免费观看| 久久午夜福利片| 中文字幕久久专区| 黄片wwwwww| 男人的好看免费观看在线视频| 亚洲av中文字字幕乱码综合| 国产主播在线观看一区二区| 在线免费十八禁| 国产精品亚洲美女久久久| 久久这里只有精品中国| 在线免费观看不下载黄p国产 | 春色校园在线视频观看| 熟妇人妻久久中文字幕3abv| 国产欧美日韩精品亚洲av| 两人在一起打扑克的视频| .国产精品久久| 男女边吃奶边做爰视频| 精品午夜福利视频在线观看一区| 国产极品精品免费视频能看的| 久久久精品欧美日韩精品| 亚洲美女搞黄在线观看 | 国产男人的电影天堂91| 精品99又大又爽又粗少妇毛片 | 两人在一起打扑克的视频| 黄色日韩在线| 欧美一区二区亚洲| 日韩大尺度精品在线看网址| 91av网一区二区| 成年人黄色毛片网站| 国产伦在线观看视频一区| 免费av毛片视频| 国产美女午夜福利| 九九在线视频观看精品| 久久精品综合一区二区三区| 级片在线观看| av在线观看视频网站免费| 丝袜美腿在线中文| 日韩在线高清观看一区二区三区 | 欧美丝袜亚洲另类 | 免费一级毛片在线播放高清视频| 欧美日韩精品成人综合77777| 18禁裸乳无遮挡免费网站照片| 免费av观看视频| 丰满的人妻完整版| 久久久国产成人精品二区| 在线观看免费视频日本深夜| 精品久久国产蜜桃| 日本在线视频免费播放| 国产三级在线视频| 1024手机看黄色片| 看十八女毛片水多多多| 免费观看人在逋| 国产精品一区二区免费欧美| 国产一区二区在线观看日韩| 国产日本99.免费观看| 国语自产精品视频在线第100页| 亚洲内射少妇av| 色av中文字幕| 最好的美女福利视频网| 亚洲欧美日韩高清专用| 一个人看的www免费观看视频| 女生性感内裤真人,穿戴方法视频| 色av中文字幕| 国产精品福利在线免费观看| 能在线免费观看的黄片| 免费观看在线日韩| 久久亚洲真实| 伦精品一区二区三区| 在线免费观看的www视频| 亚洲一区高清亚洲精品| 亚洲精品亚洲一区二区| 国产精品嫩草影院av在线观看 | 欧美日韩中文字幕国产精品一区二区三区| 级片在线观看| avwww免费| 一级黄片播放器| 国产老妇女一区| 免费大片18禁| 中文字幕高清在线视频| 久久精品影院6| 一卡2卡三卡四卡精品乱码亚洲| 欧美三级亚洲精品| 欧美成人一区二区免费高清观看| 91在线观看av| 琪琪午夜伦伦电影理论片6080| 免费黄网站久久成人精品| 欧美+亚洲+日韩+国产| 久久久久国内视频| 网址你懂的国产日韩在线| 成人特级av手机在线观看| 亚洲天堂国产精品一区在线| 国产高清激情床上av| 亚洲人成伊人成综合网2020| 国产精品综合久久久久久久免费| 日本一二三区视频观看| 亚洲精品一卡2卡三卡4卡5卡| 日韩强制内射视频| 国产又黄又爽又无遮挡在线| 国产 一区精品| 午夜精品久久久久久毛片777| 波野结衣二区三区在线| 亚洲欧美日韩卡通动漫| 国产伦一二天堂av在线观看| 亚洲精品国产成人久久av| 国产综合懂色| 国产精品久久久久久精品电影| 欧美丝袜亚洲另类 | 性欧美人与动物交配| 日韩欧美 国产精品| 男女视频在线观看网站免费| 欧美日韩瑟瑟在线播放| 亚洲av一区综合| 国产精品电影一区二区三区| 五月玫瑰六月丁香| 国产综合懂色| 一进一出好大好爽视频| 日本在线视频免费播放| 深夜a级毛片| 18禁黄网站禁片午夜丰满| 在线天堂最新版资源| 久久久久久久精品吃奶| 日本熟妇午夜| a级毛片a级免费在线| 男女那种视频在线观看| 国产精品久久久久久精品电影| 亚洲精品456在线播放app | 蜜桃亚洲精品一区二区三区| 亚洲av中文av极速乱 | av在线观看视频网站免费| 1024手机看黄色片| 日本黄大片高清| 欧美在线一区亚洲| 看免费成人av毛片| 日日干狠狠操夜夜爽| 日韩欧美在线乱码| 淫秽高清视频在线观看| 久久久久久久久中文| 亚洲国产欧洲综合997久久,| 亚洲在线自拍视频| 国产成年人精品一区二区| 直男gayav资源| 亚洲欧美日韩无卡精品| 两人在一起打扑克的视频| 国产精品综合久久久久久久免费| av女优亚洲男人天堂| 亚洲中文字幕日韩| 夜夜爽天天搞| 免费无遮挡裸体视频| 99九九线精品视频在线观看视频| 日韩欧美在线乱码| 国产成人福利小说| 免费人成在线观看视频色| 久久久久久久久久黄片| 最近最新中文字幕大全电影3| 国产不卡一卡二| 国产视频内射| 亚洲乱码一区二区免费版| 51国产日韩欧美| 精品一区二区三区人妻视频| 亚洲最大成人手机在线| 亚洲国产高清在线一区二区三| 色精品久久人妻99蜜桃| 波野结衣二区三区在线| 国模一区二区三区四区视频| 欧美日韩综合久久久久久 | 国产精品亚洲美女久久久| 综合色av麻豆| 18+在线观看网站| 亚洲精品粉嫩美女一区| 日本一本二区三区精品| 国产精品一区二区三区四区久久| 成人美女网站在线观看视频| 久久人人爽人人爽人人片va| 欧美潮喷喷水| 国产一区二区亚洲精品在线观看| 免费一级毛片在线播放高清视频| 欧美日韩中文字幕国产精品一区二区三区| 亚洲中文日韩欧美视频| 精品人妻视频免费看| 最后的刺客免费高清国语| 国产乱人伦免费视频| 男女边吃奶边做爰视频| 黄色丝袜av网址大全| 欧美日韩精品成人综合77777| 欧美中文日本在线观看视频| 自拍偷自拍亚洲精品老妇| 国内精品久久久久久久电影| 2021天堂中文幕一二区在线观| 中文字幕精品亚洲无线码一区| 悠悠久久av| 1024手机看黄色片| 日本免费一区二区三区高清不卡| 国产淫片久久久久久久久| 永久网站在线| 日韩亚洲欧美综合| 国产爱豆传媒在线观看| 三级毛片av免费| 午夜福利高清视频| 很黄的视频免费| 又粗又爽又猛毛片免费看| 成人性生交大片免费视频hd| 国内少妇人妻偷人精品xxx网站| 麻豆成人午夜福利视频| 亚洲专区中文字幕在线| 韩国av在线不卡| 波多野结衣高清无吗| 熟女人妻精品中文字幕| 男女啪啪激烈高潮av片| 亚洲成a人片在线一区二区| 99国产极品粉嫩在线观看| 免费在线观看成人毛片| 午夜影院日韩av| 99热这里只有是精品在线观看| 人人妻,人人澡人人爽秒播| bbb黄色大片| 日韩精品青青久久久久久| 少妇人妻精品综合一区二区 | 免费av不卡在线播放| 免费人成视频x8x8入口观看| 国语自产精品视频在线第100页| 毛片一级片免费看久久久久 | 男女之事视频高清在线观看| 不卡视频在线观看欧美| 性插视频无遮挡在线免费观看| 国产精品98久久久久久宅男小说| 俺也久久电影网| 国产精品久久久久久久电影| 夜夜爽天天搞| 成人特级黄色片久久久久久久| 久久国内精品自在自线图片| 我的老师免费观看完整版| 别揉我奶头~嗯~啊~动态视频| 99久久久亚洲精品蜜臀av| 熟女人妻精品中文字幕| 人人妻人人澡欧美一区二区| 亚洲精品456在线播放app | 真实男女啪啪啪动态图| 亚洲av成人av| 国产亚洲精品综合一区在线观看| 日本免费a在线| 免费看光身美女| 日韩av在线大香蕉| 国产精品亚洲一级av第二区| 在现免费观看毛片| 国产三级中文精品| 国产精品伦人一区二区| 热99在线观看视频| 十八禁国产超污无遮挡网站| 99在线人妻在线中文字幕| 国产欧美日韩精品一区二区| 18禁在线播放成人免费| 老熟妇乱子伦视频在线观看| 日本一二三区视频观看| 在线观看av片永久免费下载| 成人鲁丝片一二三区免费| 婷婷色综合大香蕉| 中文字幕人妻熟人妻熟丝袜美| 免费不卡的大黄色大毛片视频在线观看 | 麻豆一二三区av精品| 淫秽高清视频在线观看| 亚洲最大成人中文| www.www免费av| 欧美日韩黄片免| 亚洲成人久久爱视频| 日本五十路高清| 国产单亲对白刺激| 国产伦精品一区二区三区视频9| 国产av一区在线观看免费| 美女黄网站色视频| 波多野结衣巨乳人妻| 听说在线观看完整版免费高清| 成人永久免费在线观看视频| 色播亚洲综合网| 久久香蕉精品热| 精品一区二区三区视频在线| 日本与韩国留学比较| 亚洲欧美精品综合久久99| 国产蜜桃级精品一区二区三区| 国产精品一区二区免费欧美| 淫秽高清视频在线观看| 国产毛片a区久久久久| 日韩高清综合在线| 在线观看av片永久免费下载| 老司机午夜福利在线观看视频| 国内揄拍国产精品人妻在线| 九色国产91popny在线| 有码 亚洲区| 亚洲无线观看免费| 色播亚洲综合网| 亚洲电影在线观看av| 一进一出抽搐动态| 欧美日韩精品成人综合77777| 黄色丝袜av网址大全| 国产亚洲91精品色在线| 赤兔流量卡办理| 舔av片在线| 日韩精品中文字幕看吧| 亚洲欧美清纯卡通| 亚洲 国产 在线| 桃色一区二区三区在线观看| 搡女人真爽免费视频火全软件 | 久久久久免费精品人妻一区二区| 亚洲黑人精品在线| 我的老师免费观看完整版| 在线看三级毛片| 久久精品久久久久久噜噜老黄 | 噜噜噜噜噜久久久久久91| 少妇熟女aⅴ在线视频| 色精品久久人妻99蜜桃| 日本免费一区二区三区高清不卡| 如何舔出高潮| 亚洲精品色激情综合| 午夜免费男女啪啪视频观看 | 人妻久久中文字幕网| 国产精品免费一区二区三区在线| 人人妻人人看人人澡| 69人妻影院| 国产熟女欧美一区二区| 一级黄片播放器| 在线观看一区二区三区| 淫妇啪啪啪对白视频| 噜噜噜噜噜久久久久久91| 99久久九九国产精品国产免费| 亚洲四区av| 亚洲真实伦在线观看| 国产精品久久电影中文字幕| АⅤ资源中文在线天堂| 欧美+日韩+精品| 久久午夜福利片| 欧美一区二区精品小视频在线| 国产单亲对白刺激| 99热只有精品国产| 日韩在线高清观看一区二区三区 | 亚洲,欧美,日韩| 亚洲精品久久国产高清桃花| 此物有八面人人有两片| 亚洲久久久久久中文字幕| 国产高清视频在线观看网站| 久久久久久久久久成人| 春色校园在线视频观看| 自拍偷自拍亚洲精品老妇| 国产精品永久免费网站| 丰满乱子伦码专区| 婷婷亚洲欧美| 午夜爱爱视频在线播放| 在线观看一区二区三区| 可以在线观看毛片的网站| 国产黄色小视频在线观看| 九九在线视频观看精品| 麻豆精品久久久久久蜜桃| 久久精品夜夜夜夜夜久久蜜豆| 啦啦啦观看免费观看视频高清| 制服丝袜大香蕉在线| 99视频精品全部免费 在线| 99热6这里只有精品| 草草在线视频免费看| 亚洲人成网站高清观看| 亚洲av一区综合| 久久精品91蜜桃| xxxwww97欧美| 国产精品免费一区二区三区在线| 韩国av一区二区三区四区| 97热精品久久久久久| 亚洲 国产 在线| 九九爱精品视频在线观看| 男女边吃奶边做爰视频| 国产黄片美女视频| 国产精品一区www在线观看 | 99精品在免费线老司机午夜| 国产高清有码在线观看视频| 亚洲成人久久性| 99久国产av精品| 欧美精品国产亚洲| 伊人久久精品亚洲午夜| 99在线视频只有这里精品首页| 午夜福利18| 久久久久久久亚洲中文字幕| 嫩草影院新地址| 麻豆久久精品国产亚洲av| 蜜桃亚洲精品一区二区三区| 久久人人爽人人爽人人片va| 国产精品一区二区三区四区免费观看 | 久99久视频精品免费| avwww免费| 精品久久久久久久人妻蜜臀av| 天天躁日日操中文字幕| 久久热精品热| 久久精品国产鲁丝片午夜精品 | 欧美最新免费一区二区三区| 在线免费观看不下载黄p国产 | 少妇的逼水好多| 午夜福利在线观看吧| 成人国产一区最新在线观看| av黄色大香蕉| 偷拍熟女少妇极品色| 精品一区二区三区视频在线观看免费| 最新在线观看一区二区三区| 精品午夜福利在线看| 欧美日韩瑟瑟在线播放| 可以在线观看的亚洲视频| 国产伦精品一区二区三区视频9| 午夜免费男女啪啪视频观看 | 十八禁国产超污无遮挡网站| 嫩草影院新地址| 国产高潮美女av| 亚洲精品乱码久久久v下载方式| 淫秽高清视频在线观看| 国产熟女欧美一区二区| 亚洲精品456在线播放app | 亚洲欧美日韩东京热| 成人性生交大片免费视频hd| 久久99热这里只有精品18| 91在线观看av| 亚洲av免费高清在线观看| 午夜福利欧美成人| av专区在线播放| 亚洲黑人精品在线| 中文字幕av成人在线电影| 久久精品国产亚洲av天美| 亚洲无线在线观看| 亚洲专区国产一区二区| 很黄的视频免费| 高清日韩中文字幕在线| 免费av不卡在线播放| xxxwww97欧美| 波多野结衣高清作品| 婷婷亚洲欧美| 国产麻豆成人av免费视频| 国产免费av片在线观看野外av| 国产精品人妻久久久久久| 国产成人影院久久av| 国产午夜精品久久久久久一区二区三区 | videossex国产| 嫩草影视91久久| 99久久成人亚洲精品观看| 噜噜噜噜噜久久久久久91| 国产熟女欧美一区二区| 成人二区视频| 国产成人福利小说| 国产精品1区2区在线观看.| 三级毛片av免费| 久久久成人免费电影| 中文资源天堂在线| 男人的好看免费观看在线视频| 深夜a级毛片| 免费人成在线观看视频色| 国产精品久久久久久av不卡| 午夜福利在线在线| 真实男女啪啪啪动态图| 69人妻影院| 日韩人妻高清精品专区| 99久久成人亚洲精品观看| 亚洲国产欧美人成| 国产精品1区2区在线观看.|