• <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è)計
    久久久久精品人妻al黑| 天天添夜夜摸| 欧美 日韩 精品 国产| 黄色丝袜av网址大全| 亚洲国产精品一区二区三区在线| 欧美av亚洲av综合av国产av| 久久国产精品影院| 啦啦啦视频在线资源免费观看| 可以免费在线观看a视频的电影网站| 亚洲一区二区三区欧美精品| a级片在线免费高清观看视频| 午夜免费成人在线视频| 亚洲欧美一区二区三区久久| 欧美精品av麻豆av| 国产欧美日韩综合在线一区二区| 成年人免费黄色播放视频| 精品人妻1区二区| 女人爽到高潮嗷嗷叫在线视频| 欧美中文综合在线视频| 久久久久国内视频| 99re6热这里在线精品视频| 天天影视国产精品| 天天躁夜夜躁狠狠躁躁| 成年人午夜在线观看视频| 免费在线观看日本一区| 久久香蕉激情| 欧美在线一区亚洲| 777久久人妻少妇嫩草av网站| 黑人欧美特级aaaaaa片| 日本撒尿小便嘘嘘汇集6| 亚洲精华国产精华精| 精品国产乱码久久久久久男人| 女性生殖器流出的白浆| 亚洲精品国产精品久久久不卡| 色在线成人网| 亚洲va日本ⅴa欧美va伊人久久| 亚洲avbb在线观看| 日韩一卡2卡3卡4卡2021年| 高清毛片免费观看视频网站 | 久久午夜综合久久蜜桃| 狂野欧美激情性xxxx| 国产国语露脸激情在线看| 国产欧美日韩综合在线一区二区| 香蕉国产在线看| 国产av又大| 午夜精品久久久久久毛片777| 日本欧美视频一区| 亚洲自偷自拍图片 自拍| 两人在一起打扑克的视频| 人人妻人人澡人人爽人人夜夜| 身体一侧抽搐| 国产成人啪精品午夜网站| 欧美成狂野欧美在线观看| 亚洲一码二码三码区别大吗| 日韩熟女老妇一区二区性免费视频| 91精品三级在线观看| 欧美不卡视频在线免费观看 | 久久久精品国产亚洲av高清涩受| 精品国产乱码久久久久久男人| 国产精品二区激情视频| 午夜福利视频在线观看免费| 精品国产亚洲在线| 亚洲午夜理论影院| 男女高潮啪啪啪动态图| 国产单亲对白刺激| 久久性视频一级片| 欧美成狂野欧美在线观看| 99香蕉大伊视频| 大香蕉久久成人网| 精品福利永久在线观看| 日韩熟女老妇一区二区性免费视频| 成年动漫av网址| 黑人巨大精品欧美一区二区蜜桃| av免费在线观看网站| av电影中文网址| 日韩免费高清中文字幕av| 久久精品亚洲熟妇少妇任你| 在线观看舔阴道视频| 久久精品国产综合久久久| 日韩免费高清中文字幕av| 欧美日韩精品网址| 王馨瑶露胸无遮挡在线观看| 一级片'在线观看视频| 两人在一起打扑克的视频| 免费在线观看亚洲国产| x7x7x7水蜜桃| 久9热在线精品视频| 久久精品成人免费网站| 母亲3免费完整高清在线观看| 天堂中文最新版在线下载| 精品一品国产午夜福利视频| 欧美av亚洲av综合av国产av| 变态另类成人亚洲欧美熟女 | 国产男女超爽视频在线观看| 亚洲精品国产精品久久久不卡| 亚洲av第一区精品v没综合| 久久久久久久国产电影| 日日爽夜夜爽网站| 日本vs欧美在线观看视频| 高清黄色对白视频在线免费看| 日韩欧美一区视频在线观看| 搡老岳熟女国产| 日韩制服丝袜自拍偷拍| 999久久久国产精品视频| 丝瓜视频免费看黄片| 亚洲伊人色综图| 99re在线观看精品视频| 久久狼人影院| 国产淫语在线视频| 久久天堂一区二区三区四区| 日韩 欧美 亚洲 中文字幕| 亚洲熟女毛片儿| 高清欧美精品videossex| 精品一区二区三区av网在线观看| 国产亚洲欧美精品永久| 久久久精品国产亚洲av高清涩受| 国产成人影院久久av| 黄频高清免费视频| 国产精品一区二区在线观看99| 亚洲中文av在线| 美女高潮到喷水免费观看| 黑丝袜美女国产一区| 亚洲精品美女久久久久99蜜臀| 午夜福利欧美成人| 日韩制服丝袜自拍偷拍| 欧美亚洲日本最大视频资源| 水蜜桃什么品种好| 亚洲av片天天在线观看| 免费在线观看日本一区| 高清在线国产一区| 欧美日韩亚洲国产一区二区在线观看 | 丁香欧美五月| 新久久久久国产一级毛片| 国产野战对白在线观看| 老熟女久久久| 精品国产亚洲在线| 99久久99久久久精品蜜桃| 国产激情欧美一区二区| 亚洲欧美日韩高清在线视频| 久久 成人 亚洲| 久久精品亚洲熟妇少妇任你| 亚洲精品成人av观看孕妇| 亚洲国产欧美日韩在线播放| 久久国产乱子伦精品免费另类| 久久精品国产亚洲av香蕉五月 | 麻豆av在线久日| 视频在线观看一区二区三区| 中亚洲国语对白在线视频| 久久精品91无色码中文字幕| 欧美精品人与动牲交sv欧美| 亚洲黑人精品在线| 欧美色视频一区免费| 亚洲一区高清亚洲精品| 国产人伦9x9x在线观看| 精品高清国产在线一区| 欧美人与性动交α欧美软件| 久久久久久免费高清国产稀缺| 美女视频免费永久观看网站| 精品人妻1区二区| 99re在线观看精品视频| 99热只有精品国产| 天天添夜夜摸| 欧美日韩视频精品一区| 亚洲人成77777在线视频| 精品人妻在线不人妻| 丝袜在线中文字幕| 高清欧美精品videossex| 亚洲精品国产色婷婷电影| 波多野结衣一区麻豆| 性色av乱码一区二区三区2| 免费女性裸体啪啪无遮挡网站| 亚洲人成77777在线视频| 少妇裸体淫交视频免费看高清 | av片东京热男人的天堂| √禁漫天堂资源中文www| 亚洲色图综合在线观看| 免费一级毛片在线播放高清视频 | 丝瓜视频免费看黄片| 三上悠亚av全集在线观看| 国产成人精品久久二区二区免费| 国产区一区二久久| 精品免费久久久久久久清纯 | 国产99久久九九免费精品| 久久 成人 亚洲| 美国免费a级毛片| 免费日韩欧美在线观看| 一级毛片女人18水好多| 最新美女视频免费是黄的| 亚洲欧洲精品一区二区精品久久久| 色尼玛亚洲综合影院| 制服人妻中文乱码| 人人澡人人妻人| 中文字幕高清在线视频| 怎么达到女性高潮| 丁香六月欧美| 9191精品国产免费久久| 一本综合久久免费| 成人永久免费在线观看视频| 热99国产精品久久久久久7| 一本大道久久a久久精品| 午夜激情av网站| 天堂俺去俺来也www色官网| 欧美日韩亚洲高清精品| 在线观看www视频免费| 亚洲三区欧美一区| 一边摸一边抽搐一进一出视频| 变态另类成人亚洲欧美熟女 | 亚洲人成77777在线视频| 亚洲精品中文字幕一二三四区| 欧美亚洲 丝袜 人妻 在线| 操美女的视频在线观看| 法律面前人人平等表现在哪些方面| 两人在一起打扑克的视频| 久久久久国产一级毛片高清牌| 久久久久久亚洲精品国产蜜桃av| 中出人妻视频一区二区| 黑人欧美特级aaaaaa片| 午夜日韩欧美国产| 精品一区二区三区av网在线观看| 99久久99久久久精品蜜桃| 欧美国产精品va在线观看不卡| 国产亚洲av高清不卡| 99re在线观看精品视频| 日本黄色日本黄色录像| 最近最新中文字幕大全免费视频| 亚洲一区二区三区欧美精品| 一区在线观看完整版| 久久草成人影院| 中文字幕人妻熟女乱码| a在线观看视频网站| 丝袜美腿诱惑在线| 老司机影院毛片| 国产免费现黄频在线看| 亚洲av美国av| 在线观看免费视频网站a站| 久久国产精品影院| 不卡一级毛片| 三上悠亚av全集在线观看| 十八禁高潮呻吟视频| 久久精品aⅴ一区二区三区四区| 亚洲av成人av| 亚洲一区二区三区欧美精品| e午夜精品久久久久久久| 国产日韩一区二区三区精品不卡| 欧美乱码精品一区二区三区| 一二三四在线观看免费中文在| 下体分泌物呈黄色| 手机成人av网站| 九色亚洲精品在线播放| 一边摸一边做爽爽视频免费| 国产精品久久久久久精品古装| 精品少妇一区二区三区视频日本电影| 热re99久久精品国产66热6| 99国产极品粉嫩在线观看| 啦啦啦免费观看视频1| 日韩熟女老妇一区二区性免费视频| 国产av又大| 我的亚洲天堂| 色综合婷婷激情| 少妇 在线观看| 亚洲精品在线观看二区| 日本精品一区二区三区蜜桃| 精品免费久久久久久久清纯 | 精品国产乱子伦一区二区三区| 99国产精品99久久久久| 久久精品亚洲精品国产色婷小说| 悠悠久久av| 亚洲美女黄片视频| 亚洲精品一二三| 国产伦人伦偷精品视频| 91国产中文字幕| 久久草成人影院| 脱女人内裤的视频| 水蜜桃什么品种好| 婷婷精品国产亚洲av在线 | 免费不卡黄色视频| 国产aⅴ精品一区二区三区波| 香蕉久久夜色| 亚洲精品中文字幕在线视频| 国产精品久久久av美女十八| 一本一本久久a久久精品综合妖精| 丝袜美腿诱惑在线| 交换朋友夫妻互换小说| 69精品国产乱码久久久| 国产精品一区二区免费欧美| 无限看片的www在线观看| 女人高潮潮喷娇喘18禁视频| 亚洲精品粉嫩美女一区| 高清毛片免费观看视频网站 | 亚洲国产欧美一区二区综合| 精品国内亚洲2022精品成人 | 丝袜在线中文字幕| 欧美在线黄色| 亚洲精品粉嫩美女一区| 国产欧美日韩一区二区精品| 亚洲欧美日韩高清在线视频| 国产精品一区二区精品视频观看| 成人av一区二区三区在线看| 91麻豆精品激情在线观看国产 | 老汉色av国产亚洲站长工具| 国产精品成人在线| 在线播放国产精品三级| 久久久久国产精品人妻aⅴ院 | 欧美成狂野欧美在线观看| 少妇的丰满在线观看| 在线国产一区二区在线| 丝瓜视频免费看黄片| 伊人久久大香线蕉亚洲五| 精品一区二区三卡| 欧美在线黄色| 18禁国产床啪视频网站| 黄色视频不卡| 国产亚洲精品久久久久5区| 51午夜福利影视在线观看| 日本a在线网址| 午夜福利在线观看吧| 99久久国产精品久久久| 国产片内射在线| 不卡av一区二区三区| 免费黄频网站在线观看国产| 久久久水蜜桃国产精品网| 两性夫妻黄色片| 成人影院久久| 欧美精品亚洲一区二区| 日韩精品免费视频一区二区三区| 亚洲一码二码三码区别大吗| 久久久久国产一级毛片高清牌| 亚洲avbb在线观看| 淫妇啪啪啪对白视频| 亚洲av成人一区二区三| 精品福利观看| 日韩制服丝袜自拍偷拍| 日韩欧美三级三区| 久久99一区二区三区| 一进一出抽搐动态| 久久久久久久久免费视频了| 老司机在亚洲福利影院| 国产在线观看jvid| 亚洲性夜色夜夜综合| 免费在线观看黄色视频的| 精品第一国产精品| av免费在线观看网站| av在线播放免费不卡| 亚洲精品中文字幕一二三四区| 久久香蕉激情| 久久人人爽av亚洲精品天堂| 欧美日韩瑟瑟在线播放| 欧美精品高潮呻吟av久久| 国产精品九九99| 日韩大码丰满熟妇| 国产精品久久电影中文字幕 | 国产精品1区2区在线观看. | 日本黄色日本黄色录像| 91国产中文字幕| 国产高清videossex| 久久精品亚洲熟妇少妇任你| 精品人妻1区二区| 视频区欧美日本亚洲| 久99久视频精品免费| 我的亚洲天堂| 精品国产超薄肉色丝袜足j| 国产不卡av网站在线观看| 中文字幕另类日韩欧美亚洲嫩草| 欧美精品一区二区免费开放| 操出白浆在线播放| 成人18禁在线播放| 变态另类成人亚洲欧美熟女 | 亚洲精品乱久久久久久| 91九色精品人成在线观看| 国产欧美日韩一区二区三区在线| 日韩熟女老妇一区二区性免费视频| 精品一区二区三区四区五区乱码| 99国产精品99久久久久| 亚洲欧美色中文字幕在线| 午夜福利在线观看吧| 两个人看的免费小视频| 麻豆成人av在线观看| 巨乳人妻的诱惑在线观看| 中文字幕人妻丝袜制服| 亚洲第一欧美日韩一区二区三区| 欧美国产精品va在线观看不卡| 国产xxxxx性猛交| 18禁美女被吸乳视频| 中文字幕色久视频| 欧美国产精品va在线观看不卡| 精品人妻1区二区| 国产有黄有色有爽视频| 三级毛片av免费| 9191精品国产免费久久| 国产野战对白在线观看| 亚洲av日韩在线播放| √禁漫天堂资源中文www| 亚洲精品一卡2卡三卡4卡5卡| 欧美黑人欧美精品刺激| 女警被强在线播放| 搡老熟女国产l中国老女人| 丁香欧美五月| 欧美午夜高清在线| 夜夜爽天天搞| 99久久国产精品久久久| 丁香六月欧美| 不卡一级毛片| 日韩三级视频一区二区三区| 免费女性裸体啪啪无遮挡网站| 18禁黄网站禁片午夜丰满| 亚洲情色 制服丝袜| 欧美精品亚洲一区二区| 久久狼人影院| 国产精品久久久av美女十八| 五月开心婷婷网| 中出人妻视频一区二区| 亚洲五月色婷婷综合| 宅男免费午夜| av福利片在线| 午夜福利影视在线免费观看| 丰满的人妻完整版| 男人的好看免费观看在线视频 | 日本五十路高清| 五月开心婷婷网| 精品熟女少妇八av免费久了| 国产99久久九九免费精品| 99久久人妻综合| 午夜福利在线免费观看网站| 首页视频小说图片口味搜索| 99国产精品一区二区三区| 久久精品熟女亚洲av麻豆精品| 久久久久久久久久久久大奶| 国产精品一区二区免费欧美| 丝袜在线中文字幕| videosex国产| 久久久久视频综合| 国产精品99久久99久久久不卡| 在线免费观看的www视频| 国产国语露脸激情在线看| 黑人操中国人逼视频| 人妻丰满熟妇av一区二区三区 | 不卡一级毛片| 在线天堂中文资源库| 女警被强在线播放| 丝袜美足系列| 日韩成人在线观看一区二区三区| 69av精品久久久久久| 一级片'在线观看视频| 国产精品国产高清国产av | 国产亚洲欧美98| 国产精品国产高清国产av | 亚洲人成伊人成综合网2020| 午夜91福利影院| av天堂在线播放| 69av精品久久久久久| 一级毛片女人18水好多| 91成人精品电影| 18在线观看网站| 91在线观看av| 丰满的人妻完整版| bbb黄色大片| 女性被躁到高潮视频| 国产不卡av网站在线观看| av福利片在线| 亚洲 欧美一区二区三区| 免费在线观看影片大全网站| 欧美老熟妇乱子伦牲交| 男人操女人黄网站| 99久久99久久久精品蜜桃| 久久久国产欧美日韩av| 丝袜人妻中文字幕| 精品久久久久久久久久免费视频 | 午夜福利在线免费观看网站| 女人被狂操c到高潮| 国产精品一区二区精品视频观看| 一二三四社区在线视频社区8| 久久这里只有精品19| 亚洲欧美一区二区三区久久| 搡老岳熟女国产| 午夜福利在线观看吧| 久久青草综合色| 涩涩av久久男人的天堂| 亚洲一区二区三区欧美精品| 国产男女内射视频| 国产精品国产高清国产av | 黄片大片在线免费观看| 久久ye,这里只有精品| 一区二区三区激情视频| 国产亚洲精品久久久久5区| 狂野欧美激情性xxxx| www日本在线高清视频| 亚洲成人免费av在线播放| 51午夜福利影视在线观看| 老司机福利观看| 亚洲欧美激情在线| 欧美日韩中文字幕国产精品一区二区三区 | 日韩欧美免费精品| 丰满迷人的少妇在线观看| 国产亚洲精品第一综合不卡| 免费看十八禁软件| 欧美乱色亚洲激情| 中文字幕人妻丝袜一区二区| av网站在线播放免费| 天天躁夜夜躁狠狠躁躁| 久久亚洲真实| 午夜日韩欧美国产| 亚洲视频免费观看视频| 亚洲成a人片在线一区二区| 麻豆成人av在线观看| 午夜福利免费观看在线| 欧美中文综合在线视频| 12—13女人毛片做爰片一| 亚洲九九香蕉| 免费在线观看日本一区| 一区福利在线观看| 无人区码免费观看不卡| 国产精品免费视频内射| 国产免费现黄频在线看| 国产人伦9x9x在线观看| 波多野结衣av一区二区av| 免费看十八禁软件| av电影中文网址| 亚洲av美国av| 人妻久久中文字幕网| 一进一出抽搐gif免费好疼 | 日韩欧美一区二区三区在线观看 | 夫妻午夜视频| av天堂在线播放| 国产黄色免费在线视频| 天天躁日日躁夜夜躁夜夜| 亚洲,欧美精品.| 久久久久国产一级毛片高清牌| 久久久久久久久久久久大奶| 一级毛片精品| 国产蜜桃级精品一区二区三区 | av免费在线观看网站| 91成年电影在线观看| 欧美精品av麻豆av| videos熟女内射| 在线观看免费午夜福利视频| 黄色成人免费大全| 黄色视频,在线免费观看| 人人澡人人妻人| 亚洲国产毛片av蜜桃av| 久久久国产成人免费| 亚洲精品一卡2卡三卡4卡5卡| 黄色视频,在线免费观看| 黄色a级毛片大全视频| 精品人妻在线不人妻| 国产精品亚洲av一区麻豆| 欧美日韩亚洲高清精品| 免费人成视频x8x8入口观看| 精品少妇久久久久久888优播| 精品一区二区三卡| 国产成人一区二区三区免费视频网站| 美女国产高潮福利片在线看| 亚洲av电影在线进入| 免费在线观看日本一区| 亚洲第一青青草原| www.精华液| 丁香欧美五月| 最新在线观看一区二区三区| 韩国精品一区二区三区| 青草久久国产| 亚洲欧美日韩另类电影网站| 久久99一区二区三区| av不卡在线播放| 国产精品自产拍在线观看55亚洲 | 日韩欧美国产一区二区入口| 国产精华一区二区三区| 欧美精品亚洲一区二区| 交换朋友夫妻互换小说| 热99国产精品久久久久久7| 黄色视频,在线免费观看| 91老司机精品| 国产91精品成人一区二区三区| 欧美丝袜亚洲另类 | 免费在线观看黄色视频的| 国产99白浆流出| 国产黄色免费在线视频| 777久久人妻少妇嫩草av网站| 51午夜福利影视在线观看| 亚洲欧美一区二区三区久久| 啦啦啦 在线观看视频| 高潮久久久久久久久久久不卡| 亚洲精品久久成人aⅴ小说| 国产成人一区二区三区免费视频网站| 不卡一级毛片| 色精品久久人妻99蜜桃| 又黄又粗又硬又大视频| 免费高清在线观看日韩| 丰满迷人的少妇在线观看| 国产又色又爽无遮挡免费看| 777米奇影视久久| 好看av亚洲va欧美ⅴa在| 亚洲av美国av| 熟女少妇亚洲综合色aaa.| 亚洲美女黄片视频| 亚洲免费av在线视频| 亚洲少妇的诱惑av| 国产成人免费观看mmmm| 精品久久久久久,| 操美女的视频在线观看| videos熟女内射| 十分钟在线观看高清视频www| 免费看a级黄色片| 热99re8久久精品国产| 成人国产一区最新在线观看| 热99re8久久精品国产| 韩国精品一区二区三区| 日本一区二区免费在线视频| 中文字幕制服av| 黄色片一级片一级黄色片| www.自偷自拍.com| 久久久精品国产亚洲av高清涩受| 99久久人妻综合| 很黄的视频免费| 在线视频色国产色| 久久久国产一区二区| 国产主播在线观看一区二区|