• <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è)計
    h视频一区二区三区| 国产午夜精品久久久久久| 黄色a级毛片大全视频| 91国产中文字幕| kizo精华| 亚洲国产欧美日韩在线播放| 国产成人免费无遮挡视频| 亚洲成a人片在线一区二区| 亚洲精品粉嫩美女一区| 激情视频va一区二区三区| 一级毛片电影观看| 午夜福利在线观看吧| 亚洲视频免费观看视频| 日韩一区二区三区影片| 看免费av毛片| 久久久精品国产亚洲av高清涩受| 久久精品人人爽人人爽视色| 日韩制服丝袜自拍偷拍| 大片免费播放器 马上看| 精品熟女少妇八av免费久了| 日本精品一区二区三区蜜桃| 亚洲三区欧美一区| 9色porny在线观看| 日韩有码中文字幕| 亚洲av日韩在线播放| av欧美777| 天堂俺去俺来也www色官网| 国产99久久九九免费精品| 99久久国产精品久久久| 免费高清在线观看日韩| 新久久久久国产一级毛片| 日韩免费高清中文字幕av| 最新在线观看一区二区三区| 视频区欧美日本亚洲| h视频一区二区三区| 又紧又爽又黄一区二区| 国产精品一区二区在线不卡| 成年人黄色毛片网站| videosex国产| av欧美777| 丝袜在线中文字幕| 91九色精品人成在线观看| 2018国产大陆天天弄谢| 欧美精品人与动牲交sv欧美| 国产av精品麻豆| 国产精品影院久久| 咕卡用的链子| e午夜精品久久久久久久| 中文字幕人妻丝袜制服| 亚洲av美国av| 自线自在国产av| 女性被躁到高潮视频| svipshipincom国产片| av超薄肉色丝袜交足视频| 欧美黑人精品巨大| 亚洲五月色婷婷综合| 一级毛片精品| 日韩大片免费观看网站| 最近最新免费中文字幕在线| 成人免费观看视频高清| 99re在线观看精品视频| 日韩中文字幕欧美一区二区| 免费一级毛片在线播放高清视频 | 欧美黑人欧美精品刺激| 乱人伦中国视频| 黄色成人免费大全| 91九色精品人成在线观看| 男女高潮啪啪啪动态图| 中文字幕人妻熟女乱码| 狠狠狠狠99中文字幕| 18禁国产床啪视频网站| 婷婷成人精品国产| 丁香六月欧美| 精品午夜福利视频在线观看一区 | 日韩大片免费观看网站| 亚洲欧美一区二区三区黑人| av欧美777| 免费在线观看影片大全网站| 亚洲专区字幕在线| 我的亚洲天堂| 精品午夜福利视频在线观看一区 | 亚洲国产av影院在线观看| 最近最新免费中文字幕在线| 日韩视频一区二区在线观看| 一进一出好大好爽视频| av线在线观看网站| 9热在线视频观看99| 我的亚洲天堂| 国产一区二区在线观看av| 久久久久久久精品吃奶| 国产欧美日韩精品亚洲av| 日本欧美视频一区| 久久久久久久大尺度免费视频| 成年版毛片免费区| 王馨瑶露胸无遮挡在线观看| 国产高清videossex| 午夜老司机福利片| 亚洲国产精品一区二区三区在线| 日韩成人在线观看一区二区三区| 国产成人免费无遮挡视频| 亚洲精华国产精华精| 性高湖久久久久久久久免费观看| 1024香蕉在线观看| 啦啦啦免费观看视频1| 伊人久久大香线蕉亚洲五| 一级a爱视频在线免费观看| 日日摸夜夜添夜夜添小说| 大码成人一级视频| 成人手机av| 中文字幕人妻丝袜制服| 日本精品一区二区三区蜜桃| 搡老岳熟女国产| 国产高清视频在线播放一区| 亚洲一区中文字幕在线| 久久精品国产综合久久久| 久久青草综合色| 黑人欧美特级aaaaaa片| 自拍欧美九色日韩亚洲蝌蚪91| 99热国产这里只有精品6| 电影成人av| 成人免费观看视频高清| 日韩欧美一区视频在线观看| 自线自在国产av| 国产视频一区二区在线看| av不卡在线播放| 色播在线永久视频| 怎么达到女性高潮| 成人18禁在线播放| av视频免费观看在线观看| 另类精品久久| 美国免费a级毛片| 黄色丝袜av网址大全| 午夜老司机福利片| a在线观看视频网站| 亚洲av第一区精品v没综合| 五月开心婷婷网| 免费在线观看完整版高清| h视频一区二区三区| 国产精品麻豆人妻色哟哟久久| 老熟女久久久| 少妇粗大呻吟视频| 亚洲熟女毛片儿| 夜夜骑夜夜射夜夜干| 啦啦啦中文免费视频观看日本| 成人国语在线视频| 欧美人与性动交α欧美精品济南到| 亚洲精品自拍成人| 久久天堂一区二区三区四区| 涩涩av久久男人的天堂| 三上悠亚av全集在线观看| 欧美久久黑人一区二区| 99热网站在线观看| 亚洲欧美一区二区三区黑人| 欧美激情极品国产一区二区三区| av在线播放免费不卡| 国产深夜福利视频在线观看| 成人特级黄色片久久久久久久 | 国产无遮挡羞羞视频在线观看| 国产精品99久久99久久久不卡| 色老头精品视频在线观看| 欧美成人午夜精品| 一级a爱视频在线免费观看| 免费av中文字幕在线| 一区在线观看完整版| 亚洲av日韩在线播放| 五月开心婷婷网| 宅男免费午夜| 午夜福利乱码中文字幕| 欧美av亚洲av综合av国产av| 精品人妻1区二区| 韩国精品一区二区三区| 欧美 日韩 精品 国产| 美女视频免费永久观看网站| 建设人人有责人人尽责人人享有的| 一二三四在线观看免费中文在| 亚洲专区字幕在线| 国精品久久久久久国模美| 久久影院123| 多毛熟女@视频| 成人18禁在线播放| 亚洲精品国产色婷婷电影| 成人18禁在线播放| 多毛熟女@视频| 法律面前人人平等表现在哪些方面| 精品午夜福利视频在线观看一区 | av网站在线播放免费| 菩萨蛮人人尽说江南好唐韦庄| aaaaa片日本免费| 国产成人免费观看mmmm| 精品久久蜜臀av无| 精品午夜福利视频在线观看一区 | 久久人妻av系列| 久久人妻av系列| 9色porny在线观看| 久久午夜综合久久蜜桃| 成人18禁在线播放| 精品午夜福利视频在线观看一区 | 久久久国产一区二区| 欧美一级毛片孕妇| 亚洲第一欧美日韩一区二区三区 | 丰满迷人的少妇在线观看| 黑人猛操日本美女一级片| 欧美精品av麻豆av| 欧美久久黑人一区二区| xxxhd国产人妻xxx| 色播在线永久视频| 亚洲专区字幕在线| 在线 av 中文字幕| 国产老妇伦熟女老妇高清| 超碰成人久久| 久久精品亚洲精品国产色婷小说| 黑人操中国人逼视频| 欧美日韩国产mv在线观看视频| 99久久99久久久精品蜜桃| 一本久久精品| 国产精品熟女久久久久浪| 国产老妇伦熟女老妇高清| 在线十欧美十亚洲十日本专区| 一本久久精品| 国产欧美日韩综合在线一区二区| 两性夫妻黄色片| 制服诱惑二区| 最黄视频免费看| 精品国产一区二区久久| 国产精品亚洲av一区麻豆| 婷婷丁香在线五月| 黑人猛操日本美女一级片| 黑人猛操日本美女一级片| 久久久久久亚洲精品国产蜜桃av| 中亚洲国语对白在线视频| a级毛片在线看网站| 欧美日韩视频精品一区| 久久久精品国产亚洲av高清涩受| 最新在线观看一区二区三区| 可以免费在线观看a视频的电影网站| 日本vs欧美在线观看视频| 日本wwww免费看| 日韩中文字幕视频在线看片| svipshipincom国产片| 黄色a级毛片大全视频| 少妇猛男粗大的猛烈进出视频| 天天躁日日躁夜夜躁夜夜| 亚洲午夜精品一区,二区,三区| 亚洲三区欧美一区| 老汉色∧v一级毛片| 欧美亚洲日本最大视频资源| 嫩草影视91久久| 国产日韩欧美在线精品| 1024香蕉在线观看| 女人久久www免费人成看片| 亚洲九九香蕉| 天天操日日干夜夜撸| 久久精品91无色码中文字幕| 嫁个100分男人电影在线观看| 欧美日韩视频精品一区| 午夜老司机福利片| 一级片'在线观看视频| 国产一区二区 视频在线| 免费一级毛片在线播放高清视频 | 国产亚洲午夜精品一区二区久久| 国产精品一区二区免费欧美| netflix在线观看网站| 曰老女人黄片| 建设人人有责人人尽责人人享有的| 黄色 视频免费看| 涩涩av久久男人的天堂| 在线观看免费视频网站a站| 亚洲成av片中文字幕在线观看| 黑丝袜美女国产一区| 久久免费观看电影| 男人操女人黄网站| 菩萨蛮人人尽说江南好唐韦庄| 一本—道久久a久久精品蜜桃钙片| 欧美精品av麻豆av| 咕卡用的链子| 侵犯人妻中文字幕一二三四区| 每晚都被弄得嗷嗷叫到高潮| 新久久久久国产一级毛片| 久久久欧美国产精品| 日韩视频一区二区在线观看| a级毛片在线看网站| 午夜老司机福利片| 2018国产大陆天天弄谢| 99久久99久久久精品蜜桃| 麻豆成人av在线观看| 十八禁高潮呻吟视频| 欧美老熟妇乱子伦牲交| 久久国产精品大桥未久av| 麻豆av在线久日| 亚洲国产av影院在线观看| 免费久久久久久久精品成人欧美视频| 国产aⅴ精品一区二区三区波| 亚洲成人手机| 乱人伦中国视频| 成在线人永久免费视频| 亚洲成a人片在线一区二区| 亚洲自偷自拍图片 自拍| 人妻 亚洲 视频| 在线看a的网站| 热99re8久久精品国产| 成人三级做爰电影| 欧美久久黑人一区二区| 一本大道久久a久久精品| 精品国产超薄肉色丝袜足j| 在线观看一区二区三区激情| 19禁男女啪啪无遮挡网站| 亚洲国产欧美一区二区综合| 极品人妻少妇av视频| 国产精品欧美亚洲77777| 精品国内亚洲2022精品成人 | 亚洲性夜色夜夜综合| 18在线观看网站| 狂野欧美激情性xxxx| 啦啦啦中文免费视频观看日本| 国产精品99久久99久久久不卡| 亚洲av成人不卡在线观看播放网| 欧美黄色片欧美黄色片| 精品免费久久久久久久清纯 | 黄频高清免费视频| 一区在线观看完整版| 大片电影免费在线观看免费| 色视频在线一区二区三区| av不卡在线播放| 亚洲全国av大片| 一区在线观看完整版| 在线av久久热| 啦啦啦中文免费视频观看日本| 国产精品98久久久久久宅男小说| 自线自在国产av| 午夜福利,免费看| 亚洲国产欧美一区二区综合| 国产97色在线日韩免费| 日本av手机在线免费观看| 国产老妇伦熟女老妇高清| 女人被躁到高潮嗷嗷叫费观| 国产免费av片在线观看野外av| 日韩精品免费视频一区二区三区| kizo精华| 亚洲七黄色美女视频| 日本av免费视频播放| 久久久久国内视频| 久久香蕉激情| 国产精品偷伦视频观看了| 日韩中文字幕视频在线看片| 啦啦啦视频在线资源免费观看| 国产精品一区二区精品视频观看| 国产在线一区二区三区精| 亚洲成av片中文字幕在线观看| 交换朋友夫妻互换小说| 日韩大片免费观看网站| 亚洲色图综合在线观看| 国产精品1区2区在线观看. | 真人做人爱边吃奶动态| 国产野战对白在线观看| 69av精品久久久久久 | 超碰97精品在线观看| 在线观看免费视频网站a站| 黄频高清免费视频| 夜夜夜夜夜久久久久| 丁香六月欧美| 午夜视频精品福利| 91精品三级在线观看| 在线观看免费日韩欧美大片| 中文字幕av电影在线播放| 国产精品影院久久| 久久人妻福利社区极品人妻图片| 后天国语完整版免费观看| 法律面前人人平等表现在哪些方面| 一本久久精品| 老熟女久久久| 大香蕉久久网| 麻豆国产av国片精品| 久久午夜综合久久蜜桃| 午夜激情av网站| 色在线成人网| 午夜福利,免费看| 色精品久久人妻99蜜桃| 欧美成人免费av一区二区三区 | 男女免费视频国产| 99热网站在线观看| 啦啦啦视频在线资源免费观看| 精品少妇一区二区三区视频日本电影| 亚洲第一av免费看| 国产亚洲精品久久久久5区| 99久久精品国产亚洲精品| 一区二区三区国产精品乱码| 纯流量卡能插随身wifi吗| 中文字幕人妻丝袜一区二区| 一边摸一边抽搐一进一小说 | 亚洲人成电影免费在线| 高潮久久久久久久久久久不卡| 女同久久另类99精品国产91| avwww免费| 亚洲av欧美aⅴ国产| 国产男女内射视频| 欧美日韩av久久| 午夜福利乱码中文字幕| 最近最新中文字幕大全免费视频| 老司机午夜福利在线观看视频 | 亚洲色图av天堂| 国产淫语在线视频| 精品福利永久在线观看| www.精华液| 少妇 在线观看| 国产精品久久久人人做人人爽| 99国产精品一区二区蜜桃av | 香蕉国产在线看| 丰满迷人的少妇在线观看| 久久国产精品影院| 热99久久久久精品小说推荐| 中国美女看黄片| 国产高清国产精品国产三级| 日韩一区二区三区影片| 亚洲成人国产一区在线观看| 亚洲av国产av综合av卡| 色尼玛亚洲综合影院| 国产精品一区二区精品视频观看| 久久亚洲精品不卡| 亚洲av国产av综合av卡| av天堂久久9| 不卡一级毛片| 国产在线免费精品| 99久久国产精品久久久| 日日摸夜夜添夜夜添小说| 飞空精品影院首页| 女性生殖器流出的白浆| 国产精品秋霞免费鲁丝片| 99九九在线精品视频| 色在线成人网| 51午夜福利影视在线观看| 在线av久久热| 免费在线观看视频国产中文字幕亚洲| 免费少妇av软件| 纵有疾风起免费观看全集完整版| 国产亚洲午夜精品一区二区久久| 国产精品1区2区在线观看. | 亚洲精品久久午夜乱码| 国产欧美日韩一区二区三区在线| 亚洲av日韩精品久久久久久密| 色婷婷av一区二区三区视频| 香蕉久久夜色| 99re在线观看精品视频| 999久久久精品免费观看国产| 国产精品二区激情视频| 在线观看免费午夜福利视频| 咕卡用的链子| 国产又色又爽无遮挡免费看| 一区二区三区精品91| 99riav亚洲国产免费| 熟女少妇亚洲综合色aaa.| 亚洲欧美精品综合一区二区三区| 国产一卡二卡三卡精品| 亚洲一卡2卡3卡4卡5卡精品中文| 国产成人精品久久二区二区免费| 国产男靠女视频免费网站| 天天操日日干夜夜撸| 亚洲情色 制服丝袜| 天天躁夜夜躁狠狠躁躁| 免费高清在线观看日韩| 桃花免费在线播放| 精品福利观看| 日本精品一区二区三区蜜桃| 18禁黄网站禁片午夜丰满| 成人永久免费在线观看视频 | 亚洲五月婷婷丁香| 99精品在免费线老司机午夜| 精品亚洲成国产av| av超薄肉色丝袜交足视频| 看免费av毛片| 精品人妻1区二区| 69av精品久久久久久 | 丝袜喷水一区| 亚洲精华国产精华精| 久久青草综合色| 国产日韩欧美亚洲二区| 久久久久久久精品吃奶| 我的亚洲天堂| 无遮挡黄片免费观看| 中国美女看黄片| 中文字幕av电影在线播放| 91成人精品电影| 中亚洲国语对白在线视频| 亚洲 欧美一区二区三区| 色在线成人网| 国产黄色免费在线视频| 日日摸夜夜添夜夜添小说| 天天影视国产精品| 91国产中文字幕| 9191精品国产免费久久| h视频一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 99九九在线精品视频| 久热这里只有精品99| 操出白浆在线播放| 久久久久网色| 一二三四社区在线视频社区8| e午夜精品久久久久久久| 黄色视频在线播放观看不卡| 十八禁网站免费在线| 可以免费在线观看a视频的电影网站| 少妇 在线观看| 亚洲中文av在线| 色综合欧美亚洲国产小说| 麻豆乱淫一区二区| 日韩视频一区二区在线观看| 一级毛片电影观看| 久久人人97超碰香蕉20202| 亚洲七黄色美女视频| 9色porny在线观看| 少妇精品久久久久久久| 青草久久国产| 国产精品久久久久成人av| a级片在线免费高清观看视频| 视频区图区小说| 国产亚洲av高清不卡| 亚洲av电影在线进入| 国产成人啪精品午夜网站| 纯流量卡能插随身wifi吗| 欧美精品人与动牲交sv欧美| 国产免费视频播放在线视频| 久久精品熟女亚洲av麻豆精品| 一区二区av电影网| 在线亚洲精品国产二区图片欧美| 三级毛片av免费| 男女之事视频高清在线观看| 午夜91福利影院| 18禁国产床啪视频网站| 国产精品亚洲av一区麻豆| 不卡一级毛片| 欧美精品一区二区大全| 在线观看免费高清a一片| 另类亚洲欧美激情| e午夜精品久久久久久久| 日本黄色视频三级网站网址 | 亚洲av美国av| 高清av免费在线| 99久久人妻综合| 国产黄频视频在线观看| 在线十欧美十亚洲十日本专区| 日韩三级视频一区二区三区| 午夜久久久在线观看| 十八禁高潮呻吟视频| 成年人黄色毛片网站| 女性被躁到高潮视频| 亚洲av日韩精品久久久久久密| 国产片内射在线| 啦啦啦免费观看视频1| 一进一出抽搐动态| 一边摸一边做爽爽视频免费| 性色av乱码一区二区三区2| 日日爽夜夜爽网站| kizo精华| 欧美成人免费av一区二区三区 | 国产人伦9x9x在线观看| 少妇精品久久久久久久| 精品国产乱码久久久久久小说| 亚洲五月色婷婷综合| 亚洲七黄色美女视频| 大片免费播放器 马上看| 人人妻人人爽人人添夜夜欢视频| www.精华液| tube8黄色片| 国产精品秋霞免费鲁丝片| 不卡一级毛片| 亚洲av国产av综合av卡| 不卡av一区二区三区| 国产精品一区二区在线观看99| 国产成人av激情在线播放| 国产精品久久久人人做人人爽| 国产精品亚洲av一区麻豆| 男人舔女人的私密视频| 成人手机av| 黑人巨大精品欧美一区二区mp4| 亚洲欧美一区二区三区黑人| 精品国产国语对白av| 久久久水蜜桃国产精品网| 国产aⅴ精品一区二区三区波| 男女无遮挡免费网站观看| 桃花免费在线播放| 99久久人妻综合| 亚洲人成电影免费在线| 在线观看www视频免费| 男女边摸边吃奶| 正在播放国产对白刺激| 亚洲欧美色中文字幕在线| 极品人妻少妇av视频| 在线亚洲精品国产二区图片欧美| 在线观看免费日韩欧美大片| 欧美精品亚洲一区二区| 欧美在线一区亚洲| 亚洲自偷自拍图片 自拍| 久久国产精品大桥未久av| 久久 成人 亚洲| 99在线人妻在线中文字幕 | 精品久久蜜臀av无| 亚洲美女黄片视频| 桃花免费在线播放| 国产深夜福利视频在线观看| 久久香蕉激情| 美女视频免费永久观看网站| 国产黄频视频在线观看| 国产精品亚洲一级av第二区| 久久免费观看电影| 国产精品99久久99久久久不卡| 久久精品国产亚洲av高清一级| 香蕉久久夜色| 少妇裸体淫交视频免费看高清 | 久久久国产成人免费| 婷婷成人精品国产| 久久香蕉激情| 午夜视频精品福利| 久久性视频一级片| 亚洲专区字幕在线| 19禁男女啪啪无遮挡网站| 成人18禁高潮啪啪吃奶动态图|