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

    基于Nitsche方法與擬牛頓求解的二維接觸問(wèn)題等幾何分析

    2021-11-02 12:12:42胡清元沈莞薔蔣芳芳
    關(guān)鍵詞:列式算例邊界條件

    胡清元, 沈莞薔, 蔣芳芳

    (江南大學(xué) 理學(xué)院,無(wú)錫 214122)

    1 引 言

    接觸問(wèn)題是工業(yè)工程中常見(jiàn)的一類(lèi)力學(xué)問(wèn)題[1]。在接觸問(wèn)題分析方法中,解析方法精度高但適用范圍受限。有限元方法適用范圍廣,但采用的近似網(wǎng)格為接觸邊界帶來(lái)嚴(yán)重離散誤差,迭代產(chǎn)生的單元畸變可能導(dǎo)致迭代難以收斂。等幾何分析[2]方法采用非均勻有理B樣條(NURBS)為基函數(shù)建模,可精確地描述結(jié)構(gòu)邊界,適合分析接觸問(wèn)題。

    然而,當(dāng)前等幾何接觸分析常見(jiàn)方法存在一定不足。等幾何罰函數(shù)法[3]采用類(lèi)似彈簧剛度的罰系數(shù)支撐接觸面,但問(wèn)題求解受系數(shù)取值影響很大。等幾何分析中廣泛采用的Mortar方法[4]通過(guò)增加虛擬公共邊界施加接觸條件,該方法將接觸面虛擬自由度增加到系統(tǒng)總自由度中,因此總自由度數(shù)隨著接觸面演化而不斷變化。此外,等幾何配點(diǎn)法[5]也用于接觸問(wèn)題分析。半解析方法[6]將基本解與快速數(shù)值求解相結(jié)合,具有速度快和計(jì)算精度高的特點(diǎn),在等幾何分析中的應(yīng)用有待進(jìn)一步發(fā)展。

    近年來(lái),基于Nitsche方法的等幾何分析受到持續(xù)關(guān)注。Nitsche方法最早用于有限元位移邊界條件的施加[7],之后應(yīng)用到域分解問(wèn)題[8]和接觸問(wèn)題[9],以及后續(xù)的摩擦大變形接觸問(wèn)題[10]。Apostolatos等[11]綜合對(duì)比了多種域分解方法的數(shù)值性能,結(jié)果表明Nitsche方法的綜合表現(xiàn)最為穩(wěn)定。Hu等[12]首次采用Nitsche方法還原接觸邊界條件,該列式的優(yōu)勢(shì)是不增加系統(tǒng)自由度,且有完備的數(shù)學(xué)理論支撐[9]。然而,標(biāo)準(zhǔn)的Nitsche列式需指定罰系數(shù),該系數(shù)通常通過(guò)求解界面上的廣義特征值問(wèn)題得到。但對(duì)于接觸問(wèn)題,接觸面隨著迭代不斷變化,若需在每個(gè)迭代步前都求解廣義特征值問(wèn)題,將大幅度增加計(jì)算量。

    在方程求解方面,由于接觸問(wèn)題的強(qiáng)非線(xiàn)性,一般采用的Newton-Raphson方法需要對(duì)列式進(jìn)行線(xiàn)性化處理,推導(dǎo)切線(xiàn)剛度陣,并在每一步迭代中對(duì)大規(guī)模矩陣求逆。這不僅需要繁雜的推導(dǎo),也加重了迭代計(jì)算量。實(shí)際上,擬牛頓方法也可用于非線(xiàn)性問(wèn)題的求解,其核心在于采用了近似的切線(xiàn)剛度陣,并可采取BFGS逆更新格式以避免矩陣求逆。Matthies等[13]結(jié)合擬牛頓格式與步長(zhǎng)搜索,計(jì)算了材料和幾何非線(xiàn)性問(wèn)題。Laursen等[14]進(jìn)一步采用擬牛頓方法求解接觸問(wèn)題,為了保證收斂,當(dāng)?shù)环€(wěn)定時(shí)需計(jì)算切線(xiàn)剛度陣。Gabriel等[15]將第一次迭代求得的解作為第二次迭代的初始近似解。然而,Nitsche方法直接將所有可能產(chǎn)生接觸的界面納入其接觸列式,接觸面的更新會(huì)直接反映在Nitsche接觸列式當(dāng)中。因此,如何在接觸面劇烈演化和擬牛頓迭代不穩(wěn)定時(shí)進(jìn)行簡(jiǎn)單有效地修正,是有待進(jìn)一步研究的問(wèn)題。

    2 等幾何接觸問(wèn)題

    2.1 等幾何分析簡(jiǎn)述

    等幾何分析本質(zhì)上是采用NURBS為形函數(shù)的等參有限元方法,采用如下方式建模幾何場(chǎng),

    (1)

    (2)

    式中uA為控制點(diǎn)位移自由度。

    由虛位移原理,得有限元弱形式

    a(u,v)=l(v)

    (3)

    式中a(u,v)為雙線(xiàn)性項(xiàng),l(v)為線(xiàn)性項(xiàng)。為了便于推導(dǎo),式(3)忽略了位移邊界條件。

    2.2 接觸問(wèn)題與接觸條件

    不失一般性,考慮兩彈性體Ω1與Ω2發(fā)生接觸,定義接觸面分別為Γ1和Γ2。假定Γ1上點(diǎn)x1與Γ2上點(diǎn)x2發(fā)生接觸,定義接觸面外法向單位矢量為

    (4)

    定義接觸間隙函數(shù)

    g(u)=(x2-x1)·n1-[[u]]n

    (5)

    式中(x2-x1)·n1表示當(dāng)前接觸間隙,在小變形問(wèn)題中認(rèn)為其保持不變,[[u]]n=(u1-u2)·n1表示兩點(diǎn)間的相對(duì)位移。

    將接觸面記作Γ,接觸條件可寫(xiě)為

    (6)

    式中σn為其在n1方向上的分量。接觸條件1表示接觸面不能相互穿透,接觸條件2表示接觸力方向?yàn)閷⒔佑|面向外推,接觸條件3為互補(bǔ)條件。由于考慮的是無(wú)摩擦接觸,故令切向應(yīng)力σt=0。

    3 基于Nitsche方法的接觸列式

    3.1 Nitsche方法簡(jiǎn)述

    Nitsche方法的核心思想是將邊界條件以功的形式引入有限元弱形式,其列式為

    a(u,v)-Nitsche =l(v)

    (7)

    式中Nitsche部分通??杀硎緸檫吔绶e分的形式

    (8)

    式中f(σ)和f(v)分別為關(guān)于面力和位移的函數(shù),具體取法取決于施加的邊界條件類(lèi)型。

    3.2 Nitsche接觸列式

    從Nitsche統(tǒng)一列式[12]出發(fā),接觸列式可寫(xiě)為

    (9)

    式中Γ為所有可能產(chǎn)生接觸的界面,[b]R-為標(biāo)量b在R-上的投影

    (10)

    此外,γ為罰系數(shù),在文獻(xiàn)[16,17]的基礎(chǔ)上,本文采用經(jīng)驗(yàn)公式(11),

    (11)

    式中E為楊氏模量,p為形函數(shù)階次,d為問(wèn)題維度,h1和h2表示接觸面兩側(cè)的網(wǎng)格尺寸。

    在可能的接觸面Γ上,Nitsche列式通過(guò)投影算子判斷積分點(diǎn)上的函數(shù)值來(lái)還原接觸條件[9]。

    (1) 當(dāng)接觸未發(fā)生即g>0時(shí),接觸條件3要求σn=0,此時(shí)[σn+γg]R-=0。

    (2) 當(dāng)接觸發(fā)生即g=0時(shí),接觸條件2要求σn≤0,此時(shí)[σn+γg]R-=σn。

    (3) 當(dāng)g<0時(shí),違背接觸條件1,通過(guò)大數(shù)γ將接觸反力放大,從而將接觸面彈回。

    Nitsche接觸列式的物理意義是利用接觸力在虛位移v上做的虛功修正有限元弱形式。其中,接觸力主要為接觸面Γ上面力σn,以及當(dāng)接觸面發(fā)生穿透時(shí)接觸間隙g通過(guò)罰系數(shù)γ所得接觸反力。

    4 擬牛頓迭代求解格式

    4.1 擬牛頓方法與BFGS逆更新公式

    接觸問(wèn)題為強(qiáng)非線(xiàn)性問(wèn)題,接觸問(wèn)題非線(xiàn)性列式的Newton-Raphson迭代求解格式為

    KT(ui)·Δui + 1=R(ui)

    (12)

    式中KT為切線(xiàn)剛度陣,通常需要通過(guò)對(duì)列式進(jìn)行線(xiàn)性化處理得到;余量計(jì)算方式為

    R(ui)=l(v)-a(ui,v)+

    (13)

    Newton-Raphson迭代格式需要進(jìn)行線(xiàn)性化推導(dǎo),更重要的是,每次迭代都需要求解切線(xiàn)剛度陣的逆

    (14)

    這嚴(yán)重增加了計(jì)算量和分析時(shí)間。

    (15)

    并通過(guò)BFGS公式直接更新矩陣的逆,

    (16)

    4.2 迭代初始化與接觸修正

    (17)

    (18)

    其中

    (19)

    具體修正策略為,在一定迭代次數(shù)(本文設(shè)置為8次)后,當(dāng)本次余量范數(shù)大于前一次時(shí),認(rèn)為接觸面演化較為劇烈,BFGS的自修正特性不足以及時(shí)反映接觸面變化,此時(shí)需將BFGS的迭代內(nèi)核替換為式(18)所得割線(xiàn)剛度陣,從而保證加速收斂。

    5 算 例

    5.1 Taylor接觸

    Taylor接觸算例[18]是接觸問(wèn)題的分片實(shí)驗(yàn),該算例重點(diǎn)關(guān)注非協(xié)調(diào)網(wǎng)格能否均勻地傳遞接觸壓力,用于檢驗(yàn)接觸列式的正確性。如圖1所示,兩物體楊氏模量E=109Pa,泊松比υ=0.3,上方物體受均布?jí)毫=108N/m2,兩模型網(wǎng)格劃分不同。算例解析解為

    圖1 Taylor接觸算例

    ux=0.03xm,uy=-0.1ym

    σx x=0 Pa,σy y=-100 Pa,σx y=0 Pa

    接觸迭代記錄列入表1,經(jīng)過(guò)6次擬牛頓迭代后系統(tǒng)收斂。采用本文列式計(jì)算所得位移uy結(jié)果如圖2所示,結(jié)果表明,列式可以得到正確的位移響應(yīng)。應(yīng)力σy y的相對(duì)誤差如圖3所示,最大誤差為0.13%。需要說(shuō)明的是,Nitsche方法采用積分的形式將邊界條件弱施加在系統(tǒng)中,對(duì)接觸問(wèn)題,接觸力通過(guò)高斯積分點(diǎn)傳遞,因此無(wú)法嚴(yán)格地通過(guò)Taylor接觸實(shí)驗(yàn)。但隨著網(wǎng)格加密,誤差減小,本文列式可以弱通過(guò)Taylor接觸實(shí)驗(yàn)。

    圖2 位移云圖

    圖3 應(yīng)力相對(duì)誤差

    表1 Taylor接觸迭代記錄

    5.2 Hertz接觸

    Hertz接觸算例是接觸問(wèn)題中的經(jīng)典點(diǎn)接觸算例,原問(wèn)題要求兩個(gè)楊氏模量不同和半徑不同的三維圓柱體接觸。本算例將原問(wèn)題退化,考慮一個(gè)二維彈性圓面和剛性平臺(tái)(即彈性模量和半徑均為無(wú)窮大)的接觸問(wèn)題。如圖4所示,上方圓面彈性模量E=0.02×109Pa,泊松比ν=0.1,半徑r=0.2 m,承受體力gz=-1.3×106N/m3。本算例理論解為[19],接觸力壓力為2.2918×106Pa,半接觸面寬度為0.0454 m。

    圖4 Hertz接觸算例

    采用8×8網(wǎng)格計(jì)算該算例,接觸迭代記錄列入表2,其中在第13次迭代時(shí)采用了割線(xiàn)剛度陣修正,修正后的余量范數(shù)大幅度降低,加速了收斂過(guò)程。計(jì)算所得豎向位移uy云圖如圖5所示,結(jié)果表明下方剛體可以有效地支撐上方彈性體,二者之間產(chǎn)生一段平直的接觸面。在不同網(wǎng)格劃分下計(jì)算所得接觸面寬度與接觸應(yīng)力如圖6所示,其中黑色實(shí)線(xiàn)為理論解,結(jié)果表明隨著網(wǎng)格加密,本文列式可以得到較為精確的解答。

    圖5 位移云圖

    圖6 不同網(wǎng)格下接觸面寬度與接觸應(yīng)力

    表2 Hertz接觸迭代記錄

    仍采用8×8網(wǎng)格,使用單線(xiàn)程計(jì)算對(duì)比。等幾何方法計(jì)算總耗時(shí)323 s,接觸應(yīng)力最大誤差 9.91%;使用同樣的網(wǎng)格劃分,有限元計(jì)算總耗時(shí)213 s,接觸應(yīng)力最大誤差28.7%。傳統(tǒng)有限元基函數(shù)構(gòu)造簡(jiǎn)單,而等幾何分析中的基函數(shù)需遞推計(jì)算,因此等幾何分析耗時(shí)較長(zhǎng)。但值得指出的是,有限元前處理時(shí)劃分網(wǎng)格的時(shí)間并沒(méi)有統(tǒng)計(jì)在內(nèi),而等幾何分析直接采用幾何模型計(jì)算。并且,在本算例中等幾何方法僅用較為粗糙的網(wǎng)格即可精確描述接觸邊界,迭代過(guò)程穩(wěn)定,最終得到較精確的結(jié)果,而這在傳統(tǒng)有限元中是難以做到的。

    本算例的接觸表面是光滑的,對(duì)粗糙表面的接觸問(wèn)題,其接觸表面為折線(xiàn)狀。因?yàn)闊o(wú)法體現(xiàn)等幾何方法描述曲線(xiàn)邊界時(shí)的精度優(yōu)勢(shì),此時(shí)的等幾何方法作為一種高階有限元方法,其計(jì)算精度應(yīng)略?xún)?yōu)于傳統(tǒng)有限元。

    6 結(jié) 論

    本文在等幾何分析框架內(nèi),基于Nitsche方法推導(dǎo)了二維小變形無(wú)摩擦接觸問(wèn)題列式。非線(xiàn)性列式采用擬牛頓迭代格式求解,由BFGS公式直接更新近似切線(xiàn)剛度陣的逆。本文給出了一個(gè)Nitsche列式中罰系數(shù)經(jīng)驗(yàn)公式,提出了基于Nitsche界面耦合列式的擬牛頓迭代初始化方法,并采用割線(xiàn)剛度陣解決由接觸面變化導(dǎo)致的BFGS更新能力不足問(wèn)題。

    所提出的接觸分析方法具有較強(qiáng)的適用性,在粗糙網(wǎng)格下依然可以精確描述接觸邊界,列式推導(dǎo)簡(jiǎn)單,計(jì)算量小。算例表明,基于Nitsche方法的接觸列式以積分形式近似地施加接觸條件,基于擬牛頓方法和割線(xiàn)剛度陣修正的迭代過(guò)程能夠較好收斂,最終可以得到較精確的接觸面寬度和接觸面壓力解答。

    猜你喜歡
    列式算例邊界條件
    一類(lèi)帶有Stieltjes積分邊界條件的分?jǐn)?shù)階微分方程邊值問(wèn)題正解
    帶有積分邊界條件的奇異攝動(dòng)邊值問(wèn)題的漸近解
    準(zhǔn)確審題正確列式精確驗(yàn)證
    教師·下(2017年10期)2017-12-10 12:35:13
    每筐多裝多少
    基于振蕩能量的低頻振蕩分析與振蕩源定位(二)振蕩源定位方法與算例
    互補(bǔ)問(wèn)題算例分析
    基于CYMDIST的配電網(wǎng)運(yùn)行優(yōu)化技術(shù)及算例分析
    帶Robin邊界條件的2維隨機(jī)Ginzburg-Landau方程的吸引子
    燃煤PM10湍流聚并GDE方程算法及算例分析
    帶非齊次邊界條件的p—Laplacian方程正解的存在唯一性
    天美传媒精品一区二区| 免费观看在线日韩| tube8黄色片| 国产精品三级大全| 一本色道久久久久久精品综合| 交换朋友夫妻互换小说| 午夜免费男女啪啪视频观看| 欧美在线黄色| 日日啪夜夜爽| 美女主播在线视频| 亚洲美女黄色视频免费看| 边亲边吃奶的免费视频| 永久网站在线| 成年女人在线观看亚洲视频| 考比视频在线观看| 精品人妻在线不人妻| 精品亚洲成国产av| 国产色婷婷99| 少妇精品久久久久久久| 在线免费观看不下载黄p国产| 亚洲欧美日韩另类电影网站| 一区二区三区激情视频| 亚洲精品中文字幕在线视频| 高清在线视频一区二区三区| 亚洲国产毛片av蜜桃av| 午夜影院在线不卡| 天天躁狠狠躁夜夜躁狠狠躁| 一个人免费看片子| 新久久久久国产一级毛片| 一二三四在线观看免费中文在| 建设人人有责人人尽责人人享有的| 欧美成人午夜精品| 免费黄网站久久成人精品| 久久久国产欧美日韩av| 欧美成人午夜免费资源| 日本欧美视频一区| 国产成人一区二区在线| 欧美+日韩+精品| 看免费成人av毛片| 国产男女超爽视频在线观看| 夫妻性生交免费视频一级片| 永久免费av网站大全| 美女大奶头黄色视频| 免费观看a级毛片全部| 爱豆传媒免费全集在线观看| 亚洲情色 制服丝袜| 在线观看美女被高潮喷水网站| av国产精品久久久久影院| 啦啦啦在线免费观看视频4| av电影中文网址| 少妇人妻精品综合一区二区| 久久久精品区二区三区| 亚洲少妇的诱惑av| 久久久久久人妻| av线在线观看网站| 乱人伦中国视频| 免费少妇av软件| 亚洲,欧美精品.| 18禁国产床啪视频网站| 亚洲男人天堂网一区| 国产探花极品一区二区| 国产欧美亚洲国产| 亚洲,欧美,日韩| 精品人妻熟女毛片av久久网站| 国产成人精品福利久久| 老司机影院毛片| 中文天堂在线官网| 久久精品国产自在天天线| 久久久久久伊人网av| 亚洲综合色惰| 多毛熟女@视频| 啦啦啦中文免费视频观看日本| 在线 av 中文字幕| 亚洲精品第二区| 下体分泌物呈黄色| 亚洲精品中文字幕在线视频| 日本欧美视频一区| 国产成人一区二区在线| 国产精品.久久久| av.在线天堂| 天天躁夜夜躁狠狠久久av| 国产成人精品福利久久| 一级爰片在线观看| 亚洲精品国产av蜜桃| 看免费av毛片| 一本色道久久久久久精品综合| 一区二区三区精品91| 九草在线视频观看| 999精品在线视频| 久久国产亚洲av麻豆专区| 亚洲av.av天堂| 成人毛片a级毛片在线播放| 蜜桃在线观看..| 国产精品麻豆人妻色哟哟久久| 少妇人妻精品综合一区二区| 好男人视频免费观看在线| 午夜福利影视在线免费观看| 国产精品免费视频内射| 激情视频va一区二区三区| 老熟女久久久| 午夜福利一区二区在线看| 女性生殖器流出的白浆| 免费女性裸体啪啪无遮挡网站| 香蕉国产在线看| 亚洲综合色网址| 成人影院久久| 9191精品国产免费久久| 男人爽女人下面视频在线观看| 免费看av在线观看网站| 亚洲精品在线美女| 飞空精品影院首页| 欧美97在线视频| 久久久欧美国产精品| 国产探花极品一区二区| 成人毛片a级毛片在线播放| 日日撸夜夜添| 精品第一国产精品| 在线 av 中文字幕| 久热久热在线精品观看| 丝袜喷水一区| 五月开心婷婷网| 两性夫妻黄色片| 亚洲综合色惰| 制服人妻中文乱码| 人妻人人澡人人爽人人| 老司机亚洲免费影院| 在线观看www视频免费| 日本欧美视频一区| 不卡视频在线观看欧美| 在线天堂最新版资源| 午夜av观看不卡| 国产一区二区在线观看av| 亚洲精品国产色婷婷电影| 亚洲成人手机| 一本久久精品| 亚洲精品成人av观看孕妇| 999久久久国产精品视频| 国产xxxxx性猛交| 欧美最新免费一区二区三区| 欧美亚洲 丝袜 人妻 在线| 可以免费在线观看a视频的电影网站 | 最近最新中文字幕大全免费视频 | 美女午夜性视频免费| 亚洲成av片中文字幕在线观看 | 国产精品一区二区在线不卡| 91精品国产国语对白视频| 亚洲成人手机| 丰满饥渴人妻一区二区三| av在线app专区| 国产成人精品福利久久| 精品福利永久在线观看| 女人精品久久久久毛片| 看免费av毛片| 最近的中文字幕免费完整| 久久人人97超碰香蕉20202| 日韩欧美一区视频在线观看| 最近中文字幕2019免费版| 美女主播在线视频| 亚洲av福利一区| 王馨瑶露胸无遮挡在线观看| 秋霞在线观看毛片| 亚洲国产看品久久| 香蕉丝袜av| 亚洲欧洲国产日韩| 校园人妻丝袜中文字幕| 一级a爱视频在线免费观看| 国产在线一区二区三区精| 午夜福利,免费看| 午夜福利在线观看免费完整高清在| 69精品国产乱码久久久| 成人影院久久| 午夜激情av网站| 极品少妇高潮喷水抽搐| 国产精品二区激情视频| 国产精品 国内视频| 国产亚洲欧美精品永久| 黄色怎么调成土黄色| 精品少妇一区二区三区视频日本电影 | 最新中文字幕久久久久| 男女高潮啪啪啪动态图| 99九九在线精品视频| 纯流量卡能插随身wifi吗| 午夜免费男女啪啪视频观看| 精品福利永久在线观看| 国产一区亚洲一区在线观看| 久久午夜福利片| 国产精品久久久久久精品古装| 少妇的逼水好多| 夫妻性生交免费视频一级片| 一区在线观看完整版| 成人18禁高潮啪啪吃奶动态图| 国产综合精华液| 嫩草影院入口| 韩国av在线不卡| 在线观看www视频免费| 亚洲精品一区蜜桃| 久热这里只有精品99| 欧美日韩一区二区视频在线观看视频在线| 两性夫妻黄色片| 久久综合国产亚洲精品| 久久久久网色| 欧美在线黄色| 国产免费福利视频在线观看| 日本午夜av视频| 国产精品av久久久久免费| 亚洲色图综合在线观看| 人妻 亚洲 视频| 成人漫画全彩无遮挡| 国产精品一区二区在线观看99| 国产一区有黄有色的免费视频| 亚洲国产成人一精品久久久| 欧美黄色片欧美黄色片| 丝袜人妻中文字幕| 欧美日韩亚洲国产一区二区在线观看 | 美女大奶头黄色视频| 精品久久蜜臀av无| 超碰成人久久| 丁香六月天网| 十八禁网站网址无遮挡| 这个男人来自地球电影免费观看 | 亚洲国产精品一区三区| 亚洲精品成人av观看孕妇| 久久这里有精品视频免费| 母亲3免费完整高清在线观看 | 最近中文字幕高清免费大全6| 妹子高潮喷水视频| 国产伦理片在线播放av一区| 日产精品乱码卡一卡2卡三| 国产人伦9x9x在线观看 | 婷婷色麻豆天堂久久| 夫妻午夜视频| 超碰成人久久| 天天影视国产精品| 亚洲av日韩在线播放| 欧美精品国产亚洲| 99久久综合免费| 亚洲国产精品国产精品| 女人久久www免费人成看片| 成人手机av| 亚洲,一卡二卡三卡| 丝袜在线中文字幕| 亚洲视频免费观看视频| 成人毛片60女人毛片免费| 啦啦啦在线观看免费高清www| 国产精品久久久久成人av| 巨乳人妻的诱惑在线观看| 制服丝袜香蕉在线| 99久久中文字幕三级久久日本| 国产成人a∨麻豆精品| 成年av动漫网址| 日韩熟女老妇一区二区性免费视频| 亚洲av成人精品一二三区| 观看美女的网站| 亚洲成人一二三区av| 黄色配什么色好看| 日本午夜av视频| 高清不卡的av网站| 熟妇人妻不卡中文字幕| 90打野战视频偷拍视频| 国产极品粉嫩免费观看在线| 9191精品国产免费久久| 亚洲图色成人| 国产白丝娇喘喷水9色精品| 少妇人妻 视频| videossex国产| 精品亚洲成a人片在线观看| 日本-黄色视频高清免费观看| 免费久久久久久久精品成人欧美视频| a 毛片基地| 久久这里有精品视频免费| 超色免费av| 亚洲精品av麻豆狂野| 亚洲少妇的诱惑av| 精品国产露脸久久av麻豆| 国产一区二区在线观看av| 精品久久久精品久久久| 国产免费又黄又爽又色| 日韩人妻精品一区2区三区| 亚洲av福利一区| 女性被躁到高潮视频| 男女无遮挡免费网站观看| 成人手机av| 男男h啪啪无遮挡| 亚洲图色成人| 制服丝袜香蕉在线| 午夜久久久在线观看| 国产精品国产三级国产专区5o| 深夜精品福利| 美女主播在线视频| 天天躁日日躁夜夜躁夜夜| 波多野结衣av一区二区av| 交换朋友夫妻互换小说| 黑人欧美特级aaaaaa片| 在线观看www视频免费| 看十八女毛片水多多多| 成人国产麻豆网| 精品少妇内射三级| 亚洲av在线观看美女高潮| 成年美女黄网站色视频大全免费| 国产一区亚洲一区在线观看| 亚洲精品国产av成人精品| 午夜影院在线不卡| 成人毛片a级毛片在线播放| 青春草亚洲视频在线观看| 亚洲av综合色区一区| 精品国产露脸久久av麻豆| 伦理电影免费视频| 欧美亚洲日本最大视频资源| 人人妻人人爽人人添夜夜欢视频| 男人操女人黄网站| 少妇猛男粗大的猛烈进出视频| 男女午夜视频在线观看| 亚洲男人天堂网一区| 天天躁夜夜躁狠狠久久av| 久久热在线av| 久久久久久久大尺度免费视频| 精品人妻熟女毛片av久久网站| 肉色欧美久久久久久久蜜桃| 各种免费的搞黄视频| 青春草亚洲视频在线观看| 精品国产超薄肉色丝袜足j| 午夜福利,免费看| 街头女战士在线观看网站| 欧美另类一区| 中文字幕亚洲精品专区| 天堂中文最新版在线下载| 日韩av免费高清视频| 欧美日韩精品网址| 日韩在线高清观看一区二区三区| 日韩精品有码人妻一区| 丝袜脚勾引网站| 国产在线视频一区二区| 七月丁香在线播放| 99热全是精品| 午夜福利在线免费观看网站| 国产精品亚洲av一区麻豆 | www.熟女人妻精品国产| 国产成人欧美| 日韩熟女老妇一区二区性免费视频| 色吧在线观看| 麻豆乱淫一区二区| 亚洲久久久国产精品| 日韩不卡一区二区三区视频在线| 国语对白做爰xxxⅹ性视频网站| 一区二区日韩欧美中文字幕| 欧美精品一区二区免费开放| 人妻人人澡人人爽人人| 18+在线观看网站| 欧美bdsm另类| 国语对白做爰xxxⅹ性视频网站| 久热这里只有精品99| 如日韩欧美国产精品一区二区三区| 一二三四中文在线观看免费高清| 久久精品aⅴ一区二区三区四区 | 国产男女内射视频| 亚洲欧美一区二区三区国产| 亚洲精品美女久久av网站| 成人国产av品久久久| 午夜福利在线免费观看网站| 久久久久国产一级毛片高清牌| 亚洲精品成人av观看孕妇| 国产白丝娇喘喷水9色精品| 18禁裸乳无遮挡动漫免费视频| 免费黄色在线免费观看| 国产精品二区激情视频| 国产97色在线日韩免费| 国产黄色视频一区二区在线观看| 人妻一区二区av| av免费观看日本| 我要看黄色一级片免费的| 国产又色又爽无遮挡免| 日韩精品有码人妻一区| 丝袜脚勾引网站| 久久久久精品人妻al黑| 国产一区亚洲一区在线观看| 99热全是精品| 亚洲成色77777| 国产男女内射视频| 女的被弄到高潮叫床怎么办| 午夜福利在线免费观看网站| 国产精品免费大片| 国产精品久久久久久久久免| 老司机影院毛片| 又粗又硬又长又爽又黄的视频| 国产精品久久久久久精品古装| 寂寞人妻少妇视频99o| 亚洲av电影在线观看一区二区三区| 在线亚洲精品国产二区图片欧美| 成人亚洲欧美一区二区av| 男女高潮啪啪啪动态图| 伊人久久大香线蕉亚洲五| 一级片免费观看大全| 久久综合国产亚洲精品| 国产亚洲精品第一综合不卡| freevideosex欧美| 国产免费视频播放在线视频| 国产日韩欧美亚洲二区| 国产极品天堂在线| 99九九在线精品视频| 另类亚洲欧美激情| 最黄视频免费看| 成年人午夜在线观看视频| 国产探花极品一区二区| 18禁裸乳无遮挡动漫免费视频| 欧美成人午夜精品| 久久鲁丝午夜福利片| 国产精品免费大片| 成年av动漫网址| 中文天堂在线官网| 色播在线永久视频| 欧美变态另类bdsm刘玥| 十分钟在线观看高清视频www| 亚洲精品在线美女| 肉色欧美久久久久久久蜜桃| 青青草视频在线视频观看| 在线观看人妻少妇| 国产黄频视频在线观看| av电影中文网址| 久久久久精品人妻al黑| 亚洲精品国产一区二区精华液| 狠狠精品人妻久久久久久综合| 黄频高清免费视频| 亚洲成人手机| 久久久久精品人妻al黑| 亚洲欧洲日产国产| 高清在线视频一区二区三区| 只有这里有精品99| 99久久精品国产国产毛片| 大香蕉久久网| 成人影院久久| 天天影视国产精品| 国产成人免费观看mmmm| 亚洲av综合色区一区| 美女大奶头黄色视频| 人妻人人澡人人爽人人| 欧美+日韩+精品| 午夜久久久在线观看| 精品亚洲乱码少妇综合久久| 久久免费观看电影| 色播在线永久视频| 欧美黄色片欧美黄色片| 日本猛色少妇xxxxx猛交久久| 日日啪夜夜爽| 午夜av观看不卡| 97人妻天天添夜夜摸| 亚洲成人av在线免费| 天美传媒精品一区二区| 另类亚洲欧美激情| 精品久久久久久电影网| 亚洲av免费高清在线观看| 日韩欧美一区视频在线观看| 女人精品久久久久毛片| 亚洲欧美日韩另类电影网站| 亚洲欧美色中文字幕在线| 日韩中字成人| 亚洲欧美一区二区三区国产| 中文字幕人妻丝袜制服| 黄色毛片三级朝国网站| 国产精品不卡视频一区二区| 国产精品人妻久久久影院| 免费黄网站久久成人精品| 日本色播在线视频| 欧美精品亚洲一区二区| 日本免费在线观看一区| 9191精品国产免费久久| 国产精品二区激情视频| 丰满乱子伦码专区| 亚洲国产毛片av蜜桃av| 欧美bdsm另类| 有码 亚洲区| 一级毛片 在线播放| 性色av一级| 电影成人av| 亚洲一区中文字幕在线| 国产高清不卡午夜福利| 在线观看美女被高潮喷水网站| 精品视频人人做人人爽| 亚洲av.av天堂| 精品99又大又爽又粗少妇毛片| 亚洲精品,欧美精品| 亚洲国产色片| av.在线天堂| 男男h啪啪无遮挡| 亚洲精品美女久久久久99蜜臀 | av视频免费观看在线观看| 精品一区二区免费观看| 亚洲av综合色区一区| 午夜日本视频在线| 国产精品久久久久成人av| 午夜福利在线免费观看网站| 久久99热这里只频精品6学生| 精品人妻在线不人妻| 少妇的丰满在线观看| 97在线人人人人妻| 人人妻人人澡人人看| 超碰97精品在线观看| 丝袜人妻中文字幕| 超碰成人久久| 狠狠婷婷综合久久久久久88av| 99久久综合免费| www.av在线官网国产| 美女xxoo啪啪120秒动态图| 巨乳人妻的诱惑在线观看| 成人国产av品久久久| 久久久久久伊人网av| 亚洲精品中文字幕在线视频| 免费人妻精品一区二区三区视频| 99热全是精品| 国产精品 国内视频| 美女午夜性视频免费| 亚洲精品aⅴ在线观看| a级毛片黄视频| 男女下面插进去视频免费观看| av女优亚洲男人天堂| 高清黄色对白视频在线免费看| 天天躁日日躁夜夜躁夜夜| 国产精品一区二区在线观看99| 久久精品国产综合久久久| 欧美日韩综合久久久久久| 丰满少妇做爰视频| 国产有黄有色有爽视频| 久久午夜福利片| 国产精品av久久久久免费| 人妻人人澡人人爽人人| 大陆偷拍与自拍| 性色avwww在线观看| 中文字幕人妻丝袜制服| 午夜91福利影院| 黄频高清免费视频| 欧美人与性动交α欧美软件| 99国产综合亚洲精品| 丝袜美足系列| 看十八女毛片水多多多| 精品少妇内射三级| 亚洲成人av在线免费| 国语对白做爰xxxⅹ性视频网站| 精品国产乱码久久久久久男人| 好男人视频免费观看在线| 看十八女毛片水多多多| 午夜日韩欧美国产| 久久国产精品大桥未久av| 亚洲欧美一区二区三区久久| videosex国产| 亚洲成人手机| 一级毛片黄色毛片免费观看视频| 99九九在线精品视频| 肉色欧美久久久久久久蜜桃| 春色校园在线视频观看| 中文字幕另类日韩欧美亚洲嫩草| 最近手机中文字幕大全| 男女免费视频国产| 777久久人妻少妇嫩草av网站| 亚洲精品自拍成人| 日韩一区二区三区影片| 两性夫妻黄色片| 中文字幕制服av| 久久久精品国产亚洲av高清涩受| 乱人伦中国视频| 成人国语在线视频| 精品一品国产午夜福利视频| 黑人猛操日本美女一级片| 国产又色又爽无遮挡免| 青青草视频在线视频观看| 日韩,欧美,国产一区二区三区| 精品久久蜜臀av无| 一区福利在线观看| 国产精品一区二区在线不卡| 捣出白浆h1v1| 午夜91福利影院| 天堂中文最新版在线下载| 曰老女人黄片| 超色免费av| 女人被躁到高潮嗷嗷叫费观| 最新中文字幕久久久久| 新久久久久国产一级毛片| 亚洲国产欧美在线一区| 在线观看免费视频网站a站| 国产精品一区二区在线不卡| 国产在线免费精品| 一区二区三区乱码不卡18| 一级毛片黄色毛片免费观看视频| 九色亚洲精品在线播放| 只有这里有精品99| 晚上一个人看的免费电影| 国产国语露脸激情在线看| 啦啦啦啦在线视频资源| freevideosex欧美| 美女视频免费永久观看网站| 免费高清在线观看视频在线观看| 一本久久精品| 精品国产一区二区三区久久久樱花| 婷婷色综合www| 国产精品 国内视频| 久久久精品区二区三区| 视频区图区小说| 男人添女人高潮全过程视频| 久久热在线av| 久久国产精品大桥未久av| 国产精品久久久久久精品电影小说| 亚洲精品自拍成人| 欧美 日韩 精品 国产| 最近最新中文字幕免费大全7| 欧美在线黄色| 天堂8中文在线网| 岛国毛片在线播放| 在现免费观看毛片| 亚洲精品久久久久久婷婷小说| 欧美日韩亚洲高清精品| 亚洲av.av天堂| 久久精品aⅴ一区二区三区四区 | 国产乱来视频区| 看非洲黑人一级黄片| 久久久a久久爽久久v久久| 日韩av免费高清视频| 欧美人与善性xxx|