午夜1000集I99久久香蕉I三级视频网I欧美偷拍一区二区Iav噜噜噜I国产一二区视频I国产精品视频网址I久久精avI国产三级精品视频I天天操夜夜爱I91麻豆国产福利在线观看I性欧美日韩Iwww.九色91I日本a在线观看I涩avI成年人拍拍视频I第四色成人网I麻豆乱淫一区二区三区I免费手机avI香蕉视频官方网站

資料下載首頁 > 資料下載 > 單座式調節閥閥芯耦合振動

單座式調節閥閥芯耦合振動

提 供 商: 上海申弘閥門有限公司 資料大小: JPG
圖片類型: JPG 下載次數: 354 次
資料類型: JPG 瀏覽次數: 2213次
下載鏈接: 文件下載    圖片下載    
詳細介紹

                          單座式調節閥閥芯耦合振動

                        上海申弘閥門有限公司

      針對單座式調節閥閥芯-閥桿系統的流固耦合振動問題,建立了考慮閥門定位器作用的系統動態仿真模型,給出了求解閥芯-閥桿系統響應的預估-校正算法,利用ANSYS軟件對系統在固定開度與變開度情況和流開型與流閉型情況下的振動響應進行了定性分析。研究表明:對于固定開度,閥芯穩態位移偏移量以及流體力隨壓差增加以及開度減小而增大,閥芯動態位移過渡時間隨壓差增加而縮短,而流開型流向時閥芯位移響應幅度要大于流閉型流向。對于開度減小情況,采用流閉型流向時閥芯動態位移過渡時間更短,且壓差越大,閥芯動態位移過渡時間越短,而流開型流向時則相反;對于開度增大情況,閥芯動態位移過渡時間和壓差之間的關系與開度減小情況相反。調節閥在水利、電站、化工、石油、冶金等過程控制系統中起著重要作用,然而調節閥在某些工況下產生的振動往往成為引起各種事故的主要原因,振動嚴重時甚至影響系統安全平穩地運行。導致調節閥振動的主要原因是閥體內部流體流動的不穩定性,與流固耦合作用直接相關。上海申弘閥門有限公司主營閥門有:截止閥,電動截止閥,氣動截止閥,電動蝶閥,氣動蝶閥,電動球閥,氣動球閥,電動閘閥,氣動閘閥,電動調節閥,氣動調節閥,減壓閥水力控制閥真空閥門、襯膠閥門、襯氟閥門。流固耦合導致的調節閥振動可以分為兩種:一種是調節閥整體在管道上的振動;另一種是本文要研究的調節閥自身閥芯-閥桿系統振動。國內外有關閥流固耦合振動的研究較少,且大多偏重于在某些固定開度、單*向工況下分析閥體及流道內的瞬態流場及其與閥體或閥芯、閥桿的彈性振動之間的流固耦合相互作用,例如分別對調節閥、電磁閥、截止閥、換向閥以及溢流閥等不同類型的閥的動態特性與流固耦合問題進行了仿真和試驗研究。

本文以某型號單座式調節閥為對象,研究調節閥閥芯-閥桿系統在流開式與流閉式流向、固定開度與變開度工況時的流固耦合振動問題。通過建立考慮定位器作用的系統有限元流固耦合模型并進行仿真,分析流開型和流閉型調節閥閥芯-閥桿系統在固定開度和變開度條件下的動態響應。

1、閥芯-閥桿系統動力學模型
1.1、單座式調節閥結構

