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

    基于六面體網(wǎng)格的向量式有限元分析及應(yīng)用

    2018-09-06 10:00:08,
    關(guān)鍵詞:六面體質(zhì)點(diǎn)內(nèi)力

    , ,

    (1.浙江省建筑設(shè)計(jì)研究院,杭州 310006;2.浙江大學(xué) 空間結(jié)構(gòu)研究中心,杭州 310058)

    1 引 言

    實(shí)體結(jié)構(gòu)是最為通用的結(jié)構(gòu)形式,具有廣泛的工程應(yīng)用。實(shí)體單元與其他單元通過(guò)在表面連接,以模擬任意外形結(jié)構(gòu)在外界作用下的響應(yīng)。根據(jù)單元形式的不同,實(shí)體結(jié)構(gòu)單元主要有四面體和六面體兩類(lèi),前者適用性強(qiáng),而后者計(jì)算精度較高[1,2]。Pawlak等[3]推導(dǎo)了含有旋轉(zhuǎn)自由度的四面體單元形式,該單元對(duì)于具有反對(duì)稱(chēng)應(yīng)力和旋轉(zhuǎn)問(wèn)題的結(jié)構(gòu)計(jì)算具有較高精度。Tayloy等[4]提出位移不協(xié)調(diào)的非協(xié)調(diào)六面體單元形式,并通過(guò)等參變換實(shí)現(xiàn)不規(guī)則網(wǎng)格的數(shù)值模擬。Cao等[5]基于Hu-Washizu變分原理提出多變量形式的六面體單元,提高了單元抗網(wǎng)格畸變能力。李宏光等[6]引入體積坐標(biāo)方法,利用其與整體坐標(biāo)的線性關(guān)系,構(gòu)造出對(duì)網(wǎng)格畸變不敏感的六面體單元,以應(yīng)用于復(fù)雜結(jié)構(gòu)形狀。對(duì)于大位移大轉(zhuǎn)動(dòng)問(wèn)題,六面體單元通過(guò)較少的單元網(wǎng)格劃分即可獲得較為精確的結(jié)果。

    向量式有限元VFIFE(vector form intrinsic finite element)[7-9]是基于質(zhì)點(diǎn)運(yùn)動(dòng)和向量分析的新型數(shù)值方法。該法基于牛頓經(jīng)典運(yùn)動(dòng)理論,以質(zhì)點(diǎn)系的運(yùn)動(dòng)來(lái)描述結(jié)構(gòu)體系的行為,各質(zhì)點(diǎn)間的運(yùn)動(dòng)相互獨(dú)立并通過(guò)逐步計(jì)算求得整體結(jié)構(gòu)響應(yīng),而且求解時(shí)不存在迭代和剛度矩陣問(wèn)題。在變形坐標(biāo)系下,剛體位移可能造成的誤差采用逆向運(yùn)動(dòng)方式來(lái)規(guī)避。故該方法在結(jié)構(gòu)大位移大轉(zhuǎn)動(dòng)、機(jī)構(gòu)運(yùn)動(dòng)及不連續(xù)行為的數(shù)值分析領(lǐng)域具有較好的應(yīng)用[10-14]。

    本文給出了八節(jié)點(diǎn)六面體等參實(shí)體單元的向量式基本理論,采用參考面的逆向運(yùn)動(dòng)求解節(jié)點(diǎn)純變形,通過(guò)單元形函的虛功方程計(jì)算節(jié)點(diǎn)內(nèi)力,針對(duì)坐標(biāo)模式和內(nèi)力積分模式等特殊問(wèn)題提出了有效的處理方案。編制數(shù)值計(jì)算程序,并通過(guò)結(jié)構(gòu)算例驗(yàn)證理論和程序的正確性和實(shí)用性。

    2 六面體網(wǎng)格的向量式理論

    向量式有限元是由質(zhì)點(diǎn)和單元組成,結(jié)構(gòu)的質(zhì)量?jī)H分布于質(zhì)點(diǎn),本文單元為八節(jié)點(diǎn)六面體等參實(shí)體單元。求解步驟為質(zhì)點(diǎn)全位移計(jì)算、單元節(jié)點(diǎn)純變形計(jì)算和單元節(jié)點(diǎn)內(nèi)力計(jì)算,通過(guò)循環(huán)來(lái)達(dá)到逐點(diǎn)逐步求解。分析流程如圖1所示。

    圖1 分析流程圖

    Fig.1 Analysis flow chart

    2.1 質(zhì)點(diǎn)全位移

    質(zhì)點(diǎn)的運(yùn)動(dòng)滿(mǎn)足微分方程:

    (1)

    式中m為質(zhì)量,α為阻尼參數(shù),x為位移向量,F(xiàn)=fext+fint為合力向量,上標(biāo)ext和int分別表示外力和內(nèi)力。

    時(shí)間步內(nèi)的質(zhì)點(diǎn)全位移采用顯式中央差分方法數(shù)值求解運(yùn)動(dòng)方程獲得,為質(zhì)點(diǎn)新坐標(biāo)與初始坐標(biāo)的差值。

    2.2 單元節(jié)點(diǎn)純變形

    (1) 坐標(biāo)模式的轉(zhuǎn)換

    通過(guò)參考平面的逆向運(yùn)動(dòng),在質(zhì)點(diǎn)全位移中去除剛體位移,以獲得節(jié)點(diǎn)純變形。節(jié)點(diǎn)實(shí)際坐標(biāo)所組成六面體的每個(gè)面均是空間曲面四邊形,采用投影方式將每個(gè)面均投影為平面四邊形,組合后為投影六面體(即坐標(biāo)模式I),詳見(jiàn)第3.1節(jié)。

    六面體單元面i(i=1~6)的節(jié)點(diǎn)j(j=1,2,3,4)的實(shí)際坐標(biāo)向量記為xj,則平均法向量ni為

    (2)

    ni x(x-xC i)+ni y(y-yC i)+ni z(z-zC i)=0

    (3)

    (4)

    圖2 六面體單元的投影位置

    Fig.2 Projection position of the hexahedral element

    (2) 剛體運(yùn)動(dòng)的估算

    獲得投影六面體后,采用參考平面估算剛體運(yùn)動(dòng)(平移和轉(zhuǎn)動(dòng)),取投影六面體單元的abcd面為參考平面,如圖3所示。

    參考平面abcd在t0~t時(shí)段的空間運(yùn)動(dòng)包括平移、轉(zhuǎn)動(dòng)和純變形[15]。

    平移取參考面節(jié)點(diǎn)a的全位移ua,則消去平移后的質(zhì)點(diǎn)位移為

    (5)

    利用參考面估算質(zhì)點(diǎn)的逆向面外轉(zhuǎn)動(dòng)位移

    (6)

    質(zhì)點(diǎn)i(i=a~h)的逆向面內(nèi)轉(zhuǎn)動(dòng)位移為

    (7)

    質(zhì)點(diǎn)純變形是通過(guò)在質(zhì)點(diǎn)全位移中消去平移和轉(zhuǎn)動(dòng)(面外和面內(nèi))后獲得,

    (8)

    2.3 單元節(jié)點(diǎn)內(nèi)力

    (1)傳統(tǒng)形函的引入

    圖3 投影六面體單元及其參考平面

    Fig.3 Projective hexahedral element and its reference plane

    將整體坐標(biāo)系轉(zhuǎn)換至變形坐標(biāo)系,即

    (9)

    節(jié)點(diǎn)的純變形向量為

    (10)

    節(jié)點(diǎn)的純變形組合向量為

    (11)

    (12)

    (13)

    采用雅可比矩陣J將變形坐標(biāo)系轉(zhuǎn)換至局部坐標(biāo)系,參考文獻(xiàn)[16]。

    (14)

    (15)

    (2) 虛功方程的求解

    單元變形滿(mǎn)足虛功方程

    (16)

    式中fi為整體坐標(biāo)系下的節(jié)點(diǎn)內(nèi)力。

    由應(yīng)力應(yīng)變描述的變形虛功為

    (17)

    單元節(jié)點(diǎn)的內(nèi)力向量為

    (18)

    (19)

    將虛擬狀態(tài)時(shí)的節(jié)點(diǎn)內(nèi)力分量轉(zhuǎn)換至實(shí)際狀態(tài),并經(jīng)正向運(yùn)動(dòng),求得實(shí)際位置時(shí)的節(jié)點(diǎn)內(nèi)力:

    (20)

    單元作用于質(zhì)點(diǎn)的內(nèi)力finti通過(guò)節(jié)點(diǎn)內(nèi)力fi反向作用并集成獲得。

    3 關(guān)鍵問(wèn)題的處理

    3.1 坐標(biāo)模式

    (1) 節(jié)點(diǎn)坐標(biāo)模式

    單元節(jié)點(diǎn)坐標(biāo)(坐標(biāo)模式I)用于求解時(shí)間步內(nèi)的節(jié)點(diǎn)純變形和內(nèi)力,因而單元各面上節(jié)點(diǎn)均需處于同一平面。

    (2) 質(zhì)點(diǎn)坐標(biāo)模式

    3.2 內(nèi)力積分

    六面體單元內(nèi)部的應(yīng)力應(yīng)變隨位置變化,求解虛功方程時(shí)需采用數(shù)值積分計(jì)算。采用三維高斯積分方案[1]時(shí),單積分點(diǎn)和8積分點(diǎn)的相關(guān)系數(shù)列入表1,本文采用后者求解式(18),可獲得較好計(jì)算結(jié)果。

    表1 數(shù)值積分方案

    Tab.1 Numerical integration programmes

    積分點(diǎn)數(shù)焊iηiζiHijm10008.08±1/3±1/3±1/31.0,1.0,1.0,1.01.0,1.0,1.0,1.0

    (21)

    (22)

    式中k為每個(gè)坐標(biāo)方向的積分點(diǎn)數(shù)。

    4 工程結(jié)構(gòu)算例分析

    基于向量式有限元六面體等參元理論,在Matlab平臺(tái)上編制分析計(jì)算程序,實(shí)現(xiàn)向量式有限元實(shí)體結(jié)構(gòu)的計(jì)算分析。

    4.1 懸臂梁結(jié)構(gòu)靜動(dòng)力分析

    (1) 靜力分析

    懸臂梁頂面施加豎向均布靜力作用q=2.5×104kN/m2,采用斜坡-平臺(tái)方式緩慢加載,并考慮阻尼效應(yīng)來(lái)獲得收斂值。細(xì)長(zhǎng)梁自由端的撓度理論解為

    w(x)=(x2-4Lx+6L2)(qH)x2/(24EI)

    圖5為懸臂梁自由端頂部節(jié)點(diǎn)豎向位移的計(jì)算收斂過(guò)程,圖6為沿懸臂梁頂部中線各節(jié)點(diǎn)的豎向位移和von Mises應(yīng)力結(jié)果??梢钥闯觯Q向位移和von Mises應(yīng)力結(jié)果與細(xì)長(zhǎng)梁理論解及ABAQUS結(jié)果均基本一致。本文和ABAQUS的自由端節(jié)點(diǎn)豎向位移分別為0.0115 m和0.0119 m,相對(duì)理論解0.0117 m的誤差分別為-1.71%和1.71%。本文和ABAQUS的自由端節(jié)點(diǎn)von Mises應(yīng)力分別為38.06 MPa和36.4 MPa,差異約為4.36%,懸臂端節(jié)點(diǎn)由于應(yīng)力集中,結(jié)果差異相對(duì)較大,約為9.8%??梢?jiàn)本文方法在實(shí)體工程結(jié)構(gòu)靜力分析中具有較好準(zhǔn)確性。

    (2) 動(dòng)力分析

    在懸臂梁頂面施加瞬時(shí)動(dòng)力均布荷載q=2.0×104kN/m2,不考慮阻尼效應(yīng)。取t=0.01 s進(jìn)行分析,圖7和圖8分別為自由端頂部節(jié)點(diǎn)豎向位移和von Mises應(yīng)力隨時(shí)間的變化情況。

    結(jié)果表明,本文位移和應(yīng)力結(jié)果均與ABAQUS結(jié)果吻合較好,驗(yàn)證了本文方法在實(shí)體工程結(jié)構(gòu)動(dòng)力分析中的準(zhǔn)確性。圖8中兩者應(yīng)力結(jié)果的總體變化趨勢(shì)基本一致,由于端部應(yīng)力集中,局部振蕩波形有所差異。

    圖4 懸臂梁模型

    Fig.4 Cantilever-beam model

    圖5 自由端豎向位移收斂過(guò)程

    Fig.5 Convergence process of free end vertical displacement

    圖6 頂部中線位置結(jié)果

    Fig.6 Results of the top center-line

    4.2 環(huán)形套筒大變形分析

    取t=0.02 s進(jìn)行分析。圖10為環(huán)形套筒左側(cè)中心節(jié)點(diǎn)的水平位移隨時(shí)間的變化情況??梢钥闯觯疚姆椒色@得套筒大變形過(guò)程中的位移變化情況,且與ABAQUS結(jié)果基本一致。在同時(shí)采用顯式動(dòng)力分析時(shí),本文方法由于不存在剛度矩陣奇異和迭代不收斂問(wèn)題,相比ABAQUS具有更高的計(jì)算效率。

    典型的環(huán)形套筒變形和von Mises應(yīng)力云圖情況如圖11所示??梢钥闯?,套筒起初表現(xiàn)為壓扁變形(t=0.004 s),t=0.008 s時(shí),壓扁變形達(dá)到最大值;接著套筒開(kāi)始回彈,在t=0.012 s時(shí)回彈至最小變形;然后,套筒壓扁變形再次變大(t=0.015 s);依此循環(huán)往復(fù)。本文所得結(jié)構(gòu)變形和應(yīng)力分布情況均與ABAQUS計(jì)算結(jié)果吻合。

    圖7 自由端豎向位移-時(shí)間曲線

    Fig.7 Vertical displacement-time curves of the free end

    圖8 自由端von Mises應(yīng)力-時(shí)間曲線

    Fig.8 Von Mises stress-time curves of the free end

    圖9 環(huán)形套筒模型

    Fig.9 Annular sleeve model

    圖10 左側(cè)中心節(jié)點(diǎn)的位移-時(shí)間曲線

    Fig.10 Displacement-time curves of the left central node

    圖11 環(huán)形套筒壓扁變形圖

    Fig.11 Flattened deformation diagrams of annular sleeve

    5 結(jié) 論

    基于向量式有限元,給出了八節(jié)點(diǎn)六面體等參實(shí)體單元的基本公式,通過(guò)投影方式將空間曲面六面體轉(zhuǎn)換為投影六面體,采用參考面的逆向運(yùn)動(dòng)求解節(jié)點(diǎn)純變形,通過(guò)單元形函的虛功方程計(jì)算節(jié)點(diǎn)內(nèi)力。針對(duì)坐標(biāo)模式和內(nèi)力積分模式等關(guān)鍵問(wèn)題提出了有效的處理方案。

    編寫(xiě)了六面體實(shí)體單元的數(shù)值計(jì)算程序,并針對(duì)工程結(jié)構(gòu)的懸臂梁靜動(dòng)力和環(huán)形套筒大變形算例進(jìn)行分析,驗(yàn)證了本文理論和程序的正確性和實(shí)用性。

    向量式有限元適用于分析工程結(jié)構(gòu)大變形大位移行為的動(dòng)力分析,求解時(shí)不存在迭代和剛度矩陣問(wèn)題,在結(jié)構(gòu)動(dòng)力分析領(lǐng)域具有更高的效率和準(zhǔn)確性。

    猜你喜歡
    六面體質(zhì)點(diǎn)內(nèi)力
    巧用“搬運(yùn)法”解決連續(xù)質(zhì)點(diǎn)模型的做功問(wèn)題
    一個(gè)領(lǐng)導(dǎo)人的“六面體”
    孩子的生命內(nèi)力需要家長(zhǎng)去激發(fā)
    逆作法孔口邊梁內(nèi)力計(jì)算
    孩子的生命內(nèi)力需要家長(zhǎng)去激發(fā)
    一種適用于任意復(fù)雜結(jié)構(gòu)的曲六面體網(wǎng)格生成算法
    新型透空式六面體在南匯東灘促淤二期工程中的應(yīng)用
    質(zhì)點(diǎn)的直線運(yùn)動(dòng)
    質(zhì)點(diǎn)的直線運(yùn)動(dòng)
    基于六面體網(wǎng)格的水下航行體流體動(dòng)力分析
    電子制作(2017年24期)2017-02-02 07:14:27
    黄频高清免费视频| 午夜免费鲁丝| 亚洲专区国产一区二区| 午夜福利欧美成人| 国产伦人伦偷精品视频| 欧美丝袜亚洲另类 | 中文字幕最新亚洲高清| 9热在线视频观看99| 亚洲精品在线观看二区| 久久久久久久午夜电影 | 神马国产精品三级电影在线观看 | 韩国av一区二区三区四区| 色哟哟哟哟哟哟| 欧美av亚洲av综合av国产av| 久久人妻福利社区极品人妻图片| 久久这里只有精品19| 久久久久久久久免费视频了| 国产无遮挡羞羞视频在线观看| 国产主播在线观看一区二区| 一个人观看的视频www高清免费观看 | 欧美午夜高清在线| 欧美在线一区亚洲| 亚洲第一欧美日韩一区二区三区| 国产亚洲精品综合一区在线观看 | 一级作爱视频免费观看| 国产精品电影一区二区三区| 999久久久国产精品视频| 99精品欧美一区二区三区四区| 精品福利永久在线观看| 精品久久久久久成人av| 国产精品免费视频内射| 男女下面进入的视频免费午夜 | 亚洲狠狠婷婷综合久久图片| 午夜精品久久久久久毛片777| www.熟女人妻精品国产| 麻豆成人av在线观看| 午夜a级毛片| 国产伦人伦偷精品视频| 色综合欧美亚洲国产小说| 一本综合久久免费| 女人高潮潮喷娇喘18禁视频| 99国产综合亚洲精品| 人人妻人人澡人人看| 欧美中文综合在线视频| 亚洲一区中文字幕在线| 男女做爰动态图高潮gif福利片 | 免费在线观看黄色视频的| 50天的宝宝边吃奶边哭怎么回事| 久久婷婷成人综合色麻豆| 亚洲精品一二三| 一级片免费观看大全| 国产亚洲欧美精品永久| 手机成人av网站| 搡老岳熟女国产| 国产有黄有色有爽视频| 成人三级做爰电影| 国产免费男女视频| 男女午夜视频在线观看| 久久精品人人爽人人爽视色| aaaaa片日本免费| 色播在线永久视频| 99re在线观看精品视频| 亚洲一码二码三码区别大吗| 欧美一级毛片孕妇| 一级黄色大片毛片| 亚洲熟妇中文字幕五十中出 | 自线自在国产av| 国产亚洲欧美在线一区二区| e午夜精品久久久久久久| 在线国产一区二区在线| 韩国精品一区二区三区| 视频区欧美日本亚洲| 亚洲人成电影免费在线| 国产激情欧美一区二区| 亚洲少妇的诱惑av| 亚洲一码二码三码区别大吗| 制服人妻中文乱码| 纯流量卡能插随身wifi吗| 久久久水蜜桃国产精品网| 真人一进一出gif抽搐免费| 91老司机精品| 韩国精品一区二区三区| 亚洲精品中文字幕一二三四区| 成人国语在线视频| 欧美激情久久久久久爽电影 | 性色av乱码一区二区三区2| 国产精品亚洲av一区麻豆| 午夜福利影视在线免费观看| 一区福利在线观看| 51午夜福利影视在线观看| 久久草成人影院| 老司机亚洲免费影院| 亚洲国产看品久久| 国产精华一区二区三区| 免费观看精品视频网站| 精品人妻1区二区| 久久久久久大精品| 成人国语在线视频| 97人妻天天添夜夜摸| 精品人妻在线不人妻| 视频在线观看一区二区三区| 国产人伦9x9x在线观看| 纯流量卡能插随身wifi吗| 国产成人精品在线电影| 桃色一区二区三区在线观看| 国产乱人伦免费视频| 日韩欧美一区视频在线观看| 中文字幕人妻熟女乱码| 成人永久免费在线观看视频| 亚洲欧美日韩无卡精品| 久久久久九九精品影院| 国产在线精品亚洲第一网站| 欧美亚洲日本最大视频资源| av在线天堂中文字幕 | 色哟哟哟哟哟哟| 美女高潮到喷水免费观看| aaaaa片日本免费| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲五月婷婷丁香| 精品国内亚洲2022精品成人| 国产野战对白在线观看| 亚洲七黄色美女视频| 久久久国产一区二区| 国产精品一区二区精品视频观看| 人人妻人人爽人人添夜夜欢视频| 日韩一卡2卡3卡4卡2021年| 又黄又粗又硬又大视频| 欧美老熟妇乱子伦牲交| 女警被强在线播放| 高清欧美精品videossex| 精品国产一区二区三区四区第35| 在线观看免费日韩欧美大片| 国产熟女xx| 丁香欧美五月| www.www免费av| 久久久久久久久免费视频了| 美女大奶头视频| av网站免费在线观看视频| 久久久久精品国产欧美久久久| 国产无遮挡羞羞视频在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 国产黄色免费在线视频| 50天的宝宝边吃奶边哭怎么回事| 午夜免费成人在线视频| 色精品久久人妻99蜜桃| 黄色毛片三级朝国网站| 国产亚洲精品第一综合不卡| 久久久久九九精品影院| 国产亚洲精品综合一区在线观看 | 天堂影院成人在线观看| 神马国产精品三级电影在线观看 | 丰满饥渴人妻一区二区三| 国产一区二区三区综合在线观看| 亚洲成a人片在线一区二区| 在线十欧美十亚洲十日本专区| 黄色女人牲交| 日本wwww免费看| 欧美国产精品va在线观看不卡| 欧美日韩黄片免| 亚洲中文av在线| 国产97色在线日韩免费| 亚洲精品国产区一区二| 日韩视频一区二区在线观看| 国产精品久久电影中文字幕| 妹子高潮喷水视频| 国产av精品麻豆| 在线av久久热| 麻豆成人av在线观看| 亚洲,欧美精品.| 久9热在线精品视频| 黄色毛片三级朝国网站| 女人高潮潮喷娇喘18禁视频| 久久天躁狠狠躁夜夜2o2o| 天堂影院成人在线观看| 999精品在线视频| 国产午夜精品久久久久久| 国产欧美日韩综合在线一区二区| 欧美最黄视频在线播放免费 | 国产有黄有色有爽视频| 国产av一区在线观看免费| 欧美成人午夜精品| 亚洲av电影在线进入| 在线av久久热| 动漫黄色视频在线观看| 黄色视频,在线免费观看| 成年女人毛片免费观看观看9| 一边摸一边抽搐一进一出视频| 两个人免费观看高清视频| 久久性视频一级片| 免费一级毛片在线播放高清视频 | 在线永久观看黄色视频| 真人一进一出gif抽搐免费| 亚洲av片天天在线观看| 在线视频色国产色| 亚洲精华国产精华精| 夜夜躁狠狠躁天天躁| 亚洲专区中文字幕在线| 啪啪无遮挡十八禁网站| 国产av一区在线观看免费| 亚洲 欧美 日韩 在线 免费| 日本黄色日本黄色录像| 欧美黄色淫秽网站| 日韩中文字幕欧美一区二区| 90打野战视频偷拍视频| 亚洲一区高清亚洲精品| 黄网站色视频无遮挡免费观看| 老汉色av国产亚洲站长工具| 精品欧美一区二区三区在线| 亚洲精品中文字幕在线视频| 少妇裸体淫交视频免费看高清 | 成人影院久久| 国产三级在线视频| 久久精品亚洲av国产电影网| 91在线观看av| 久久精品国产清高在天天线| 国产区一区二久久| 中文亚洲av片在线观看爽| 五月开心婷婷网| 首页视频小说图片口味搜索| 久久香蕉精品热| av在线天堂中文字幕 | 国产精品香港三级国产av潘金莲| 亚洲一卡2卡3卡4卡5卡精品中文| 啦啦啦在线免费观看视频4| 久久久久亚洲av毛片大全| 性色av乱码一区二区三区2| 很黄的视频免费| 亚洲国产精品合色在线| 国产亚洲欧美在线一区二区| 欧美久久黑人一区二区| 在线视频色国产色| 国产精品秋霞免费鲁丝片| 国产精品98久久久久久宅男小说| 黄频高清免费视频| 免费在线观看日本一区| 少妇被粗大的猛进出69影院| 精品久久蜜臀av无| 国产高清视频在线播放一区| 亚洲五月婷婷丁香| 九色亚洲精品在线播放| 中文字幕人妻丝袜制服| 久久国产亚洲av麻豆专区| 国产精品香港三级国产av潘金莲| e午夜精品久久久久久久| 精品国产亚洲在线| 视频在线观看一区二区三区| 又黄又爽又免费观看的视频| 国内毛片毛片毛片毛片毛片| 91精品国产国语对白视频| 国产免费av片在线观看野外av| 无人区码免费观看不卡| 亚洲人成网站在线播放欧美日韩| 人成视频在线观看免费观看| 又黄又粗又硬又大视频| 国产三级黄色录像| 人人澡人人妻人| 国产真人三级小视频在线观看| 一级作爱视频免费观看| 变态另类成人亚洲欧美熟女 | 亚洲一区高清亚洲精品| 12—13女人毛片做爰片一| 亚洲黑人精品在线| 精品国产国语对白av| 人人妻人人爽人人添夜夜欢视频| 亚洲精品国产色婷婷电影| 狂野欧美激情性xxxx| av电影中文网址| 国产成年人精品一区二区 | 88av欧美| 免费av中文字幕在线| 久久人人精品亚洲av| 真人一进一出gif抽搐免费| 手机成人av网站| 成在线人永久免费视频| 中文字幕高清在线视频| 亚洲精品美女久久av网站| 午夜成年电影在线免费观看| 亚洲成人免费av在线播放| 波多野结衣一区麻豆| 国产精品98久久久久久宅男小说| 亚洲国产精品一区二区三区在线| 韩国av一区二区三区四区| 久久精品影院6| 欧洲精品卡2卡3卡4卡5卡区| 在线观看一区二区三区激情| 色婷婷av一区二区三区视频| 亚洲 国产 在线| 青草久久国产| 亚洲成人免费电影在线观看| 久久久久久久久久久久大奶| 757午夜福利合集在线观看| 国产精品爽爽va在线观看网站 | 亚洲色图 男人天堂 中文字幕| 手机成人av网站| 亚洲九九香蕉| 韩国av一区二区三区四区| 成年人免费黄色播放视频| 国产片内射在线| 亚洲国产精品sss在线观看 | 日韩免费av在线播放| 99riav亚洲国产免费| 不卡av一区二区三区| 国产亚洲欧美在线一区二区| 亚洲在线自拍视频| 国产成年人精品一区二区 | 国产午夜精品久久久久久| 最近最新中文字幕大全电影3 | 国产aⅴ精品一区二区三区波| 欧美在线黄色| 少妇的丰满在线观看| 神马国产精品三级电影在线观看 | 亚洲成人精品中文字幕电影 | 国产一级毛片七仙女欲春2| 91久久精品国产一区二区成人| 搡老熟女国产l中国老女人| 床上黄色一级片| 又黄又爽又刺激的免费视频.| 日本 av在线| 国产成人欧美在线观看| 国产野战对白在线观看| 欧美3d第一页| 亚洲国产精品sss在线观看| 免费在线观看成人毛片| 久久欧美精品欧美久久欧美| 午夜激情福利司机影院| 国产三级中文精品| 美女 人体艺术 gogo| 99久久久亚洲精品蜜臀av| 久久久久九九精品影院| 夜夜看夜夜爽夜夜摸| 男女视频在线观看网站免费| 一个人看的www免费观看视频| 精品免费久久久久久久清纯| 熟妇人妻久久中文字幕3abv| 男女之事视频高清在线观看| 一本综合久久免费| 国产精品精品国产色婷婷| 欧美潮喷喷水| 精品久久国产蜜桃| 国产在线男女| 国产一区二区三区在线臀色熟女| 国产亚洲精品av在线| 久久精品综合一区二区三区| 在线国产一区二区在线| 最近视频中文字幕2019在线8| a级一级毛片免费在线观看| 国产欧美日韩精品亚洲av| 中文字幕人妻熟人妻熟丝袜美| 搡女人真爽免费视频火全软件 | 久久久久久久精品吃奶| 伊人久久精品亚洲午夜| 亚洲欧美日韩卡通动漫| 国产白丝娇喘喷水9色精品| 天美传媒精品一区二区| 国内精品一区二区在线观看| 久久久久久久精品吃奶| 国产精品1区2区在线观看.| 一区二区三区四区激情视频 | 国产亚洲av嫩草精品影院| 麻豆成人av在线观看| 亚洲精品影视一区二区三区av| 人妻制服诱惑在线中文字幕| 69av精品久久久久久| 亚洲va日本ⅴa欧美va伊人久久| 一级毛片久久久久久久久女| 国产午夜精品论理片| 露出奶头的视频| 人人妻人人看人人澡| 国产精品日韩av在线免费观看| 日韩大尺度精品在线看网址| 女生性感内裤真人,穿戴方法视频| 脱女人内裤的视频| 亚洲精品粉嫩美女一区| 色综合婷婷激情| 亚洲av电影在线进入| 免费av观看视频| 嫩草影院精品99| 亚洲午夜理论影院| 99国产综合亚洲精品| 亚洲一区二区三区不卡视频| 亚洲最大成人av| 91在线观看av| 51午夜福利影视在线观看| 日本熟妇午夜| 亚洲国产精品成人综合色| 久久久久国产精品人妻aⅴ院| 一区二区三区四区激情视频 | 欧美色视频一区免费| 999久久久精品免费观看国产| 亚洲最大成人手机在线| 日韩欧美 国产精品| 99精品久久久久人妻精品| 午夜福利高清视频| 97超级碰碰碰精品色视频在线观看| 神马国产精品三级电影在线观看| 久久久国产成人精品二区| 91av网一区二区| 午夜福利成人在线免费观看| 久久这里只有精品中国| av在线老鸭窝| АⅤ资源中文在线天堂| 午夜日韩欧美国产| 亚洲av第一区精品v没综合| 亚洲国产精品合色在线| 中国美女看黄片| 夜夜看夜夜爽夜夜摸| 日韩亚洲欧美综合| 久久热精品热| 国产精品亚洲av一区麻豆| 美女免费视频网站| 国产av一区在线观看免费| 色吧在线观看| 啦啦啦观看免费观看视频高清| 国产av不卡久久| a级一级毛片免费在线观看| 一个人免费在线观看的高清视频| 成人永久免费在线观看视频| 夜夜爽天天搞| 中文亚洲av片在线观看爽| 午夜精品在线福利| 波多野结衣高清无吗| 99久久精品国产亚洲精品| 久久国产乱子伦精品免费另类| 亚洲av不卡在线观看| 久久精品国产亚洲av香蕉五月| 天堂影院成人在线观看| 床上黄色一级片| 日韩欧美在线乱码| av视频在线观看入口| .国产精品久久| 青草久久国产| 男人舔奶头视频| 女人被狂操c到高潮| 一a级毛片在线观看| 99久久成人亚洲精品观看| 深夜精品福利| 国产av麻豆久久久久久久| 99热这里只有精品一区| 最新中文字幕久久久久| 亚洲第一电影网av| 天堂动漫精品| 国产精品人妻久久久久久| 深爱激情五月婷婷| 婷婷精品国产亚洲av| 国产精华一区二区三区| 欧美一区二区精品小视频在线| 免费在线观看成人毛片| 久久精品久久久久久噜噜老黄 | 国产精品久久久久久久电影| 一进一出抽搐动态| 国产av在哪里看| 欧美成人免费av一区二区三区| 在线观看一区二区三区| 国产av一区在线观看免费| 亚洲国产欧洲综合997久久,| 在线播放无遮挡| 性色av乱码一区二区三区2| 激情在线观看视频在线高清| a在线观看视频网站| 91久久精品电影网| 成人鲁丝片一二三区免费| 久久亚洲精品不卡| 国产探花在线观看一区二区| 露出奶头的视频| 又黄又爽又免费观看的视频| 禁无遮挡网站| 日本五十路高清| 中文字幕免费在线视频6| 在线免费观看的www视频| 麻豆一二三区av精品| 欧美一区二区亚洲| 亚洲专区国产一区二区| 日韩有码中文字幕| 国产真实伦视频高清在线观看 | 午夜久久久久精精品| 久久婷婷人人爽人人干人人爱| 很黄的视频免费| 精品国产亚洲在线| 男女视频在线观看网站免费| 色综合婷婷激情| 十八禁网站免费在线| 尤物成人国产欧美一区二区三区| 蜜桃久久精品国产亚洲av| 91av网一区二区| 国产精品久久电影中文字幕| eeuss影院久久| 蜜桃亚洲精品一区二区三区| 久久人人精品亚洲av| 又爽又黄无遮挡网站| 日韩成人在线观看一区二区三区| 非洲黑人性xxxx精品又粗又长| 日本成人三级电影网站| 亚洲av成人不卡在线观看播放网| 亚洲精品粉嫩美女一区| 午夜精品一区二区三区免费看| 在线免费观看的www视频| 老女人水多毛片| 日本 av在线| 亚洲av电影不卡..在线观看| 亚洲精品日韩av片在线观看| 夜夜夜夜夜久久久久| 中文字幕av成人在线电影| 欧美xxxx性猛交bbbb| 欧美三级亚洲精品| 久久午夜福利片| 久久久久国产精品人妻aⅴ院| 黄色视频,在线免费观看| 日韩欧美在线乱码| 露出奶头的视频| 看黄色毛片网站| 久久精品影院6| 亚洲天堂国产精品一区在线| 亚洲美女视频黄频| 国产成人欧美在线观看| 少妇人妻一区二区三区视频| 国产单亲对白刺激| 欧美在线黄色| 中文字幕免费在线视频6| 少妇丰满av| 色综合欧美亚洲国产小说| eeuss影院久久| 亚洲avbb在线观看| 精品国产三级普通话版| 久久伊人香网站| 国产伦精品一区二区三区视频9| 欧美极品一区二区三区四区| 男女做爰动态图高潮gif福利片| 蜜桃久久精品国产亚洲av| 免费av毛片视频| 国产在视频线在精品| 久久久久久久久久黄片| 啦啦啦观看免费观看视频高清| www.999成人在线观看| 中出人妻视频一区二区| av专区在线播放| 欧美日韩瑟瑟在线播放| 久久6这里有精品| 久久香蕉精品热| 青草久久国产| 免费无遮挡裸体视频| 99riav亚洲国产免费| 999久久久精品免费观看国产| 精品一区二区三区视频在线| 九九在线视频观看精品| 国产亚洲av嫩草精品影院| 12—13女人毛片做爰片一| 久久午夜亚洲精品久久| 久久国产精品影院| 最近最新免费中文字幕在线| 老熟妇乱子伦视频在线观看| 精品久久久久久久末码| 真实男女啪啪啪动态图| 91久久精品电影网| 97热精品久久久久久| 99久久九九国产精品国产免费| 一个人免费在线观看电影| 中文字幕高清在线视频| 99在线人妻在线中文字幕| 看黄色毛片网站| 俄罗斯特黄特色一大片| 国产伦一二天堂av在线观看| 搡老妇女老女人老熟妇| 国语自产精品视频在线第100页| 毛片女人毛片| 亚洲国产高清在线一区二区三| 别揉我奶头~嗯~啊~动态视频| 午夜免费成人在线视频| 中文字幕免费在线视频6| 久久久久久久久久黄片| 色综合婷婷激情| 欧美zozozo另类| 国产亚洲精品久久久久久毛片| 国产精品嫩草影院av在线观看 | 国产精品久久电影中文字幕| 亚洲一区二区三区色噜噜| 国产精品自产拍在线观看55亚洲| 中文资源天堂在线| 97超视频在线观看视频| 色视频www国产| 精品久久久久久久久亚洲 | 午夜精品一区二区三区免费看| 久久精品综合一区二区三区| 一本综合久久免费| 狠狠狠狠99中文字幕| 在线观看免费视频日本深夜| 99在线视频只有这里精品首页| 91字幕亚洲| 如何舔出高潮| 深夜精品福利| 午夜免费男女啪啪视频观看 | 特级一级黄色大片| 天美传媒精品一区二区| 亚洲色图av天堂| 亚洲aⅴ乱码一区二区在线播放| 亚洲第一区二区三区不卡| 成年人黄色毛片网站| av天堂在线播放| 国产久久久一区二区三区| 亚洲av成人av| 亚洲午夜理论影院| 99国产精品一区二区蜜桃av| 日本黄色视频三级网站网址| 久久九九热精品免费| 波多野结衣高清无吗| 亚洲av电影不卡..在线观看| 国产精品,欧美在线| 亚洲在线自拍视频| 国产精品爽爽va在线观看网站| 久久久久国产精品人妻aⅴ院| 俺也久久电影网| 精品人妻偷拍中文字幕| 亚洲欧美精品综合久久99| 欧美黑人欧美精品刺激| netflix在线观看网站| 69av精品久久久久久|