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

    基于亞區(qū)結(jié)構(gòu)的人體二尖瓣的建模與仿真

    2015-10-13 01:24:38趙曉佳黃曉陽蘇茂龍
    關(guān)鍵詞:乳頭肌亞區(qū)瓣葉

    趙曉佳,宋 曉,黃曉陽*,蘇茂龍,鐘 琪

    (1.廈門大學(xué)信息科學(xué)與技術(shù)學(xué)院,福建 廈門361005;2.南陽師范學(xué)院計算機與信息技術(shù)學(xué)院,河南南陽473061;3.廈門大學(xué)附屬中山醫(yī)院超聲醫(yī)學(xué)科,福建廈門361004)

    基于亞區(qū)結(jié)構(gòu)的人體二尖瓣的建模與仿真

    趙曉佳1,宋 曉2,黃曉陽1*,蘇茂龍3,鐘 琪1

    (1.廈門大學(xué)信息科學(xué)與技術(shù)學(xué)院,福建 廈門361005;2.南陽師范學(xué)院計算機與信息技術(shù)學(xué)院,河南南陽473061;3.廈門大學(xué)附屬中山醫(yī)院超聲醫(yī)學(xué)科,福建廈門361004)

    通過建立精細(xì)的人體二尖瓣有限元模型,仿真前后瓣葉厚度不均勻的二尖瓣模型在心臟收縮初期的動態(tài)閉合過程,與前后瓣葉厚度均勻的二尖瓣模型的仿真結(jié)果相比較.仿真結(jié)果表明:兩種模型的瓣葉應(yīng)力分布不均勻,最大應(yīng)力值都出現(xiàn)在后瓣葉亞區(qū)與亞區(qū)之間的裂口位置,與臨床二尖瓣撕裂常發(fā)生位置一致,表明了仿真方法以及結(jié)果的正確性,瓣葉厚度不均勻模型中瓣葉應(yīng)力水平相對較低,更符合實際.最后根據(jù)臨床數(shù)據(jù),設(shè)置厚度不均勻模型中瓣環(huán)與乳頭肌的二維運動,再次進行仿真.仿真結(jié)果表明該精細(xì)模型瓣葉的閉合度比之前瓣環(huán)與乳頭肌固定不動模型的閉合度高,瓣葉應(yīng)力分布相對較均勻,更接近實際生理情況,為后續(xù)的病理等的仿真奠定了基礎(chǔ).

    二尖瓣;有限元模型;應(yīng)力分布;最大應(yīng)力

    二尖瓣是人體心臟的重要瓣膜組織之一,它位于左心房與左心室之間,由二尖瓣環(huán)、兩個瓣葉、腿索、乳頭肌等結(jié)構(gòu)組成[1].二尖瓣的瓣葉分為較大的前瓣葉和較小的后瓣葉,后瓣葉可細(xì)化為3個較窄的亞區(qū)結(jié)構(gòu)[2]:P1、P2、P3.二尖瓣的瓣葉與瓣環(huán)、腿索相連,腿索又與乳頭肌相連.根據(jù)Lam等[3]的分類方式將腿索分為3級:1級腿索連接于瓣葉的游離緣,2級腿索連接于瓣葉的粗糙部,3級腿索從室壁連接到后瓣葉的基底部.二尖瓣維持正常的生理職能,需要各個組成結(jié)構(gòu)的完整、乳頭肌的正常收縮及腿索的牽引[4].

    二尖瓣的生物力學(xué)是一個復(fù)雜的問題:結(jié)構(gòu)體的大變形,瓣膜組織的各向異性的非線性彈性行為以及心動周期內(nèi)的血液動力學(xué)的載荷等.建模仿真是分析瓣膜功能有效性的技術(shù)方法.隨著科技的發(fā)展、仿真技術(shù)的不斷完善,人們可以對瓣膜功能有更加精確的了解,可以選擇更適合的病理治療方式.通過對人體二尖瓣進行建模仿真,可以更好地認(rèn)識二尖瓣的力學(xué)機理,對后期分析二尖瓣關(guān)閉不全等病理結(jié)構(gòu)以及指導(dǎo)瓣膜修復(fù)術(shù)具有重要的意義.

    目前關(guān)于二尖瓣的研究分為結(jié)構(gòu)模型和流固耦合模型兩類.結(jié)構(gòu)模型是指模型中只有二尖瓣的結(jié)構(gòu)組件,不考慮血液;流固耦合模型是指二尖瓣的結(jié)構(gòu)與血液流體域相互作用的模型,它主要研究二尖瓣結(jié)構(gòu)在流體作用下的各種行為以及結(jié)構(gòu)對流體域的影響.其中結(jié)構(gòu)模型是流固耦合模型的基礎(chǔ),主要研究的學(xué)者有Kunzelman團隊[5-7]、Lim團隊[8]等.Kunzelman等[5]建立了最原始的人體二尖瓣模型,模擬了二尖瓣的閉合過程.Lim等[8]建立厚度均勻的膜結(jié)構(gòu)瓣葉,模擬正常二尖瓣工作機理,但該模型中瓣葉自由邊上無尖端,忽略了瓣葉形態(tài)的局部細(xì)節(jié).Kunzelman等[7]、羅云等[9]將二尖瓣的后瓣葉簡化為1個瓣尖,成功模擬出瓣葉的開閉過程,但并未根據(jù)實際生理情況將后瓣葉分為3個亞區(qū).這些研究都沒有模擬瓣環(huán)與乳頭肌的二維運動過程.本文研究的目的,就是為了比較更加精細(xì)的模型對仿真結(jié)果的影響.

    本文根據(jù)人體二尖瓣的實際生理結(jié)構(gòu),以二尖瓣的亞區(qū)結(jié)構(gòu)特征,建立比較精細(xì)的二尖瓣三維有限元模型:后瓣葉分為3個亞區(qū)結(jié)構(gòu),腿索分為游離邊緣腿索,粗糙帶以及基底部腿索.二尖瓣瓣膜厚度分為兩種情況:1)前后瓣葉厚度相同的厚度均勻模型; 2)前后瓣葉厚度不相同的厚度不均勻模型.兩種模型分別模擬其瓣葉在心臟收縮初期心房心室壓力差的作用下的閉合過程,并比較了兩種仿真模型的差異.在瓣葉厚度不均勻的模型基礎(chǔ)上,按照臨床數(shù)據(jù)及相關(guān)描述,設(shè)置瓣環(huán)與乳頭肌的二維運動過程,更加精細(xì)地模擬二尖瓣的閉合過程.

    1 方 法

    采用ANSXS建立二尖瓣有限元模型,并用LSDXNA求解器進行求解,得到二尖瓣模型的應(yīng)力隨時間的變化的結(jié)果,用LS-PREPOST后處理器生成可視化結(jié)果.

    1.1 二尖瓣幾何模型

    二尖瓣模型由瓣環(huán)、瓣葉、腿索和乳頭肌組成.本文根據(jù)文獻[3,10-11]建立的模型如圖1所示:xy平面的瓣環(huán)呈D字,乳頭肌簡化為2個固定的點,瓣葉從瓣環(huán)向z軸正方向延伸,后瓣葉細(xì)分為3個亞區(qū)結(jié)構(gòu).與瓣葉游離邊緣相連的腿索有24根,與后瓣葉基底部相連的腿索有2根,與粗糙帶相連的腿索分支有18根.

    圖1 人體二尖瓣幾何模型Fig.1 Geometric model of human mitral valve

    1.2 物理性質(zhì)

    二尖瓣的所有組織均設(shè)置為各向同性的彈性材料.瓣葉密度為1 040 kg/m3[10],泊松比設(shè)置為0.45,彈性模量設(shè)置為4 MPa[12];腿索的密度為1 040 kg/ m3[10],其應(yīng)力-應(yīng)變曲線如圖2所示[13].

    1.3 網(wǎng)格劃分

    在ANSXS中對二尖瓣模型進行網(wǎng)格劃分,采用自由網(wǎng)格劃分,共劃分為4 396個Belytschko殼單元;腿索橫截面積均設(shè)置為0.4 mm2[14],共劃分為760個桁架單元.模型中瓣膜的厚度分2種情況討論:1)瓣葉厚度均勻的模型:前后瓣葉厚度均設(shè)置為1.3 mm; 2)前后瓣葉厚度不相同的厚度不均勻模型:根據(jù)文獻[15-16],前瓣葉厚度設(shè)置為1.32 mm,后瓣葉厚度設(shè)置為1.26 mm.兩種情況下的網(wǎng)格劃分方式一樣,僅前后瓣葉厚度不同.

    圖2 二尖瓣腿索材料應(yīng)力-應(yīng)變曲線Fig.2 Stress-strain behavior of mitral valve chordae tendineae

    1.4 邊界與負(fù)載條件

    基于二尖瓣的面對稱結(jié)構(gòu),在對稱面上節(jié)點x軸方向位移為0.二尖瓣瓣環(huán)與乳頭肌在整個心動周期內(nèi)的位置相對固定,因此假設(shè)瓣環(huán)和乳頭肌固定不動,瓣環(huán)上節(jié)點的位移為0.

    左心室與左心房在單個心動周期內(nèi)壓力值的臨床數(shù)據(jù)[17]如圖3所示,由于二尖瓣在心臟收縮初期壓力斜坡(early systolic pressure ramp,ESPR)區(qū),瓣膜載荷會隨著時間動態(tài)增長至最大值,此時的二尖瓣較容易出現(xiàn)應(yīng)力集中現(xiàn)象,因此著重研究ESPR區(qū)瓣膜的應(yīng)力分布情況.通過測量圖3上同一時刻心室與心房的壓強,然后求得它們之間的差值即左心室房壓強差,這個左心室房壓強差就是施加在二尖瓣上的動態(tài)載荷,如圖4所示.

    圖3 心室心房壓力曲線Fig.3 Atrial and ventricular pressure profiles

    2 結(jié) 果

    2.1 二尖瓣瓣葉閉合

    分別模擬兩種模型在ESPR區(qū)二尖瓣瓣葉的閉合情況,仿真結(jié)果如圖5所示.

    由圖5可知,2個模型中的瓣葉均在左房室壓差以及腿索的共同作用下慢慢對合,最終實現(xiàn)前后瓣葉的接觸.兩個模型中瓣環(huán)的形狀沒有發(fā)生變化,模型的閉合程度基本一樣,無顯著差異.

    圖4 二尖瓣載荷曲線Fig.4 Transvalvular pressure load applied to the mitral valve

    圖5 二尖瓣瓣葉閉合情況Fig.5 Configuration of the valve at closure

    2.2 二尖瓣瓣葉應(yīng)力分布

    分別對二尖瓣在ESPR區(qū)應(yīng)力分布進行仿真,對應(yīng)的最大von Mises應(yīng)力分布如圖6所示.

    如圖所示,在ESPR區(qū)400 ms時刻,左房室壓差趨于穩(wěn)定,約16 k Pa,兩個模型的二尖瓣von Mises應(yīng)力均達(dá)到最大,最大應(yīng)力值均出現(xiàn)在亞區(qū)的裂口位置,與臨床上瓣葉容易撕裂的位置一致[18],與鐘琪等[19]的結(jié)果一致.表1列出了兩種模型的仿真結(jié)果對比情況.兩個模型應(yīng)力分布都不太均勻,均出現(xiàn)了應(yīng)力集中現(xiàn)象,兩個模型除了后瓣葉亞區(qū)間應(yīng)力最大外,前后瓣葉以及瓣葉與瓣環(huán)連合處的von Mises應(yīng)力也都比較大,但厚度不均勻模型的總體von Mises應(yīng)力水平比厚度均勻模型的von Mises應(yīng)力水平低,更不容易出現(xiàn)撕裂的狀況,與實際生理相符.

    表1 2種模型仿真結(jié)果對比Tab.1 Comparison of the results between two models

    2.3 瓣環(huán)和乳頭肌的運動

    根據(jù)文獻[20-23]對瓣環(huán)運動的描述,瓣環(huán)在xy平面的移動可以描述為:1)前瓣環(huán)長度增加11%,但是瓣環(huán)總長度縮短13%;2)瓣環(huán)橫向長度(即x軸長度)增加5%.

    根據(jù)文獻[16]描述,以及整個ESPR區(qū),瓣環(huán)與乳頭肌的位置是相對不動,乳頭肌沿x軸朝著前后瓣葉連合處移動0.85 mm,沿y軸朝著前瓣葉方向移動1.10 mm.

    在厚度不均勻模型基礎(chǔ)上,仿真瓣環(huán)以及乳頭肌在xy平面的移動以及瓣葉的閉合過程,仿真結(jié)果如圖7所示.

    圖6 400 ms時刻瓣葉最大von Mises應(yīng)力分布圖Fig.6 Maximum von Mises stress distribution on the valve leaflets at 400 ms

    如圖7所示,ESPR區(qū)初期,二尖瓣的von Mises應(yīng)力比較低,左心房室的壓力差約為0 kPa;在300 ms時刻,左心房室的壓差約為10 kPa,二尖瓣的von Mises應(yīng)力分布在0.042~0.253 MPa之間;在375 ms時刻,左心房室壓差約為16 kPa,二尖瓣的von Mises應(yīng)力主要在0.047~0.280 MPa,腿索與后瓣葉的連接處及前后瓣葉連合處出現(xiàn)應(yīng)力集中現(xiàn)象;在425 ms時刻,左心房室壓差約為16 kPa,二尖瓣的von Mises應(yīng)力分布在0.048~0.291 MPa,基底腿索與后瓣葉的連接處及前后瓣葉連合處出現(xiàn)應(yīng)力集中現(xiàn)象.總體看,二尖瓣瓣環(huán)在整個ESPR區(qū)的von Mises應(yīng)力分布相對比較均勻,最大von Mises應(yīng)力集中在腿索與瓣葉的連接處以及前后瓣葉之間的連合處,最大應(yīng)力值為0.485 MPa.在二尖瓣瓣葉的閉合過程中,二尖瓣瓣環(huán)的形態(tài)發(fā)生了很大的變化,由原來的D字形變成了一個類似橢圓形,瓣環(huán)在x軸方向長度變長,y軸方向長度縮短.除了形態(tài)的變化,瓣葉的閉合度也明顯發(fā)生了變化,與原先厚度不均勻模型相比,瓣葉的閉合程度明顯提高,與實際生理情況符合,如表2所示.

    圖7 不同時刻二尖瓣瓣葉von Mises應(yīng)力分布圖Fig.7 von Mises stress distribution on the valve leaflets at different time

    表2 瓣環(huán)與乳頭肌的運動對仿真結(jié)果的影響Tab.2 Impact of annulus and papillary muscles movement to the simulation results

    3 結(jié) 論

    為了研究比較精細(xì)模型對仿真結(jié)果的影響,本文參照實際生理結(jié)構(gòu),將二尖瓣的后瓣葉細(xì)分為3個亞區(qū),建立了瓣膜厚度分布均勻和分布不均勻2種模型,對二尖瓣瓣葉在ESPR區(qū)的閉合過程進行了仿真.仿真結(jié)果表明,在整個ESPR區(qū),瓣膜應(yīng)力最大值出現(xiàn)在400 ms左右,位置在P2亞區(qū)與P1、P3亞區(qū)的裂口位置,此時心房心室的壓力差最大,約16 k Pa.厚度不均勻模型的應(yīng)力水平比較低,不易出現(xiàn)撕裂現(xiàn)象,更接近實際的生理情況.在此基礎(chǔ)上,本文根據(jù)實際生理情況,設(shè)置瓣環(huán)與乳頭肌的運動,再次仿真瓣葉的閉合過程,仿真結(jié)果顯示前后瓣葉的連合處及腿索與后瓣葉連接處應(yīng)力較大,仿真的結(jié)果與臨床二尖瓣撕裂研究[18]相符合,同時與Kunzelman等[14]的模擬結(jié)果接近,前瓣葉的最大應(yīng)力值略小一些,約為0.300 MPa,后瓣葉的應(yīng)力值均約為0.200 MPa.瓣葉的閉合度更高,更符合實際生理情況.本文研究結(jié)果表明,建立瓣膜厚度不均勻模型,并在仿真過程中設(shè)置瓣環(huán)與乳頭肌的運動,仿真模型及過程更符合人體真實二尖瓣情況,仿真結(jié)果也更接近臨床數(shù)據(jù),更適合后續(xù)二尖瓣病理情況的仿真與分析.

    瓣葉實際材料屬性是非線性各向異性的,本文瓣葉簡化為線性各向同性,以后的研究將細(xì)化瓣葉屬性.血液的流動對二尖瓣的行為產(chǎn)生影響,二尖瓣瓣膜在血液影響下的流固耦合運動情況將是后續(xù)值得深入研究的一個方向.

    [1] 馮敏,黃曉陽,蘇茂龍,等.基于雙彈性腔的左心循環(huán)系統(tǒng)的建模與仿真[J].廈門大學(xué)學(xué)報:自然科學(xué)版,2012,51 (4):676-681.

    [2] Dominik J,Zacek P.Heart valve surgery:an illustrated guide[M].Berlin:Spring-Verlag,2010:10-11.

    [3] Lam J H,Ranganathan N,Wigle E D,et al.Morphology of the human mitral valve:Ⅰ.Chordae tendinese:a new classification[J].Circulation,1970,41(3):449-458.

    [4] 黃化清.二尖瓣腿索斷裂[J].臨床心血管病雜志,1985,1 (1/2):80-83.

    [5] Kunzelman K S,Cochran R P,Chuong C,et al.Finite element analysis of the mitral valve[J].J Heart Valve Dis, 1993,2:326-340.

    [6] Kunzelman K S,Cochran R P,Verrier E D,et al.Anatomic basis for mitral valve modeling[J].J Heart Valve Dis, 1994,3:491-496.

    [7] Kunzelman K S,Einstein D R,Cochran R P.Fluid-structure interaction models of the mitral valve:function in normal and pathological states[J].Philosophical Transactions of the Royal Society B:Biological Sciences,2007, 362(1484):1393-1406.

    [8] Lim K H,Xeo J H,Duran C.Three-dimensional asymmetrical modeling of the mitral valve:a finite element study with dynamic boundaries[J].J Heart Valve Dis, 2005,14:386-392.

    [9] 羅云.二尖瓣有限元建模與生物力學(xué)分析[D].廈門:廈門大學(xué),2014.

    [10] Lau K D,Diaz V,Scambler P,et al.Mitral valve dynamics in structural and fluid-structure interaction models [J].Medical Engineering&Physics,2010,32(9): 1057-1064.

    [11] Ranganathan N,Lam J H C,Wigle E D,et al.Morphology of the human mitral valve:II.The valve leaflets[J]. Circulation,1970,41:459-467.

    [12] Pan F D,Donzella G,Fucci C,et al.Structural effects of an innovative surgical technique to repair heart valve defects[J].J Biomech,2005,38:2460-2471.

    [13] Kunzelman K S,Cochran K P.Mechanical properties of basal and marginal mitral valve chordae tendineae[J]. ASAIO Trans,1990,36(3):405-408.

    [14] Kunzelman K S,Cochran R P,Murphree S S,et al.Differential collagen distribution in the mitral valve and its influence on biomechanical behaviour[J].J Heart Valve Dis,1993,2:236-244.

    [15] 鐘琪.人體心臟二尖瓣有限元建模與仿真[D].廈門:廈門大學(xué),2014.

    [16] Stevanella M,Votta E,Redaelli A.Mitral valve finite element modeling:implicatioons of tissues′nonlinear response and annular motion[J].Journal of Bionechanical Engineering,2009,131:1-9.

    [17] Levick J R.An introduction to cardiovascular physiology [M].London:Hodder Arnold Publication,2003.

    [18] 盧曉芳,呂清,謝明星,等.實時三維超聲心動圖評價二尖瓣裂的初步研究[J].中國醫(yī)學(xué)影像技術(shù),2006,22 (4):490-492.

    [19] 鐘琪,曾文華,黃曉陽,等.人體二尖瓣建模及生物力學(xué)分析[J].醫(yī)用生物力學(xué),2014,29(1):53-59.

    [20] Lansac E,Lim K H,Shomura X,et al.Dynamic balance of the aortomitral junction[J].J Thorac Cardiovasc Surg,2002,123(5):911-918.

    [21] Ormiston J A,Shah P M,Tei C,et al.Size and motion of the mitral valve annulus in man.I.A two-dimensional echocardiographic method and findings in normal subjects[J].Circulation,1981,64(1):113-120.

    [22] Fyrenius A,Engvall J,Janerot-Sjoberg B.Major and minor axes of the normal mitral annulus[J].J Heart Valve Dis,2001,10(2):146-152.

    [23] Gorman III J H,Jackson B M,Enomoto X,et al.The effect of regional ischemia on mitral valve annular saddle shape[J].Ann Thorac Surg,2004,77(2):544-548.

    Modeling and Simulation Methods for Mitral Valve Based on the Mitral Valve Sub-region

    ZHAO Xiao-jia1,SONG Xiao2,HUANG Xiao-yang1*,SU Mao-long3,ZHONG Qi1
    (1.School of Information Science and Engineering,Xiamen University,Xiamen 361005,China;2.College of Computer and Information Technology,Nanyang Normal University,Nanyang 473061,China;3 Ultrasound Department, Zhongshan Hospital,Xiamen University,Xiamen 361004,China)

    A fine model of the human mitral valve with non-uniform thickness has been established.The valve closure has been simulated using the finite element method during early systolic period,compared with the uniform one.Simulation results show that the stress is unevenly distributed on the leaflets,that the mitral valve with non-uniform has lower stress level than the uniform one.Furthermore,the maximum principal stresses occurs at the clefts between the three scallops in the posterior leaflet,indicating positions of the mitral valve cleft in clinic. The result is consistent with the physiology.It suggests that the method and the results are correct.Finally,based on clinical data,the motion of annulus and papillary muscles has been simulated,showing higher closure and relatively uniform stress distribution,closer to the actual physiological conditons.It laid the foundation for the simulation of subsequent pathology.

    the mitral valve;finite element model;stress distribution;maximum principal stress

    10.6043/j.issn.0438-0479.2015.02.021

    TP 391

    A

    0438-0479(2015)02-0276-05

    2014-08-28 錄用日期:2014-11-19

    國家自然科學(xué)基金(61102137,61271336,61327001);福建省自然科學(xué)基金(2014J01438)

    *通信作者:xyhuang@xmu.edu.cn

    趙曉佳,宋曉,黃曉陽,等.基于亞區(qū)結(jié)構(gòu)的人體二尖瓣的建模與仿真[J].廈門大學(xué)學(xué)報:自然科學(xué)版,2015,54(2): 276-280.

    :Zhao Xiaojia,Song Xiao,Huang Xiaoyang,et al.Research on modeling and simulation methods for mitral valve based on the mitral valve sub-region[J].Journal of Xiamen University:Natural Science,2015,54(2):276-280.(in Chinese)

    猜你喜歡
    乳頭肌亞區(qū)瓣葉
    開在夾縫里的花
    詩歌月刊(2023年9期)2023-09-28 07:40:58
    起源于左室乳頭肌的室性心律失常的研究進展
    右心室乳頭肌的解剖學(xué)特征及臨床意義
    京津風(fēng)沙源區(qū)生態(tài)保護與建設(shè)工程對防風(fēng)固沙服務(wù)功能的影響
    淺析福建深部高溫巖體地震異常響應(yīng)
    阿爾茨海默病前扣帶回亞區(qū)體積與認(rèn)知損傷相關(guān)性
    左心室乳頭肌起源室性心律失常的電生理特征及導(dǎo)管消融
    Circulation:TAVI術(shù)后晚期瓣葉血栓常見,是否需要抗凝治療?
    基于海馬亞區(qū)的阿爾茨海默病磁共振結(jié)構(gòu)和功能連接研究
    超聲心動圖診斷乳頭肌肥厚型心肌病1例
    欧美人与性动交α欧美精品济南到| 又大又爽又粗| 亚洲视频免费观看视频| 91成人精品电影| 欧美日韩黄片免| 黄色成人免费大全| 欧美日韩乱码在线| 18禁观看日本| 午夜福利高清视频| 亚洲 欧美一区二区三区| 麻豆国产av国片精品| 亚洲九九香蕉| 亚洲精品国产精品久久久不卡| 午夜福利欧美成人| 亚洲精品一区av在线观看| 神马国产精品三级电影在线观看 | 久久久水蜜桃国产精品网| 午夜免费成人在线视频| 99国产精品一区二区三区| 夜夜躁狠狠躁天天躁| 成熟少妇高潮喷水视频| 国产av又大| 男女之事视频高清在线观看| 亚洲熟妇中文字幕五十中出| 亚洲熟妇中文字幕五十中出| 日韩欧美三级三区| 国产蜜桃级精品一区二区三区| 亚洲国产精品999在线| 99久久久亚洲精品蜜臀av| 美国免费a级毛片| 亚洲中文日韩欧美视频| 国产精品 欧美亚洲| 色尼玛亚洲综合影院| 亚洲第一欧美日韩一区二区三区| 欧美丝袜亚洲另类 | 亚洲人成网站在线播放欧美日韩| 欧美精品亚洲一区二区| 亚洲一区二区三区不卡视频| 一区二区三区国产精品乱码| 久久久久久久久中文| 精品卡一卡二卡四卡免费| 免费观看精品视频网站| 天堂影院成人在线观看| 亚洲欧洲精品一区二区精品久久久| 免费在线观看黄色视频的| 欧美成人性av电影在线观看| 啦啦啦 在线观看视频| 日韩欧美一区视频在线观看| 午夜影院日韩av| 黑人巨大精品欧美一区二区蜜桃| 搞女人的毛片| 国产精品乱码一区二三区的特点 | 欧美老熟妇乱子伦牲交| 久久精品国产清高在天天线| 午夜免费观看网址| 美女免费视频网站| 伊人久久大香线蕉亚洲五| 黑人巨大精品欧美一区二区mp4| 女人爽到高潮嗷嗷叫在线视频| 久久久久久久久中文| 不卡av一区二区三区| 免费在线观看视频国产中文字幕亚洲| 少妇熟女aⅴ在线视频| 亚洲色图av天堂| 国产视频一区二区在线看| 午夜两性在线视频| 久久久久久免费高清国产稀缺| 亚洲三区欧美一区| 亚洲欧美日韩另类电影网站| 欧美国产精品va在线观看不卡| 美国免费a级毛片| 色播在线永久视频| 国产亚洲精品综合一区在线观看 | 色综合亚洲欧美另类图片| 一级a爱片免费观看的视频| 丁香欧美五月| 亚洲美女黄片视频| 美女免费视频网站| 99re在线观看精品视频| 日本欧美视频一区| avwww免费| 一区二区日韩欧美中文字幕| 男女做爰动态图高潮gif福利片 | 国产午夜精品久久久久久| 国产又爽黄色视频| 搡老岳熟女国产| 巨乳人妻的诱惑在线观看| 国产xxxxx性猛交| 日韩欧美一区视频在线观看| 国产私拍福利视频在线观看| 一二三四社区在线视频社区8| 美女高潮喷水抽搐中文字幕| av有码第一页| 久久久国产精品麻豆| 香蕉丝袜av| 人成视频在线观看免费观看| 一边摸一边抽搐一进一小说| 一边摸一边做爽爽视频免费| 国产乱人伦免费视频| 在线观看免费视频日本深夜| 老汉色∧v一级毛片| 国产欧美日韩一区二区三区在线| 亚洲av成人不卡在线观看播放网| 此物有八面人人有两片| 久久久国产欧美日韩av| 中文字幕人成人乱码亚洲影| 在线观看免费视频网站a站| 亚洲成人国产一区在线观看| 一区二区日韩欧美中文字幕| 亚洲精品美女久久久久99蜜臀| 久久久水蜜桃国产精品网| 91av网站免费观看| 夜夜爽天天搞| 久久性视频一级片| 美女 人体艺术 gogo| 亚洲成国产人片在线观看| 精品久久久久久成人av| 国产av一区在线观看免费| 欧美最黄视频在线播放免费| 久久 成人 亚洲| 亚洲国产精品合色在线| 777久久人妻少妇嫩草av网站| 亚洲三区欧美一区| 一区二区三区高清视频在线| 97碰自拍视频| 亚洲第一av免费看| 久久精品国产清高在天天线| 国产成年人精品一区二区| 大陆偷拍与自拍| 亚洲一码二码三码区别大吗| 嫁个100分男人电影在线观看| 青草久久国产| 国产精品秋霞免费鲁丝片| 亚洲欧美精品综合久久99| 啦啦啦观看免费观看视频高清 | 国产成年人精品一区二区| 俄罗斯特黄特色一大片| 1024视频免费在线观看| 中文字幕人妻熟女乱码| 狂野欧美激情性xxxx| 9色porny在线观看| 国产99久久九九免费精品| 手机成人av网站| 身体一侧抽搐| 亚洲国产欧美日韩在线播放| 日韩欧美国产一区二区入口| 日韩欧美一区视频在线观看| 19禁男女啪啪无遮挡网站| 狠狠狠狠99中文字幕| 视频在线观看一区二区三区| 美女 人体艺术 gogo| 欧美日韩亚洲国产一区二区在线观看| 我的亚洲天堂| 咕卡用的链子| 亚洲专区国产一区二区| 色哟哟哟哟哟哟| 欧美亚洲日本最大视频资源| 亚洲精品在线美女| 亚洲精品在线观看二区| 在线观看免费视频网站a站| 国产成年人精品一区二区| 91av网站免费观看| 久久中文字幕一级| 成人国产综合亚洲| 精品少妇一区二区三区视频日本电影| 日韩国内少妇激情av| 亚洲 欧美 日韩 在线 免费| 日日摸夜夜添夜夜添小说| 精品久久久久久久久久免费视频| 黑人操中国人逼视频| 精品福利观看| 午夜福利高清视频| 久久精品aⅴ一区二区三区四区| 亚洲精品美女久久久久99蜜臀| 亚洲性夜色夜夜综合| 18禁裸乳无遮挡免费网站照片 | 欧美黑人精品巨大| 久久精品国产亚洲av香蕉五月| 精品人妻在线不人妻| 亚洲男人天堂网一区| 久久人妻熟女aⅴ| 又大又爽又粗| 久久这里只有精品19| 一卡2卡三卡四卡精品乱码亚洲| 高清在线国产一区| 在线国产一区二区在线| 亚洲av片天天在线观看| 欧美精品亚洲一区二区| 此物有八面人人有两片| 99re在线观看精品视频| 午夜免费观看网址| 午夜福利免费观看在线| 香蕉丝袜av| 国产欧美日韩一区二区精品| 中文字幕色久视频| av天堂久久9| 咕卡用的链子| 麻豆成人av在线观看| 亚洲国产高清在线一区二区三 | svipshipincom国产片| 亚洲avbb在线观看| 啦啦啦韩国在线观看视频| 国产免费男女视频| 中文亚洲av片在线观看爽| 日韩免费av在线播放| 亚洲精品久久成人aⅴ小说| 日韩欧美三级三区| 日韩视频一区二区在线观看| 一本久久中文字幕| 女人精品久久久久毛片| 最新美女视频免费是黄的| 亚洲专区字幕在线| 欧美成人午夜精品| 久久人人精品亚洲av| 国产蜜桃级精品一区二区三区| 亚洲狠狠婷婷综合久久图片| 十八禁人妻一区二区| 99久久国产精品久久久| 色av中文字幕| 一级毛片精品| av免费在线观看网站| 精品人妻1区二区| 国产欧美日韩精品亚洲av| www.www免费av| 日本撒尿小便嘘嘘汇集6| 黑人欧美特级aaaaaa片| 亚洲色图av天堂| 人人妻,人人澡人人爽秒播| 国产主播在线观看一区二区| 欧美日韩中文字幕国产精品一区二区三区 | 国产成人欧美| 亚洲成人精品中文字幕电影| av视频免费观看在线观看| 叶爱在线成人免费视频播放| 欧美日韩福利视频一区二区| 午夜老司机福利片| 成人亚洲精品av一区二区| 人人澡人人妻人| 欧美绝顶高潮抽搐喷水| 亚洲国产欧美日韩在线播放| 国产又色又爽无遮挡免费看| or卡值多少钱| 后天国语完整版免费观看| www.精华液| 欧美国产精品va在线观看不卡| 好男人电影高清在线观看| 变态另类丝袜制服| 夜夜爽天天搞| 99久久综合精品五月天人人| 亚洲欧美日韩高清在线视频| 老熟妇乱子伦视频在线观看| 精品久久久久久,| 亚洲狠狠婷婷综合久久图片| 在线观看免费日韩欧美大片| 国产单亲对白刺激| 精品久久久久久久人妻蜜臀av | 久久婷婷成人综合色麻豆| av欧美777| 久久人妻福利社区极品人妻图片| 国产av又大| 亚洲男人的天堂狠狠| 色婷婷久久久亚洲欧美| 村上凉子中文字幕在线| 成人亚洲精品一区在线观看| 香蕉国产在线看| 国产aⅴ精品一区二区三区波| 精品少妇一区二区三区视频日本电影| 国产精品98久久久久久宅男小说| 亚洲精品在线美女| 久久久久久免费高清国产稀缺| 亚洲精品久久成人aⅴ小说| 亚洲中文av在线| 99riav亚洲国产免费| 欧美中文综合在线视频| 岛国在线观看网站| 国产在线观看jvid| 亚洲九九香蕉| av片东京热男人的天堂| 国产成人免费无遮挡视频| 欧美亚洲日本最大视频资源| 久久天堂一区二区三区四区| 亚洲五月天丁香| 亚洲人成伊人成综合网2020| 亚洲色图综合在线观看| 又紧又爽又黄一区二区| 一区二区日韩欧美中文字幕| 久久精品国产亚洲av高清一级| x7x7x7水蜜桃| 黑人欧美特级aaaaaa片| 国产精品99久久99久久久不卡| 国产又爽黄色视频| 国产精品爽爽va在线观看网站 | 久久精品aⅴ一区二区三区四区| 免费观看人在逋| 国产精品秋霞免费鲁丝片| 啦啦啦观看免费观看视频高清 | 亚洲国产精品成人综合色| 精品一区二区三区四区五区乱码| 久久人人爽av亚洲精品天堂| 国产精品免费视频内射| 久久精品国产清高在天天线| 高清毛片免费观看视频网站| 日本三级黄在线观看| 欧美不卡视频在线免费观看 | 日本vs欧美在线观看视频| 看片在线看免费视频| 老熟妇仑乱视频hdxx| 免费在线观看完整版高清| 亚洲黑人精品在线| 国产成人欧美在线观看| 手机成人av网站| 午夜a级毛片| 99国产极品粉嫩在线观看| 午夜福利一区二区在线看| 国产成人精品无人区| 好男人在线观看高清免费视频 | 老司机午夜福利在线观看视频| 国产亚洲精品一区二区www| 午夜福利高清视频| 国产亚洲欧美在线一区二区| 国产区一区二久久| 搡老熟女国产l中国老女人| 亚洲伊人色综图| 亚洲精品国产精品久久久不卡| 操出白浆在线播放| 在线国产一区二区在线| 乱人伦中国视频| 亚洲狠狠婷婷综合久久图片| 非洲黑人性xxxx精品又粗又长| 国产伦一二天堂av在线观看| 亚洲欧美日韩无卡精品| 男女之事视频高清在线观看| 中文字幕色久视频| 日韩精品中文字幕看吧| 国产精品久久久久久精品电影 | 叶爱在线成人免费视频播放| 国产亚洲欧美98| 国产熟女xx| 亚洲精品一区av在线观看| 90打野战视频偷拍视频| 亚洲人成伊人成综合网2020| 午夜福利影视在线免费观看| 国产亚洲精品第一综合不卡| 极品教师在线免费播放| 两人在一起打扑克的视频| 国产伦一二天堂av在线观看| 亚洲欧美日韩无卡精品| 电影成人av| 两个人看的免费小视频| 国产一区二区三区视频了| 亚洲中文av在线| 久久精品国产99精品国产亚洲性色 | 国产av一区在线观看免费| 在线观看66精品国产| 亚洲久久久国产精品| 麻豆国产av国片精品| 黄色毛片三级朝国网站| 波多野结衣高清无吗| 国产三级黄色录像| 亚洲人成电影免费在线| 99精品在免费线老司机午夜| 91字幕亚洲| 欧美日韩中文字幕国产精品一区二区三区 | 久久久久久大精品| 国产成人欧美在线观看| 久久精品亚洲精品国产色婷小说| 日本 av在线| 成人手机av| 韩国av一区二区三区四区| 久久国产精品男人的天堂亚洲| 免费高清在线观看日韩| 国产97色在线日韩免费| 中文字幕精品免费在线观看视频| 欧美在线黄色| 国产精华一区二区三区| 看免费av毛片| 老熟妇仑乱视频hdxx| а√天堂www在线а√下载| 国产高清videossex| 久久久久久免费高清国产稀缺| 亚洲av日韩精品久久久久久密| 国产麻豆69| 一二三四在线观看免费中文在| 亚洲久久久国产精品| 性少妇av在线| 老熟妇乱子伦视频在线观看| 99精品欧美一区二区三区四区| 亚洲人成电影免费在线| 超碰成人久久| 久久精品国产清高在天天线| 麻豆国产av国片精品| 搞女人的毛片| 亚洲av成人av| 校园春色视频在线观看| 黄色成人免费大全| 亚洲成av片中文字幕在线观看| bbb黄色大片| 成人国产综合亚洲| 男女床上黄色一级片免费看| 两人在一起打扑克的视频| 精品少妇一区二区三区视频日本电影| 久久久久久久午夜电影| 久久久久国产一级毛片高清牌| 国产一级毛片七仙女欲春2 | 99国产精品免费福利视频| 十分钟在线观看高清视频www| 一二三四在线观看免费中文在| videosex国产| 婷婷丁香在线五月| 国产亚洲精品久久久久久毛片| 欧美日韩亚洲综合一区二区三区_| 欧美黄色淫秽网站| 亚洲欧美精品综合一区二区三区| 99在线人妻在线中文字幕| 老司机午夜福利在线观看视频| 国产午夜福利久久久久久| 成在线人永久免费视频| 丰满人妻熟妇乱又伦精品不卡| 久久久久久人人人人人| 日本免费a在线| 久久久久精品国产欧美久久久| 日本免费一区二区三区高清不卡 | 中文字幕精品免费在线观看视频| 国产精品免费一区二区三区在线| av福利片在线| 亚洲一区二区三区不卡视频| 欧美激情 高清一区二区三区| 亚洲免费av在线视频| 在线观看66精品国产| 免费高清在线观看日韩| 老司机福利观看| 一区二区三区高清视频在线| 两个人视频免费观看高清| 黄色a级毛片大全视频| 亚洲国产欧美日韩在线播放| 日本免费a在线| 脱女人内裤的视频| 欧美av亚洲av综合av国产av| 国产精品99久久99久久久不卡| 久久久久亚洲av毛片大全| 亚洲国产精品成人综合色| 亚洲色图综合在线观看| 在线观看免费视频网站a站| 亚洲一卡2卡3卡4卡5卡精品中文| 又黄又爽又免费观看的视频| 欧美激情极品国产一区二区三区| 日本三级黄在线观看| 午夜免费成人在线视频| 亚洲五月天丁香| 成人精品一区二区免费| 午夜福利一区二区在线看| e午夜精品久久久久久久| 老汉色∧v一级毛片| a级毛片在线看网站| 咕卡用的链子| 亚洲熟妇中文字幕五十中出| 色av中文字幕| 亚洲av成人不卡在线观看播放网| 给我免费播放毛片高清在线观看| 欧美乱码精品一区二区三区| 亚洲,欧美精品.| 亚洲av电影不卡..在线观看| 日韩欧美国产一区二区入口| 一区二区三区精品91| 亚洲最大成人中文| 亚洲无线在线观看| 日韩欧美在线二视频| 国产男靠女视频免费网站| 亚洲色图 男人天堂 中文字幕| 岛国在线观看网站| 露出奶头的视频| 国产片内射在线| 中文字幕人成人乱码亚洲影| 亚洲欧美激情综合另类| 国产午夜福利久久久久久| av网站免费在线观看视频| netflix在线观看网站| 欧美老熟妇乱子伦牲交| 久久久国产成人精品二区| 91字幕亚洲| 国产av又大| 亚洲欧美日韩另类电影网站| 满18在线观看网站| 黑人操中国人逼视频| 亚洲熟妇熟女久久| 国产激情久久老熟女| 日韩成人在线观看一区二区三区| 欧美日韩福利视频一区二区| 久久香蕉精品热| 久久精品国产亚洲av高清一级| 亚洲国产欧美网| 女生性感内裤真人,穿戴方法视频| 亚洲天堂国产精品一区在线| 中文字幕人成人乱码亚洲影| 欧美色欧美亚洲另类二区 | 中亚洲国语对白在线视频| 超碰成人久久| 叶爱在线成人免费视频播放| 黄色毛片三级朝国网站| 亚洲,欧美精品.| 美女扒开内裤让男人捅视频| 亚洲少妇的诱惑av| 亚洲人成电影观看| 又黄又粗又硬又大视频| 国产精品日韩av在线免费观看 | 午夜a级毛片| 国产亚洲欧美98| 亚洲美女黄片视频| 97超级碰碰碰精品色视频在线观看| 亚洲国产欧美日韩在线播放| 亚洲精品久久国产高清桃花| 女性被躁到高潮视频| 精品国产乱码久久久久久男人| 麻豆一二三区av精品| 欧美黑人欧美精品刺激| av天堂久久9| 精品人妻在线不人妻| 色av中文字幕| 亚洲国产毛片av蜜桃av| 无人区码免费观看不卡| 久久婷婷成人综合色麻豆| 日韩高清综合在线| 无限看片的www在线观看| 精品久久久久久久毛片微露脸| 99国产精品99久久久久| 亚洲av成人av| xxx96com| 国产私拍福利视频在线观看| 午夜久久久在线观看| 欧美在线黄色| 久久久国产成人精品二区| 国产av在哪里看| 精品国产一区二区久久| 国产成人欧美| 两性夫妻黄色片| 亚洲成人免费电影在线观看| 一个人观看的视频www高清免费观看 | 久久久久久久精品吃奶| 人人澡人人妻人| 在线免费观看的www视频| 伦理电影免费视频| 午夜福利一区二区在线看| 熟女少妇亚洲综合色aaa.| 人人妻人人澡欧美一区二区 | 亚洲精品一卡2卡三卡4卡5卡| 久久中文字幕一级| 精品久久久精品久久久| 国产精品自产拍在线观看55亚洲| 国产成人免费无遮挡视频| 成人免费观看视频高清| 精品久久久久久久人妻蜜臀av | 午夜精品在线福利| 亚洲 欧美 日韩 在线 免费| 少妇裸体淫交视频免费看高清 | 久久久久久久久中文| 老汉色av国产亚洲站长工具| 成人欧美大片| 久久久久久久午夜电影| 热re99久久国产66热| 女性生殖器流出的白浆| 成人精品一区二区免费| 嫩草影院精品99| 午夜日韩欧美国产| 国产一区二区三区在线臀色熟女| 亚洲av片天天在线观看| 免费高清在线观看日韩| 咕卡用的链子| 久久人妻熟女aⅴ| 免费在线观看视频国产中文字幕亚洲| 亚洲无线在线观看| 亚洲一区二区三区不卡视频| 叶爱在线成人免费视频播放| 熟妇人妻久久中文字幕3abv| 日本免费a在线| 色老头精品视频在线观看| 高清黄色对白视频在线免费看| 日韩国内少妇激情av| 国产成人精品无人区| 美女国产高潮福利片在线看| 国产精品久久久久久人妻精品电影| 欧美不卡视频在线免费观看 | av免费在线观看网站| 啪啪无遮挡十八禁网站| 九色亚洲精品在线播放| 视频在线观看一区二区三区| 一边摸一边做爽爽视频免费| 午夜免费激情av| 午夜两性在线视频| 好男人在线观看高清免费视频 | 国产乱人伦免费视频| 日本 欧美在线| 亚洲三区欧美一区| 国产av又大| 亚洲全国av大片| 两人在一起打扑克的视频| 久久国产乱子伦精品免费另类| 国产97色在线日韩免费| 搡老岳熟女国产| 亚洲男人的天堂狠狠| 1024香蕉在线观看| 夜夜躁狠狠躁天天躁| 国产一区二区在线av高清观看| 窝窝影院91人妻| 国产成人精品在线电影| 亚洲熟妇熟女久久| 国产精品野战在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 91大片在线观看| 国产一区二区三区综合在线观看| 一级毛片女人18水好多| 日韩三级视频一区二区三区| 久久 成人 亚洲| 成人亚洲精品一区在线观看| 成人永久免费在线观看视频| 制服丝袜大香蕉在线|