單座式調節閥整體結構如圖1所示,主要由執行機構和調節機構組成。執行機構將控制器輸出電流信號轉換為調節閥閥桿的直線位移;調節機構通過閥桿與執行機構相連,并將位移信號轉換為閥芯和閥座之間流通面積的變化,從而改變流經閥體的流體運動狀態。為了使閥芯和閥座之間的開度穩定在某個特定位置以保證實現控制器所要求的目標狀態,調節閥設有閥門定位器附件,組成以閥桿實際位移為測量信號、以控制器電流輸出為設定信號的反饋控制系統。電子式電動單座調節閥依靠流經閥內介質自身的壓力、溫度作為能源驅動閥門自動工作,不需要外接電源和二次儀表。這種自力式調節閥都利用閥輸出端的反饋信號(壓力、壓差、溫度)通過信號管傳遞到執行機構驅動閥瓣改變閥門的開度,達到調節壓力、流量、溫度的目的。電子式電動單座調節閥由PS系列和3810系列直行程電動執行機構和低流阻直通單座閥組成。(也可配其它品牌直行程執行機構)。電動執行機構為電子式一體化結構,內有伺服放大器,輸入控制信號(4-20mADC或1-5VDC)及電源即可控制閥門開度,達到對壓力、流量、液位、溫度等工參數的調節。具有動作靈敏、連線簡單、流量大、體積小、調節精度高等特點。控制精度和性能比DKZ型有明顯提高。單座閥適用于對泄露量求嚴格、閥前后壓差低及有一定粘度的場合。廣泛應用于電力、冶金、化工、石油、環保、輕工、設備等行業。電子式電動單座調節閥是一種壓力平衡式的調節閥。接受調節儀表來的直流電流信號,改變被調介質流量,使被控工藝參數保持在給定值。廣泛應用于電力、冶金、化工、石油、輕紡、制藥、造紙等工業部門的生產自動化控制。
電子式電動單座調節閥的特點是:
1、可調節范圍大,固有可調比為50,流量特性有直線和等百分比。
2、閥體按流體力學原理設計的等截面低流阻流道,額定流量系數增大30%。
3、調節切斷型采用軟密封結構閥芯,達VI級泄漏標準(零泄漏)。
4、伺服放大器采用深度動態負反饋,可提高自動調節精度。
5、電動操作器有多種形式,可適用于4~20mA。DC或0~10mA.DC。
6、波紋管密封型調節閥對移動的閥桿形成*的密封,堵絕流體外漏。
7、電子式電動單座調節閥可直接由電流信號控制閥門開度,無需伺服放大器。

單座式調節閥結構示意圖
圖1 單座式調節閥結構示意圖

1.彈簧2.膜片3.氣室4.位移刻度5.閥桿6.填料7.閥芯8.閥座9.閥體10.調節機構11.執行機構1.2 閥芯-閥桿系統動力學模型

閥芯-閥桿系統是調節閥執行-調節機構的主體,主要包括薄膜膜片、推桿、彈簧、填料、閥桿、閥芯和閥座等,如圖1所示。

對于特定調節閥結構,其調節性能主要取決于閥芯-閥桿相對閥座的運動。閥芯-閥桿與氣室彈簧、填料、流體等組成一個流固耦合動力學系統。在建立閥芯-閥桿動力學模型之前作如下假設:調節閥內部的流體(水)是不可壓縮的;在調節閥工作過程中,水沒有熱交換;閥桿與填料之間的摩擦力認為是一種粘性阻尼力;整個閥體為剛性體,忽略其彈性變形;閥芯-閥桿只沿軸線運動,不考慮橫向運動;單座閥的執行機構為氣動正作用執行機構。

根據以上假設以及單座閥的結構圖,可以將閥芯-閥桿系統簡化為一個單自由度質量-彈簧-阻尼系統,如圖2所示。閥芯位移坐標原點設在閥門全關的位置,取向上為正方向。

單座式調節閥結構示意圖
圖2 閥芯-閥桿系統動力學模型

根據圖2所示,單座閥閥芯-閥桿系統動力學方程為

單座式調節閥閥芯-閥桿系統流固耦合振動研究 (1) 式中
單座式調節閥閥芯-閥桿系統流固耦合振動研究
———閥芯t時刻的位移、速度及加速度

m———閥芯-閥桿系統的總質量

k———彈簧剛度系數

c———等效粘性阻尼系數

Fc(t)———執行機構作用在閥芯上的控制力

Fl(t)———流體不平衡力

g———重力加速度

式(1)中流體不平衡力Fl(t)(以下簡稱流體力)是指行程調節機構中流體作用在閥芯上的軸向合力,其大小與調節閥兩端的壓差、閥芯大小和形狀、調節閥類型、閥上游壓力、流體流向及流體物理特性等因素有關。對于流開式流向(圖1中流體從左側流入,右側流出),流體力為

單座式調節閥閥芯-閥桿系統流固耦合振動研究 (2)
式中 p1———閥前壓力 p2———閥后壓力

Δp———前后壓差 ds———閥芯直徑

dz———閥桿直徑

對于流閉式流向(圖1中流體從右側流入,左側流出),流體力為

單座式調節閥閥芯-閥桿系統流固耦合振動研究 (3)
對于動態問題,無論流開式還是流閉式流向,閥芯受到的實際流體力Fl(t)由閥芯上下分布的壓力(通過預估-校正有限元流固耦合計算)對閥芯表面面積積分求得。為了使求解過程簡化,可以將流體力近似認為流體壓力在閥芯等效橫截面積上作用的結果。

式(1)中控制力Fc(t)是指執行機構作用在調節機構閥芯-閥桿系統上的軸向力,其大小與電-氣閥門定位器結構參數、控制器輸出電流信號產生的電磁力以及閥芯位移有關。閥門定位器是調節閥的主要附件。它將閥桿位移信號作為輸入的反饋測量信號,以控制器輸出信號作為設定信號,將二者進行比較,當兩者有偏差時,改變其到執行機構的輸出信號,使執行機構動作,從而建立了閥桿位移信號與控制器輸出信號之間的一一對應關系,即

