• <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例
    亚洲欧美一区二区三区国产| 日韩 亚洲 欧美在线| 中文字幕av电影在线播放| av免费观看日本| 亚洲av欧美aⅴ国产| 插逼视频在线观看| 伦理电影大哥的女人| 中国美白少妇内射xxxbb| 国产欧美日韩一区二区三区在线 | 久久久久久久亚洲中文字幕| 丝袜喷水一区| 黑人巨大精品欧美一区二区蜜桃 | 中文字幕精品免费在线观看视频 | 亚洲国产欧美在线一区| 久久久久久伊人网av| 汤姆久久久久久久影院中文字幕| 天天操日日干夜夜撸| 国产免费现黄频在线看| 亚洲经典国产精华液单| 国产精品国产av在线观看| 精品人妻在线不人妻| 久久久久久人妻| 亚洲在久久综合| 久久精品熟女亚洲av麻豆精品| 国产在线一区二区三区精| 婷婷色av中文字幕| 97精品久久久久久久久久精品| 男人爽女人下面视频在线观看| 少妇猛男粗大的猛烈进出视频| 久久久久久久精品精品| 午夜久久久在线观看| 国产一区亚洲一区在线观看| 少妇丰满av| 美女国产视频在线观看| 九九久久精品国产亚洲av麻豆| 超色免费av| 一区在线观看完整版| 亚洲一级一片aⅴ在线观看| 黑人猛操日本美女一级片| 久久精品久久精品一区二区三区| 久久青草综合色| 最近手机中文字幕大全| 国产精品久久久久成人av| 国产黄色视频一区二区在线观看| 国产精品久久久久久久久免| 91精品国产九色| 亚洲欧美色中文字幕在线| 国产淫语在线视频| 伊人久久精品亚洲午夜| 国产免费福利视频在线观看| 国产精品成人在线| 啦啦啦在线观看免费高清www| 一级毛片电影观看| 国产成人a∨麻豆精品| 插阴视频在线观看视频| 少妇 在线观看| 国产69精品久久久久777片| 伊人久久国产一区二区| 99久久精品一区二区三区| 女人久久www免费人成看片| 国产精品 国内视频| 人人妻人人澡人人看| 亚洲国产欧美在线一区| 青春草亚洲视频在线观看| 人妻人人澡人人爽人人| 欧美3d第一页| 国产亚洲av片在线观看秒播厂| 人妻 亚洲 视频| 欧美人与善性xxx| 中文字幕人妻丝袜制服| 午夜日本视频在线| av免费观看日本| 日本猛色少妇xxxxx猛交久久| 久久99热这里只频精品6学生| 黄色毛片三级朝国网站| 国产亚洲午夜精品一区二区久久| xxxhd国产人妻xxx| 女性被躁到高潮视频| 男男h啪啪无遮挡| 亚洲四区av| 视频中文字幕在线观看| 王馨瑶露胸无遮挡在线观看| 欧美精品国产亚洲| 最近最新中文字幕免费大全7| 国产精品久久久久久久电影| 国产精品不卡视频一区二区| 亚洲精品国产av蜜桃| 国产成人aa在线观看| 国产精品久久久久久精品古装| 乱码一卡2卡4卡精品| 国产毛片在线视频| 夜夜看夜夜爽夜夜摸| 欧美日韩在线观看h| 一本色道久久久久久精品综合| 午夜免费鲁丝| 精品一区二区三区视频在线| 亚洲不卡免费看| 午夜视频国产福利| 大陆偷拍与自拍| 欧美日韩视频精品一区| 色吧在线观看| 国产成人av激情在线播放 | 黄色欧美视频在线观看| 伊人久久精品亚洲午夜| 午夜日本视频在线| 亚洲综合色网址| 91久久精品国产一区二区成人| 少妇 在线观看| 亚洲少妇的诱惑av| 亚洲欧美成人精品一区二区| 欧美xxⅹ黑人| 国产伦精品一区二区三区视频9| 人人妻人人添人人爽欧美一区卜| 国产精品国产三级国产专区5o| 精品一区二区三卡| 国产熟女欧美一区二区| 成人国产麻豆网| 中文字幕最新亚洲高清| 亚洲,一卡二卡三卡| 我要看黄色一级片免费的| 午夜精品国产一区二区电影| 99久久精品一区二区三区| 全区人妻精品视频| 日韩伦理黄色片| 成人漫画全彩无遮挡| 亚洲av免费高清在线观看| 18禁动态无遮挡网站| 少妇被粗大的猛进出69影院 | 免费看av在线观看网站| 成人影院久久| 久久影院123| 妹子高潮喷水视频| 日本-黄色视频高清免费观看| 色哟哟·www| 国产片特级美女逼逼视频| 午夜免费男女啪啪视频观看| 日韩不卡一区二区三区视频在线| 国产精品熟女久久久久浪| h视频一区二区三区| 大香蕉久久成人网| 能在线免费看毛片的网站| 美女中出高潮动态图| 婷婷色麻豆天堂久久| 在线观看美女被高潮喷水网站| 久久久久久人妻| 夜夜看夜夜爽夜夜摸| 97在线视频观看| 蜜桃国产av成人99| 亚洲欧美精品自产自拍| 久久精品久久精品一区二区三区| 亚洲精品成人av观看孕妇| 午夜激情福利司机影院| 久久99热6这里只有精品| 女的被弄到高潮叫床怎么办| 丁香六月天网| kizo精华| 亚洲欧美精品自产自拍| www.色视频.com| 人人澡人人妻人| 久久国内精品自在自线图片| 久久精品久久久久久久性| 免费看光身美女| 搡女人真爽免费视频火全软件| 日韩亚洲欧美综合| 看免费成人av毛片| 十分钟在线观看高清视频www| 最近2019中文字幕mv第一页| 亚洲av二区三区四区| 免费不卡的大黄色大毛片视频在线观看| 下体分泌物呈黄色| 成人亚洲精品一区在线观看| 中文字幕最新亚洲高清| 丝瓜视频免费看黄片| 中文字幕人妻丝袜制服| 免费观看性生交大片5| 又黄又爽又刺激的免费视频.| 亚洲激情五月婷婷啪啪| 亚洲成人手机| 国产精品人妻久久久久久| 精品国产乱码久久久久久小说| 日本欧美国产在线视频| 国产探花极品一区二区| 久久久国产一区二区| 蜜桃国产av成人99| 美女大奶头黄色视频| av在线播放精品| 亚洲av成人精品一区久久| 午夜日本视频在线| 久久国内精品自在自线图片| 国产探花极品一区二区| xxx大片免费视频| 麻豆乱淫一区二区| 精品熟女少妇av免费看| 啦啦啦视频在线资源免费观看| 久久久a久久爽久久v久久| 男女啪啪激烈高潮av片| 国产成人a∨麻豆精品| av专区在线播放| 午夜av观看不卡| 美女cb高潮喷水在线观看| 国产又色又爽无遮挡免| 18禁裸乳无遮挡动漫免费视频| 韩国av在线不卡| 色94色欧美一区二区| 九色成人免费人妻av| 丝袜在线中文字幕| 纵有疾风起免费观看全集完整版| 精品卡一卡二卡四卡免费| 两个人免费观看高清视频| 日韩精品有码人妻一区| a 毛片基地| 岛国毛片在线播放| 一级二级三级毛片免费看| 国产亚洲精品第一综合不卡 | 国产男女内射视频| 中文字幕制服av| 亚洲av日韩在线播放| 久久精品久久精品一区二区三区| 少妇的逼水好多| 日韩制服骚丝袜av| 亚洲av综合色区一区| 在线看a的网站| 一级毛片电影观看| 天堂俺去俺来也www色官网| 精品人妻一区二区三区麻豆| 一级二级三级毛片免费看| 纯流量卡能插随身wifi吗| 亚洲欧美精品自产自拍| 男女国产视频网站| 特大巨黑吊av在线直播| videosex国产| 亚洲精品乱久久久久久| 国产乱来视频区| 少妇高潮的动态图| 国产成人91sexporn| 人妻人人澡人人爽人人| 亚洲精品久久成人aⅴ小说 | 精品国产乱码久久久久久小说| 国产亚洲最大av| 久久久久久久久久久丰满| 免费高清在线观看视频在线观看| 久久久久网色| 久久久久久久久久人人人人人人| 久久97久久精品| 人妻 亚洲 视频| 美女国产视频在线观看| 国产爽快片一区二区三区| 国产永久视频网站| 欧美 亚洲 国产 日韩一| 国产欧美亚洲国产| 人人澡人人妻人| 国产色婷婷99| 一级毛片 在线播放| 欧美亚洲日本最大视频资源| 国产极品天堂在线| 男人爽女人下面视频在线观看| 亚洲av二区三区四区| 18禁在线播放成人免费| 赤兔流量卡办理| 亚洲一级一片aⅴ在线观看| 香蕉精品网在线| 欧美精品一区二区大全| 亚洲国产最新在线播放| 看十八女毛片水多多多| 青春草亚洲视频在线观看| 久久国产精品男人的天堂亚洲 | 高清不卡的av网站| 久久av网站| 国产精品免费大片| 国产精品麻豆人妻色哟哟久久| 国产熟女午夜一区二区三区 | 人人妻人人澡人人看| 久久99热6这里只有精品| 午夜福利网站1000一区二区三区| av不卡在线播放| 国产精品一区二区三区四区免费观看| 国产成人av激情在线播放 | 最后的刺客免费高清国语| 国产成人一区二区在线| 有码 亚洲区| 嘟嘟电影网在线观看| 三级国产精品欧美在线观看| 国内精品宾馆在线| 免费观看在线日韩| 最近手机中文字幕大全| 精品国产露脸久久av麻豆| 亚洲欧美日韩卡通动漫| 欧美人与善性xxx| 久久久精品94久久精品| 欧美97在线视频| 极品人妻少妇av视频| 丝袜在线中文字幕| 亚洲天堂av无毛| 一边摸一边做爽爽视频免费| 亚洲精品日韩在线中文字幕| 久久久久国产精品人妻一区二区| 天天操日日干夜夜撸| 最近手机中文字幕大全| 夜夜看夜夜爽夜夜摸| av卡一久久| 不卡视频在线观看欧美| 在线观看美女被高潮喷水网站| av视频免费观看在线观看| 在线观看三级黄色| 国产色婷婷99| 国产精品国产三级国产av玫瑰| 亚洲精品乱码久久久v下载方式| 国产一区二区三区av在线| 一级黄片播放器| 亚洲av欧美aⅴ国产| 99久久精品国产国产毛片| 女性生殖器流出的白浆| 91在线精品国自产拍蜜月| 尾随美女入室| 天堂8中文在线网| 久久久精品免费免费高清| 国产精品人妻久久久影院| 亚洲av在线观看美女高潮| 一区在线观看完整版| 天堂俺去俺来也www色官网| a级毛片在线看网站| 久久久久网色| 男的添女的下面高潮视频| 国产精品熟女久久久久浪| 九九在线视频观看精品| 午夜福利影视在线免费观看| 肉色欧美久久久久久久蜜桃| 中文字幕精品免费在线观看视频 | 亚洲综合色惰| 亚洲国产精品国产精品| 国产伦精品一区二区三区视频9| 18禁动态无遮挡网站| 制服诱惑二区| 999精品在线视频| 日本与韩国留学比较| 狂野欧美激情性bbbbbb| 久久免费观看电影| 最近的中文字幕免费完整| 久久99一区二区三区| 国语对白做爰xxxⅹ性视频网站| 成人免费观看视频高清| 嫩草影院入口| 人体艺术视频欧美日本| 国产精品久久久久成人av| 美女脱内裤让男人舔精品视频| 亚洲婷婷狠狠爱综合网| 精品少妇黑人巨大在线播放| 国产av国产精品国产| 国产成人freesex在线| 亚洲国产精品一区三区| 91精品一卡2卡3卡4卡| 丰满乱子伦码专区| 中国美白少妇内射xxxbb| 久久久久人妻精品一区果冻| 国产成人一区二区在线| av国产精品久久久久影院| 高清视频免费观看一区二区| 在线观看美女被高潮喷水网站| 午夜福利,免费看| 亚洲情色 制服丝袜| 熟女av电影| 久久影院123| 在线播放无遮挡| 国产视频首页在线观看| 熟女av电影| 成人18禁高潮啪啪吃奶动态图 | 精品久久久久久电影网| 久久精品久久久久久久性| 久久狼人影院| 日韩免费高清中文字幕av| 麻豆乱淫一区二区| 在线精品无人区一区二区三| 日韩av免费高清视频| 国产爽快片一区二区三区| 免费观看在线日韩| 国产成人免费无遮挡视频| 少妇的逼水好多| 国产一区二区在线观看日韩| 曰老女人黄片| 一个人免费看片子| 一级,二级,三级黄色视频| 亚洲国产av新网站| 日韩 亚洲 欧美在线| 免费观看在线日韩| 一边亲一边摸免费视频| 国产一级毛片在线| 亚洲欧美精品自产自拍| 国产成人免费观看mmmm| 亚洲四区av| 狂野欧美白嫩少妇大欣赏| a级片在线免费高清观看视频| 亚洲怡红院男人天堂| 麻豆乱淫一区二区| 边亲边吃奶的免费视频| 婷婷色综合www| kizo精华| 国产精品三级大全| 中文字幕制服av| 亚洲av二区三区四区| 国产精品女同一区二区软件| 青春草视频在线免费观看| 菩萨蛮人人尽说江南好唐韦庄| 久久精品国产亚洲网站| 国产淫语在线视频| 日韩av不卡免费在线播放| 王馨瑶露胸无遮挡在线观看| 久久av网站| 中文精品一卡2卡3卡4更新| 亚洲精品乱码久久久v下载方式| 亚洲精品自拍成人| 18禁在线播放成人免费| 飞空精品影院首页| 免费观看在线日韩| 夜夜骑夜夜射夜夜干| 久久女婷五月综合色啪小说| 中文欧美无线码| 99热全是精品| 亚洲av日韩在线播放| 久久99热这里只频精品6学生| videos熟女内射| 亚洲欧美一区二区三区黑人 | 日本黄色日本黄色录像| 多毛熟女@视频| 国产色爽女视频免费观看| av福利片在线| 麻豆精品久久久久久蜜桃| 亚洲精品亚洲一区二区| 国产精品久久久久久精品古装| 日韩伦理黄色片| 亚洲少妇的诱惑av| 大码成人一级视频| av专区在线播放| 丰满迷人的少妇在线观看| 91精品三级在线观看| 黄色怎么调成土黄色| 久久精品人人爽人人爽视色| 婷婷色综合www| 亚洲av欧美aⅴ国产| 在线免费观看不下载黄p国产| 街头女战士在线观看网站| 中文字幕免费在线视频6| 国产成人免费观看mmmm| 人妻 亚洲 视频| 青青草视频在线视频观看| 国产在线视频一区二区| av国产久精品久网站免费入址| 看十八女毛片水多多多| 免费观看的影片在线观看| 国产成人精品无人区| 狂野欧美白嫩少妇大欣赏| 亚洲色图 男人天堂 中文字幕 | 亚洲图色成人| 国产国拍精品亚洲av在线观看| 欧美变态另类bdsm刘玥| 亚洲欧美日韩卡通动漫| 久久人妻熟女aⅴ| 日韩亚洲欧美综合| 国产黄频视频在线观看| 亚洲av在线观看美女高潮| 免费不卡的大黄色大毛片视频在线观看| 精品一品国产午夜福利视频| 美女xxoo啪啪120秒动态图| 性高湖久久久久久久久免费观看| 在现免费观看毛片| 人妻夜夜爽99麻豆av| 97精品久久久久久久久久精品| 国产黄色免费在线视频| 香蕉精品网在线| 亚洲国产av新网站| av专区在线播放| 午夜免费鲁丝| 国产高清不卡午夜福利| videosex国产| 午夜福利影视在线免费观看| 久久免费观看电影| 国产欧美日韩一区二区三区在线 | 久久毛片免费看一区二区三区| 精品少妇久久久久久888优播| 欧美日韩综合久久久久久| 亚洲国产精品成人久久小说| 国产一区二区三区av在线| 伊人久久国产一区二区| 黑人猛操日本美女一级片| 日本黄色日本黄色录像| 国产亚洲精品久久久com| 最新的欧美精品一区二区| 狠狠精品人妻久久久久久综合| 亚洲欧美成人综合另类久久久| 18禁动态无遮挡网站| 最黄视频免费看| 国产精品三级大全| 如何舔出高潮| 超色免费av| 亚洲欧美日韩另类电影网站| 成人毛片a级毛片在线播放| 国产在视频线精品| 亚洲成色77777| 十分钟在线观看高清视频www| 国产日韩欧美在线精品| 亚洲欧洲日产国产| 少妇被粗大的猛进出69影院 | 成年美女黄网站色视频大全免费 | 高清午夜精品一区二区三区| 久久国产亚洲av麻豆专区| 久久鲁丝午夜福利片| 亚洲欧洲精品一区二区精品久久久 | 国产成人精品婷婷| 乱人伦中国视频| 午夜av观看不卡| 特大巨黑吊av在线直播| 看非洲黑人一级黄片| 老熟女久久久| 久久久久久伊人网av| 国产精品麻豆人妻色哟哟久久| 一级片'在线观看视频| 久久狼人影院| 2018国产大陆天天弄谢| 中文字幕久久专区| 国产精品欧美亚洲77777| 97超碰精品成人国产| 欧美国产精品一级二级三级| 精品人妻在线不人妻| 国产视频首页在线观看| 久久久欧美国产精品| 国产精品免费大片| 一级爰片在线观看| 国产国拍精品亚洲av在线观看| 日韩欧美一区视频在线观看| xxx大片免费视频| 欧美97在线视频| 日本vs欧美在线观看视频| 精品人妻一区二区三区麻豆| 黑人高潮一二区| 青春草视频在线免费观看| 精品久久国产蜜桃| 97在线人人人人妻| 亚洲国产精品专区欧美| 久久久a久久爽久久v久久| 视频区图区小说| 亚洲人成网站在线播| 日韩中字成人| 欧美日韩亚洲高清精品| 亚洲精品日韩av片在线观看| 中文字幕免费在线视频6| 美女国产视频在线观看| 黑人巨大精品欧美一区二区蜜桃 | 欧美人与性动交α欧美精品济南到 | 18禁观看日本| 美女国产视频在线观看| tube8黄色片| 精品久久久久久久久av| 国产淫语在线视频| 大码成人一级视频| 亚洲丝袜综合中文字幕| 久久久精品94久久精品| 91精品一卡2卡3卡4卡| xxxhd国产人妻xxx| 免费高清在线观看日韩| av卡一久久| 免费观看av网站的网址| 九色亚洲精品在线播放| 久久久午夜欧美精品| 另类精品久久| 日韩伦理黄色片| 在线亚洲精品国产二区图片欧美 | 久久久久久久大尺度免费视频| 九九久久精品国产亚洲av麻豆| 久久久久久久亚洲中文字幕| 一级毛片电影观看| 看免费成人av毛片| av电影中文网址| 天堂中文最新版在线下载| 日韩 亚洲 欧美在线| av播播在线观看一区| 成年av动漫网址| 亚洲一级一片aⅴ在线观看| 亚洲精品久久成人aⅴ小说 | 下体分泌物呈黄色| 免费人妻精品一区二区三区视频| 中文字幕制服av| 日韩大片免费观看网站| 秋霞伦理黄片| 午夜免费鲁丝| 国产精品一区二区在线观看99| 女性被躁到高潮视频| 中文欧美无线码| 自线自在国产av| 亚洲综合色惰| 成年av动漫网址| 十分钟在线观看高清视频www| 日韩视频在线欧美| 伦理电影免费视频| 看免费成人av毛片| 欧美 亚洲 国产 日韩一| 日韩一本色道免费dvd| 五月天丁香电影| 欧美精品一区二区大全| 亚洲国产欧美日韩在线播放| 亚洲怡红院男人天堂| 男女无遮挡免费网站观看| 青青草视频在线视频观看| 亚洲人成77777在线视频| 亚洲四区av| 亚洲精品乱码久久久久久按摩| 欧美成人精品欧美一级黄| 国产精品秋霞免费鲁丝片| 国产免费又黄又爽又色| 在线观看国产h片| 女性生殖器流出的白浆| 免费观看的影片在线观看| 日韩人妻高清精品专区| 99热这里只有是精品在线观看| 久久精品国产自在天天线| 精品亚洲乱码少妇综合久久| 亚洲综合色惰|