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

    航空發(fā)動(dòng)機(jī)管路流固耦合振動(dòng)的固有頻率分析

    2017-11-10 09:24:08李占營(yíng)王建軍邱明星
    航空發(fā)動(dòng)機(jī) 2017年1期
    關(guān)鍵詞:固支輸液管外徑

    李占營(yíng),王建軍,邱明星

    (1.北京航空航天大學(xué)能源與動(dòng)力工程學(xué)院,北京100191;2.中國(guó)航發(fā)沈陽(yáng)發(fā)動(dòng)機(jī)研究所,沈陽(yáng)110015)

    航空發(fā)動(dòng)機(jī)管路流固耦合振動(dòng)的固有頻率分析

    李占營(yíng)1,王建軍1,邱明星2

    (1.北京航空航天大學(xué)能源與動(dòng)力工程學(xué)院,北京100191;2.中國(guó)航發(fā)沈陽(yáng)發(fā)動(dòng)機(jī)研究所,沈陽(yáng)110015)

    為研究流體哥氏力和管路參數(shù)等因素對(duì)航空發(fā)動(dòng)機(jī)管路固有振動(dòng)頻率的影響規(guī)律,采用G alerkin方法建立了管路流固耦合數(shù)學(xué)模型,并通過(guò)復(fù)特征值分析得到了系統(tǒng)的固有頻率。通過(guò)將采用G alerkin方法的計(jì)算結(jié)果與試驗(yàn)測(cè)試數(shù)據(jù)進(jìn)行比較,驗(yàn)證了G alerkin方法的正確性;給出了試驗(yàn)管路的臨界流速,并研究了流體哥氏力和管路截面尺寸對(duì)系統(tǒng)固有頻率的影響。結(jié)果表明:試驗(yàn)燃油管路實(shí)際流速遠(yuǎn)小于發(fā)生屈曲失穩(wěn)的臨界流速;哥氏力對(duì)不銹鋼和鈦合金2種管材燃油管路固有頻率的影響很小;相同壁厚管路,外徑越小,流固耦合對(duì)固有頻率的影響越大。

    輸液管路;流固耦合;G alerkin方法;固有頻率;哥氏力;航空發(fā)動(dòng)機(jī)

    0 引言

    航空發(fā)動(dòng)機(jī)管路系統(tǒng)用于發(fā)動(dòng)機(jī)部件之間及其與飛機(jī)間的流體輸送。管路系統(tǒng)振動(dòng)會(huì)引發(fā)管路故障,影響到管路系統(tǒng)甚至發(fā)動(dòng)機(jī)安全可靠地運(yùn)行。該振動(dòng)屬于典型的輸液管路振動(dòng),而作為振動(dòng)的前沿課題,國(guó)內(nèi)外很多學(xué)者[1-3]采用不同方法對(duì)輸液管路的動(dòng)力學(xué)行為進(jìn)行了深入研究。Holmes[4]利用Lyapunov直接法研究了兩端固支輸液直管在自激勵(lì)作用下全局動(dòng)態(tài)行為,表明該類管路不會(huì)發(fā)生顫振失穩(wěn);張智勇[5]采用傳遞矩陣法研究了流固耦合對(duì)輸液管路固有頻率的影響;金基鐸[6]采用平均法研究了兩端固支輸液直管在脈動(dòng)內(nèi)流作用下的參數(shù)共振問(wèn)題;Tornabene[7]采用廣義微分求積法研究了輸液直管的臨界流速。

    在固有振動(dòng)特性研究方面,楊曉東等[8]以復(fù)模態(tài)分析方法為基準(zhǔn),研究了Galerkin截?cái)喾椒▽?duì)兩端鉸支輸液管路固有頻率分析的精確性;齊歡歡等[9]采用Galerkin離散和復(fù)模態(tài)2種方法計(jì)算了不同流速下懸臂輸液管的固有頻率;楊超等[10]采用特征線法研究了管道結(jié)構(gòu)阻尼、管材泊松比、管道壁厚以及管道材料等參數(shù)對(duì)蓄水池-管道-閥門系統(tǒng)聲彈耦合振動(dòng)特性的影響;Lin等[11]采用復(fù)模態(tài)法研究了基礎(chǔ)激勵(lì)下航空發(fā)動(dòng)機(jī)液壓直管的流固耦合振動(dòng),得出了流體壓力、流速和軸向力的增加,會(huì)引起管路固有頻率降低的結(jié)論;Huang等[12]基于Ferrari法,推導(dǎo)出兩端支承輸液管路固有頻率方程,并研究了兩端支承輸液直管的固有頻率和穩(wěn)定性;Wang等[13]采用有限元法基于ADINA建立輸液管路3維流固耦合動(dòng)力學(xué)模型,計(jì)算了直管和曲管的固有頻率,并與解析法進(jìn)行了對(duì)比。

    從現(xiàn)有的研究工作可見,輸液管路固有振動(dòng)特性的研究主要集中在研究管路的臨界流速,以及流體壓力、流速等因素對(duì)固有頻率的影響,還未發(fā)現(xiàn)由流體引起的哥氏力對(duì)固有頻率影響的研究。

    本文采用Galerkin方法建立輸液管路系統(tǒng)的流固耦合動(dòng)力學(xué)模型,計(jì)算并試驗(yàn)驗(yàn)證燃油管路流固耦合固有頻率,研究其臨界流速及哥氏力和管路外徑對(duì)輸液管路固有頻率的影響規(guī)律。

    1 管路系統(tǒng)流固耦合模型

    1.1 研究對(duì)象

    研究對(duì)象為兩端固支的輸液管路系統(tǒng),如圖1所示。模型兩端固定約束,流體從左端流入、右端流出。管路的長(zhǎng)度和抗彎剛度分別為L(zhǎng)和EI,其中E為彈性模量、I為管路截面慣性矩。管內(nèi)流體的平均流速為U。

    1.2 運(yùn)動(dòng)方程及其離散特性

    假定管內(nèi)為無(wú)黏不可壓縮的穩(wěn)定流動(dòng)流體,忽略重力、阻尼力、管路軸向載荷和管內(nèi)流體壓力,并對(duì)管路采用小變形假設(shè),由牛頓力學(xué)原理,得到管路微段的運(yùn)動(dòng)微分方程[14]

    式中:w、x、t、mp、mf分別為管路橫向位移、管路橫截面所處的位置、時(shí)間、管路線密度、管內(nèi)流體線密度,各項(xiàng)的物理意義分別為彈性恢復(fù)力、離心力、哥氏力和慣性力。

    兩端固支邊界條件為

    引入

    得到無(wú)量綱方程

    其中

    則無(wú)量綱邊界條件為

    考慮定常流,采用4階Galerkin方法對(duì)式(4)進(jìn)行離散,即令

    式中:φ(rξ)、q(rτ)分別為相同邊界條件梁的無(wú)量綱振型函數(shù)、離散系統(tǒng)廣義坐標(biāo)。

    將式(7)代入式(4),左乘 φ(rξ),并在[0,1]區(qū)間內(nèi)積分,同時(shí)利用正弦函數(shù)的正交性(即δsr是 Kronecker’s delta函數(shù))及,得到離散后的系統(tǒng)微分方程

    式中:s=1,2,…,4;λr為兩端固支梁無(wú)量綱特征值;bsr=

    式(8)可以寫成如下矩陣形式

    式(9)可以寫成標(biāo)準(zhǔn)形式

    由于哥氏力矩陣B的存在,使得系統(tǒng)阻尼矩陣[C]為反對(duì)稱矩陣,因此只能采用復(fù)特征值方法求解式(10)的特征值。首先引入狀態(tài)方程,將式(10)特征值問(wèn)題轉(zhuǎn)換為標(biāo)準(zhǔn)特征值問(wèn)題,再進(jìn)行求解。

    2 輸液管路的固有頻率試驗(yàn)

    對(duì)兩端固支燃油管路進(jìn)行固有頻率測(cè)試,將試驗(yàn)測(cè)試結(jié)果與采用Galerkin方法的計(jì)算結(jié)果進(jìn)行對(duì)比。

    燃油鈦合金直管的固有頻率測(cè)量采用錘擊法[15]。試驗(yàn)設(shè)備包括加速度計(jì)、力錘、INV303系列智能信號(hào)采集處理分析儀、計(jì)算機(jī)及DASP大容量數(shù)據(jù)采集與信號(hào)處理分析軟件。錘擊法對(duì)燃油管路進(jìn)行振動(dòng)測(cè)試如圖2所示。

    試驗(yàn)件管長(zhǎng)1000 mm,厚1 mm,外徑為14 mm。采用臺(tái)鉗夾緊管路兩端,近似實(shí)現(xiàn)固支約束。約束位置點(diǎn)距離導(dǎo)管兩端100 mm。

    數(shù)值計(jì)算采用Galerkin方法進(jìn)行建模,模型中管路材料的彈性模量、泊松比和密度分別取96 GPa、0.39 和 4470 kg/m3,燃油密度取其20℃時(shí)的密度值802.4 kg/m3。兩端固支空管的前2階固有頻率見表1。

    從表中可見,燃油管路的前2階固有頻率數(shù)值計(jì)算結(jié)果與試驗(yàn)結(jié)果比較吻合。相對(duì)誤差在3%以內(nèi),因此Galerkin方法的計(jì)算精度得到了驗(yàn)證,計(jì)算值與試驗(yàn)值的差異原因在于臺(tái)鉗夾緊不能實(shí)現(xiàn)完全固支邊界條件。

    不同流速下兩端固支燃油管路的第1階固有頻率見表2。

    表1 兩端固支空管的固有頻率

    表2 不同流速兩端固支燃油管路的固有頻率

    從表中可見,計(jì)算結(jié)果與試驗(yàn)結(jié)果比較吻合,隨著流速的增加第1階固有頻率均減??;燃油流速在2.62 m/s以下時(shí),管路的第1階固有頻率減小的幅值很小,這是因?yàn)榇藭r(shí)的無(wú)量綱流速僅為0.0692;計(jì)算結(jié)果與試驗(yàn)結(jié)果的相對(duì)誤差為4.26%~5.59%,差異的原因在于試驗(yàn)不能實(shí)現(xiàn)計(jì)算中假設(shè)的完全固支邊界條件。

    3 臨界流速及固有頻率分析

    3.1 臨界流速分析

    采用4階Galerkin離散方法建立試驗(yàn)管路的數(shù)學(xué)模型,計(jì)算其特征值問(wèn)題,得到試驗(yàn)管路固有頻率隨流速變化的規(guī)律,如圖3所示。

    從圖中可見,隨著流速的增加,管路的剛度減小,固有頻率降低,這與現(xiàn)有文獻(xiàn)的研究結(jié)論一致[14];試驗(yàn)管路的臨界流速為150 m/s,實(shí)際發(fā)動(dòng)機(jī)燃油管路流速遠(yuǎn)小于這一流速,因此不會(huì)發(fā)生屈曲失穩(wěn)。

    3.2 哥氏力對(duì)固有頻率的影響

    之后,通過(guò)計(jì)算研究哥氏力對(duì)管路系統(tǒng)固有頻率的影響規(guī)律。燃油流速在100m/s之內(nèi),是否考慮哥氏力時(shí)管路系統(tǒng)的固有頻率及相對(duì)誤差分別如圖4、5所示。

    從圖5中可見,燃油流速在100 m/s以內(nèi)時(shí),忽略哥氏力引入誤差隨著流速的增加而增加,但最大誤差不超過(guò)1.5E-3。這說(shuō)明哥氏力對(duì)鈦合金燃油管路前2階固有頻率的影響幾乎可以忽略。這是因?yàn)橛扇加土魉僖鸬母缡狭?duì)阻尼矩陣的改變很小,從而使管路系統(tǒng)的阻尼固有頻率和無(wú)阻尼固有頻率幾乎相同。

    相同截面尺寸的不銹鋼管路是否考慮哥氏力時(shí)的固有頻率及相對(duì)誤差分別如圖6、7所示。

    從圖6、7中可見,與鈦合金燃油管路的規(guī)律相同,即隨著流速增加,管路剛度減小,固有頻率降低;忽略哥氏力得到的固有頻率最大誤差不超過(guò)5E-4。另外,與鈦合金管路相比,不銹鋼質(zhì)量增加對(duì)固有頻率的影響小于截面剛度增加對(duì)固有頻率的影響,因此,不銹鋼管路的固有頻率大于相同截面尺寸鈦合金管路的。

    3.3 管路外徑對(duì)固有頻率的影響

    壁厚為1 mm,外徑分別為20、14和8 mm的兩端固支鈦合金燃油管路第1階固有頻率隨流速的變化規(guī)律如圖8所示。在不同外徑下,流速對(duì)鈦合金管路固有頻率影響的相對(duì)變化情況如圖9所示。

    從圖8、9中可見,相同壁厚管路在相同流速下:

    (1)隨著外徑的增大,鈦合金燃油管路的固有頻率升高。這是因?yàn)橥ㄟ^(guò)增大管路外徑,系統(tǒng)剛度增大,進(jìn)而固有頻率提高。

    (2)隨著外徑的增大,相對(duì)于流速為零時(shí)的固有頻率變化率減小,這說(shuō)明流固耦合對(duì)于直徑小的管路影響更大。例如,在流速為100 m/s時(shí),外徑為8 mm的管路固有頻率是流速為零時(shí)的55%。

    壁厚為1 mm,外徑分別為20、14和8 mm的兩端固支不銹鋼燃油管路第1階固有頻率隨流速的變化規(guī)律如圖10所示。在不同外徑下,流速對(duì)不銹鋼管路固有頻率影響的相對(duì)變化情況如圖11所示。

    從圖10、11中可見,與鈦合金管路相同,相同壁厚管路在相同流速下:

    (1)隨著外徑的增大,不銹鋼燃油管路的固有頻率升高。

    (2)隨著外徑的增大,相對(duì)于流速為零時(shí)的固有頻率變化率減小,這說(shuō)明流固耦合對(duì)于直徑小的管路影響更大。例如,在流速為100 m/s時(shí),外徑為8 mm的不銹鋼管路固有頻率是流速為零時(shí)的82%。

    另外,通過(guò)對(duì)比可以發(fā)現(xiàn),相同截面特性的不銹鋼和鈦合金管路,流固耦合對(duì)鈦合金管路的影響更大。

    4 結(jié)論

    (1)通過(guò)Galerkin方法截取適當(dāng)?shù)哪B(tài),可建立1個(gè)能反映輸液管路流固耦合動(dòng)力學(xué)特性的低階模型,用于計(jì)算其固有振動(dòng)特性。

    (2)航空發(fā)動(dòng)機(jī)鈦合金燃油管路的實(shí)際流動(dòng)速度遠(yuǎn)小于其失穩(wěn)臨界流速,因此其不會(huì)發(fā)生屈曲失穩(wěn)。

    (3)進(jìn)行航空發(fā)動(dòng)機(jī)管路系統(tǒng)振動(dòng)特性分析時(shí),可忽略哥氏力對(duì)不銹鋼和鈦合金2種管材固有頻率的影響,即其阻尼固有頻率和無(wú)阻尼固有頻率相差很小,可以直接使用對(duì)稱矩陣特征值求解器求其固有振動(dòng)特性。

    (4)相同壁厚的燃油管路,流固耦合對(duì)外徑小的管路固有頻率影響更大;相同截面尺寸的燃油管路,流固耦合對(duì)鈦合金管路固有頻率的影響大于對(duì)不銹鋼管路固有頻率的影響。

    [1]Paidoussis M P,Li G X.Pipes conveying fluid:a model dynamical problem[J].Journal of Fluids and Structures,1993,7(2):137-204.

    [2]任建亭,姜節(jié)勝.輸流管道系統(tǒng)振動(dòng)研究進(jìn)展 [J].力學(xué)進(jìn)展,2003,33(3):313-324.REN Jianting,JIANG Jiesheng.Advances and trends on vibration of pipes conveying fluid [J].Advances in Mechanics,2003,33(3):313-324.(in Chinese)

    [3]Ibrahim R A.Overview of mechanics of pipes conveying fluids—part I:fundamental studies[J].Journal of Pressure Vessel Technology,2010,132(3):034001-034032.

    [4]Holmes P J.Pipes supported at both ends cannot flutter[J].Journal of Applied Mechanics,1978,45(3):619-622.

    [5]張智勇,沈榮瀛,王強(qiáng).充液管道系統(tǒng)的模態(tài)分析 [J].固體力學(xué)學(xué)報(bào),2001,22(2):143-149.ZHANG Zhiyong,SHEN Rongying,WANG Qiang.The modal analysis of the liquid-filled pipe system[J].ACTA Mechanica Solida Sinica,2001,22(2):143-149.(in Chinese)

    [6]金基鐸,宋志勇,楊曉東.兩端固定輸流管道的穩(wěn)定性和參數(shù)共振[J].振動(dòng)工程學(xué)報(bào),2004,17(2):190-195.JIN Jiduo,SONG Zhiyong,YANG Xiaodong.Stability and parametric resonances of a clamped-clamped pipe conveying fluid[J].Journal of Vibration Engineer,2004,17(2):190-195.(in Chinese)

    [7]Tornabene F,Marzani A,Viola E,et al.Critical flow speeds of pipes conveying fluid by the generalized differential quadrature method[J].Journal of Theoretical and Applied Mechanics 2010,3(3):121-138.

    [8]楊曉東,金基鐸.輸流管道流-固耦合振動(dòng)的固有頻率分析 [J].振動(dòng)與沖擊,2008,27(3):80-81.YANG Xiaodong,JIN Jiduo.Comparison of Galerkin method and complex mode method in natural frequency analysis of tube conveying fluid[J].Journal of Vibration and Shock,2008,27(3):80-81.(in Chinese)

    [9]齊歡歡,徐鑒.Galerkin模態(tài)截?cái)鄬?duì)計(jì)算懸臂輸液管道固有頻率的影響[J].振動(dòng)與沖擊,2011,30(1):148-151.QI Huanhuan,XU jian.Effect of Galerkin modal truction on natural frequency analysis of a cantilevered pipe conveying fluid[J].Journal of Vibration and Shock,2011,30(1):148-151.(in Chinese)

    [10]楊超,范士娟.管材參數(shù)對(duì)輸液管流固耦合振動(dòng)的影響[J].振動(dòng)與沖擊,2011,30(7):210-213.YANG Chao,FAN shijuan.Influence of pipe parameters on fluidstructure coupled vibration of a fluid-conveying pipe[J].Journal of Vibration and Shock,2011,30(7):210-213.(in Chinese)

    [11]Lin J Z,Qin L,Zhou E T,et al.Fluid-structure interaction vibration of hydraulic pipe system[J].Advanced Engineering Forum,2011(2):822-827.

    [12]Huang Y M,Ge S,Wu W,et al.A direct method of natural frequency analysis on pipeline conveying fluid with both ends supported[J].Nuclear Engineering&Design,2012,253(12):12-22.

    [13]Wang L,Gan J,Ni Q.Natural frequency analysis of fluid-conveying pipes in the ADINA system[J].Journal of Physics Conference,2013,448(1):12014-12020(7).

    [14]Paidoussis M P.Fluid-structure interactions:slender structures and axial flow[M].London:Academic Press Limited,1998:89-90.

    [15]Ewins D J.Modal testing:theory and practice[M].Letchworth:Research Studies Press,1986:156-160.

    Analysis for Natural Frequencies of Pipe Conveying Fluid Considering Fluid-Structure Interaction

    LI Zhan-ying1,WANG Jian-jun1,QIU Ming-xing2
    (1.School of Energy and Power Engineering,Beihang University,Beijing 100191,China;2.AECC Shenyang Engine Research Institute,Shenyang 110015,China)

    In order to analyze the effect of Coriolis force and pipe parameters on the natural frequencies,the fluid-structure coupling model of pipe conveying fluid was established using the Galerkin method,and the natural frequencies of the pipe were obtained by the complex-eigenvalue method.Comparing with the experimental data,the analytical results were verified,which indicated that the analytical model is rational.The critical velocity of the test pipe conveying fuel was obtained,and the influence of Coriolis force and the section characteristics to the natural frequencies were investigated.It is found that the fuel velocity of the test pipe is far below the critical velocity.The Coriolis force has little effect on the natural frequencies of both the stainless steel and the titanium alloy pipe.For the pipe in the same thickness,the fluid-structure interaction has much effect on the natural frequency as the out diameter of the pipe decrease.

    pipe conveying fluid;fluid-structure interaction;Galerkin method;natural frequency;Coriolis force;aeroengine

    V 233

    A

    10.13477/j.cnki.aeroengine.2017.01.012

    2016-05-18

    李占營(yíng)(1983),男,在讀博士研究生,研究方向?yàn)楹娇瞻l(fā)動(dòng)機(jī)強(qiáng)度、振動(dòng)和可靠性;E-mail:lizhanying@buaa.edu.cn。

    李占營(yíng),王建軍,邱明星.航空發(fā)動(dòng)機(jī)管路流固耦合振動(dòng)的固有頻率分析[J].航空發(fā)動(dòng)機(jī),2017,43(1):66-70.LIZhanying,WANGJianjun,QIU Mingxing.Analysisfornaturalfrequenciesofpipeconveyingfluidconsideringfluid-structureinteraction[J].Aeroengine,2017,43(1):66-70.

    (編輯:栗樞)

    猜你喜歡
    固支輸液管外徑
    固支結(jié)構(gòu)對(duì)壓阻式壓力傳感器的輸出影響研究
    ◆敷設(shè)線纜用材料
    敷設(shè)線纜用材料
    三角形邊長(zhǎng)與內(nèi)外徑之間的三個(gè)不等式
    一次性防脫落輸液管針頭的設(shè)計(jì)
    基于復(fù)合材料的山體滑坡固支模型及配套傳感器布置方法申請(qǐng)?zhí)? 202010514582.0
    傳感器世界(2021年2期)2021-03-27 19:36:46
    一角點(diǎn)支撐對(duì)面兩邊固支正交各向異性矩形薄板彎曲問(wèn)題的辛疊加解
    一種便于連接的塑料管
    約束分布對(duì)混合邊界薄板固有頻率影響的研究★
    山西建筑(2018年6期)2018-03-22 07:14:15
    輸液管插入法治療牛非泡沫性瘤胃脹氣
    国产av麻豆久久久久久久| 爱豆传媒免费全集在线观看| 岛国在线免费视频观看| 婷婷色av中文字幕| 欧美日韩国产亚洲二区| 欧美又色又爽又黄视频| 男女边吃奶边做爰视频| 亚洲国产精品成人久久小说 | 寂寞人妻少妇视频99o| 国产精品av视频在线免费观看| 色哟哟·www| 国产精品久久视频播放| 在线播放国产精品三级| 99视频精品全部免费 在线| 亚洲乱码一区二区免费版| 国产一区二区在线av高清观看| 一边摸一边抽搐一进一小说| 亚洲美女搞黄在线观看| 在线观看66精品国产| 黄片wwwwww| 国产精品女同一区二区软件| a级毛色黄片| 国产成人精品一,二区 | 韩国av在线不卡| 丝袜喷水一区| 国内精品美女久久久久久| 亚洲av男天堂| 毛片一级片免费看久久久久| 高清在线视频一区二区三区 | 国产成人午夜福利电影在线观看| 欧美一区二区亚洲| 精品午夜福利在线看| 亚洲乱码一区二区免费版| 久久99热这里只有精品18| 在线观看一区二区三区| 国产极品天堂在线| 久久国产乱子免费精品| 岛国在线免费视频观看| 国产在视频线在精品| 亚洲国产精品合色在线| 草草在线视频免费看| 69av精品久久久久久| 男人狂女人下面高潮的视频| 日本熟妇午夜| 少妇丰满av| 99九九线精品视频在线观看视频| 午夜激情福利司机影院| 国产伦理片在线播放av一区 | 听说在线观看完整版免费高清| 亚洲欧美中文字幕日韩二区| 亚洲精华国产精华液的使用体验 | 国产v大片淫在线免费观看| 国产午夜精品论理片| 蜜臀久久99精品久久宅男| 久久久精品94久久精品| 久久精品国产亚洲网站| 不卡视频在线观看欧美| 搡女人真爽免费视频火全软件| 两个人视频免费观看高清| 国产一级毛片在线| 少妇高潮的动态图| 国产伦精品一区二区三区视频9| 嫩草影院新地址| 在线免费观看不下载黄p国产| 99久久精品一区二区三区| 97超视频在线观看视频| 亚洲精品国产av成人精品| 99久久成人亚洲精品观看| 国产黄片美女视频| 日韩三级伦理在线观看| 51国产日韩欧美| 成人午夜精彩视频在线观看| 亚洲自拍偷在线| 偷拍熟女少妇极品色| 色视频www国产| 在线播放无遮挡| av天堂在线播放| 在线观看美女被高潮喷水网站| 夜夜夜夜夜久久久久| av在线亚洲专区| 黄片无遮挡物在线观看| 校园人妻丝袜中文字幕| 国产乱人偷精品视频| 精品不卡国产一区二区三区| 99在线人妻在线中文字幕| av女优亚洲男人天堂| 国产精品电影一区二区三区| 久久人人精品亚洲av| 青青草视频在线视频观看| 免费看av在线观看网站| av天堂中文字幕网| 国产国拍精品亚洲av在线观看| 精品欧美国产一区二区三| 免费观看精品视频网站| 自拍偷自拍亚洲精品老妇| 综合色丁香网| 老司机福利观看| 久久久久久久久中文| 99视频精品全部免费 在线| 日韩精品青青久久久久久| 少妇被粗大猛烈的视频| 国产精品精品国产色婷婷| 日韩av在线大香蕉| 亚洲aⅴ乱码一区二区在线播放| 欧美成人a在线观看| 国产真实乱freesex| 国产大屁股一区二区在线视频| 国产熟女欧美一区二区| 久久久久久久午夜电影| 欧美成人精品欧美一级黄| 人人妻人人澡人人爽人人夜夜 | 成年版毛片免费区| 一边亲一边摸免费视频| 性欧美人与动物交配| 国内精品一区二区在线观看| 亚州av有码| 国产成人aa在线观看| 97超碰精品成人国产| 亚洲乱码一区二区免费版| 两个人视频免费观看高清| 尤物成人国产欧美一区二区三区| 欧美一区二区精品小视频在线| 国国产精品蜜臀av免费| 一卡2卡三卡四卡精品乱码亚洲| 国产白丝娇喘喷水9色精品| 丰满的人妻完整版| 禁无遮挡网站| 国产精品不卡视频一区二区| 亚洲欧美成人综合另类久久久 | 国产又黄又爽又无遮挡在线| 日本免费一区二区三区高清不卡| а√天堂www在线а√下载| 内地一区二区视频在线| 国产高潮美女av| 免费看美女性在线毛片视频| 日韩欧美国产在线观看| 神马国产精品三级电影在线观看| 日韩av在线大香蕉| 亚洲在线自拍视频| 久久久午夜欧美精品| 校园人妻丝袜中文字幕| 亚洲三级黄色毛片| 村上凉子中文字幕在线| 男人的好看免费观看在线视频| 日韩视频在线欧美| 18禁黄网站禁片免费观看直播| 女人十人毛片免费观看3o分钟| 色尼玛亚洲综合影院| 永久网站在线| 久久6这里有精品| 级片在线观看| 日本五十路高清| 国产成人a∨麻豆精品| 国产 一区 欧美 日韩| 22中文网久久字幕| 久久韩国三级中文字幕| 亚洲av中文字字幕乱码综合| 丰满的人妻完整版| 久久精品久久久久久久性| 国产精品.久久久| 丰满人妻一区二区三区视频av| 精品一区二区三区视频在线| 成人高潮视频无遮挡免费网站| 九九久久精品国产亚洲av麻豆| 午夜福利视频1000在线观看| 麻豆一二三区av精品| 亚洲国产精品sss在线观看| 一边摸一边抽搐一进一小说| 亚洲,欧美,日韩| 久久这里只有精品中国| 99热6这里只有精品| 亚洲欧美日韩东京热| av免费在线看不卡| 国产老妇伦熟女老妇高清| 高清毛片免费看| 成人特级黄色片久久久久久久| 亚洲成人精品中文字幕电影| 久久午夜福利片| 一级毛片aaaaaa免费看小| 村上凉子中文字幕在线| 久久久久网色| 看十八女毛片水多多多| 国产精品爽爽va在线观看网站| 欧美成人免费av一区二区三区| 在线观看66精品国产| 久久精品综合一区二区三区| 精品久久久久久久久久久久久| 边亲边吃奶的免费视频| 欧美成人a在线观看| 亚洲丝袜综合中文字幕| 成人三级黄色视频| 人妻少妇偷人精品九色| 变态另类成人亚洲欧美熟女| 99九九线精品视频在线观看视频| 亚洲三级黄色毛片| 色综合亚洲欧美另类图片| 人妻少妇偷人精品九色| h日本视频在线播放| 哪里可以看免费的av片| 精品日产1卡2卡| 国产午夜福利久久久久久| 最近视频中文字幕2019在线8| 亚洲久久久久久中文字幕| av在线蜜桃| 一个人免费在线观看电影| 老师上课跳d突然被开到最大视频| 午夜福利高清视频| 午夜老司机福利剧场| 国产精品永久免费网站| 亚洲激情五月婷婷啪啪| 在现免费观看毛片| av在线亚洲专区| 丰满人妻一区二区三区视频av| 国产视频内射| 成人毛片60女人毛片免费| 国产高清有码在线观看视频| 国产在线男女| 91麻豆精品激情在线观看国产| 精品无人区乱码1区二区| 成人亚洲欧美一区二区av| 中文字幕熟女人妻在线| 日韩大尺度精品在线看网址| 亚洲成人中文字幕在线播放| 国产人妻一区二区三区在| 岛国毛片在线播放| 深爱激情五月婷婷| 亚洲丝袜综合中文字幕| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 人人妻人人看人人澡| 欧美日韩精品成人综合77777| 国内精品宾馆在线| 成人亚洲欧美一区二区av| 一级二级三级毛片免费看| 中文在线观看免费www的网站| 日韩精品青青久久久久久| videossex国产| 99久久精品一区二区三区| 九九久久精品国产亚洲av麻豆| 精品无人区乱码1区二区| 国产亚洲精品久久久久久毛片| 听说在线观看完整版免费高清| 久久久a久久爽久久v久久| 日本与韩国留学比较| 老司机影院成人| 日韩,欧美,国产一区二区三区 | 熟女人妻精品中文字幕| 日本-黄色视频高清免费观看| 久久综合国产亚洲精品| 国产精品电影一区二区三区| 亚洲精品色激情综合| 中出人妻视频一区二区| 欧美日本亚洲视频在线播放| 久99久视频精品免费| 网址你懂的国产日韩在线| 欧美一区二区国产精品久久精品| 国产69精品久久久久777片| 欧美人与善性xxx| 久久精品91蜜桃| 99久久精品热视频| 搡老妇女老女人老熟妇| 蜜桃久久精品国产亚洲av| 日本成人三级电影网站| 亚洲国产欧美在线一区| 久久精品夜色国产| 欧美+日韩+精品| 91午夜精品亚洲一区二区三区| 最近视频中文字幕2019在线8| 久久久a久久爽久久v久久| www日本黄色视频网| 乱系列少妇在线播放| .国产精品久久| 听说在线观看完整版免费高清| 成人美女网站在线观看视频| 免费av毛片视频| 久久久国产成人精品二区| 我的老师免费观看完整版| 99久久无色码亚洲精品果冻| www.av在线官网国产| 久久久精品94久久精品| 青春草视频在线免费观看| 成人国产麻豆网| 99热全是精品| 日韩人妻高清精品专区| 好男人在线观看高清免费视频| 少妇高潮的动态图| 成人午夜高清在线视频| 久久久久久久久久成人| 免费看美女性在线毛片视频| 日日摸夜夜添夜夜添av毛片| 激情 狠狠 欧美| 国产探花在线观看一区二区| 99热只有精品国产| 国产又黄又爽又无遮挡在线| 日日撸夜夜添| 国内久久婷婷六月综合欲色啪| 日韩欧美精品v在线| 简卡轻食公司| av在线蜜桃| 亚洲精品久久国产高清桃花| 国产黄片美女视频| 国产黄片视频在线免费观看| 亚洲性久久影院| 变态另类丝袜制服| 欧美一级a爱片免费观看看| 国内揄拍国产精品人妻在线| 如何舔出高潮| 3wmmmm亚洲av在线观看| 久久精品国产99精品国产亚洲性色| 大型黄色视频在线免费观看| 国产真实乱freesex| 亚洲人成网站在线播放欧美日韩| 一区二区三区免费毛片| 国产黄片美女视频| 天天躁日日操中文字幕| 精品少妇黑人巨大在线播放 | 亚洲人与动物交配视频| 又粗又爽又猛毛片免费看| 国内精品宾馆在线| 精品一区二区三区人妻视频| 国产亚洲av嫩草精品影院| 国产黄片美女视频| 亚洲五月天丁香| 亚洲欧美日韩东京热| 国产精品精品国产色婷婷| 国产在线男女| or卡值多少钱| 国产成人午夜福利电影在线观看| 国产视频首页在线观看| 精品国内亚洲2022精品成人| 有码 亚洲区| 成人综合一区亚洲| 99在线人妻在线中文字幕| 精品久久久久久久久久免费视频| 国产精品一及| 国产大屁股一区二区在线视频| 成人国产麻豆网| 成熟少妇高潮喷水视频| eeuss影院久久| 国产精品伦人一区二区| 精品久久久久久久久亚洲| 色视频www国产| 亚洲精品久久久久久婷婷小说 | 亚洲精品456在线播放app| 人人妻人人澡欧美一区二区| 久久久久性生活片| 国产极品天堂在线| 成人鲁丝片一二三区免费| 久久99精品国语久久久| 久久99热6这里只有精品| 亚洲激情五月婷婷啪啪| 日本撒尿小便嘘嘘汇集6| 神马国产精品三级电影在线观看| 国产日本99.免费观看| 中文字幕av成人在线电影| 乱系列少妇在线播放| 日韩欧美 国产精品| 给我免费播放毛片高清在线观看| 国产精品久久久久久久电影| 午夜福利视频1000在线观看| 1000部很黄的大片| 国产黄色视频一区二区在线观看 | 国产国拍精品亚洲av在线观看| 精品久久久久久久久久免费视频| 亚洲国产欧美在线一区| 99久久精品一区二区三区| 国产精品三级大全| 欧美变态另类bdsm刘玥| 淫秽高清视频在线观看| 欧美+日韩+精品| 少妇裸体淫交视频免费看高清| 97超碰精品成人国产| 激情 狠狠 欧美| 在线播放无遮挡| 欧美性猛交╳xxx乱大交人| 18+在线观看网站| 一级毛片电影观看 | 欧美激情久久久久久爽电影| 亚洲aⅴ乱码一区二区在线播放| 99九九线精品视频在线观看视频| 亚洲国产日韩欧美精品在线观看| 成人漫画全彩无遮挡| 国产av麻豆久久久久久久| 国产亚洲91精品色在线| 亚洲人成网站在线观看播放| 乱人视频在线观看| 亚洲欧美日韩东京热| 两个人视频免费观看高清| 黄色一级大片看看| 欧美成人a在线观看| av天堂中文字幕网| 毛片女人毛片| 欧美精品一区二区大全| 熟女人妻精品中文字幕| 国产三级在线视频| 欧美极品一区二区三区四区| 少妇被粗大猛烈的视频| 国产在视频线在精品| 在线免费观看的www视频| 欧美成人精品欧美一级黄| 成人毛片a级毛片在线播放| 1024手机看黄色片| 国产女主播在线喷水免费视频网站 | 亚洲第一区二区三区不卡| 午夜福利视频1000在线观看| 22中文网久久字幕| 国产69精品久久久久777片| av国产免费在线观看| 精品日产1卡2卡| 男女做爰动态图高潮gif福利片| 亚洲图色成人| 97超视频在线观看视频| 看片在线看免费视频| 亚洲最大成人av| 日韩成人伦理影院| 最近中文字幕高清免费大全6| 九草在线视频观看| 国产成人影院久久av| 日本黄大片高清| av女优亚洲男人天堂| 亚洲精品自拍成人| 免费看av在线观看网站| 久久精品国产亚洲av香蕉五月| 国产在视频线在精品| 国产精品人妻久久久影院| 国产激情偷乱视频一区二区| 观看美女的网站| 看免费成人av毛片| 成人无遮挡网站| 欧美三级亚洲精品| 午夜激情福利司机影院| 精华霜和精华液先用哪个| 欧美极品一区二区三区四区| 亚洲av.av天堂| 国产在线精品亚洲第一网站| 国产一区二区亚洲精品在线观看| 国内精品久久久久精免费| 精品久久国产蜜桃| 精品少妇黑人巨大在线播放 | 免费黄网站久久成人精品| 美女黄网站色视频| 中文字幕免费在线视频6| 国产亚洲5aaaaa淫片| 亚洲五月天丁香| 尾随美女入室| 欧美性猛交╳xxx乱大交人| 精品久久久久久久人妻蜜臀av| 97超视频在线观看视频| 欧美高清性xxxxhd video| 天天躁夜夜躁狠狠久久av| 欧美激情在线99| 不卡视频在线观看欧美| 天美传媒精品一区二区| 亚洲欧美日韩东京热| videossex国产| 国产精品一区二区性色av| 韩国av在线不卡| 久久精品国产自在天天线| 国产精品久久久久久久久免| 美女黄网站色视频| 国产精品人妻久久久久久| 日本免费a在线| 亚洲色图av天堂| 97热精品久久久久久| 寂寞人妻少妇视频99o| 秋霞在线观看毛片| 3wmmmm亚洲av在线观看| 亚洲在线自拍视频| 欧美性感艳星| 人人妻人人澡人人爽人人夜夜 | 亚洲精品色激情综合| 久久久久久久久中文| 成人午夜高清在线视频| 亚洲中文字幕一区二区三区有码在线看| 国产精品,欧美在线| 免费观看精品视频网站| 少妇裸体淫交视频免费看高清| 伦理电影大哥的女人| 给我免费播放毛片高清在线观看| 久久欧美精品欧美久久欧美| 高清毛片免费看| 国产av不卡久久| 亚洲欧美成人精品一区二区| 精品久久久久久久末码| 国产女主播在线喷水免费视频网站 | 麻豆国产97在线/欧美| 久久久久久久久久成人| 久久久久九九精品影院| 又粗又爽又猛毛片免费看| 亚洲av不卡在线观看| 亚洲国产精品sss在线观看| 婷婷色av中文字幕| 深夜a级毛片| 国产中年淑女户外野战色| 亚洲第一区二区三区不卡| 久久6这里有精品| 日韩一区二区视频免费看| 久久99热这里只有精品18| 成人亚洲精品av一区二区| 亚洲欧美精品自产自拍| 日韩欧美精品免费久久| www.av在线官网国产| 青青草视频在线视频观看| 日韩欧美在线乱码| 亚州av有码| 别揉我奶头 嗯啊视频| 哪个播放器可以免费观看大片| 国产日本99.免费观看| 三级国产精品欧美在线观看| 18禁裸乳无遮挡免费网站照片| 麻豆成人av视频| 深爱激情五月婷婷| 欧美人与善性xxx| 搡老妇女老女人老熟妇| 亚洲欧美日韩东京热| 亚洲av熟女| 看免费成人av毛片| 免费搜索国产男女视频| 亚洲av二区三区四区| 岛国在线免费视频观看| 少妇高潮的动态图| 成人午夜高清在线视频| 我的老师免费观看完整版| 亚洲精品456在线播放app| av免费在线看不卡| 真实男女啪啪啪动态图| 爱豆传媒免费全集在线观看| 国产在线精品亚洲第一网站| 久久精品夜夜夜夜夜久久蜜豆| 国产69精品久久久久777片| 日韩一区二区三区影片| 身体一侧抽搐| 老司机福利观看| 一边亲一边摸免费视频| 色综合色国产| 简卡轻食公司| 亚洲美女视频黄频| 99久久精品一区二区三区| 久久韩国三级中文字幕| 欧美+亚洲+日韩+国产| 亚洲国产精品久久男人天堂| a级毛片a级免费在线| 网址你懂的国产日韩在线| 最近2019中文字幕mv第一页| 日韩精品有码人妻一区| 国产精品电影一区二区三区| 夜夜爽天天搞| 此物有八面人人有两片| 午夜久久久久精精品| 亚洲精品国产成人久久av| 国产高清不卡午夜福利| 一级黄片播放器| 国产精品一及| 亚洲在久久综合| 麻豆国产av国片精品| 午夜a级毛片| 中文字幕av在线有码专区| 在线观看午夜福利视频| 欧美又色又爽又黄视频| 丝袜喷水一区| 日日干狠狠操夜夜爽| 国产精品综合久久久久久久免费| 成人亚洲精品av一区二区| 国产一区二区在线av高清观看| 国产精品国产高清国产av| 国产麻豆成人av免费视频| 国产成人91sexporn| 欧美高清成人免费视频www| 一本久久精品| 免费观看a级毛片全部| 欧美一级a爱片免费观看看| 亚洲欧美日韩卡通动漫| 国产三级中文精品| 成人特级黄色片久久久久久久| 秋霞在线观看毛片| 在线免费观看的www视频| 欧美+日韩+精品| 久久久国产成人免费| 久久亚洲精品不卡| 97超碰精品成人国产| 美女cb高潮喷水在线观看| 亚洲欧美日韩无卡精品| 欧美潮喷喷水| 成人无遮挡网站| 欧美日韩在线观看h| 免费人成在线观看视频色| 美女内射精品一级片tv| 亚洲av中文av极速乱| 男女那种视频在线观看| 女人十人毛片免费观看3o分钟| 久久精品人妻少妇| 老女人水多毛片| 观看美女的网站| 午夜福利高清视频| 午夜久久久久精精品| 最近中文字幕高清免费大全6| 22中文网久久字幕| 国内久久婷婷六月综合欲色啪| 综合色av麻豆| av在线亚洲专区| 亚洲国产精品成人久久小说 | 国产成人精品一,二区 | 成人毛片60女人毛片免费| 日本三级黄在线观看| 免费观看人在逋| 日本免费一区二区三区高清不卡| 精品人妻一区二区三区麻豆| 日本爱情动作片www.在线观看| 中文字幕av成人在线电影| 内地一区二区视频在线| 欧美一区二区亚洲| ponron亚洲| 又爽又黄a免费视频| 中国国产av一级| 美女脱内裤让男人舔精品视频 | 国产爱豆传媒在线观看|