單座式調節閥閥芯-閥桿系統流固耦合振動研究 (4)
式中 p、q、T———與定位器結構性能有關的參數

N(X0)———與目標位移X0成比例的電磁力

因此,調節閥閥芯-閥桿系統動力學方程可寫為

單座式調節閥閥芯-閥桿系統流固耦合振動研究 (5)
2、閥芯-閥桿系統流固耦合問題求解
由于前面建立的調節閥閥芯-閥桿系統動力學方程中的流體力Fl(t)無法用準確的解析表達式表示出來,需要通過有限元流固耦合方法計算得到。本文采用ANSYS12.0中的CFD模塊對調節閥內部流場進行分析計算,首先建立內部流場的幾何模型,然后對不同工況下的調節閥流場進行計算,求得流體力。調節閥內部流場幾何模型可以通過三維建模軟件SolidWorks來建立,計算流場時需要對ANSYS中的CFD模塊的流固耦合方法進行改進。

2.1、流固耦合分析典型步驟

通過ANSYS軟件物理環境方法對閥芯-閥桿系統進行流固耦合分析。閥芯-閥桿會在流體力作用下發生位移,位移的大小將改變閥芯移動壁面邊界,從而顯著影響流場的形狀。流固耦合分析通過在結構分析中得到閥芯移動壁面位移,用于流場分析。具體流固耦合分析的典型步驟如下:

(1)創建整個幾何模型:包括流體區域和調節閥閥體結構區域。

(2)創建流體物理環境:給流體區域賦予單元類型,還要確定迭代次數,激活湍流模型,施加邊界條件。

(3)創建結構物理環境:清除在流體物理環境中設定的信息,準備定義結構物理環境。轉換單元類型并設定單元選項,將流體區域單元設定為NULL,將結構區域賦予單元類型,施加結構邊界條件,定義合適的載荷步和求解選項,然后寫入結構物理環境文件。

(4)流體/結構求解循環:在本系統中,入口的速度作為總體收斂的準則。當兩次Flotran求解的入口速度差值足夠小時,求解結束。初始Flotran分析設置的迭代次數應當多一些,以利于較好地收斂。隨后的流體分析由于是在前一次流體分析基礎上重啟動,因此,迭代次數可以少一些。結構分析同樣也需要重啟動。對于非線性分析,節點必須在重啟動以前恢復到初始位置。本文計算一共執行了5次耦合迭代,第1次Flotran分析共迭代100次,以得到較高的收斂精度,隨后的4次Flotran分析各迭代40次,即可滿足精度要求。一共迭代260次,取前250次數據即前0.25s的位移響應。

2.2、閥芯-閥桿系統流固耦合模型預估-校正求解算法

由于調節閥閥芯與流場接觸面為移動壁面,ANSYS軟件的典型流固耦合方法難以滿足計算要求,本文采用預估-校正算法[8]來解決這種情況下的流固耦合問題。根據閥芯-閥桿系統動力學方程(5)和預估-校正算法,可得出閥芯-閥桿系統的預測步和校正步計算表達式,分別為

預測步

單座式調節閥閥芯-閥桿系統流固耦合振動研究 (6)
單座式調節閥閥芯-閥桿系統流固耦合振動研究 (7)
單座式調節閥閥芯-閥桿系統流固耦合振動研究 (8)
第r步校正步

單座式調節閥閥芯-閥桿系統流固耦合振動研究 (9)
(10)
單座式調節閥閥芯-閥桿系統流固耦合振動研究
2.3、閥芯-閥桿系統流固耦合有限元模型建立

2.3.1、調節閥內部流場幾何模型創建

調節閥的內部流場就是調節閥內部充滿流體后所占的空間。假設閥體、閥芯是剛體,流場的邊界是閥體、閥芯與流體的耦合面,閥體與流體之間的耦合面是不動的,而閥芯與流體之間的耦合面是移動的(包括固定開度下的振動和變開度時的移動)。由于假設閥芯是剛體,所以在做流場分析時,需要將閥芯部分挖去,對于不同的開度,只需要將對應閥芯位移的不同位置挖去即可,由于閥桿對流場的影響比較小,所以在建立的模型中忽略閥桿。調節閥的內部流場模型如圖3所示。

調節閥內部流場模型
圖3 調節閥內部流場模型

2.3.2、網格劃分

本文采用自由網格劃分方法。在劃分網格之前首先要設置單元類型,并為實體模型分配單元屬性,本文中的流場模型采用的是Fluid142單元。閥芯與流體之間的耦合面流場變化比較大,采用zui大為0.002的網格劃分,而其他部位閥體與流體之間的耦合面,采用zui大為0.004的網格劃分,根據劃分結果,共有單元102306個,節點20531個。網格劃分結果如圖4所示。

