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

    幾何非線性分析的高效高階無網(wǎng)格法

    2020-12-23 02:16:42陳嵩濤段慶林馬今偉
    計算力學(xué)學(xué)報 2020年6期
    關(guān)鍵詞:網(wǎng)格法構(gòu)形軸承座

    陳嵩濤, 段慶林, 馬今偉

    (大連理工大學(xué) 工業(yè)裝備結(jié)構(gòu)分析國家重點實驗室,大連116023)

    1 引 言

    許多工程問題的力學(xué)分析需考慮結(jié)構(gòu)構(gòu)形的變化,如橡膠材料的大變形等。這一類問題的剛體轉(zhuǎn)動和應(yīng)變很劇烈,小變形假設(shè)失效,需要采用幾何非線性模型進行計算。

    傳統(tǒng)有限元法的形函數(shù)依賴于網(wǎng)格單元,在處理幾何非線性問題時容易因網(wǎng)格畸變而導(dǎo)致計算失敗。與此不同,無網(wǎng)格法[1]的形函數(shù)僅依賴于離散點,無需網(wǎng)格單元,因而處理大變形問題更具優(yōu)勢。目前,無網(wǎng)格法已在幾何非線性的數(shù)值分析中得到廣泛應(yīng)用[2-4]。

    與線性分析不同,求解幾何非線性問題時參考構(gòu)形的選取關(guān)系到分析效率與結(jié)果的精確性,有限元通常選擇相鄰構(gòu)形作為參考構(gòu)形,即UL法,相比于有限元,無網(wǎng)格法求解形函數(shù)導(dǎo)數(shù)過程較為復(fù)雜,使用UL法不斷更新構(gòu)形,必然會降低計算效率。針對該問題,需要建立既能及時跟蹤變形又不產(chǎn)生極大計算代價的構(gòu)形更新方式。Léger等[5]發(fā)展了一種介于TL法和UL法之間的構(gòu)形更新方法,即在每個載荷步迭代平衡后將此構(gòu)形作為新的參考構(gòu)形,這種方法比較適用于形函數(shù)復(fù)雜、更新構(gòu)形代價較大的無網(wǎng)格法。

    由于無網(wǎng)格形函數(shù)的豐富性,Galerkin弱形式的準確積分要求較多的區(qū)域積分點,嚴重降低了計算效率。無網(wǎng)格領(lǐng)域已發(fā)展了一些高效積分方法,如Wang等[6]的子域積分、Beissel等[7]的節(jié)點積分和Duan等[8]的穩(wěn)定應(yīng)力點積分等。其中,Chen等[9]發(fā)展的基于應(yīng)變光順的穩(wěn)定相容節(jié)點積分SCNI(Stabilized Conforming Nodal Integration)方法極具吸引力,其突出優(yōu)點是不引入任何人工參數(shù),且能精確通過線性分片試驗。該方法已成功應(yīng)用于幾何非線性分析[10,11],在高度不規(guī)則計算網(wǎng)格下仍展現(xiàn)出良好的計算穩(wěn)定性和效率。Wang等[12]進一步將該方法拓展到板單元的大變形問題,成功克服了剪切自鎖。

    Puso等[13]發(fā)現(xiàn)SCNI仍然可能導(dǎo)致邊界附近鋸齒形振蕩模式的出現(xiàn)。Duan等[14]也指出,即使對于二次無網(wǎng)格近似,SCNI也只能通過線性分片試驗,而不能通過二次分片試驗。本質(zhì)上,SCNI是一點積分方案,只能反映出常數(shù)應(yīng)變場。對于應(yīng)變分布可能較為復(fù)雜的非線性問題,高階方法更能夠精確地再現(xiàn)變形過程,避免體積自鎖。Ortiz等[15]指出,在一些大應(yīng)變算例中低階方法魯棒性較差,對這一類問題發(fā)展高效的高階方法是必要的。Duan等[14]針對線彈性小變形問題發(fā)展了一個三點積分方案,積分點上的形函數(shù)導(dǎo)數(shù)由其與形函數(shù)之間散度定理的離散形式計算。該方法可精確通過二次分片試驗,因而稱之為二階一致三點積分方法QC3(Quadratically Consistent 3-point integration method)。相應(yīng)的,SCNI可稱為線性一致一點積分方法LC1(Linearly Consistent 1-point integration method)。目前,QC3積分方案已成功應(yīng)用于不可壓縮固體[15]、斷裂相場[16]和斯托克斯流[17]等問題,但這些研究仍主要集中在線性分析方面。

    本文把QC3積分方案拓展到幾何非線性分析,并進一步引入Léger等[5]的更新構(gòu)形方法,發(fā)展能夠高效分析大變形問題的高階無網(wǎng)格法,并通過數(shù)值結(jié)果考察和驗證方法的計算精度與效率。

    2 幾何非線性分析列式

    如圖1所示,考慮當前構(gòu)形為Ωn + 1,坐標為xn + 1,自然邊界為Γn + 1,本質(zhì)邊界為Γun + 1,Γun + 1∪Γtn + 1=Γn + 1,Γun + 1∩Γtn + 1= ?。與標準的TL和UL方法不同,本文選取第n個載荷步迭代平衡后的構(gòu)形Ωn作為參考構(gòu)形,該構(gòu)形上的各物理量均標注下標n,如圖1的Γtn和Γun。相應(yīng)地,初始構(gòu)形Ω0的各物理量均標注下標0。圖1的Fn和fn + 1分別表示Ωn相對于Ω0以及Ωn + 1相對于Ωn的變形梯度。

    當前構(gòu)形下,平衡方程及邊界條件為

    (1)

    (2)

    σi jnj=tn + 1ionΓtn + 1

    (3)

    (4)

    式中左端為內(nèi)力虛功Wint,右端為外力虛功Wext。采用本文選取的參考構(gòu)形Ωn,式(4)可表示為

    (5)

    式中bn和tn為定義在參考構(gòu)形下的體力和面力載荷,Sn和En分別為

    (6)

    (7)

    圖1 計算構(gòu)形

    式中Jn= |Fn|,I為二階單位張量,S和E分別為參考初始構(gòu)形的第二Piola-Kirchhoff應(yīng)力(PK2應(yīng)力)和格林應(yīng)變。顯然,Sn和En分別為前推到參考構(gòu)形Ωn的PK2應(yīng)力和格林應(yīng)變[5]。

    采用節(jié)點的無網(wǎng)格形函數(shù)建立位移近似u=NIuI,并代入式(5),則該式左端可表示為

    (8)

    (9)

    (10)

    (11)

    (12)

    (13)

    (14)

    (15)

    (16)

    式中m為節(jié)點數(shù)目,un=xn + 1-xn為參考構(gòu)形到當前構(gòu)形的位移增量。須著重指出的是,本文方法采用Ωn作為計算參考構(gòu)形,因而式(12,14,16)的空間導(dǎo)數(shù)均為對xn的導(dǎo)數(shù),這是與參考初始構(gòu)形的傳統(tǒng)TL法以及參考相鄰構(gòu)形的傳統(tǒng)UL法的核心區(qū)別。

    非線性方程的求解采用牛頓-拉普森迭代算法,計算流程簡述如下。

    (17)

    (3) 建立切線剛度陣K,并通過罰函數(shù)法施加本質(zhì)邊界條件求解增量平衡方程

    (18)

    (19)

    3 二階一致三點積分格式

    由以上列式可以看到,剛度陣的建立需要計算形函數(shù)導(dǎo)數(shù)的區(qū)域積分。由于無網(wǎng)格形函數(shù)是非多項式的有理函數(shù),精確計算該區(qū)域積分需要較多的積分點,故計算效率下降。針對該問題,本文引入Duan等[14]針對線彈性小變形問題發(fā)展的三點積分格式QC3來減少積分點數(shù)目。

    如圖2所示,QC3方法采用背景三角形單元進行積分計算,每個單元僅使用三個積分點(五角星),邊界積分在每條邊上使用兩個一維高斯點(三角形)。QC3方法的核心思想是導(dǎo)數(shù)修正,即積分點上形函數(shù)的空間導(dǎo)數(shù)不采用標準形式,而是由散度定理(20)確定

    (20)

    式中ΩSn為參考構(gòu)形下的背景積分子域,ΓSn為其邊界,q(x) =p,x(x)∪p,y(x),p(x)為MLS形函數(shù)的基底函數(shù)向量。本文采用二次近似,即p(x) =[1xyx2xyy2]T,則q(x) = [1xy]T,因而

    圖2 QC3積分方法

    式(20)實際上表示了三個方程,剛好可以用于確定三個積分點上形函數(shù)的空間導(dǎo)數(shù)。采用數(shù)值積分,式(20)的離散形式為

    (21)

    4 數(shù)值算例

    采用三個算例考察和驗證本文發(fā)展的無網(wǎng)格方法處理幾何非線性問題的有效性。作為對比,對有限元(線性三節(jié)點三角形單元,標識為FEM)、LC1以及積分點上直接使用標準導(dǎo)數(shù)(而非修正導(dǎo)數(shù))的無網(wǎng)格法(標識為ST16,即每個背景三角形單元使用16個積分點)也進行了程序?qū)崿F(xiàn)。

    4.1 淺拱

    如圖3所示,兩端固支弧形淺拱在頂部中心點A處受向下的集中載荷作用[10]。淺拱中性軸半徑為100 mm,拱厚為2 mm,弧對應(yīng)圓心角為 28.06°,楊氏模量E=4.8×103N/mm2,泊松比υ=0.0。無網(wǎng)格離散采用105個節(jié)點,分20個加載步,每步載荷為1 N,各種方法得到的A點處的載荷-位移響應(yīng)如圖4所示??梢钥闯觯疚牡腝C3方法以及ST16方法的結(jié)果與解析解吻合良好,但是FEM和LC1方法的計算結(jié)果明顯偏離了解析解。應(yīng)說明的是,F(xiàn)EM的結(jié)果通過加密網(wǎng)格和改變構(gòu)形更新方式(如采用UL法)可以得到顯著改善。

    圖3 受集中載荷的淺拱

    圖4 淺拱A點處的載荷-位移響應(yīng)

    QC3與ST16的計算精度相當,但QC3更為高效,表1比較了三種無網(wǎng)格方法在四種計算節(jié)點配置下消耗的CPU時間。顯然,QC3的計算時間大約僅為ST16的1/4。

    4.2 懸臂梁

    如圖5所示,懸臂梁一端固支,另一端部中心A處受向上的集中載荷作用[10],梁長為20 mm,厚度為1 mm,楊氏模量E=4.8×103N/mm2,泊松比υ=0.0。采用405個節(jié)點離散,分100個加載步加載,每步載荷為0.1 N,QC3方法得到的變形過程如圖6所示,載荷-位移響應(yīng)的比較如圖7所示。QC3和ST16的結(jié)果與解析解吻合,LC1的精度相對較差,F(xiàn)EM則明顯偏離解析解。

    圖8比較了各方法得到的應(yīng)力場。FEM和LC1的應(yīng)力場均出現(xiàn)了明顯的虛假數(shù)值振蕩。本文QC3方法使用的積分點明顯少于ST16,但兩種方法得到了類似的光滑無振蕩的應(yīng)力場。

    表1 淺拱算例的CPU時間比較

    圖5 受集中載荷的懸臂梁

    圖6 QC3方法得到的不同載荷下變形后的懸臂梁

    4.3 軸承座

    如圖9所示,軸承座尺寸為10 cm×10 cm,大孔位于正中心,四角小孔對稱分布,h=1 cm,r=0.5 cm,R=4.5 cm。楊氏模量E=7.9×104N/cm2,泊松比υ=0.33。軸承座的上下表面分別施加向下和向上的不隨構(gòu)形變化的均勻面載荷。無網(wǎng)格計算采用了3830個節(jié)點,如圖9所示,分50個加載步,前25步每步加載2 N直至達到極值載荷50 N,后25步每步卸載2 N直至完全卸載。

    圖7 懸臂梁A點處的載荷-位移響應(yīng)

    圖8 懸臂梁的σx x場

    圖9 軸承座的幾何模型和計算節(jié)點

    圖10為P=20 N時軸承座的y向正應(yīng)力場。以Abaqus結(jié)果為參照,QC3得到的應(yīng)力極值較為準確,而LC1明顯偏離參考解。圖11展示了QC3方法得到的軸承座加卸載的變形過程。在完全卸載后軸承座與初始構(gòu)形吻合良好,無虛假殘余變形,這進一步驗證了本方法計算復(fù)雜結(jié)構(gòu)大變形的能力。

    圖10 軸承座的σy y場

    圖11 QC3得到的軸承座的變形過程

    5 結(jié) 論

    本文將QC3方法拓展到了幾何非線性分析,并采用上一載荷步的收斂構(gòu)形作為計算參考構(gòu)形,針對超彈性材料建立了大變形分析的高階無網(wǎng)格法。與標準的高階無網(wǎng)格法(ST16)相比,本文方法大幅度減少了積分點數(shù)目,節(jié)省了大量的計算時間。與一點積分的高階無網(wǎng)格法(LC1)相比,本文方法具有高得多的計算精度。因此,本文方法有效提高了高階無網(wǎng)格法幾何非線性分析的計算效率,是超彈性材料大變形問題值得考慮的計算方案。

    猜你喜歡
    網(wǎng)格法構(gòu)形軸承座
    調(diào)相機軸承座振動優(yōu)化改進方案研究
    基于ANSYS Workbench軟件在軸承座模態(tài)分析中的應(yīng)用
    防爆電機(2022年1期)2022-02-16 01:14:00
    雙星跟飛立體成像的構(gòu)形保持控制
    雷擊條件下接地系統(tǒng)的分布參數(shù)
    通有構(gòu)形的特征多項式
    角接觸球軸承的優(yōu)化設(shè)計算法
    基于遺傳算法的機器人路徑規(guī)劃研究
    基于GIS的植物葉片信息測量研究
    425輕量化橋軸承座工藝改進
    對一個幾何構(gòu)形的探究
    麻豆成人av在线观看| 在线永久观看黄色视频| 妹子高潮喷水视频| 免费无遮挡裸体视频| 欧美乱色亚洲激情| 亚洲精品av麻豆狂野| 高潮久久久久久久久久久不卡| 操美女的视频在线观看| 好男人在线观看高清免费视频 | 他把我摸到了高潮在线观看| 日本欧美视频一区| 乱人伦中国视频| 欧美一区二区精品小视频在线| 亚洲一区二区三区色噜噜| 777久久人妻少妇嫩草av网站| 国产伦一二天堂av在线观看| 99精品在免费线老司机午夜| 午夜久久久在线观看| 大码成人一级视频| 国产精品美女特级片免费视频播放器 | 成人18禁在线播放| 校园春色视频在线观看| 黄片小视频在线播放| 韩国精品一区二区三区| 欧美亚洲日本最大视频资源| 国产一区二区三区视频了| 国产精品亚洲一级av第二区| av免费在线观看网站| 亚洲色图av天堂| 在线观看免费视频日本深夜| 久久精品亚洲熟妇少妇任你| 国产精品亚洲美女久久久| 国产亚洲精品av在线| 国产精品野战在线观看| 黑人欧美特级aaaaaa片| 精品第一国产精品| 久久久久精品国产欧美久久久| 亚洲男人的天堂狠狠| 精品日产1卡2卡| 大码成人一级视频| 国产精品久久电影中文字幕| 国产精品亚洲av一区麻豆| 大码成人一级视频| 亚洲av电影不卡..在线观看| 首页视频小说图片口味搜索| 在线观看免费午夜福利视频| 国产色视频综合| 国产亚洲av嫩草精品影院| 亚洲精品久久国产高清桃花| 久久久久九九精品影院| 成人手机av| 两性夫妻黄色片| 香蕉国产在线看| 国内毛片毛片毛片毛片毛片| 中文字幕人妻丝袜一区二区| 悠悠久久av| 激情视频va一区二区三区| 久久精品91蜜桃| 亚洲熟女毛片儿| 少妇 在线观看| 国产亚洲精品久久久久久毛片| 久久精品亚洲熟妇少妇任你| 久久久久国产精品人妻aⅴ院| 激情在线观看视频在线高清| 又紧又爽又黄一区二区| 精品少妇一区二区三区视频日本电影| 欧美国产精品va在线观看不卡| 女性被躁到高潮视频| 老熟妇乱子伦视频在线观看| 精品福利观看| 色婷婷久久久亚洲欧美| 精品福利观看| 日韩中文字幕欧美一区二区| 妹子高潮喷水视频| 亚洲成人免费电影在线观看| 国产午夜福利久久久久久| 亚洲一区二区三区色噜噜| 亚洲国产欧美日韩在线播放| 日韩免费av在线播放| 国产av在哪里看| 波多野结衣av一区二区av| 一a级毛片在线观看| 日韩欧美一区二区三区在线观看| 男女下面进入的视频免费午夜 | 国产亚洲精品av在线| 91在线观看av| 777久久人妻少妇嫩草av网站| 欧美日韩黄片免| 午夜福利高清视频| 韩国精品一区二区三区| 亚洲第一av免费看| 亚洲成人国产一区在线观看| 亚洲精品粉嫩美女一区| 搡老妇女老女人老熟妇| 亚洲aⅴ乱码一区二区在线播放 | 国产精品免费视频内射| 麻豆成人av在线观看| 一区二区三区激情视频| 亚洲精品国产区一区二| 成人国产一区最新在线观看| 亚洲少妇的诱惑av| 亚洲精品一区av在线观看| 午夜a级毛片| 国产av一区二区精品久久| 国产私拍福利视频在线观看| 国产一区二区在线av高清观看| 久久中文字幕一级| 免费搜索国产男女视频| 欧美色视频一区免费| bbb黄色大片| 国产97色在线日韩免费| 天堂影院成人在线观看| 一区二区三区国产精品乱码| 手机成人av网站| 免费在线观看影片大全网站| 久久欧美精品欧美久久欧美| 99re在线观看精品视频| 亚洲精品久久成人aⅴ小说| 亚洲熟妇中文字幕五十中出| or卡值多少钱| 法律面前人人平等表现在哪些方面| 久久这里只有精品19| 在线av久久热| 高清在线国产一区| 看片在线看免费视频| 亚洲国产日韩欧美精品在线观看 | 热99re8久久精品国产| 成人亚洲精品av一区二区| 伦理电影免费视频| 久久 成人 亚洲| 在线十欧美十亚洲十日本专区| 亚洲精品在线美女| 亚洲少妇的诱惑av| 色综合亚洲欧美另类图片| 欧美日本中文国产一区发布| 90打野战视频偷拍视频| 色播亚洲综合网| 一边摸一边抽搐一进一小说| 午夜精品在线福利| 欧美一级毛片孕妇| 国产激情欧美一区二区| 成人国产综合亚洲| 色哟哟哟哟哟哟| 天天躁狠狠躁夜夜躁狠狠躁| 国产高清激情床上av| 久久精品人人爽人人爽视色| 日韩高清综合在线| 在线av久久热| 50天的宝宝边吃奶边哭怎么回事| 国产亚洲精品av在线| 国产视频一区二区在线看| 免费看a级黄色片| 99热只有精品国产| 免费看美女性在线毛片视频| 国产亚洲av嫩草精品影院| 性欧美人与动物交配| 午夜久久久久精精品| 欧美乱色亚洲激情| 亚洲人成网站在线播放欧美日韩| 久久欧美精品欧美久久欧美| 他把我摸到了高潮在线观看| 午夜福利免费观看在线| 久久久久久久精品吃奶| 欧美人与性动交α欧美精品济南到| 亚洲第一欧美日韩一区二区三区| 可以在线观看毛片的网站| 一级a爱视频在线免费观看| 18禁观看日本| 国产精品久久久人人做人人爽| 久久久精品欧美日韩精品| 午夜福利视频1000在线观看 | 亚洲成人久久性| ponron亚洲| 久久婷婷人人爽人人干人人爱 | 久久久久久久久免费视频了| 亚洲精品av麻豆狂野| 精品福利观看| 久久久国产欧美日韩av| 俄罗斯特黄特色一大片| 国产精品美女特级片免费视频播放器 | 亚洲,欧美精品.| 露出奶头的视频| 国产97色在线日韩免费| 韩国av一区二区三区四区| 精品电影一区二区在线| 国产一区二区三区视频了| 高潮久久久久久久久久久不卡| 韩国av一区二区三区四区| 日韩免费av在线播放| 在线永久观看黄色视频| 男女之事视频高清在线观看| 久99久视频精品免费| 国产日韩一区二区三区精品不卡| 亚洲av成人一区二区三| 夜夜看夜夜爽夜夜摸| 成人三级黄色视频| 琪琪午夜伦伦电影理论片6080| 又大又爽又粗| 色综合欧美亚洲国产小说| 精品国产国语对白av| 悠悠久久av| 他把我摸到了高潮在线观看| 18禁美女被吸乳视频| 精品国产一区二区久久| 夜夜躁狠狠躁天天躁| 亚洲专区国产一区二区| 国产亚洲精品第一综合不卡| 亚洲久久久国产精品| 免费观看人在逋| 欧美色视频一区免费| 国产片内射在线| 成人国产综合亚洲| 久久久国产欧美日韩av| 国产精品自产拍在线观看55亚洲| 后天国语完整版免费观看| √禁漫天堂资源中文www| 亚洲一区二区三区不卡视频| 久久久久亚洲av毛片大全| 国产精品二区激情视频| 国产成人精品无人区| 人妻久久中文字幕网| 午夜福利影视在线免费观看| 伊人久久大香线蕉亚洲五| 欧美黑人欧美精品刺激| 久久精品成人免费网站| 国产激情久久老熟女| 丝袜美足系列| 欧美亚洲日本最大视频资源| 97人妻精品一区二区三区麻豆 | 亚洲色图 男人天堂 中文字幕| 欧美丝袜亚洲另类 | 国产精品精品国产色婷婷| 视频在线观看一区二区三区| 久久草成人影院| 国产成+人综合+亚洲专区| 欧美+亚洲+日韩+国产| 久久精品影院6| 给我免费播放毛片高清在线观看| 欧美色欧美亚洲另类二区 | 亚洲狠狠婷婷综合久久图片| 人人妻人人澡人人看| 真人做人爱边吃奶动态| 丝袜在线中文字幕| 国产91精品成人一区二区三区| 视频在线观看一区二区三区| 高清在线国产一区| 精品久久久精品久久久| 此物有八面人人有两片| 成人亚洲精品一区在线观看| av福利片在线| 亚洲 欧美一区二区三区| 老汉色∧v一级毛片| 精品乱码久久久久久99久播| 免费搜索国产男女视频| 久久草成人影院| 欧美日本视频| 久久人人97超碰香蕉20202| 亚洲成av片中文字幕在线观看| 亚洲精品久久国产高清桃花| 法律面前人人平等表现在哪些方面| 成人国产综合亚洲| 久久国产精品男人的天堂亚洲| 最近最新中文字幕大全免费视频| 欧美中文日本在线观看视频| 91在线观看av| 精品人妻1区二区| 香蕉丝袜av| 婷婷丁香在线五月| 在线观看舔阴道视频| 精品人妻1区二区| 国产精品av久久久久免费| 性少妇av在线| 精品国产国语对白av| 18禁观看日本| or卡值多少钱| 久久久久九九精品影院| 亚洲欧洲精品一区二区精品久久久| 亚洲五月色婷婷综合| 久久精品aⅴ一区二区三区四区| 免费观看人在逋| 亚洲成国产人片在线观看| 欧美最黄视频在线播放免费| 欧美+亚洲+日韩+国产| 他把我摸到了高潮在线观看| 国产精品久久久久久人妻精品电影| 美女扒开内裤让男人捅视频| 无限看片的www在线观看| av在线播放免费不卡| 色精品久久人妻99蜜桃| av天堂久久9| 身体一侧抽搐| 久久久久久国产a免费观看| 久久这里只有精品19| 熟妇人妻久久中文字幕3abv| 国产xxxxx性猛交| av超薄肉色丝袜交足视频| 制服人妻中文乱码| 麻豆久久精品国产亚洲av| 国产成年人精品一区二区| 99香蕉大伊视频| 精品卡一卡二卡四卡免费| 精品人妻1区二区| 亚洲av熟女| 国产高清videossex| 亚洲精品在线美女| 99国产精品免费福利视频| 高清毛片免费观看视频网站| 国产日韩一区二区三区精品不卡| 亚洲一区二区三区色噜噜| 久久久久久大精品| 亚洲色图 男人天堂 中文字幕| 99国产精品99久久久久| 免费看a级黄色片| 久久人妻熟女aⅴ| 欧洲精品卡2卡3卡4卡5卡区| 欧美国产精品va在线观看不卡| 中文字幕人成人乱码亚洲影| 欧美在线一区亚洲| 免费久久久久久久精品成人欧美视频| 欧美性长视频在线观看| 国产又色又爽无遮挡免费看| 久久天躁狠狠躁夜夜2o2o| 欧美最黄视频在线播放免费| 一区二区三区高清视频在线| 色综合站精品国产| 免费不卡黄色视频| 女性被躁到高潮视频| 一进一出好大好爽视频| 午夜免费激情av| 中文字幕av电影在线播放| 亚洲精品一区av在线观看| 亚洲精品中文字幕一二三四区| 成人18禁在线播放| 999精品在线视频| 国产一区二区三区在线臀色熟女| 成人亚洲精品av一区二区| 亚洲国产高清在线一区二区三 | 亚洲成人免费电影在线观看| 在线天堂中文资源库| 女人精品久久久久毛片| 激情在线观看视频在线高清| 国产一卡二卡三卡精品| 波多野结衣一区麻豆| 国内精品久久久久精免费| 国产在线观看jvid| 国产成人影院久久av| 一级毛片高清免费大全| 丁香欧美五月| 十八禁人妻一区二区| 热re99久久国产66热| 久久国产乱子伦精品免费另类| 久久精品成人免费网站| 窝窝影院91人妻| 69av精品久久久久久| 法律面前人人平等表现在哪些方面| 久久久久久免费高清国产稀缺| 国产精品免费一区二区三区在线| 午夜精品久久久久久毛片777| 咕卡用的链子| 亚洲国产精品999在线| av超薄肉色丝袜交足视频| 国产精品电影一区二区三区| 亚洲av成人不卡在线观看播放网| 性色av乱码一区二区三区2| 最近最新中文字幕大全电影3 | 亚洲欧美日韩高清在线视频| 国产午夜福利久久久久久| 久久久国产精品麻豆| 精品日产1卡2卡| 精品国产超薄肉色丝袜足j| 黑丝袜美女国产一区| 日本三级黄在线观看| 亚洲成人久久性| 搡老妇女老女人老熟妇| 国产片内射在线| 国产99白浆流出| 亚洲中文av在线| 欧美日韩亚洲国产一区二区在线观看| 在线十欧美十亚洲十日本专区| 精品国产美女av久久久久小说| 这个男人来自地球电影免费观看| 欧美大码av| 亚洲 国产 在线| 99国产精品一区二区蜜桃av| 日韩欧美国产在线观看| 黄色 视频免费看| 国内毛片毛片毛片毛片毛片| svipshipincom国产片| 12—13女人毛片做爰片一| 亚洲欧美一区二区三区黑人| 久久人妻av系列| 亚洲伊人色综图| 18禁裸乳无遮挡免费网站照片 | a级毛片在线看网站| 50天的宝宝边吃奶边哭怎么回事| 日韩精品青青久久久久久| a级毛片在线看网站| 午夜日韩欧美国产| 他把我摸到了高潮在线观看| 亚洲人成伊人成综合网2020| 麻豆成人av在线观看| 免费高清在线观看日韩| 法律面前人人平等表现在哪些方面| 亚洲中文字幕一区二区三区有码在线看 | 超碰成人久久| 久久中文看片网| 成人国语在线视频| 夜夜看夜夜爽夜夜摸| 久久久国产精品麻豆| 国产高清videossex| 国产野战对白在线观看| 电影成人av| 中文字幕高清在线视频| 啦啦啦免费观看视频1| 国产亚洲av嫩草精品影院| 亚洲一区二区三区不卡视频| 国产精品,欧美在线| 国产单亲对白刺激| 999久久久精品免费观看国产| 国产亚洲精品久久久久5区| 人人澡人人妻人| 此物有八面人人有两片| 亚洲狠狠婷婷综合久久图片| 丰满人妻熟妇乱又伦精品不卡| 亚洲国产欧美日韩在线播放| 国产伦人伦偷精品视频| 成人18禁高潮啪啪吃奶动态图| av天堂在线播放| 中国美女看黄片| 亚洲一区二区三区不卡视频| 亚洲成人国产一区在线观看| 国产精品秋霞免费鲁丝片| 精品少妇一区二区三区视频日本电影| 国产av精品麻豆| 91大片在线观看| 免费少妇av软件| 欧美人与性动交α欧美精品济南到| 亚洲国产精品合色在线| 亚洲熟妇中文字幕五十中出| 在线观看免费午夜福利视频| 在线播放国产精品三级| 国产精品一区二区免费欧美| 在线观看舔阴道视频| 亚洲第一青青草原| 久久久久国产一级毛片高清牌| 亚洲国产欧美日韩在线播放| 久久久久久久久久久久大奶| 午夜日韩欧美国产| 99精品欧美一区二区三区四区| 在线av久久热| 亚洲中文字幕日韩| 国产高清有码在线观看视频 | 99国产极品粉嫩在线观看| 亚洲五月婷婷丁香| 国产免费男女视频| 天堂动漫精品| 久久伊人香网站| 亚洲精品一区av在线观看| 亚洲黑人精品在线| 久久国产亚洲av麻豆专区| 日本 av在线| 精品无人区乱码1区二区| 国产视频一区二区在线看| а√天堂www在线а√下载| 中文字幕av电影在线播放| 国产高清激情床上av| 操出白浆在线播放| 又紧又爽又黄一区二区| 天天添夜夜摸| 中文字幕久久专区| 久热这里只有精品99| 国产成人一区二区三区免费视频网站| 美女国产高潮福利片在线看| 88av欧美| www.精华液| 此物有八面人人有两片| 黄色视频,在线免费观看| 国产精品一区二区免费欧美| 国产av一区二区精品久久| 久99久视频精品免费| 一级黄色大片毛片| 亚洲av熟女| 午夜日韩欧美国产| АⅤ资源中文在线天堂| 国产xxxxx性猛交| 男人操女人黄网站| 国产精品亚洲一级av第二区| 国产成人精品无人区| 9191精品国产免费久久| 啦啦啦 在线观看视频| 制服人妻中文乱码| 日韩欧美一区视频在线观看| 国产亚洲欧美98| 精品电影一区二区在线| 午夜免费激情av| 久热爱精品视频在线9| 欧美日韩亚洲综合一区二区三区_| 啦啦啦观看免费观看视频高清 | 色av中文字幕| 12—13女人毛片做爰片一| 久久久久九九精品影院| 麻豆成人av在线观看| 亚洲欧美一区二区三区黑人| 9191精品国产免费久久| 亚洲精品国产精品久久久不卡| 制服人妻中文乱码| 男女下面进入的视频免费午夜 | av视频在线观看入口| 亚洲 欧美一区二区三区| 亚洲av成人不卡在线观看播放网| 99在线视频只有这里精品首页| 国产精品免费视频内射| 亚洲aⅴ乱码一区二区在线播放 | 18禁美女被吸乳视频| 亚洲专区字幕在线| 亚洲狠狠婷婷综合久久图片| 一区二区三区国产精品乱码| 亚洲久久久国产精品| 亚洲天堂国产精品一区在线| 婷婷六月久久综合丁香| 国产亚洲欧美精品永久| 一级a爱片免费观看的视频| 在线av久久热| 亚洲自拍偷在线| 欧美色欧美亚洲另类二区 | 好男人在线观看高清免费视频 | 纯流量卡能插随身wifi吗| 国产精华一区二区三区| 亚洲熟女毛片儿| 嫁个100分男人电影在线观看| 老司机在亚洲福利影院| 男女下面插进去视频免费观看| 欧美激情久久久久久爽电影 | 咕卡用的链子| 一个人免费在线观看的高清视频| 久久国产亚洲av麻豆专区| 国产亚洲欧美在线一区二区| 精品一区二区三区av网在线观看| 国产私拍福利视频在线观看| 多毛熟女@视频| 黄色女人牲交| 国产片内射在线| 午夜福利一区二区在线看| 国内久久婷婷六月综合欲色啪| 精品人妻在线不人妻| 午夜久久久在线观看| 午夜a级毛片| 亚洲avbb在线观看| 欧美在线黄色| 校园春色视频在线观看| 国产又爽黄色视频| 国产亚洲欧美精品永久| 在线观看免费视频网站a站| 国内久久婷婷六月综合欲色啪| 国产av一区二区精品久久| svipshipincom国产片| 免费av毛片视频| 久久人妻av系列| 中文字幕人妻丝袜一区二区| 50天的宝宝边吃奶边哭怎么回事| 国产精品一区二区三区四区久久 | 欧美色视频一区免费| 精品一区二区三区视频在线观看免费| 亚洲国产精品成人综合色| 国产亚洲精品一区二区www| 欧美日韩亚洲综合一区二区三区_| 在线观看舔阴道视频| 免费高清在线观看日韩| 欧美日本视频| 亚洲七黄色美女视频| 手机成人av网站| 69av精品久久久久久| 国产精品98久久久久久宅男小说| 国产精品久久久久久人妻精品电影| 亚洲午夜精品一区,二区,三区| 两个人视频免费观看高清| 91九色精品人成在线观看| 嫩草影院精品99| 国产伦一二天堂av在线观看| 午夜亚洲福利在线播放| 在线十欧美十亚洲十日本专区| 狠狠狠狠99中文字幕| 久久精品亚洲熟妇少妇任你| 亚洲国产精品成人综合色| 欧美 亚洲 国产 日韩一| 搡老妇女老女人老熟妇| 国产精品乱码一区二三区的特点 | 啦啦啦韩国在线观看视频| 亚洲精品美女久久av网站| 成人欧美大片| 亚洲av日韩精品久久久久久密| 一本综合久久免费| cao死你这个sao货| 欧美成人午夜精品| 99riav亚洲国产免费| 亚洲精品国产区一区二| 18禁黄网站禁片午夜丰满| av天堂在线播放| 欧美成人免费av一区二区三区| 男女做爰动态图高潮gif福利片 | 精品国产超薄肉色丝袜足j| 欧美成狂野欧美在线观看| 女人被躁到高潮嗷嗷叫费观| 操出白浆在线播放| 国产亚洲精品第一综合不卡| 中文字幕色久视频| 一a级毛片在线观看| 国产高清videossex| 长腿黑丝高跟| 亚洲国产日韩欧美精品在线观看 | 亚洲成人精品中文字幕电影| 午夜福利成人在线免费观看| 国产av一区二区精品久久| 国产日韩一区二区三区精品不卡| 99热只有精品国产|