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

    頸動脈分叉的非穩(wěn)態(tài)數(shù)值模擬分析

    2018-11-16 09:37:32王汝良胡霖霖郭金興張凱旋張臘梅戴志穎陳廣新
    軟件 2018年10期
    關(guān)鍵詞:剪切應(yīng)力心動渦流

    王汝良,胡霖霖,郭金興,張凱旋,張臘梅,戴志穎,陳廣新*

    ?

    頸動脈分叉的非穩(wěn)態(tài)數(shù)值模擬分析

    王汝良1,胡霖霖2,郭金興3,張凱旋2,張臘梅2,戴志穎2,陳廣新2*

    (1. 牡丹江醫(yī)學(xué)院附屬紅旗醫(yī)院影像科,黑龍江 牡丹江 157011;2. 牡丹江醫(yī)學(xué)院附屬紅旗醫(yī)院介入科, 黑龍江 牡丹江 157011;3. 牡丹江醫(yī)學(xué)院醫(yī)學(xué)影像學(xué)院,黑龍江 牡丹江 157011)

    本課題旨在研究頸動脈分叉心動周期內(nèi)的血液動力學(xué)參數(shù),如頸動脈的壁切應(yīng)力、壓力、血流速度及渦流變化及各力學(xué)參數(shù)對頸動脈斑塊形成的影響。選取患者的頸動脈CT血管成像數(shù)據(jù)構(gòu)建有限元模型,通過計(jì)算流體力學(xué)軟件仿真計(jì)算。在設(shè)定了初始條件的情況下,獲得了頸動脈的血液動力學(xué)參數(shù),得出了頸動脈壁面剪切應(yīng)力及壁面壓力的的血流規(guī)律,并對頸動脈的壁切應(yīng)力、壓力、血流速度及渦流變化進(jìn)行了詳細(xì)的分析。頸動脈竇部位的壁面剪切應(yīng)力最小,血管壁的壁面壓力存在“負(fù)壓”效應(yīng),負(fù)壓效應(yīng)導(dǎo)致血流速度趨緩、血管直徑減小、血流量減小、血流阻力增加。在脈動周期一定時(shí)刻血流在頸動脈的分叉附近會形成渦流,渦流區(qū)內(nèi)血液低速流動且存在回流及二次流等復(fù)雜血流,實(shí)驗(yàn)結(jié)果符合血液的流體力學(xué)運(yùn)動方程的理論分析。

    頸動脈;有限元模型;渦流;血流動力學(xué)分析

    0 引言

    心腦血管疾病嚴(yán)重威脅人類健康,給社會和家庭造成了巨大的生活與財(cái)務(wù)負(fù)擔(dān)[1],是老年人致死的主要病因之一,危害性僅次于惡性腫瘤。我國腦卒發(fā)病率的增長速度高達(dá)8.7%,據(jù)文獻(xiàn)記載,10~20%的腦卒中是由于血管內(nèi)粥樣硬化的斑塊脫落所引起的[2,3],研究發(fā)現(xiàn),動脈硬化及血管內(nèi)膜的增厚,血液運(yùn)輸通道變狹窄,血液流速增加,對管壁的沖擊增大,管壁破損后形成斑塊,最終斑塊 脫落會形成血栓,而頸動脈處斑塊脫落可導(dǎo)致腦卒中[4,5]。病理解剖發(fā)現(xiàn)斑塊具有病灶選擇性,好發(fā)于頸動脈竇及頸動脈分叉等部位。管壁壓力、血流速度、回流、渦流、二次流、震蕩剪切因子(Oscillatory shear factor,OSI)及低平均壁面剪切應(yīng)力(Low average wall shear stress,TAWSS)是斑塊的血流動力學(xué)成因。粒子滯留時(shí)間(Particle retention time,RRT)的長短將會直接影響動脈粥樣硬化的發(fā)生和發(fā)展[6,10]。Malek等人對血管的壁面切應(yīng)力進(jìn)行了詳細(xì)的研究,研究結(jié)果顯示壁面切應(yīng)力的臨界值低于0.4 Pa時(shí)會發(fā)生動脈粥樣硬化[8,9]。因此,研究頸動脈的血液流體力學(xué)因素具有重要的臨床意義。

    本實(shí)驗(yàn)利用實(shí)際病人的頸動脈CT血管成像(CT angiography,CTA)數(shù)據(jù)構(gòu)建有限元模型,通過計(jì)算流體力學(xué)(Computational fluid mechanics,CFD)軟件進(jìn)行仿真計(jì)算,獲得頸動脈在心動周期內(nèi)的血液動力學(xué)參數(shù),分析頸動脈的壁切應(yīng)力、壓力、血流速度及渦流變化情況。

    1 材料與方法

    1.1 資料

    圖像數(shù)據(jù):數(shù)據(jù)來源于一男性64歲患者,選取80張斷層厚度為0.5 mm,圖像矩陣為512×512的頸部CTA斷層掃描Dicom原始數(shù)據(jù)。儀器應(yīng)用東芝64排螺旋CT,CT掃描參數(shù):電壓120 V,電流250 mA。

    圖像后處理工作站:戴爾Precision T7810,Xeon E5-2609 v3 CPU、nVIDIA Quadro k2200圖像顯卡;圖像后處理及計(jì)算流體力學(xué)分析軟件:Mimics 20.0、3-matic 12.0、Ansys 16.0。

    1.2 方法

    1.2.1 建立有限元模型

    將80張Dicom格式CTA數(shù)據(jù)文件導(dǎo)入Mimics中應(yīng)用閾值分割、動態(tài)區(qū)域增長等方法,初步重建建頸動脈三維模型。將重建后的三維模型導(dǎo)入3-matic軟件中進(jìn)行光滑處理,去除細(xì)小分支,切平出口、入口面,最后完成頸動脈三維幾何模型(圖1),該模型包括頸總動脈(Common carotid artery,CCA)、頸內(nèi)動脈(Internal carotid artery,ICA)、頸外動脈(External carotid artery,ECA)。將該模型以stl格式導(dǎo)入Ansys icem cfd劃分網(wǎng)格,邊界層實(shí)施五層棱柱加密以保證計(jì)算的精度。網(wǎng)格類型為四面體,網(wǎng)格單元節(jié)點(diǎn)共計(jì)218524個(gè),網(wǎng)格單元為121872個(gè)。將劃分之后的網(wǎng)格導(dǎo)入Ansys fluent 16.0中轉(zhuǎn)換為多面體網(wǎng)格,可以進(jìn)一步提升計(jì)算精度。圖1中A、B、C、D、E、F為選定的位置觀測點(diǎn)。

    圖1 頸動脈分叉幾何模型

    1.2.2 血流控制方程

    本研究假定血液為不可壓縮、粘性的牛頓流體,血管壁為剛性壁,控制方程為:

    式中為速度矢量,p為流場壓力,ρ為血流密度,μ為動力粘度。本文中,μ = 0.00356 Pa×s,ρ = 1060 kg/m3。

    1.2.3 邊界條件

    本文研究采用瞬態(tài)計(jì)算,入口速度為隨時(shí)間變化的脈動曲線(圖2)。出口壓力設(shè)定為0,血管壁無滑動,各速度分量為0,不考慮重力影響。計(jì)算指定時(shí)間步長為0.01 s,為獲得穩(wěn)定周期結(jié)果,共進(jìn)行三個(gè)心動周期計(jì)算,對第三個(gè)周期的計(jì)算結(jié)果進(jìn)行分析。

    圖2 入口的血流波形

    2 結(jié)果分析

    2.1 壁面剪切應(yīng)力(WSS)分析

    研究結(jié)果表明[1,3],動脈粥樣硬化一般發(fā)生在血管壁面的剪切應(yīng)力低于1 Pa的區(qū)域。圖3是位置觀測點(diǎn)A、B、C、D的壁面切應(yīng)力曲線圖。圖3結(jié)果顯示,頸動脈竇處(A點(diǎn))的壁面剪切應(yīng)力震蕩范圍為0~0.12 Pa,始終保持在較低的波動范圍。頸動脈B點(diǎn)的壁面剪切應(yīng)力在整個(gè)心動周期內(nèi),始終遠(yuǎn)高于A點(diǎn)。頸動脈D點(diǎn)的壁面剪切應(yīng)力在整個(gè)心動周期內(nèi),始終遠(yuǎn)高于C點(diǎn)。B、C、D、E各點(diǎn)的壁切應(yīng)力變化趨勢與A點(diǎn)保持一致,并隨時(shí)間震蕩。圖4是一個(gè)心動周期內(nèi)典型時(shí)刻壁面剪切力的變化情況,圖4表明,頸動脈竇部位的壁面剪切應(yīng)力分布始終是整個(gè)頸動脈區(qū)域最小的,而臨床研究已經(jīng)證實(shí),動脈粥樣硬化易發(fā)生在低壁面切應(yīng)力處[7]。

    圖3 A、B、C、D各點(diǎn)壁面剪切應(yīng)力

    圖4 不同時(shí)刻剪切應(yīng)力變化情況

    2.2 壁面壓力分析

    圖5為頸動脈B、E、F壁面壓強(qiáng)隨時(shí)間的變化曲線,圖6為整個(gè)頸動脈各個(gè)典型時(shí)刻的壓強(qiáng)對時(shí)間的變化曲線。實(shí)驗(yàn)結(jié)果顯示,在時(shí)間段=0.08 s到=0.3 s(心動收縮減速期),血管壁包括B、E、F點(diǎn)的壁面壓力存在負(fù)值現(xiàn)象即“負(fù)壓”效應(yīng),與入口血流速度降低正相關(guān)。在低壁面壓力作用下,這種“負(fù)壓”效應(yīng)導(dǎo)致血流速度趨緩、血管直徑減小、血流量減小、血流阻力增加及腦部供血不足,繼而引發(fā)缺血性腦卒中,或誘發(fā)局部血管痙攣,致使腦缺血進(jìn)一步加重。此外,血流速度慢將導(dǎo)致血液中的脂質(zhì)附著血管壁時(shí)間增長,可促使粥樣斑塊的進(jìn)一步發(fā)展[12],當(dāng)初始條件完全設(shè)定的情況下該實(shí)驗(yàn)結(jié)論與血液的流體力學(xué)運(yùn)動方程的理論分析和求解完全吻合。

    圖5 B、E、F各點(diǎn)壁面壓力變化情況

    圖6 各典型時(shí)刻壓力變化情況

    2.3 血流分析

    頸動脈血流線圖結(jié)果顯示(圖5),在心動收縮加速期,頸總動脈(CCA)、頸內(nèi)動脈(ICA)、頸外動脈(ECA)的血流速度相差不多,方向與軸向平行。在=0.08 s時(shí),入口速度達(dá)到峰值,CCA、ICA、ECA的血流速度都達(dá)到最高值。在ICA的A、B(分叉處)、C點(diǎn)附近區(qū)域速度較低,血流停滯面積較小,血細(xì)胞與管壁之間存在碰撞和摩擦。在心動收縮期內(nèi)血流速度劇烈變化,對頸動脈結(jié)構(gòu)和功能會產(chǎn)生影響,有利于血液中脂質(zhì)等成分沉積,為斑塊形成提供了條件[12]。在脈動收縮期末(=0.3 s),頸動脈整體血流速度最低;在舒張期,血流速度保持管壁產(chǎn)生作用,尤其是頸動脈竇附近區(qū)域、分叉部位,更有利于斑塊的形成、發(fā)展[14,15]。

    圖7是不同典型時(shí)刻的血流流線圖。=0.06 s,=0.08 s為心動周期中的收縮上升期,=0.09s,= 0.1 s,=0.12s,=0.16 s,=0.21 s,=0.3 s為心動周期的收縮下降期。其余時(shí)刻為心動周期的舒張期。從收縮下降期開始,渦流逐漸形成,在0.3s時(shí)刻渦流最大,渦流位置發(fā)生在頸動脈竇部附近。從舒張期開始,渦流強(qiáng)度逐漸降低。研究結(jié)果顯示,渦流區(qū)內(nèi)血流速度小且存在回流及二次流等復(fù)雜血流,血液中脂質(zhì)及纖維蛋白易沉積于此,有利于斑塊的形成。此外,血流速度大小、方向的頻繁改變對斑塊的形成與發(fā)展也有促進(jìn)作用[13,16]。

    3 結(jié)論

    本課題對頸動脈在心動周期內(nèi)的血液動力學(xué)參數(shù),如頸動脈的壁切應(yīng)力、壓力、血流速度及渦流變化進(jìn)行了研究,同時(shí)探討了各力學(xué)參數(shù)對頸動脈斑塊形成的影響,實(shí)驗(yàn)結(jié)果與血液的流體力學(xué)運(yùn)動方程的理論分析和求解完全吻合。研究結(jié)果證實(shí),頸動脈竇部位的壁面剪切應(yīng)力分布是頸動脈區(qū)域最小的,動脈粥樣硬化易發(fā)生在低壁面切應(yīng)力處。血管壁的壁面壓力存在負(fù)值現(xiàn)象即“負(fù)壓”效應(yīng)。這種“負(fù)壓”效應(yīng)導(dǎo)致血流速度趨緩、血管直徑減小、血流量減小、血流阻力增加及腦部供血不足,繼而引發(fā)缺血性腦卒中,或誘發(fā)局部血管痙攣。在心動收縮期內(nèi)血流速度劇烈變化,對頸動脈結(jié)構(gòu)和功能會產(chǎn)生影響,為斑塊形成提供了條件。研究結(jié)果顯示,渦流區(qū)內(nèi)血流速度小且存在回流及二次流等復(fù)雜血流,血流速度大小、方向的頻繁改變對斑塊的形成與發(fā)展有一定的促進(jìn)作用。

    圖7 不同時(shí)刻流線圖

    [1] 袁煒, 陳忠利. 頸動脈血液動力學(xué)的數(shù)值模擬[J]. 中國組織工程研究. 2014, 18(42): 6785-6786.

    [2] Razavi A, Shirani E, Sadeghi MR. Numerical simulation of blood pulsatile flow in a stenosed carotid artery using different rheological models[J]. J Biomech. 2011; 44(11): 2021-2030.

    [3] 張耀楠, 李艷, 康雁. 基于CT影像的頸動脈分叉血流動力學(xué)特性分析[J]. 東北大學(xué)學(xué)報(bào)(自然科學(xué)版), 2014, 35(3): 356-357.

    [4] 周志尊, 王佳美, 胡明成, 陳廣新, 劉陽, 董默, 孫強(qiáng), 孫鵬, 劉佳維. 頸動脈三維重建及臨床應(yīng)用的研究[J]. 軟件, 2017, 38(8): 32-35.

    [5] Bai-Nan X, Fu-Yu W, Lei L, Xiao-Jun Z, Hai-Yue J. Hemodynamics model of fluid-solid interaction in internal carotid artery aneurysms[J]. Neurosurg Rev, 2011, 34(1): 39-47.

    [6] Karmonik C, Yen C, Diaz O, Klucznik R, Grossman RG, Benndorf G. Temporal variations of wall shear stress parameters in intracranial aneurysms-importance of patient-specific inflow waveforms for CFD calculations[J]. Acta neurochirurgica, 2010, 152(8): 1391-8.

    [7] Watanabe T, Isoda H, Takehara Y, Terada M, Naito T, Kosugi T, Onishi Y, Tanoi C, Izumi T. Hemodynamic vascular biomarkers for initiation of paraclinoid internalcarotidarteryaneurysms using patient-specific computational fluid dynamic simulation based on magnetic resonance imaging[J]. Neuroradiology, 2018, 60(5): 545-555.

    [8] Xu B, Zhong H, Duan S. Modeling of internalcarotidarteryaneurysm and blood flow simulation[J]. Technology and health care: official journal of the European Society for Engineering and Medicine, 2015, 23(1): 43-8.

    [9] Ren G, Cao X, Wang D, Jiang P, Li Y, Pei B. hemodynamic numerical simulationofcarotidartery aneurysm before and after surgery based on CT date[J]. Journal of biomedical engineering, 2014, 31(2): 341-6.

    [10] Wang F, Xu B, Sun Z, Wu C, Zhang X. Wall shear stress in intracranial aneurysms and adjacent arteries. Neural regeneration research[J]. 2013, 8(11): 1007-15.

    [11] Lal BK, Beach KW, Sumner DS. Intracranial collateralization determineshemodynamicforces forcarotid plaque disruption. Journal of vascular surgery, 2011, 54(5): 1461-71.

    [12] De Wilde D, Trachet B, Debusschere N, Iannaccone F, Swillens A, Degroote J, Vierendeels J, De Meyer GRY, Segers P. Assessment of shear stress related parameters in thecarotidbifurcation using mouse-specific FSI simulations[J]. Journal of biomechanics, 2016, 49(11): 2135-2142.

    [13] Saho T, Onishi H. Evaluation of the impact ofcarotidarterybifurcationangle on hemodynamics by use of computational fluid dynamics: asimulationand volunteer study[J]. Radio-logical physics and technology, 2016, 9(2): 277-85.

    [14] Gharahi H, Zambrano BA, Zhu DC, DeMarco JK, Baek S. Computational fluid dynamicsimulationof humancarotidarterybifurcationbased on anatomy and volumetric blood flow rate measured with magnetic resonance imaging[J]. International journal of advances in engineering sciences and applied math-ematics, 2016, 8(1): 40-60.

    [15] Younis HF, Kaazempur-Mofrad MR, Chan RC, Isasi AG, Hinton DP, Chau AH, Kim LA, Kamm RD. Hemodynamics and wall mechanics in humancarotidbifurcationand its conseq-uences for atherogenesis: investigation of inter-individual vari-ation[J]. Bio-m-echanics and modeling in mechanobiology 2004, 3(1): 17-32.

    [16] Younis HF, Kaazempur-Mofrad MR, Chung C, Chan RC, Kamm RD. Computational analysis of the effects of exercise on hemodynamics in thecarotidbifurcation[J]. Annals of biomed-ical engineering, 2003, 31(8): 995-1006.

    Carotid artery Bifurcation Unsteady State Numerical Simulation Analysis

    WANG Ru-liang1, HU Lin-lin2, GUO Jin-xing3, ZHANG Kai-xuan2, ZHANG La-mei2, DAI Zhi-ying2, CHEN Guang-xin2*

    (1. Department of radiology, Hongqi Hospital Affiliated to Mudanjiang medical University, Mudanjiang Heilongjiang 157011; 2. Medical image college of Mudanjiang medical university,Mudanjiang Heilongjiang 157011; 3. Department of intervention, Hongqi Hospital Affiliated to Mudanjiang medical University, Mudanjiang Heilongjiang 157011)

    The purpose of this project is to study carotid artery bifurcation hemodynamic parameters in the carotid cardiac cycle, such as wall shear stress, pressure, blood flow velocity and eddy current velocity and the influence of various mechanical parameters on the formation of carotid plaque. A finite element model was constructed based on patient’s CT angiographic data of the carotid artery. Through the simulation calculation of fluid mechanics software, the flow line was designed in the experimental to observe the blood turbulence.Results: Under the initial conditions, the hemodynamic parameters of the carotid artery were obtained, the shear stress of the carotid wall and the blood flow pattern were obtained, and the wall shear stress, pressure, blood flow velocity and eddy flow changes of the carotid artery were analyzed in detail. The wall shear stress of the carotid sinus is the smallest, and the wall pressure of the vascular wall has the effect of "negative pressure", which slow down the blood flow rate, decrease blood vessel diameter, and increase blood flow resistance. When the velocity reaches a specific value,the blood flow forms vortex near the carotid bifurcation. Eddy blood flow velocity in this area is small and there is reflux and secondary flow. The experimental results conform to the fluid mechanics equations of the blood of theoretical analysis.

    Carotid artery; Finite element model; Eddy current; Hemodynamic analysis

    TP319

    A

    10.3969/j.issn.1003-6970.2018.10.008

    黑龍江省自然基金項(xiàng)目(No. H2015082);黑龍江省教育廳項(xiàng)目(No.12541846);黑龍江省大學(xué)生創(chuàng)業(yè)項(xiàng)目(No.201646)

    王汝良(1969-),教授、主任醫(yī)師,研究方向:影像與核醫(yī)學(xué);胡霖霖(1993-),碩士研究生,研究方向:醫(yī)學(xué)圖像處理;郭金興(1984-),主管護(hù)師,研究方向:介入治療;張凱旋(1993-),碩士研究生,研究方向:醫(yī)學(xué)圖像處理;張臘梅(1993-),碩士研究生,研究方向:超聲診斷;戴志穎(1996-),學(xué)生,研究方向:醫(yī)學(xué)圖像處理。

    陳廣新(1978-),講師,研究方向:人工智能與醫(yī)學(xué)圖像處理。牡丹江醫(yī)學(xué)院醫(yī)學(xué)影像學(xué)院醫(yī)學(xué)圖像處理教研室。

    王汝良,胡霖霖,郭金興,等. 頸動脈分叉的非穩(wěn)態(tài)數(shù)值模擬分析[J]. 軟件,2018,39(10):36-41

    猜你喜歡
    剪切應(yīng)力心動渦流
    您了解心動過緩嗎
    基于CFD仿真分析的各缸渦流比一致性研究
    是心動啊,糟糕,躲不掉雪中的清華路156號
    心瓣瓣膜區(qū)流場中湍流剪切應(yīng)力對瓣膜損害的研究進(jìn)展
    渦流傳感器有限元仿真的研究與實(shí)施
    剪切應(yīng)力對聚乳酸結(jié)晶性能的影響
    中國塑料(2016年6期)2016-06-27 06:34:24
    上一次令編輯們心動的事
    Coco薇(2015年5期)2016-03-29 23:44:49
    “嗔哥”心動一日美食
    Coco薇(2015年3期)2015-12-24 03:26:08
    動脈粥樣硬化病變進(jìn)程中血管細(xì)胞自噬的改變及低剪切應(yīng)力對血管內(nèi)皮細(xì)胞自噬的影響*
    硫化氫在低剪切應(yīng)力導(dǎo)致內(nèi)皮細(xì)胞自噬障礙中的作用
    亚洲最大成人手机在线| 亚洲精品久久久久久婷婷小说 | 色噜噜av男人的天堂激情| 亚洲av成人av| 午夜爱爱视频在线播放| 又黄又爽又刺激的免费视频.| 久久鲁丝午夜福利片| 国产亚洲5aaaaa淫片| 精品午夜福利在线看| 99久久无色码亚洲精品果冻| 99久国产av精品| 亚洲欧美成人精品一区二区| 日本五十路高清| 亚洲av成人av| 久久久久九九精品影院| 亚洲av成人精品一区久久| 伦精品一区二区三区| 欧美精品一区二区大全| 久久久久久久午夜电影| 久久精品国产亚洲网站| av在线亚洲专区| 男女啪啪激烈高潮av片| 亚洲18禁久久av| 最后的刺客免费高清国语| 亚洲人成网站在线播| 美女被艹到高潮喷水动态| 午夜免费男女啪啪视频观看| 老女人水多毛片| 日本三级黄在线观看| 婷婷色av中文字幕| 尤物成人国产欧美一区二区三区| 中国美白少妇内射xxxbb| 在线a可以看的网站| 看片在线看免费视频| 韩国av在线不卡| 禁无遮挡网站| 在现免费观看毛片| 亚洲最大成人手机在线| 午夜精品在线福利| 熟女人妻精品中文字幕| 在线免费观看不下载黄p国产| 美女cb高潮喷水在线观看| 日韩欧美精品免费久久| 免费观看的影片在线观看| 亚洲最大成人手机在线| 国产精品久久久久久精品电影小说 | 亚洲欧美一区二区三区国产| 精品国内亚洲2022精品成人| 欧美三级亚洲精品| 天堂√8在线中文| 久久久色成人| 看片在线看免费视频| ponron亚洲| 中文精品一卡2卡3卡4更新| 国语对白做爰xxxⅹ性视频网站| 亚洲av电影不卡..在线观看| 美女脱内裤让男人舔精品视频| 搡女人真爽免费视频火全软件| 免费大片18禁| 99热这里只有是精品在线观看| 成人午夜精彩视频在线观看| 中文天堂在线官网| 久久久久久国产a免费观看| 成人性生交大片免费视频hd| 国产又色又爽无遮挡免| 91午夜精品亚洲一区二区三区| 一边亲一边摸免费视频| 日本欧美国产在线视频| 国产激情偷乱视频一区二区| 精品人妻偷拍中文字幕| 在线观看美女被高潮喷水网站| 久久精品国产亚洲网站| 国产毛片a区久久久久| 国产精品无大码| 少妇裸体淫交视频免费看高清| a级毛色黄片| 免费av观看视频| 亚洲激情五月婷婷啪啪| 狂野欧美激情性xxxx在线观看| 能在线免费观看的黄片| 欧美日本视频| 男人舔女人下体高潮全视频| 最后的刺客免费高清国语| 中文字幕人妻熟人妻熟丝袜美| 黄色日韩在线| 精品午夜福利在线看| 亚洲欧美日韩高清专用| 波多野结衣巨乳人妻| 国产免费福利视频在线观看| 中文精品一卡2卡3卡4更新| 久久久久网色| 在线播放无遮挡| 亚洲精品,欧美精品| 青春草亚洲视频在线观看| 国产乱来视频区| 国产精品不卡视频一区二区| 五月伊人婷婷丁香| 插阴视频在线观看视频| 美女被艹到高潮喷水动态| 天天躁夜夜躁狠狠久久av| 国产精品三级大全| 国产美女午夜福利| 国产探花极品一区二区| 搡女人真爽免费视频火全软件| 亚洲国产最新在线播放| 一级毛片aaaaaa免费看小| 22中文网久久字幕| 最近视频中文字幕2019在线8| 韩国av在线不卡| 高清视频免费观看一区二区 | 中国国产av一级| 青春草国产在线视频| 中国美白少妇内射xxxbb| 村上凉子中文字幕在线| 99久久无色码亚洲精品果冻| 成年版毛片免费区| 热99在线观看视频| 一区二区三区免费毛片| 国产日韩欧美在线精品| 中文字幕av成人在线电影| 久久久久久大精品| 亚洲国产成人一精品久久久| 全区人妻精品视频| 99在线人妻在线中文字幕| 久99久视频精品免费| 秋霞伦理黄片| av国产免费在线观看| av视频在线观看入口| 久久精品国产亚洲av天美| 人体艺术视频欧美日本| 日韩三级伦理在线观看| 人妻系列 视频| 搞女人的毛片| 久久久久精品久久久久真实原创| 大香蕉97超碰在线| 欧美+日韩+精品| 最新中文字幕久久久久| 男人的好看免费观看在线视频| 夜夜爽夜夜爽视频| 久久国产乱子免费精品| 欧美成人a在线观看| 丝袜喷水一区| 日本五十路高清| 超碰av人人做人人爽久久| 亚洲av不卡在线观看| 免费av不卡在线播放| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲伊人久久精品综合 | .国产精品久久| 又粗又硬又长又爽又黄的视频| 午夜激情欧美在线| 老司机影院毛片| 五月玫瑰六月丁香| 欧美丝袜亚洲另类| 波野结衣二区三区在线| 精品久久久久久久久av| 99热这里只有是精品在线观看| 久久久久国产网址| 国产精品一区二区三区四区久久| 熟女人妻精品中文字幕| 久久鲁丝午夜福利片| www.色视频.com| 禁无遮挡网站| 91狼人影院| eeuss影院久久| videossex国产| 亚洲三级黄色毛片| 2021天堂中文幕一二区在线观| 色播亚洲综合网| 久久久a久久爽久久v久久| 人妻夜夜爽99麻豆av| 综合色av麻豆| 久久人人爽人人片av| 欧美成人午夜免费资源| 亚洲一级一片aⅴ在线观看| 亚洲av一区综合| 97在线视频观看| 一区二区三区乱码不卡18| 国产亚洲一区二区精品| videos熟女内射| 国产午夜精品论理片| 伦理电影大哥的女人| 亚洲av日韩在线播放| 深夜a级毛片| 我的女老师完整版在线观看| 国产高清国产精品国产三级 | 日韩亚洲欧美综合| 精品99又大又爽又粗少妇毛片| 国产淫语在线视频| 精品国产一区二区三区久久久樱花 | 色综合色国产| 久久久久久久午夜电影| 久久久a久久爽久久v久久| 亚洲激情五月婷婷啪啪| 国模一区二区三区四区视频| 亚洲性久久影院| 欧美日韩精品成人综合77777| 一本一本综合久久| 亚洲国产精品合色在线| 少妇的逼好多水| 国产精品美女特级片免费视频播放器| 国产淫语在线视频| 纵有疾风起免费观看全集完整版 | 国产三级在线视频| 国产精品一及| 久久久久九九精品影院| videos熟女内射| av在线蜜桃| 五月玫瑰六月丁香| 国产老妇伦熟女老妇高清| 国产精品麻豆人妻色哟哟久久 | 三级男女做爰猛烈吃奶摸视频| 精品国产一区二区三区久久久樱花 | 观看免费一级毛片| 97超碰精品成人国产| 午夜福利高清视频| 亚洲美女视频黄频| 校园人妻丝袜中文字幕| 99热这里只有是精品在线观看| 免费人成在线观看视频色| 日本黄色视频三级网站网址| 色噜噜av男人的天堂激情| 免费观看人在逋| 99热这里只有是精品在线观看| 中国美白少妇内射xxxbb| 波多野结衣高清无吗| 国产一区二区在线av高清观看| 一本一本综合久久| 午夜久久久久精精品| 国产精品.久久久| 久久精品国产自在天天线| 99在线人妻在线中文字幕| 国产片特级美女逼逼视频| 久久人妻av系列| 伦精品一区二区三区| av在线观看视频网站免费| 久久人人爽人人爽人人片va| 久久综合国产亚洲精品| 国产精品一二三区在线看| 成人高潮视频无遮挡免费网站| 亚洲伊人久久精品综合 | 国产色爽女视频免费观看| 蜜臀久久99精品久久宅男| 日韩欧美精品免费久久| 国产精品一区二区三区四区免费观看| 国产视频首页在线观看| 高清在线视频一区二区三区 | 免费观看人在逋| 亚洲在线自拍视频| 日韩制服骚丝袜av| 一卡2卡三卡四卡精品乱码亚洲| 成人毛片a级毛片在线播放| 日日摸夜夜添夜夜爱| 久久这里有精品视频免费| 精品一区二区免费观看| 三级国产精品片| 一级黄色大片毛片| 少妇裸体淫交视频免费看高清| 亚洲在线自拍视频| 亚洲国产精品专区欧美| 老司机福利观看| 国产三级中文精品| 男人和女人高潮做爰伦理| 亚洲最大成人中文| 舔av片在线| 日本一本二区三区精品| 成人国产麻豆网| 少妇熟女欧美另类| 麻豆久久精品国产亚洲av| 日日干狠狠操夜夜爽| 欧美一级a爱片免费观看看| 亚洲精品成人久久久久久| 简卡轻食公司| 国产一级毛片七仙女欲春2| 视频中文字幕在线观看| 国产精品人妻久久久久久| 色综合亚洲欧美另类图片| 建设人人有责人人尽责人人享有的 | 在线观看一区二区三区| 麻豆精品久久久久久蜜桃| 国产熟女欧美一区二区| 天天躁日日操中文字幕| 亚洲欧美日韩东京热| 精品久久国产蜜桃| 日本一本二区三区精品| 亚洲国产精品专区欧美| 国产久久久一区二区三区| 日本五十路高清| 老司机影院成人| 色综合色国产| 久久99热这里只频精品6学生 | 久久久久久伊人网av| 大香蕉97超碰在线| 美女脱内裤让男人舔精品视频| 91在线精品国自产拍蜜月| 亚洲三级黄色毛片| 国产成人精品婷婷| 天堂影院成人在线观看| 少妇丰满av| 麻豆精品久久久久久蜜桃| 亚洲av日韩在线播放| 卡戴珊不雅视频在线播放| 国产精品永久免费网站| 日韩av不卡免费在线播放| 久久韩国三级中文字幕| 国产日韩欧美在线精品| av国产免费在线观看| 青春草国产在线视频| 1000部很黄的大片| 国产亚洲av嫩草精品影院| 成人高潮视频无遮挡免费网站| 国产精品爽爽va在线观看网站| 精品久久久久久成人av| 日本黄大片高清| 精品国产三级普通话版| 校园人妻丝袜中文字幕| 国产麻豆成人av免费视频| 久久精品国产亚洲av涩爱| 久久精品国产自在天天线| 一区二区三区乱码不卡18| 中文字幕制服av| 夫妻性生交免费视频一级片| 亚洲怡红院男人天堂| 一边摸一边抽搐一进一小说| 九九爱精品视频在线观看| 免费看光身美女| 日本wwww免费看| 欧美zozozo另类| 99久久人妻综合| av天堂中文字幕网| 一区二区三区高清视频在线| 国产精品人妻久久久影院| 亚洲成人久久爱视频| 一个人看视频在线观看www免费| 久久婷婷人人爽人人干人人爱| 大香蕉97超碰在线| 国产成人午夜福利电影在线观看| 国产三级在线视频| 搡老妇女老女人老熟妇| 男女那种视频在线观看| 国产 一区精品| 亚洲最大成人手机在线| 国产精品电影一区二区三区| 亚洲国产欧洲综合997久久,| 国产日韩欧美在线精品| 长腿黑丝高跟| 精品久久久噜噜| 亚洲综合色惰| videos熟女内射| 亚洲欧美成人综合另类久久久 | www.色视频.com| 精品久久久久久成人av| 欧美日韩精品成人综合77777| 国产亚洲一区二区精品| 人人妻人人澡欧美一区二区| 午夜激情欧美在线| 婷婷色麻豆天堂久久 | 看非洲黑人一级黄片| 大香蕉97超碰在线| 九九在线视频观看精品| 国产免费视频播放在线视频 | 黄片wwwwww| 免费看av在线观看网站| 中文欧美无线码| 99热全是精品| 免费看a级黄色片| 国产视频内射| 亚洲真实伦在线观看| 午夜激情欧美在线| 2022亚洲国产成人精品| 久久人人爽人人片av| 久久久午夜欧美精品| 亚洲欧美一区二区三区国产| 成人毛片a级毛片在线播放| 一本一本综合久久| 边亲边吃奶的免费视频| 99在线视频只有这里精品首页| 国产av码专区亚洲av| 精品久久久久久久末码| av天堂中文字幕网| 18禁裸乳无遮挡免费网站照片| 一个人看的www免费观看视频| 天天躁日日操中文字幕| 色网站视频免费| 最近视频中文字幕2019在线8| 国产高清三级在线| 91久久精品国产一区二区成人| 热99re8久久精品国产| 精品人妻一区二区三区麻豆| 亚洲精品,欧美精品| 久久精品综合一区二区三区| 中文字幕精品亚洲无线码一区| 熟妇人妻久久中文字幕3abv| 黄片wwwwww| 国产片特级美女逼逼视频| 搞女人的毛片| 亚洲av男天堂| 国产黄色视频一区二区在线观看 | 日日撸夜夜添| 国产亚洲一区二区精品| 在线播放无遮挡| 日韩大片免费观看网站 | 麻豆成人av视频| 久久久久久伊人网av| 桃色一区二区三区在线观看| 夫妻性生交免费视频一级片| 国产精品一及| 干丝袜人妻中文字幕| 婷婷色麻豆天堂久久 | 久久这里只有精品中国| av黄色大香蕉| 国产黄a三级三级三级人| 联通29元200g的流量卡| 男女国产视频网站| 国产精品不卡视频一区二区| 一区二区三区四区激情视频| 欧美最新免费一区二区三区| 亚洲va在线va天堂va国产| 精品国产一区二区三区久久久樱花 | 精品99又大又爽又粗少妇毛片| 男女下面进入的视频免费午夜| 老司机福利观看| 热99在线观看视频| 久久久精品大字幕| 麻豆国产97在线/欧美| 亚洲欧美成人精品一区二区| 久久6这里有精品| 久久久久久久久久久免费av| 午夜免费激情av| 国产成人精品久久久久久| 国产真实乱freesex| 日韩av不卡免费在线播放| 99热全是精品| 欧美日韩在线观看h| 久久精品国产99精品国产亚洲性色| 亚洲精品成人久久久久久| 成年av动漫网址| 亚洲怡红院男人天堂| 国产亚洲av片在线观看秒播厂 | 成人漫画全彩无遮挡| 亚洲在久久综合| 69av精品久久久久久| 久久久久久久久久久免费av| 禁无遮挡网站| 男人舔奶头视频| 欧美丝袜亚洲另类| 黄片wwwwww| 久久精品国产鲁丝片午夜精品| 国产精品国产三级国产av玫瑰| 免费观看精品视频网站| 久久久精品欧美日韩精品| 国产精品乱码一区二三区的特点| 一区二区三区免费毛片| 国产精品一区二区性色av| 最近2019中文字幕mv第一页| 人人妻人人澡欧美一区二区| 久久精品国产亚洲av天美| 日韩大片免费观看网站 | 国产精品一及| 啦啦啦韩国在线观看视频| 日韩强制内射视频| 亚洲欧美日韩东京热| 精品久久久久久久末码| 亚洲av日韩在线播放| 最后的刺客免费高清国语| 岛国毛片在线播放| 六月丁香七月| 麻豆久久精品国产亚洲av| 最后的刺客免费高清国语| 亚洲av二区三区四区| 亚洲精品影视一区二区三区av| 日本黄色视频三级网站网址| 最近的中文字幕免费完整| 91aial.com中文字幕在线观看| 色尼玛亚洲综合影院| 亚洲精品色激情综合| 欧美变态另类bdsm刘玥| 免费观看a级毛片全部| 日本一二三区视频观看| 久久久久久九九精品二区国产| 99热6这里只有精品| 国产一区二区在线观看日韩| 久久精品国产亚洲av天美| 精品午夜福利在线看| 日日摸夜夜添夜夜添av毛片| 色播亚洲综合网| 免费观看在线日韩| 欧美xxxx性猛交bbbb| 久久国产乱子免费精品| 国产在线一区二区三区精 | 亚洲美女搞黄在线观看| 爱豆传媒免费全集在线观看| 中文资源天堂在线| av卡一久久| 国产午夜福利久久久久久| 久久久久性生活片| 亚洲欧洲国产日韩| 国产高清三级在线| 中国国产av一级| 黄色日韩在线| 狂野欧美白嫩少妇大欣赏| 国产午夜精品久久久久久一区二区三区| 国产精品永久免费网站| 麻豆国产97在线/欧美| 又黄又爽又刺激的免费视频.| 国产精品久久久久久精品电影| av又黄又爽大尺度在线免费看 | 欧美极品一区二区三区四区| 91午夜精品亚洲一区二区三区| 国产片特级美女逼逼视频| 在线a可以看的网站| 天天躁日日操中文字幕| 久久久久久久久久成人| 亚洲av一区综合| 国产亚洲精品av在线| 天堂√8在线中文| 国产成人a区在线观看| 美女国产视频在线观看| 99久久人妻综合| 亚洲av成人精品一二三区| 亚洲人成网站在线播| 久久精品国产自在天天线| 国产v大片淫在线免费观看| 亚洲国产精品sss在线观看| 久热久热在线精品观看| 久久国产乱子免费精品| 国产精品蜜桃在线观看| 亚洲精品成人久久久久久| 国产91av在线免费观看| 久久精品国产自在天天线| 欧美成人午夜免费资源| 日韩精品有码人妻一区| 欧美高清性xxxxhd video| 少妇丰满av| 亚洲国产成人一精品久久久| 亚洲一级一片aⅴ在线观看| 亚洲无线观看免费| 一本久久精品| 久久99精品国语久久久| 久久久a久久爽久久v久久| 久久综合国产亚洲精品| 99视频精品全部免费 在线| 99热全是精品| 美女黄网站色视频| 国产精品不卡视频一区二区| 久久久久九九精品影院| 国产一区二区三区av在线| 国产又黄又爽又无遮挡在线| 美女xxoo啪啪120秒动态图| 欧美区成人在线视频| 精品无人区乱码1区二区| 亚洲综合色惰| 91精品一卡2卡3卡4卡| 亚洲精品影视一区二区三区av| 在线观看av片永久免费下载| 中文字幕亚洲精品专区| 亚洲内射少妇av| 国产毛片a区久久久久| 乱系列少妇在线播放| 26uuu在线亚洲综合色| 精品国内亚洲2022精品成人| 一个人观看的视频www高清免费观看| 纵有疾风起免费观看全集完整版 | 黄片无遮挡物在线观看| 亚洲天堂国产精品一区在线| 免费观看人在逋| 天堂中文最新版在线下载 | av在线天堂中文字幕| 国产av码专区亚洲av| 99热这里只有是精品在线观看| 日韩大片免费观看网站 | 日本猛色少妇xxxxx猛交久久| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲最大成人手机在线| 欧美bdsm另类| 成人美女网站在线观看视频| 狠狠狠狠99中文字幕| 精品一区二区三区人妻视频| 国产单亲对白刺激| 国产成人精品婷婷| 国产精品国产高清国产av| 国产淫片久久久久久久久| 在线观看一区二区三区| 成人午夜精彩视频在线观看| 性色avwww在线观看| 美女大奶头视频| av在线老鸭窝| av视频在线观看入口| 欧美一区二区精品小视频在线| 蜜桃久久精品国产亚洲av| 淫秽高清视频在线观看| 精品国产露脸久久av麻豆 | 高清在线视频一区二区三区 | 国产精品野战在线观看| 美女内射精品一级片tv| 日产精品乱码卡一卡2卡三| 一级爰片在线观看| 国产黄色小视频在线观看| 亚洲精品aⅴ在线观看| 成人亚洲精品av一区二区| 一级毛片电影观看 | av视频在线观看入口| 22中文网久久字幕| 国产av一区在线观看免费| 国产高清国产精品国产三级 | 久久精品久久久久久久性| 日韩欧美三级三区| 日韩,欧美,国产一区二区三区 | 久久精品国产自在天天线| 精品久久久久久久人妻蜜臀av| 国产免费视频播放在线视频 | 青青草视频在线视频观看| 国产精品一区二区性色av| 国产一级毛片在线| a级毛片免费高清观看在线播放|