流場模型網格劃分
圖4 流場模型網格劃分

2.3.3、載荷施加

在網格劃分之后,要對模型施加載荷以及邊界條件,本文中對模型施加的具體條件如下:

(1)入口:定義壓力為0.2~2.1MPa。

(2)出口:定義壓力為0.1MPa。

(3)固定壁面條件:所有與閥體接觸的耦合面上流體流動速度為零,相對位移為零,即除了進出口端面和閥芯周圍的耦合面外,其余邊界上速度和位移均為零。

(4)移動壁面條件:與閥芯接觸的耦合面上的流體速度和位移應該與閥芯的運動速度和位移一致,即當閥芯靜止時,耦合面上的流體速度和位移均為零;當閥芯運動時,耦合面上每個節點的任何時刻的速度和位移大小都與閥芯上對應節點的速度和位移大小相同。由于本文只考慮閥芯軸向運動,所以以上提到的位移和速度均指軸向位移和軸向速度,其余方向的位移和速度都為零。

(5)不考慮溫度變化,即不涉及熱交換。

3、閥芯-閥桿系統流固耦合動態響應分析
通過求解閥芯-閥桿系統流固耦合模型,仿真計算調節閥在不同邊界條件下閥芯受到的動態流體力以及閥芯-閥桿系統的動態位移響應,分析定位器對調節閥閥芯振動的影響,考察閥內介質流向不同時,閥芯在固定開度和變開度的振動情況,找出閥芯受到的動態力與位移、壓差之間的關系。

3.1、閥門定位器對閥芯-閥桿系統動態響應的影響

假設閥芯在目標位移15mm處,分別在進出口壓差為0.1、0.5、1.0、2.0MPa下工作,介質流向為流開狀態,初始相對位移為-0.1mm,初始速度和加速度均為零,初始流體力為50N。分別對比分析帶有閥門定位器和沒有閥門定位器時閥芯-閥桿系統的閥芯位移和流體力響應,如圖5、6所示。

單座閥閥芯的振動響應對比
圖5 單座閥閥芯的振動響應對比

(a)壓差0.1MPa(b)壓差0.5MPa(c)壓差1.0MPa(d)壓差2.0MPa

由圖5可得:帶有定位器比無定位器時閥芯振動波動小得多,而且趨近目標平衡位置的時間更短。在較大壓差(0.5~2.0MPa)下,無論是否帶有定位器,壓差越大,閥芯振動偏離目標位置越遠,并且振幅越大,而帶有定位器的調節閥工作狀態更為穩定、。

由圖6可以看出:隨著壓差增大,無定位器的閥芯受到的流體力波動也隨之增大,而帶有定位器的閥芯受到的流體力變化則一直保持平穩,并且逐漸從小于前者的均值過渡到大于其均值。

單座閥閥芯受到的流體力對比
圖6 單座閥閥芯受到的流體力對比

3.2、單座閥在固定開度下的動態特性分析

假設調節閥閥門全關位置為坐標原點,取向上為正方向,總行程為0(全關位置)~25mm(全開位置),進出口壓差變化范圍為0.1~2.0MPa,介質流向分別為流開和流閉。考察閥芯在目標位移(開度)作自由振動時閥芯的振動情況。

3.2.1、不同流向對單座閥動態特性的影響

在某個目標位移下,閥芯在受到初始運動條件的擾動時,會在平衡位置作自由振動。假設閥芯分別在不同目標位移(3、5、10、15、20、25mm)下和不同進出口壓差(0.1、0.5、1.0、2.0MPa)下工作,初始相對位移為-0.1mm,初始速度和加速度均為0,不同工況下初始流體力不同,考慮定位器作用,得出介質流向不同時單座閥閥芯的位移與閥芯受到的流體力隨時間變化歷程。圖7為某固定開度(目標位移為10mm)時閥芯位移與閥芯所受流體力的時間響應歷程。圖8為位移、壓差與流體力的三維關系圖和二維關系圖(為便于比較,圖中流閉時的流體力均取正值表示)。

 固定開度不同壓差下的振動響應
圖7 固定開度不同壓差下的振動響應

(a)位移響應(b)流體力響應

動態流體力與位移、壓差之間的關系
圖8 動態流體力與位移、壓差之間的關系

(a)流開狀態三維表示(b)流閉狀態三維表示(c)位移為參數(d)壓差為參數

由圖8可以看出:

(1)在不同的目標位移下,無論介質流向為流開或者流閉,流體力總是隨壓差增大而線性增長。

