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

    基于B樣條插值法的柔性矩形薄板的動力學(xué)分析

    2019-11-27 03:21王勇章定國范紀華
    振動工程學(xué)報 2019年5期
    關(guān)鍵詞:固有頻率

    王勇 章定國 范紀華

    摘要: 采用B樣條插值方法研究柔性矩形薄板的動力學(xué)特性??紤]薄板的面外變形、面內(nèi)變形以及面外變形引起的面內(nèi)變形,利用B樣條插值方法對柔性薄板的變形場進行離散,以拉格朗日方程為基礎(chǔ)推導(dǎo)出作大范圍運動柔性薄板的動力學(xué)方程,并運用MATLAB軟件對薄板動力學(xué)仿真問題進行編程。通過動力學(xué)仿真,對比分析了B樣條插值法、假設(shè)模態(tài)法以及有限元法的仿真結(jié)果,驗證了B樣條插值方法的正確性,并表明B樣條插值法在處理柔性薄板的大變形問題的計算精度上具有優(yōu)良性能和推廣潛力。

    關(guān)鍵詞: 多體動力學(xué); 柔性矩形薄板; B樣條插值; 固有頻率

    中圖分類號: O313.7 ?文獻標志碼: A ?文章編號: 1004-4523(2019)05-0811-11

    DOI:10.16385/j.cnki.issn.1004-4523.2019.05.009

    引 言

    柔性體的變形場離散問題是柔性多體系統(tǒng)動力學(xué)中的基本問題[1]。目前假設(shè)模態(tài)法、有限段法和有限元法是三種使用較為廣泛的變形場離散方法。假設(shè)模態(tài)法最大的優(yōu)點是計算效率高[2]。但當工程結(jié)構(gòu)為非規(guī)則的梁式或板式結(jié)構(gòu)的時候,假設(shè)模態(tài)法便有了一定的局限性。一是對于復(fù)雜的柔性體結(jié)構(gòu)很難求出振型函數(shù),二是當工程模型發(fā)生改變的時候,為了滿足精度要求需要疊加更多的模態(tài),這樣降低了計算效率。有限段法在梁式構(gòu)件的離散上應(yīng)用較廣[1]。有限段法精度存在問題的原因在于彈簧的總剛度是由每一小段彈簧剛度向節(jié)點等效移植得到的,而且隨著系統(tǒng)自由度數(shù)增大到一定程度,數(shù)值計算難度也會增大。有限元法的缺點是精度分析需要耗費大量計算資源。因此,對于多體系統(tǒng)動力學(xué)問題仍需要關(guān)注和探索新的變形場離散方法。

    文獻[3]采用樣條有限點法研究了功能梯度壓電梁板的靜力問題,并與已有文獻資料對比證明了該方法的合理性和有效性。文獻[4]提出一個基于B樣條函數(shù)、高階剪切變形理論及變分原理分析復(fù)合材料板殼的樣條有限點法,建立了區(qū)別于有限元法和有限段法的針對靜力學(xué)、動力學(xué)、熱效應(yīng)等分析的新計算格式。文獻[5]利用無網(wǎng)格配點法計算了復(fù)合材料層合板自由振動的固有頻率和振型,文獻還研究了薄板樣條徑向基函數(shù)中形狀參數(shù)的選取對計算結(jié)果的收斂性。文獻[6]將板的撓度和剪切應(yīng)變作為場變量的基礎(chǔ),在此基礎(chǔ)上考慮面內(nèi)位移,采用樣條無網(wǎng)格法建立了熱環(huán)境下厚/薄壓電功能梯度板動力分析新的計算格式。文獻[7-8]給出了B 樣條插值方法與絕對節(jié)點坐標有限元法之間的關(guān)系,從而達成了計算機輔助設(shè)計(Computer Aided Design, CAD)中的幾何造型與計算機輔助分析(Computer Aided Analysis, CAA)中變形場描述的統(tǒng)一。在此基礎(chǔ)上文獻[9-10]進一步研究了B樣條插值、Bezier曲線等計算幾何方法和絕對節(jié)點坐標有限元法之間的轉(zhuǎn)化,通過仿真實例研究表明,將CAD 中幾何造型法和CAA中變形場描述相統(tǒng)一既可以達到CAD和CAA之間的等幾何分析,又可以提高CAA 分析精度和效率。

    本文以B樣條插值法作為離散方法,有別于傳統(tǒng)的假設(shè)模態(tài)、有限元等方法,研究其在柔性薄板動力學(xué)建模和分析上的性能,一定程度上拓展了柔性多體變形場離散方法。旋轉(zhuǎn)結(jié)構(gòu)為實際工程中重要的研究對象[11-12],因此本文對做大范圍旋轉(zhuǎn)運動的柔性矩形薄板的動力學(xué)問題進行研究,在計及面外變形、面內(nèi)變形以及面外變形引起的面內(nèi)變形的情況下建立動力學(xué)模型,利用B樣條插值方法對柔性薄板的變形位移場進行離散,以拉格朗日方程為基礎(chǔ)推導(dǎo)出作大范圍運動柔性矩形薄板的動力學(xué)方程,并運用MATLAB軟件對薄板動力學(xué)仿真問題進行編程。進行動力學(xué)仿真,對比分析了B樣條插值法、假設(shè)模態(tài)法以及有限元法的仿真結(jié)果,驗證了B樣條插值方法的正確性,并表明B樣條插值法在處理柔性薄板的大變形問題的計算精度上具有優(yōu)良性能和推廣潛力。

    圖3-6中圖(b)為圖(a)的30-40 s的放大圖,也即為恒定角速度下的響應(yīng)圖。從圖中可以看出,B樣條插值法的計算結(jié)果同假設(shè)模態(tài)法以及有限元法吻合度很高,這說明B樣條插值法的確能夠用于離散矩形薄板位移場。通過比較圖3(a)與圖4(a)、圖5(a)可以看出,當Ω取0.2 Hz時,零次近似模型與一次近似模型的計算結(jié)果比較接近,而當Ω取0.4 Hz和0.6 Hz時,一次近似模型的計算結(jié)果明顯小于零次近似模型,進一步當Ω取0.8 Hz時,零次近似模型的計算結(jié)果發(fā)散,與實際情況不符;而一次近似模型的計算結(jié)果收斂,與實際相符合。證明了附加剛度項在很大程度上影響了系統(tǒng)的動力學(xué)性態(tài), 說明了薄板作高速大范圍運動時,附加剛度項不容忽略[13]。圖3(b)、圖4(b)中零次近似模型振幅比一次近似模型要小,而在圖5(b)中零次近似模型振幅要大于一次近似模型,這主要是兩種模型在30 s時的狀態(tài)量不同導(dǎo)致的。圖3(b)、圖4(b)中一次近似模型在30 s時變形量大于零次近似模型,導(dǎo)致在之后30-40 s內(nèi)加速度不變的情況下變形量大于零次近似模型,圖5(b)情況剛好相反。從振動頻率角度分析,薄板作大范圍旋轉(zhuǎn)運動的角速度為定值時,f一次近似模型>f零次近似模型。且當角速度逐漸增大時,兩者的差值也隨之增大。

    對以上仿真數(shù)據(jù)研究發(fā)現(xiàn),薄板的面外變形量很小。接下來的仿真中取薄板的彈性模量分別為E=5×1010 N/m2,E=3×1010 N/m2和E=1010 N/m2,薄板其他參數(shù)保持不變。薄板仍繞y軸作定軸轉(zhuǎn)動,旋轉(zhuǎn)規(guī)律采用式(55)的形式,Ω=5 Hz,仿真中取T=30 s,總仿真時間t=40 s。圖7,8和9表示薄板外側(cè)角點z方向變形量(面外變形)在三種離散方法下的計算結(jié)果,仿真模型均為考慮耦合變形量的一次近似模型。

    圖7(b),8(b)和9(b)分別為圖7(a),8(a)和9(a)中30-40 s內(nèi)薄板外側(cè)角點在z方向的變形放大圖,也即為恒定角速度薄板振動變形圖。從圖7,8和9可以看出,薄板最大面外變形隨薄板彈性模量的減小而增大,這是因為減小彈性模量導(dǎo)致薄板的柔性增大,在同樣角速度條件下薄板的變形更大。三種不同彈性模量條件下B樣條插值法與有限元法的仿真結(jié)果吻合度很高,而假設(shè)模態(tài)法同有限元法存在一定的誤差,且誤差主要體現(xiàn)在薄板面外變形量的衰減以及恒定角速度時的振動過程。由上可以說明,當薄板面外變形較大時,假設(shè)模態(tài)法的精度會降低,而B樣條插值法的精度滿足要求。

    為了進一步說明大變形條件下假設(shè)模態(tài)法的局限性,將彈性模量減小為E=4×108 N/m2,仿真周期T=16 s,仿真總時間為20 s。圖10表示Ω=0.8 Hz時薄板外側(cè)角點z方向變形量。圖中柔性薄板的最大面外變形超過了4 m,屬于大變形,B樣條插值法同有限元法的結(jié)果仍然基本一致,而假設(shè)模態(tài)法的仿真結(jié)果誤差很大。說明源于結(jié)構(gòu)力學(xué)中固有振型的假設(shè)模態(tài)法,僅適用于小變形情況,不能處理大變形問題,而B樣條插值法和有限元法適用于大變形問題。

    從表1可以看出薄板的固有頻率隨的增大而增大。圖13為采用B樣條插值離散得到的薄板橫向振動前4階固有頻率隨角速度的變化曲線??梢钥闯霰“宓牡?、第2階固有頻率隨角速度增大而增大的速度較為緩慢;第3、第4階固有頻率隨角速度增大的速度先快后慢。

    從表5的數(shù)據(jù)可以看出,當角速度的范圍在0-20 rad/s時,采用B樣條插值方法計算得到的薄板的前3階固有頻率值相對有限元法的計算結(jié)果誤差小于1%,而采用假設(shè)模態(tài)法計算得到的前3階固有頻率值相對有限元法的計算結(jié)果誤差存在大于1%的情形。說明當角速度小于20 rad/s時,采用B樣條插值方法計算柔性矩形薄板的低階固有頻率能夠滿足精度要求。

    4 結(jié) 論

    (1) 對于剛?cè)狁詈闲D(zhuǎn)矩形薄板的動力學(xué)問題,B樣條插值法能很好地離散薄板的變形場,且當轉(zhuǎn)軸位于薄板平面內(nèi)時,該方法同樣適用。

    (2) 當柔性薄板具有較大面外變形時,B樣條插值法的計算精度比假設(shè)模態(tài)法要高,因此B樣條插值法在柔性板大變形動力學(xué)響應(yīng)問題上具有推廣潛力。

    (3)當轉(zhuǎn)速較低時,運用B樣條插值法計算得到的薄板的前3階固有頻率同有限元法的結(jié)果很相近,且誤差小于假設(shè)模態(tài)法同有限元法的誤差,能夠滿足精度要求。

    參考文獻:

    [1] 范紀華,章定國.旋轉(zhuǎn)懸臂梁動力學(xué)的B樣條插值方法[J].機械工程學(xué)報,2012,48(23): 59-64.

    Fan Jihua,Zhang Dingguo. B-spline interpolation method for the dynamics of rotating cantilever beam[J]. Journal of Mechanical Engineering,2012,48(23): 59-64.

    [2] 吳勝寶,章定國,康 新. 剛體-微梁系統(tǒng)的動力學(xué)特性[J]. 機械工程學(xué)報,2010,46(3):76-82.

    Wu Shengbao,Zhang Dingguo,Kang Xin.Dynamic properties of hub-microbeam system[J]. Journal of Mechanical Engineering,2010,46(3):76-82.

    [3] 沈德洲. 功能梯度壓電梁板的樣條有限點分析[D].合肥:合肥工業(yè)大學(xué),2008.

    Shen Dezhou. Spline finite-point method for analysis of functionally graded piezoelectric beam and plate[D]. Hefei:Hefei University of Technology,2008.

    [4] 秦 榮. 復(fù)合材料板殼分析的樣條有限點法[J].工程力學(xué),2001,18(1): 14-22.

    Qin Rong. Spline finite point method for analysis of laminated composite plates and shells[J]. Engineering Mechanics,2001,18(1): 14-22.

    [5] 郭向群,楊 康,項 松,等. 復(fù)合材料層合板自由振動的薄板樣條無網(wǎng)格解[J].材料導(dǎo)報,2016,(S2):166-169.

    Guo Xiangqun, Yang Kang, Xiang Song, et al. Meshless solution based on the thin plate spline for free vibration of laminated composite plates[J]. Materials Review,2016,(S2):166-169.

    [6] 李雙蓓,吳 海,莫春美. 基于樣條無網(wǎng)格法的厚/薄壓電功能梯度板動力分析[J]. 振動與沖擊,2017,36(7):116-122.

    Li Shuangbei, Wu Hai, Mo Chunmei. Dynamic analysis of thick/thin piezoelectric FGM plates with spline meshless method[J].Journal of Vibration and Shock,2017,36(7):116-122.

    [7] SANBORN G G, SHABANA A A. On the integration of computer aided design and analysis using the finite element absolute nodal coordinate formulation[J].Multibody Syst. Dyn.,2009,22:181-197.

    [8] SANBORN G G,SHABANA A A. A rational finite element method based on the absolute modal coordinate formulation[J]. Nonlinear Dyn.,2009,58:565-572.

    [9] LAN Peng,SHABANA A A. Integration of B-spline geometry and ANCF finite element analysis[J]. Nonlinear Dyn., 2010,61:193-206.

    [10] LAN Peng, SHABANA A A. Rational finite elements and flexible body dynamics[J]. Journal of Vibration and Acoustics,2010,132:041007-1-9.

    [11] SALERNO A, OSTROVSKAYA S, ANGELES J. Dynamics of a novel rolling robot-analysis and simulation[J]. Chinese Journal of Mechanical Engineering,2004,17(S1):13-17.

    [12] HU Hong, SHI Hongyan, BEN M R. Active vibration suppression via linearizing hysteresis of piezoceramic actuators[J]. Chinese Journal of Mechanical Engineering,2007,20(4):84-87.

    [13] 吳勝寶. 作大范圍運動柔性梁和柔性薄板剛?cè)狁詈蟿恿W(xué)建模與仿真[D].南京:南京理工大學(xué),2009.

    Wu Shengbao. Modeling and simulation of rigid-flexible coupling dynamics of flexible beams and flexible thin plates with large overall motion[D]. Nanjing:Nanjing University of Science and Technology, 2009.

    [14] Yoo H H,Chung J. Dynamics of rectangular plates undergoing prescribed overall motion[J]. Journal of Sound and Vibration, 2001,239(1):123-137.

    [15] 劉又午,王建明,張大鈞,等.作大范圍運動的矩形板動力分析[J]. 振動與沖擊,1998,17(1):38-43.

    Liu Youwu, Wang Jianmin, Zhang Dajun,et al. Dynamic analysis of rectangnlar plates undergoing large overall motion[J]. Journal of Vibration and Shock, 1998,17(1):38-43.

    [16] 趙飛云,洪嘉振,劉錦陽,等.高速旋轉(zhuǎn)柔性矩形薄板的動力學(xué)建模和近似算法[J]. 振動工程學(xué)報,2006,19(3):416-421.

    Zhao Feiyun, Hong Jiazhen, Liu Jinyang, et al. Dynamic modeling and modal truncation approach for a high-speed rotating thin elastic rectangular plate[J]. Journal of Vibration Engineering, 2006,19(3):416-421.

    [17] 秦 榮.樣條無網(wǎng)格法[M].北京:科學(xué)出版社,2012:1-29,52-81,101-139,141-160,164-222.

    Qin Rong. Spline Meshless Method[M]. Beijing:Science Press,2012:1-29,52-81,101-139,141-160,164-222.

    [18] 蔡國平,洪嘉振. 考慮附加質(zhì)量的中心剛體-柔性懸臂梁系統(tǒng)的動力特性研究[J]. 機械工程學(xué)報,2005,41(2):33-40.

    Cai Guoping,Hong Jiazhen. Dynamics study of hub-beam system with tip mass[J]. Journal of Mechanical Engineering,2005,41(2):33-40.

    [19] 鄭 彤,章定國,廖連芳,等. 航空發(fā)動機葉片剛?cè)狁詈蟿恿W(xué)分析[J].機械工程學(xué)報,2014,50(23): 42-49.

    Zheng Tong,Zhang Dingguo,Liao Lianfang, et al. Rigid-flexible coupling dynamic analysis of aero-engine blades[J]. Journal of Mechanical Engineering, 2014,50(23): 42-49.

    [20] 王建明.柔性體剛?cè)狁詈蟿恿W(xué)建模理論及動力剛化有限元算法研究[R].上海:上海交通大學(xué)博士后研究工作報告,1999.

    Wang Jianming. Research on flexible body with rigid-flexible coupling dynamics modeling theory and dynamics rigid finite element method[R]. Shanghai:Shanghai Jiao Tong University Postdoctoral Research Report,1999.

    Abstract: The dynamic behavior of flexible rectangular thin plates is investigated by the B-spline interpolation method. The off-plane deformation, in-plane deformation and the in-plane shortening caused by the off-plane deformation of the thin plate are all considered. A dynamic model is established and the B-spline interpolation method is used to discretize the deformation field of the flexible thin plate. The rigid-flexible coupling dynamic equations of the flexible thin plate with large overall motion are established via employing the second kind of Lagrange′s equation, and the MATLAB software is used to program for the dynamic simulation problems of the flexible thin plate. To validate the method, dynamic simulations are carried out. The simulation results of the B-spline interpolation method are compared with the ones of the assumed mode method and the finite element method, which demonstrate that the B-spline interpolation method has the potential of popularization in the calculation accuracy to deal with the discretization of the deformation field of the flexible thin plate.

    Key words: multibody dynamics; flexible rectangular thin plate; B-spline interpolation method; natural frequency

    作者簡介: 王 勇(1992-),男,碩士研究生。電話:15195758184;E-mail:1342506994@qq.com

    通訊作者: 章定國(1967-),男,博士,教授,博士生導(dǎo)師。E-mail:zhangdg419@mail.njust.edu.cn

    猜你喜歡
    固有頻率
    梯形和三角形波紋夾芯板的聲振特性研究
    基于SolidWorks某型自卸車車架的模態(tài)分析
    基于舒適度的壓型鋼板混凝土樓蓋板厚取值研究
    某調(diào)速型液力偶合器泵輪的模態(tài)分析
    TMD等強度懸臂梁實際應(yīng)用時的頻率精確分析
    預(yù)應(yīng)力作用下橋式起重機變截面主梁振動分析
    振動壓路機跳振的原因及解決方案
    模態(tài)分析在浮置板軌道參數(shù)優(yōu)化中的應(yīng)用
    某SUV車架多目標拓撲優(yōu)化設(shè)計
    久久午夜亚洲精品久久| 欧美高清成人免费视频www| 中出人妻视频一区二区| 一级黄片播放器| 麻豆成人午夜福利视频| 国产亚洲精品久久久com| 国产精品久久视频播放| 99久久无色码亚洲精品果冻| 亚洲精品影视一区二区三区av| 18禁在线播放成人免费| 极品教师在线免费播放| 日本a在线网址| 国产黄色小视频在线观看| 免费观看人在逋| 女人十人毛片免费观看3o分钟| 久久久久久久精品吃奶| 少妇猛男粗大的猛烈进出视频 | 一区二区三区免费毛片| 给我免费播放毛片高清在线观看| 别揉我奶头~嗯~啊~动态视频| 少妇人妻精品综合一区二区 | 国产成人一区二区在线| 真人做人爱边吃奶动态| 欧美激情国产日韩精品一区| 国产精品一区二区三区四区免费观看 | 国国产精品蜜臀av免费| 日本精品一区二区三区蜜桃| 一级黄色大片毛片| 变态另类成人亚洲欧美熟女| 国产一级毛片七仙女欲春2| 此物有八面人人有两片| 成人三级黄色视频| 精品国产三级普通话版| www.色视频.com| 久久久久久国产a免费观看| 偷拍熟女少妇极品色| 免费人成视频x8x8入口观看| 久久国内精品自在自线图片| 三级毛片av免费| 老司机福利观看| 国产色爽女视频免费观看| 99热精品在线国产| 国产激情偷乱视频一区二区| 不卡视频在线观看欧美| 搞女人的毛片| 日本成人三级电影网站| 国产三级中文精品| 少妇丰满av| 久久久久久久久久黄片| 我要看日韩黄色一级片| 免费搜索国产男女视频| 国产亚洲精品久久久com| 三级毛片av免费| 亚洲avbb在线观看| 国产色爽女视频免费观看| 国产精品人妻久久久久久| 尾随美女入室| 精品福利观看| 国产高潮美女av| 国产三级在线视频| 国产免费av片在线观看野外av| 亚洲熟妇中文字幕五十中出| 国内毛片毛片毛片毛片毛片| 午夜福利在线观看吧| 午夜精品在线福利| 亚洲人成伊人成综合网2020| 婷婷丁香在线五月| 观看免费一级毛片| 赤兔流量卡办理| 欧美日韩国产亚洲二区| 老熟妇乱子伦视频在线观看| 麻豆成人av在线观看| 日韩一本色道免费dvd| 国产 一区精品| 亚州av有码| 日本免费a在线| 亚洲国产精品成人综合色| 可以在线观看的亚洲视频| 色视频www国产| 国产综合懂色| 国产午夜精品论理片| 欧美日韩瑟瑟在线播放| 国产在线精品亚洲第一网站| 国产精品永久免费网站| 色av中文字幕| 日韩欧美在线二视频| 成人精品一区二区免费| 特级一级黄色大片| 一区二区三区激情视频| 国产精品福利在线免费观看| 精品一区二区免费观看| a级毛片免费高清观看在线播放| 在线观看午夜福利视频| 三级毛片av免费| 亚洲第一电影网av| 午夜免费男女啪啪视频观看 | 成人国产综合亚洲| 国产一级毛片七仙女欲春2| 少妇人妻一区二区三区视频| 91麻豆精品激情在线观看国产| 丰满乱子伦码专区| 性插视频无遮挡在线免费观看| 亚洲成人久久爱视频| 国产久久久一区二区三区| 中文字幕高清在线视频| 国产精品1区2区在线观看.| ponron亚洲| 国内揄拍国产精品人妻在线| 99在线人妻在线中文字幕| 欧美成人a在线观看| 国产麻豆成人av免费视频| 午夜福利成人在线免费观看| 日韩一区二区视频免费看| 成年版毛片免费区| 久久久午夜欧美精品| 99热只有精品国产| 日韩av在线大香蕉| 国产精品美女特级片免费视频播放器| 亚洲欧美精品综合久久99| 级片在线观看| 99久久精品一区二区三区| 熟女电影av网| 久久国产乱子免费精品| 亚洲av美国av| 一级毛片久久久久久久久女| 午夜老司机福利剧场| 久久精品国产99精品国产亚洲性色| 九色国产91popny在线| 国产精品乱码一区二三区的特点| 成人特级黄色片久久久久久久| 亚洲欧美日韩无卡精品| 国产精品人妻久久久影院| 91在线精品国自产拍蜜月| 美女 人体艺术 gogo| 高清在线国产一区| 国产伦在线观看视频一区| 日韩欧美 国产精品| 美女高潮的动态| 午夜精品一区二区三区免费看| .国产精品久久| 亚洲自拍偷在线| 国产国拍精品亚洲av在线观看| av视频在线观看入口| 少妇人妻精品综合一区二区 | 免费看a级黄色片| 国产精品国产三级国产av玫瑰| 国产黄色小视频在线观看| 三级毛片av免费| 久久国产精品人妻蜜桃| 日韩精品中文字幕看吧| 日韩欧美国产一区二区入口| 国产精品久久久久久精品电影| 国产精品女同一区二区软件 | 动漫黄色视频在线观看| 婷婷色综合大香蕉| 十八禁网站免费在线| 亚洲精品久久国产高清桃花| 蜜桃亚洲精品一区二区三区| 看片在线看免费视频| av.在线天堂| 国产主播在线观看一区二区| 最近视频中文字幕2019在线8| 毛片女人毛片| 日韩欧美在线二视频| 国产成人一区二区在线| 亚洲成人久久性| 高清日韩中文字幕在线| 一区二区三区四区激情视频 | 最好的美女福利视频网| 久久久久久久久久黄片| 日日干狠狠操夜夜爽| 美女 人体艺术 gogo| 精品一区二区三区av网在线观看| 黄色一级大片看看| 国产午夜精品久久久久久一区二区三区 | 内地一区二区视频在线| a级一级毛片免费在线观看| av福利片在线观看| 久久人人精品亚洲av| 亚洲av成人av| 午夜福利在线观看吧| 欧美中文日本在线观看视频| 精品久久久噜噜| 亚洲最大成人中文| 又粗又爽又猛毛片免费看| 热99在线观看视频| 极品教师在线免费播放| 99在线视频只有这里精品首页| 色尼玛亚洲综合影院| 欧美+日韩+精品| 欧美一区二区精品小视频在线| 国产91精品成人一区二区三区| 在线观看美女被高潮喷水网站| 乱系列少妇在线播放| 国产精品永久免费网站| 91在线观看av| 国产精品一区二区性色av| 午夜影院日韩av| 亚洲av美国av| 中国美女看黄片| 欧美3d第一页| 99精品在免费线老司机午夜| 伦理电影大哥的女人| 自拍偷自拍亚洲精品老妇| 免费高清视频大片| 亚洲五月天丁香| 99精品久久久久人妻精品| 乱码一卡2卡4卡精品| .国产精品久久| 亚洲成人中文字幕在线播放| 九九热线精品视视频播放| 久久久久九九精品影院| 国内久久婷婷六月综合欲色啪| 少妇熟女aⅴ在线视频| 色在线成人网| 亚洲无线在线观看| 亚洲最大成人手机在线| 亚洲第一电影网av| 不卡视频在线观看欧美| 亚洲乱码一区二区免费版| 亚洲精品国产成人久久av| 男人和女人高潮做爰伦理| 精品日产1卡2卡| 尾随美女入室| 日韩精品青青久久久久久| 久久中文看片网| 国产成人影院久久av| 亚洲国产精品久久男人天堂| 不卡视频在线观看欧美| av在线观看视频网站免费| 九色成人免费人妻av| 毛片女人毛片| 中文字幕免费在线视频6| 国产伦在线观看视频一区| 嫩草影视91久久| 床上黄色一级片| 日本一二三区视频观看| 色尼玛亚洲综合影院| 久久久久国产精品人妻aⅴ院| 国产av在哪里看| 九九热线精品视视频播放| 乱系列少妇在线播放| 又黄又爽又刺激的免费视频.| 在线观看免费视频日本深夜| 欧美bdsm另类| 69av精品久久久久久| 久久精品国产99精品国产亚洲性色| 两性午夜刺激爽爽歪歪视频在线观看| 村上凉子中文字幕在线| 国产精品女同一区二区软件 | 黄色日韩在线| 日韩精品有码人妻一区| av在线天堂中文字幕| 久久亚洲精品不卡| 久久久久久伊人网av| 精品久久久久久久末码| 久久久久精品国产欧美久久久| 白带黄色成豆腐渣| 草草在线视频免费看| 国产精品久久视频播放| 亚洲国产精品久久男人天堂| 黄色日韩在线| 日韩精品有码人妻一区| av视频在线观看入口| 欧美极品一区二区三区四区| 欧美激情久久久久久爽电影| 美女免费视频网站| 国产高清有码在线观看视频| 中出人妻视频一区二区| 天天一区二区日本电影三级| 久久精品久久久久久噜噜老黄 | 久久久色成人| 人妻丰满熟妇av一区二区三区| 在线观看免费视频日本深夜| 久久久久久伊人网av| 男人的好看免费观看在线视频| 午夜久久久久精精品| 国产精品亚洲一级av第二区| 黄色女人牲交| 午夜老司机福利剧场| 久久人妻av系列| 欧美区成人在线视频| 99热网站在线观看| 性插视频无遮挡在线免费观看| 亚洲男人的天堂狠狠| 国产精品99久久久久久久久| 自拍偷自拍亚洲精品老妇| 日本免费a在线| 色av中文字幕| 欧美日韩综合久久久久久 | 99久久无色码亚洲精品果冻| bbb黄色大片| 亚洲中文日韩欧美视频| 亚洲精品在线观看二区| 一级毛片久久久久久久久女| 国产视频内射| 国产精品野战在线观看| 国产一区二区三区视频了| 欧美日韩精品成人综合77777| 中出人妻视频一区二区| 联通29元200g的流量卡| 又紧又爽又黄一区二区| 美女大奶头视频| 亚洲专区中文字幕在线| 免费不卡的大黄色大毛片视频在线观看 | 久久久久久大精品| 精品久久久久久久末码| 又爽又黄无遮挡网站| 久久久久免费精品人妻一区二区| 国产高清视频在线播放一区| 亚洲精品色激情综合| 精品人妻熟女av久视频| av在线蜜桃| 最近最新免费中文字幕在线| 黄色一级大片看看| 一区二区三区四区激情视频 | 天堂影院成人在线观看| 国产精品av视频在线免费观看| 两人在一起打扑克的视频| ponron亚洲| 久久午夜福利片| 搞女人的毛片| 舔av片在线| 久久久久久久亚洲中文字幕| 亚洲精品成人久久久久久| 久久精品国产亚洲网站| 国产爱豆传媒在线观看| 1024手机看黄色片| 永久网站在线| 美女cb高潮喷水在线观看| 免费无遮挡裸体视频| 天美传媒精品一区二区| 国产精品99久久久久久久久| 偷拍熟女少妇极品色| 精品久久久久久久人妻蜜臀av| 亚洲精品色激情综合| 日韩欧美国产在线观看| 久久这里只有精品中国| 亚洲国产精品久久男人天堂| 精品国产三级普通话版| 内地一区二区视频在线| 男女做爰动态图高潮gif福利片| 国产v大片淫在线免费观看| 91午夜精品亚洲一区二区三区 | 少妇人妻一区二区三区视频| 成人二区视频| 欧美日韩精品成人综合77777| 又爽又黄无遮挡网站| 亚洲精品一区av在线观看| 熟女人妻精品中文字幕| 91午夜精品亚洲一区二区三区 | 美女大奶头视频| 欧美3d第一页| 亚洲成人精品中文字幕电影| 亚洲成a人片在线一区二区| 亚洲 国产 在线| 精品一区二区三区视频在线| ponron亚洲| 久99久视频精品免费| 国产精品一区二区性色av| 欧美成人一区二区免费高清观看| 国产高清不卡午夜福利| 禁无遮挡网站| 不卡一级毛片| 国产真实乱freesex| 精品一区二区三区人妻视频| 精品一区二区三区视频在线| 亚洲 国产 在线| bbb黄色大片| 一进一出抽搐gif免费好疼| 精品欧美国产一区二区三| av福利片在线观看| 午夜老司机福利剧场| 国产高清不卡午夜福利| 国产亚洲av嫩草精品影院| 熟女电影av网| 日本欧美国产在线视频| 我的老师免费观看完整版| 俺也久久电影网| 变态另类成人亚洲欧美熟女| 久久天躁狠狠躁夜夜2o2o| 欧美最新免费一区二区三区| 春色校园在线视频观看| 国产熟女欧美一区二区| 国内精品宾馆在线| 精品人妻偷拍中文字幕| 亚洲欧美激情综合另类| 国产三级在线视频| 国产精品伦人一区二区| 欧美bdsm另类| 日韩欧美在线二视频| www日本黄色视频网| 成人鲁丝片一二三区免费| 亚洲精品一区av在线观看| 免费看av在线观看网站| 久久国产乱子免费精品| 免费黄网站久久成人精品| 午夜福利欧美成人| 国产精品永久免费网站| 久久精品国产清高在天天线| 午夜精品一区二区三区免费看| 日韩欧美一区二区三区在线观看| 十八禁国产超污无遮挡网站| 亚洲熟妇熟女久久| 午夜日韩欧美国产| 免费大片18禁| 黄色配什么色好看| 嫩草影院入口| 亚洲成人免费电影在线观看| 成年人黄色毛片网站| 亚洲一区二区三区色噜噜| 日韩亚洲欧美综合| 国产三级中文精品| 啦啦啦啦在线视频资源| 午夜a级毛片| 九九爱精品视频在线观看| 国产蜜桃级精品一区二区三区| 亚洲第一电影网av| 五月玫瑰六月丁香| 天堂av国产一区二区熟女人妻| 亚洲无线在线观看| 特级一级黄色大片| 亚洲aⅴ乱码一区二区在线播放| 少妇高潮的动态图| 最近最新免费中文字幕在线| 日本黄大片高清| 国产精品一区www在线观看 | 亚洲精品乱码久久久v下载方式| 香蕉av资源在线| 国产在线精品亚洲第一网站| 欧美日韩中文字幕国产精品一区二区三区| 亚洲,欧美,日韩| 极品教师在线免费播放| 欧美精品啪啪一区二区三区| 99riav亚洲国产免费| 动漫黄色视频在线观看| 中文字幕精品亚洲无线码一区| 欧美日本亚洲视频在线播放| 成熟少妇高潮喷水视频| 久久久久久久久大av| 男女那种视频在线观看| 国国产精品蜜臀av免费| 国产午夜精品论理片| 91在线精品国自产拍蜜月| 精品福利观看| 日本欧美国产在线视频| 日韩欧美国产一区二区入口| netflix在线观看网站| 丰满乱子伦码专区| 两性午夜刺激爽爽歪歪视频在线观看| 91久久精品国产一区二区成人| 91精品国产九色| 国内毛片毛片毛片毛片毛片| 婷婷丁香在线五月| 国语自产精品视频在线第100页| 我的女老师完整版在线观看| 麻豆成人午夜福利视频| 中国美白少妇内射xxxbb| 成人国产麻豆网| 久久久久久久久大av| 成人av在线播放网站| 欧美+日韩+精品| 热99在线观看视频| 久久九九热精品免费| 日韩欧美在线乱码| avwww免费| 久久九九热精品免费| 免费av观看视频| avwww免费| 午夜精品久久久久久毛片777| 亚洲国产日韩欧美精品在线观看| 午夜福利欧美成人| 亚洲av免费高清在线观看| 18+在线观看网站| www日本黄色视频网| 成年人黄色毛片网站| 久久久久免费精品人妻一区二区| 热99在线观看视频| 人妻少妇偷人精品九色| 日韩,欧美,国产一区二区三区 | 免费观看在线日韩| or卡值多少钱| 国产av一区在线观看免费| 午夜福利高清视频| 日日干狠狠操夜夜爽| 久久久久久久亚洲中文字幕| 日韩欧美国产一区二区入口| 久久午夜福利片| 欧美高清成人免费视频www| 久久99热这里只有精品18| 99热6这里只有精品| 国产高清视频在线播放一区| 性色avwww在线观看| 久久久久久久久中文| 黄色女人牲交| 在线免费观看不下载黄p国产 | 制服丝袜大香蕉在线| 亚洲在线自拍视频| 国产一区二区三区av在线 | 国产老妇女一区| or卡值多少钱| 久久国产精品人妻蜜桃| 国产一区二区三区av在线 | av在线蜜桃| 久久久国产成人免费| 两个人的视频大全免费| 久久九九热精品免费| 麻豆久久精品国产亚洲av| 99在线视频只有这里精品首页| 国产精品综合久久久久久久免费| 99精品久久久久人妻精品| 在线观看av片永久免费下载| 三级毛片av免费| 自拍偷自拍亚洲精品老妇| 国产单亲对白刺激| 久久精品久久久久久噜噜老黄 | 尾随美女入室| 乱人视频在线观看| 中文字幕高清在线视频| 精品久久国产蜜桃| 窝窝影院91人妻| 国产精品人妻久久久影院| 亚洲人成网站高清观看| 久久久久久伊人网av| 最新在线观看一区二区三区| 国产精品久久久久久久电影| 欧美日韩乱码在线| 大又大粗又爽又黄少妇毛片口| 久久亚洲精品不卡| 国产日本99.免费观看| 免费观看在线日韩| 夜夜爽天天搞| 中国美白少妇内射xxxbb| 久久精品影院6| 国产在视频线在精品| 国产一区二区三区av在线 | 亚洲国产高清在线一区二区三| 欧美成人性av电影在线观看| 成人美女网站在线观看视频| avwww免费| 成人二区视频| 嫩草影视91久久| 午夜精品在线福利| 亚洲色图av天堂| 午夜亚洲福利在线播放| 国产av麻豆久久久久久久| 九九热线精品视视频播放| 在线观看舔阴道视频| 美女xxoo啪啪120秒动态图| a级一级毛片免费在线观看| 国国产精品蜜臀av免费| 亚洲av中文字字幕乱码综合| 日韩强制内射视频| 色5月婷婷丁香| 岛国在线免费视频观看| 九九在线视频观看精品| 麻豆国产av国片精品| 日韩一本色道免费dvd| 欧美色欧美亚洲另类二区| 色在线成人网| 久久精品国产亚洲av天美| 亚洲av免费在线观看| 亚洲人成网站高清观看| 国产精品三级大全| 中文资源天堂在线| 一个人免费在线观看电影| 99久久九九国产精品国产免费| 成人无遮挡网站| 麻豆久久精品国产亚洲av| 亚洲人成伊人成综合网2020| 一a级毛片在线观看| 又黄又爽又刺激的免费视频.| 天堂√8在线中文| 看黄色毛片网站| 69人妻影院| 麻豆成人av在线观看| 在现免费观看毛片| av天堂在线播放| 人妻久久中文字幕网| 久久精品久久久久久噜噜老黄 | 亚洲第一区二区三区不卡| 国产精品美女特级片免费视频播放器| 免费看a级黄色片| 亚洲不卡免费看| 亚洲国产欧洲综合997久久,| 搡老熟女国产l中国老女人| 亚洲黑人精品在线| 不卡一级毛片| 国产视频内射| 婷婷六月久久综合丁香| 国产av麻豆久久久久久久| 久久6这里有精品| 国产一级毛片七仙女欲春2| 国产视频内射| 国产精品久久久久久精品电影| 可以在线观看毛片的网站| 最近中文字幕高清免费大全6 | 国产一区二区亚洲精品在线观看| 伦精品一区二区三区| 久久天躁狠狠躁夜夜2o2o| 22中文网久久字幕| 亚洲av不卡在线观看| 少妇丰满av| 成人综合一区亚洲| 亚洲七黄色美女视频| 少妇高潮的动态图| 国产精品一区www在线观看 | 级片在线观看| 亚洲最大成人av| 一级黄色大片毛片| 三级毛片av免费| 亚洲综合色惰| 免费av观看视频| 丝袜美腿在线中文| 午夜福利在线观看免费完整高清在 |