(2)在不同的進出口壓差下,無論介質流向為流開或者流閉,流體力隨著位移的增大(開度增大)而減小。

(3)在較小壓差(0.1MPa)下,流開狀態下的流體力總是小于流閉狀態下的流體力;在較大壓差(0.5~2.0MPa)下,在閥門接近全開和全關位置時,流開狀態下的流體力小于流閉狀態下的流體力,而閥門在中間位置時,流開狀態下的流體力大于流閉狀態下的流體力。

3.2.2、不同流向對閥芯位移偏移量的影響

閥芯在不同壓差下趨向目標開度的過程中,實際工作開度相對設置目標開度對應的閥芯位移之間有一定的偏移量。根據不同壓差、位移下閥芯運動仿真,可以得到不同介質流向時閥芯偏移量與壓差、位移之間的關系,如圖9所示(為便于比較,圖中流閉時的偏移量均取正值表示)。

偏移量與壓差、位移之間的關系
圖9 偏移量與壓差、位移之間的關系

(a)流開狀態(b)流閉狀態

由圖9可以看出:

(1)在較小壓差(0.1MPa)下,流開狀態下閥芯偏移量要比流閉狀態下的偏移量大;在較大壓差(0.5~2.0MPa)下,流閉狀態下的偏移量總體要比流開狀態下的偏移量大;且各種工況下的偏移量都在允許的范圍(3%~5%)內,能夠滿足定位精度。

(2)在同一個目標位置處,壓差越大,偏移量越大。

(3)在同一壓差下,目標位移越大,偏移量越小。

3.3、單座閥在變開度下的動態特性分析

閥芯在受到流體力和控制力作用時,會從初始位置向目標位置運動,并zui終保持在目標位置附近作振動。假設閥芯運動過程為20~15mm(開度變小)和5~10mm(開度變大),進出口壓差分別為0.1、0.5、1.0、2.0MPa,初始相對位移為-0.1mm,初始速度和加速度均為0,不同工況下初始流體力不同,并考慮定位器作用,得出介質流向不同時單座閥閥芯的位移與閥芯受到的流體力隨時間變化歷程,如圖10、11所示。

 開度為20~15mm時在不同壓差下的閥芯振動位移和流體力
圖10 開度為20~15mm時在不同壓差下的閥芯振動位移和流體力

(a)位移響應(b)流體力響應

由圖10、11可以看出:

(1)閥芯位移從初始位置以振蕩衰減的方式到達目標位置,并在目標位置附近作自由振動。

(2)閥芯在向下運動減小開度過程(20~15mm)中,采用流閉流向比采用流開流向,閥芯到達目標位置所用時間更短;當介質為流閉流向時,壓差越大,閥芯達到目標位置所用時間越短,而介質為流開流向時,壓差越大,閥芯達到目標位置所用時間越長。

(3)閥芯在向上運動增大開度過程(5~10mm)中,采用流開流向比采用流閉流向,閥芯到達目標位置所用時間更短;當介質為流開流向時,壓差越大,閥芯到達目標位置所用時間越短,而介質為流閉流向時,壓差越大,閥芯達到目標位置所用時間越長。

開度為5~10mm時在不同壓差下的閥芯振動位移和流體力
圖11 開度為5~10mm時在不同壓差下的閥芯振動位移和流體力

(a)位移響應(b)流體力響應

(4)閥芯在向下運動(開度減小)過程中,閥芯受到的流體力從初始位置以振蕩衰減的方式到達平衡位置,并在平衡位置附近振動,且到達平衡位置所用時間與位移到達目標位置所用時間相同;介質為流開流向時,壓差越大,閥芯受到的流體力達到平衡位置所用時間越長;介質為流閉流向時,壓差越大,閥芯受到的流體力達到平衡位置所用時間越短;總體上看,在同一壓差下,閥芯在流開流向時受到的流體力要比流閉流向時受到的流體力大。

(5)調節閥在向上運動(開度增大)過程中,閥芯受到的流體力始終處于平衡位置附近,并且在同一壓差下,閥芯在流閉流向時受到的流體力比流開流向時受到的流體力大。

4、結論
(1)建立了一個考慮定位器作用的調節閥閥芯-閥桿系統流固耦合動力學模型,并利用ANSYS軟件對固定開度和變開度情況、流體流向為流開和流閉情況進行了動態響應仿真分析。

(2)研究結果表明,閥門定位器具有減小閥芯位移振動響應和閥芯位移偏移量的作用,選擇合適的定位器參數可以保證調節閥動態性能處于良好狀態。對于固定開度情況,無論流閉型還是流開型流向,閥芯穩態位移偏移量以及流體力隨壓差增加和目標位移減小(開度減小)而增大,閥芯動態位移過渡時間隨壓差增加而縮短,而流開型流向時閥芯位移響應幅度要大于流閉型流向。對于開度減小情況,采用流閉流向相對流開流向時閥芯動態位移過渡時間更短,采用流閉流向時,壓差越大,閥芯動態位移過渡時間越短,而流開流向時則相反;對于開度增大情況,閥芯動態位移過渡時間和壓差之間的關系與開度減小情況相反。與本文相關的論文有:礦山電動插板閥

主站蜘蛛池模板: 成年人三级网站| 麻豆免费下载| 午夜福利三级理论电影| 波多野结衣在线观看一区| 国产原创剧情| 国产成人福利| 欧洲高潮三级做爰| 日韩丝袜另类精品av二区| 色窝窝无码一区二区三区| 色诱久久av| 成年人看的黄色| 欧美在线免费观看一区| 国产精品免费一区二区三区都可以| 超碰蜜桃| 黄页网站在线看| 日韩av一| 人人妻人人澡人人爽精品日本| 91新视频| 粉色视频网站| 6080午夜伦理| 色在线免费观看| 91精品大片| 欧美日韩一卡二卡三卡| 91免费看片网站| 国产又爽又黄视频| 欧洲成人午夜精品无码区久久| 亚洲精品婷婷| 2015成人永久免费视频| 人与拘一级a毛片| 成人h视频在线| 国产成人无码一二三区视频| 伊人剧场| 华人永久免费视频| 欧美日本在线| 毛片中文字幕| 激情xxx| 亚洲午夜免费| 日韩一级一级| 四虎精品| 91成人高清| 国产精品网址| 操干网| 在线观看黄色片| 色盈盈影院| av免费观看大全| 污av| 99操| 日韩av爽爽爽久久久久久| 黑人巨大精品欧美一区二区免费| 黄污视频在线观看| 久欧美| 少妇高潮一区二区三区99小说| 狠狠干中文字幕| 欧美日韩在线综合| 日韩精品黄色片| 国产精品成人一区二区三区电影毛片 | 四虎免费看黄| 尹人综合网| 日女人网站| 一级黄色电影片| 香蕉网站视频| 国产大片一区二区三区| 久久亚洲精品视频| 杨幂一区二区国产精品| 欧美猛交免费| 天天色综合小说| 久久久久久久极品内射| 97理伦| 天天鲁| 日韩福利片| 亚洲精品丝袜久久久久久| 日韩精品在线播放| 欧美日韩精品亚洲精品| 亚洲做受高潮| 国产日韩精品在线| 夜夜cao| wwwxxx黄色| 日本免费观看视频| www.久久久| 蜜桃视频在线播放| 国产二区视频| 少妇出轨精品中出一区二区| 激情久久久久久| 被各种性器调教到哭vk| 亚洲午夜精品在线观看| 免费日韩网站| 在线99热| 日本福利在线| 午夜性片| 青青视频二区| 亚洲欧美日韩综合在线| 精品国产欧美一区二区三区成人| 色狠狠干| 国产一区二区三区精品在线观看| 欧美二三区| 黄色大片免费网站| 日本成人免费电影| 欧美激情 亚洲| 日韩美女一区| 欧美成人黑人xx视频免费观看| 无码人妻丰满熟妇啪啪网站| 亚洲天堂av网站| 国产精品5| 小视频在线观看| 另类激情亚洲| 奇米成人| 四川黄色一级片| 国产传媒一区二区三区| 4438全国成人免费| 农夫色综合| 色眯眯影院| 麻豆av一区二区三区久久| 91成人免费看片| av免费中文字幕| 中国少妇精品| 亚洲网站免费观看| 黄色在线观看国产| 精品久久久久久久久久久久久久久久久久 | 少妇高潮一区二区三区99小说| 依依综合网| 中国丰满熟妇xxxx性| 亚洲精品aⅴ| 视频在线观看一区| 插插插综合视频| 天天舔天天爱| 91玉足脚交白嫩脚丫| 久久精品激情| 一级黄色大片免费看| 青娱乐最新官网| 免费成人深夜夜行p站| 亚洲综合久久久| 全黄一级片| www.九九精品| 夜夜天天| yjizz国产| 伊人久久久久久久久久| 影音先锋日韩av| 女人色极品影院| 深夜福利成人| 欧美黑白大战| 玖玖久久| 大奶子在线| 天天操天天射天天射| 亚洲欧洲中文日韩久久av乱码| 亚洲最大的网站| 久久不射电影网| 牛牛免费视频| 五月婷婷综合在线| 国产最新在线视频| 日韩欧美中| 夜色综合网| 日韩毛片免费观看| xxxxx麻豆| 成人欧美一区二区三区黑人免费| 国产一区二区影院| 成人三级图片| 黄色18网站| av黄色在线观看| 久久精品这里| 夜av| 国产精品剧情av| 国产欧美日韩综合精品一区| 天天综合天天综合| 视频在线观看你懂的| 96视频在线免费观看| 国产在线欧美| 在线观看成人网| 久久国产三级| 黄网址在线观看| 在线观看av不卡| 日韩欧美亚洲精品高清国产| 国产一区在线视频| 在线免费观看黄色av| 日韩在线视屏| 日本japanese丰满白浆| 日韩欧美成人免费电影| 国语free性xxxxx按摩| 华丽的外出在线观看| 亚洲第一精品网站| 亚洲国产视频在线| 欧美日韩精品电影| 国产日产久久欧美清爽| 无码无套少妇毛多18pxxxx| 美女被男人插| 日躁夜躁狠狠躁2001| 天天综合天天做天天综合| 亚天堂| 中文字幕制服丝袜在线| 色天堂在线视频| 四虎在线免费视频| 丝袜中文字幕| 91久久国产露脸精品国产护士| avwww在线| 日批免费看| 淫视频在线观看| 欧美日韩六区| 凸凹人妻人人澡人人添| 91免费网站视频| 丰满少妇大力进入av亚洲| 怡春院欧美| 极品在线观看| 日韩欧美大陆| 国产无套精品| 色噜噜综合| 中文字幕12页| 不卡一区二区在线观看| 亚洲精品一区在线观看| 五月婷婷久久综合| 欧美精品在线一区| 久久高清无码视频| 特级大毛片| 一级肉体全黄裸片| 国产高清sp| 久爱视频在线观看| 欧美a∨| 久久国产99| 琪琪色综合网| 中文字幕区| 精品国产乱码久久久久久88av | 欧美一区二区三区四区五区六区| 欧洲久久精品| 好了av在线| 天天干天天操天天透| av网站导航| 男女av片| 国产精品自拍网| 好看的国产精彩视频| 亚洲综合欧美日韩| 日韩一级精品| 成人深夜网站| japan捆绑紧缚bdsm在线视频| 精产国品一二三产品蜜桃| 乱lun合集男男高h| 欧美毛片网站| 韩国19禁床戏大尺度片| 在线观看亚洲成人| 风韵多水的老熟妇| 夫妻毛片| 日韩av直播| 91免费视频播放| 男人插入女人下面视频| 欧美一区二三区| 真实人妻互换毛片视频| 美女脱衣服一干二净| 在线视频观看国产| 免费日本在线| 日韩在线中文| 免费涩涩网站| 最新的av网站| 欧美性一区| 调教女m荡骚贱淫故事| 三级黄色免费| 国产思思| 亚洲美女自拍| 亚洲人 女学生 打屁股 得到| 免费九九视频| 先锋影音av在线资源| 福利影院在线| 日日噜噜夜夜狠狠久久波多野| 国产一级二级三级精品| 极品超粉嫩尤物69xx| 五月婷婷六月色| 国产精品乱码久久久久久久久| 亚洲区一| 国产精品二区视频| 日本少妇喷水视频| 夜夜爽网址| 色悠悠av| 男女裸体无遮挡做爰| www.亚洲视频.com| 欧美日韩亚洲国产一区| jizz 亚洲大全| 好爽又高潮了毛片| 奇米影视四色7777| 麻豆小视频| 国产欧美日韩三级| 欧美精品黑人猛交高潮| 亚洲精品视频二区| www射我里面在线观看| 久久国产精品一区二区三区| 日韩免费黄| 久久橹| 欧美精品一级在线观看| 中文字幕在线观看91| 狠狠爱综合网| 老版水浒传83版免费播放| 一区二区三区精品在线| 超碰在线3| 国内精品国产成人国产三级| 国产视频高清在线观看| 小蝌蚪视频色| 久久国产精品久久久| 欧美aaaaa| 久久综合88| 自拍偷拍视频亚洲| 色哟哟精品一区| 国产毛片精品| 蜜桃色av| 99自拍| 男人天堂视频在线观看| 少妇在线视频| aaa淫| sm在线观看| 四虎影视成人精品| 日本一区二区在线播放| 亚洲色图av在线| 91网站观看| 白石茉莉奈中文字幕在| 国产激情自拍| 婷婷激情社区| 看片网站在线观看| 国产成人av一区| 新香蕉视频| 制服.丝袜.亚洲.另类.中文| 啪啪拍网站| 成人快手免费看片| 茄子爱啪啪| 欧美高清hd18日本| 亚洲加勒比久久88色综合| 国产精品视频伦| aa一级片| 欧美精品一区在线发布| 一区二区三区伦理片| 91福利影院| 欧美黄色a级片| 宇都宫紫苑在线播放| 丰满少妇大力进入av亚洲| 黄色福利网| 中文字幕在线看| 97在线免费| 国产黄色录相| 亚洲综合色在线观看| 国产又黄又爽无遮挡在线观看| 免费三级a| 40一50一60老女人毛片| jzzjzzjzz亚洲成熟少妇| 在线婷婷| 欧美激情免费视频| 麻豆爱爱| 新疆毛片| 欧美成人免费高清视频| 91一起草| 国产91大片| 亚洲天堂avav| 日本高清网站| 精品久久一| 在线观看视频h| 暖暖视频日本| 亚洲性图视频| 精品免费国产| 做视频免费观看网站| 超碰在97| 日本高清www午色夜在线视频 | 国产高清免费在线播放| 免费看成人毛片| 淫视频网站| 国产精品乱码一区二区 | 久久久久久久久久国产精品| 兔费看少妇性l交大片免费 | 亚洲第一自拍| 日本v视频| 能看的黄色网址| 韩国av永久免费| 欧美区二区三区| 另一种灿烂生活| 免费激情视频网站| 国产视频一区二区三| 好看的国产精品| 久久性片| 男人在线视频| 97香蕉超级碰碰久久免费直播间| 丝袜 亚洲 另类 欧美 重口| 一级绝黄| 色爱av综合网| 不卡日本视频| 99啦porny丨首页入口| 狠狠涩| 亚洲一区二区色图| 91视频精品| 玖玖爱免费视频| 欧美成人播放| 夜夜免费视频| 午夜激情网| 亚洲伦理自拍| 午夜激情免费| 久久久久亚洲av无码专区喷水| 午夜看片在线观看| 黄色片视频免费| 久久男人的天堂| 日本黄色免费视频| 亚洲免费黄色网| 久久怡红院| 自拍愉拍| 亚洲精品视| 日韩Av无码精品| 短篇山村男同肉耽h| 精品99999| 红桃成人网| 日韩欧美123| 老司机午夜福利视频| 日本不卡1| 高清毛片aaaaaaaaa片| 麻豆性生活| www.97视频| 无遮挡免费网站| 成人国产三级| 67194国产| 香蕉a| 国产精品成人av电影| 亚洲你懂的| 亚洲欧美日本在线观看| 日本va视频| 中老年妇女性色视频| 无码人妻aⅴ一区二区三区| 超碰女优| 激情小说欧美色图| 亚洲激情欧美另类| 成人免费区| 亚洲国产第一页| 一级黄色欧美| 黄频网站在线观看| 天天舔夜夜| 黄色片视频免费| 亚洲精品视频导航| 国产中文字幕一区二区三区| 午夜一区| 国产精品一二三区视频出来一| xxxx18日本| 天堂免费在线视频| 高清不卡av| 中文字幕色| 在线有码视频| 日韩二级片| 久久久精品久久| 16一17女人毛片| 夜夜夜噜噜噜| 波多野结衣视频播放| 久久精品h| 男人天堂亚洲| 国产四虎在线| 91操操| 欧美激情一区二区在线| heyzo精品少妇| 国产精品视频网站| 天天做天天爱天天操| 91精品国产欧美一区二区成人| 国产乡下妇女三片| 精品人妻一区二区免费视频| 久久精品av| 大陆明星乱淫(高h)小说| 女性向片在线观看| 波多野结衣精品视频| 欧美老司机| 日韩 欧美| 国产九色sp调教91| 美女超碰在线| 蜜臀久久精品久久久用户群体| av导航福利| 男人av的天堂| 午夜手机福利| 成人做爰100| 亚洲五月婷婷| 欧美做受69| 亚洲一品道| av在线影视| 99精品欧美| 秋霞av在线| 色在线高清| 成人久久久电影| 久久综合99| 麻豆tv在线观看| 91蝌蚪视频在线| 久热国产精品| 国产曰肥老太婆无遮挡| 欧美丰满美乳xxx高潮www| 狠狠躁夜夜躁人人爽视频| 人人澡超碰碰| 欧美激情16p| 国产精品久久久午夜夜伦鲁鲁| 日本熟妇毛耸耸xxxxxx| 日韩亚洲视频| 亚洲最大av网站| 狂野欧美性猛交blacked| 高清av在线| 综合色自拍| 国产综合精品在线| 久热这里只有| 国产成人精品a视频一区| 国产免国产免费| 三浦理惠子av在线播放| 国产黑丝在线看| 黄污视频在线播放| 国产精品久久久久久免费| 99色综合| 18久久久| 欧洲一区二区三区在线| 亚洲国产精品久久久久久6q| 一区欧美| 91手机在线视频| 少妇淫片免费看大片动漫版app|