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

    醫(yī)學有限元仿真實驗系統(tǒng)研究

    2020-11-30 09:02:22金山馮雨欣王玥涵曲佳欣董寧周成鈺徐會敏宋海南
    軟件 2020年9期
    關鍵詞:生物力學有限元

    金山 馮雨欣 王玥涵 曲佳欣 董寧 周成鈺 徐會敏 宋海南

    摘? 要: 將有限元分析應用于醫(yī)學研究,可以解決很多物理學在醫(yī)學應用上需要解決卻難以解決的復雜問題,為醫(yī)學研究及臨床治療提供理論指導和科學依據。在醫(yī)學院校開設醫(yī)學有限元仿真實驗課程對學生學習、教師教學科研、學校與附屬醫(yī)院合作有著十分重要的意義。

    關鍵詞: 有限元;CFD;計算機建模與仿真;生物力學

    中圖分類號: TP311.52? ? 文獻標識碼: A? ? DOI:10.3969/j.issn.1003-6970.2020.09.020

    本文著錄格式:金山,馮雨欣,王玥涵,等. 醫(yī)學有限元仿真實驗系統(tǒng)研究[J]. 軟件,2020,41(09):7275

    【Abstract】: To solve the complex human health problems, computer biomechanical modeling and simulation experimental system can be used to model and make quantitative analysis of human tissues and organs, and establish a precise and standardized non-invasive detection and analysis technology. Based on the clinical case (image), using fluid mechanics and solid mechanics theory, and computer numerical simulation technology to solve human tissues and make organs quantitative analysis, and formulate individualized treatment scheme, accomplish biomechanical-design.Through the biomechanical modeling and simulation experiment system, the theory instruction to the clinic and the medical college clinical can be accomplished.

    【Key words】: Finite element; CFD; Simulation; Modeling and simulation; Biomechanics

    0? 引言

    隨著計算機技術的快速發(fā)展,應用有限元方法的計算機建模仿真方法(數值計算方法)與理論分析方法、實驗方法并稱為科學研究的三大主要方法[1]。人體處于力學環(huán)境之中,人體各系統(tǒng),如循環(huán)系統(tǒng)、運動系統(tǒng)、消化系統(tǒng)、呼吸系統(tǒng)和泌尿系統(tǒng)等生理活動均受力學因素影響。隨著臨床影像技術等生物醫(yī)學技術的發(fā)展,基于有限元方法的計算機生物力學建模與仿真在人體發(fā)育、生長和疾病發(fā)病機制以及個體化防治中發(fā)揮越來越重要的作用,將生物醫(yī)學基礎研究與力學數值模擬的定量化研究有機結合,體現(xiàn)學科交叉,為醫(yī)學研究及臨床治療提供理論指導和科學依據。在醫(yī)學院校開設基于有限元分析的醫(yī)學仿真實驗課程,對醫(yī)學生的學習、教師的教學科研、學校與醫(yī)院的合作有重要的意義。

    1? 基于有限元技術的醫(yī)學仿真實驗的優(yōu)勢

    有限元方法最早應用于骨科研究,開始于脊柱生物力學。幾十年來其在解決生物力學問題上得到了廣泛應用。醫(yī)學有限元仿真實驗是一個交叉研究領域,包含了生物力學、醫(yī)學圖像處理、臨床解剖學、虛擬現(xiàn)實以及有限元技術的多種學科技術[2]?;谟邢拊ǖ尼t(yī)學仿真實驗具有以下幾個方面的優(yōu)勢:相較于現(xiàn)場實驗,基于有限元方法的醫(yī)學仿真實驗可檢測全場數據、能夠實現(xiàn)理想條件的優(yōu)勢,現(xiàn)場實驗數據精度受儀器、測量方法、測量環(huán)境影響較大[3];利用有限元法進行的模擬實驗具有實驗時間短,費用少,可模擬復雜條件、力學性能全面,可重復性好。此外,可根據需要構建模型,對模型根據實驗條件仿真,可在不同實驗條件下模擬人體器官變形、應力/應變分布、內部能量變化、極限破壞分析等。

    有限元分析適用于生物醫(yī)學材料、內科、胸外科、腦外科、婦科、兒科、眼科、神經科等臨床研究與醫(yī)療器械設計,應用領域廣。在醫(yī)學院校開設醫(yī)學有限元仿真實驗教學,可以為臨床、康復、生物醫(yī)學工程等專業(yè)學生提供便捷的實驗平臺,尤其是為“醫(yī)工結合”創(chuàng)造了橋梁。此外,在仿真實驗的基礎上,結合生物醫(yī)學的理論研究,為臨床科室的診斷、預測上提供科學的指導,促進醫(yī)學院校與臨床醫(yī)院的教學結合。

    2? 醫(yī)學有限元實驗內容

    醫(yī)學有限元實驗內容是基于醫(yī)學院校的本科、研究生培養(yǎng)方案的指導下,結合相應的課程來建設。構建分層次、模塊化、共享的模擬仿真體系。目前,醫(yī)學有限元建模仿真實驗主要集中在血管系統(tǒng)、肌骨系統(tǒng)、視覺系統(tǒng)等與力學因素密切相關的人體組織和器官,緊密聯(lián)系臨床問題,以臨床病例(影像)為基礎,應用流體力學和固體力學理論、系統(tǒng)生物信息與控制理論,結合先進的流場和應力場測試和醫(yī)學影像技術、宏觀與微觀結合,動物實驗與力學模型及數值模擬相結合,對相關組織進行建模與定量分析,從而建立精確規(guī)范的無創(chuàng)檢測和分析技術,以及進行個體化治療方案的生物力學設計。

    2.1? 血管系統(tǒng)建模與仿真系統(tǒng)

    血管系統(tǒng)建模與仿真主要包括兩大部分:血管系統(tǒng)的生物力學建模與仿真;局部血管及植介入器械的生物力學建模與仿真。血管系統(tǒng)建模與仿真模塊的主要方法是借助于計算流體力學技術(computational fluid dynamics,CFD)的血流動力學數值模擬仿真,分析血液循環(huán)系統(tǒng)的血流動力學因素,評估、設計心血管植入/介入器械。動脈粥樣硬化多發(fā)生在人體動脈特殊的位置,如冠狀動脈、頸動脈分支、主動脈弓、腎動脈分支等。這些部位都伴有彎曲和分支,具有復雜的流體動力學問題,目前大量的實驗研究和仿真研究集中在這些特殊的幾何部位和特征部位[3-5]:

    (1)介入治療與外科手術仿真評價;

    (2)瓣膜動力學仿真;

    (3)動脈狹窄仿真;

    (4)動脈瘤(腹主動脈瘤和腦動脈瘤)仿真;

    (5)冠狀動脈分支仿真;

    (6)主動脈分支仿真:

    (7)頸動脈分支仿真:

    (8)人體主動脈仿真。

    血流動力學仿真分析需要用到下列技術:

    計算流體力學分析技術:控制流體運動的方程被稱為納維-斯托克斯(Navier-Stokes,N-S)方程。這些方程是根據質量守恒、動量守恒及能量守恒推導出來的。

    連續(xù)方程: 式中,u、v、w為速度在三個坐標上的分量;

    離散化技術:計算機仿真需要將物理流動域離散為網格形式的計算域。當前在使用中的離散化思想有:有限差分、有限體積、有限元法。

    時間積分:廣泛使用的時間相關技術主要分兩類:顯式和隱式。

    計算機建模與仿真技術:基于人體真實生理病理參數和醫(yī)學影像信息進行詳細的個體化建模。如借助高精度的CT掃描技術,可以對冠狀動脈等血管進行準確的三維模型重建,設置準確的邊界條件和材料物理屬性,利用計算機對模型進行準確的力學模擬[6-8]。如圖1所示,利用MIMICS軟件(MATERIALIST公司,比利時)對基于CT掃描獲得的冠狀動脈影像數據進行三維重建,獲得三維模型;如圖2所示,利用ANYS? FLUENT(ANSYS公司,美國)進行冠狀動脈數值模擬結果(速度流線圖)。

    2.2? 骨肌系統(tǒng)生物力學建模與仿真

    骨肌系統(tǒng)建模與仿真實驗思路為:緊密圍繞臨床問題,以臨床病例(影像)為基礎,應用固體力學理論、系統(tǒng)生物信息和控制理論,結合先進的應力場測試和醫(yī)學影像技術,宏觀與微觀相結合,對相關組織和器官進行建模與定量分析,從而建立精確規(guī)范的無創(chuàng)檢測和分析技術,以及進行個體化治療方案的生物力學設計[9-10]。當前,醫(yī)學影像技術對骨折、骨質疏松等諸多骨科疾病診斷非常重要,基于醫(yī)療成像設備獲取的二維圖像序列構建骨組織的三維幾何模型,在醫(yī)學診斷、手術規(guī)劃、醫(yī)學教學等方面有很高的應用價值。骨科植入體的構型優(yōu)化、介入宿主之間的力學關系、介入體的生物力學評價方面是骨肌系統(tǒng)建模仿真實驗的最主要的領域[11-12]。基于CAD(Computer Aided Design)的骨肌系統(tǒng)精細化處理方法,完成一個植介入體的仿真分析需要以下幾個方面:

    圖像獲?。豪门R床的成像設備如CT、MRI等獲取受試者的數據;

    圖像的三維重建:應用三維建模軟件構建骨肌系統(tǒng)的三維模型;

    模型的完善與CAD造型:應用三維逆向工程軟件進行完善模型,應用三維CAD軟件構建植入體、裝配組合體;

    有限元仿真計算:應用有限元仿真軟件進行仿真計算。圖3為肱骨遠端骨折切開復位固定術仿真分析應力云圖。

    骨肌系統(tǒng)本構關系:基于不同組織對力學刺激的響應有所差異,在相同的應力狀態(tài)下會呈現(xiàn)出不同的應變,本構關系體現(xiàn)其應力與應變的關系,體現(xiàn)了組織固有的材料特性。在骨肌系統(tǒng)建模與仿真中,常用的本構關系有線彈性模型,超彈性模型和粘彈性模型。表1為常用組織的線彈性本構參數[13-15];

    (1)生物力學模型的控制方程:(1)平衡方程:,其中為應力張量,為應力張量的散度,為單位體積上的外力矢量;

    (2)幾何方程:,其中為Cauchy應變張量,為位移矢量,與分別為位移的左右梯度;

    (3)本構方程:或W=W(),左式為應力關于應變的函數,右式為應變能密度關于應變的函數。

    目前骨肌系統(tǒng)的建模與仿真主要包括:下肢生物力學建模與仿真:髖關節(jié)建模與仿真;膝關節(jié)損傷建模與仿真;足踝損傷建模與仿真。脊柱生物力學建模與仿真:頸椎建模與仿真;腰椎建模與仿真。頭部生物力學與仿真:頭頸動力學仿真;眼球損傷生物力學仿真;口腔生物力學仿真;基于Mico-CT的骨微觀力學仿真。多剛體動力學建模與仿真:多剛體系統(tǒng)假人建模與仿真;胸腰段脊柱沖擊響應動力學研究;骨重建過程的建模預仿真。

    3? 醫(yī)學有限元仿真實驗軟件

    3.1? 醫(yī)學圖像處理軟件

    MIMICS:MIMICS是Materialise公司開發(fā)的一款具有交互式醫(yī)學影像控制系統(tǒng)[16]。具有高度模塊化的3D醫(yī)學圖像生成及編輯處理軟件,它支持導入各種掃描的醫(yī)學圖像數據,如CT、MRI等。通過面繪制法對數據進行三維重建并進行編輯,然后輸出為通用CAD、FEA(Finite Element Analysis,F(xiàn)EA),快速成型格式。MIMICS包括圖像導入模塊、圖像分割模塊、圖像可視化模塊、圖像配準模塊、圖像測量模塊等[17]。

    Simpleware:Simpleware是由英國Simpleware Ltd公司開發(fā)的一套實現(xiàn)三維圖像到CAD轉化、快速成型和有限元建模的集成化軟件。目前Simpleware軟件已廣泛應用于逆向工程、材料工程、生物力學工程、有限元分析等多工業(yè)、多學科領域。Simpleware軟件包括ScanIP、ScanFE、ScanCAD三大部分。ScanIP為圖像處理軟件、ScanFE網格生成模塊、ScanCAD為CAD整合模塊。

    3.2? 逆向工程軟件

    Geomagic:Geomagic是結合了三維點云、三角網格編輯功能以及CAD造型設計功能的三維逆向工程軟件。該軟件的主要特點是支持多種醫(yī)學成像設備的文件格式的讀取和轉換。在醫(yī)學建模與仿真中,用戶需要根據人體的掃描數據建立感興趣區(qū)的模型[18-20],因此,Geomagic逆向工程軟件在醫(yī)學仿真中應用十分廣泛。

    Rapidform:Rapidform軟件在生物醫(yī)學領域,Rapidform的優(yōu)勢在于其較強大的三角面片處理功能。該軟件可通過處理掃描點云數據生成NURBS曲面,便于在CAD軟件進行設計和改進、在FEA軟件中進行仿真和分析。

    3.3? 正向工程軟件

    matic:3-matic是比利時Materilise公司開發(fā)的一款數字化CAD軟件。該軟件支持網格劃分、逆向工程?;赟TL格式的文件處理模式,3-matic對人體組織模型編輯功能更為強大,效率極高,極大地保證了模型結構的真實性。3-matic被稱為醫(yī)學CAD軟件。

    3.4? 醫(yī)學有限元常用CAD軟件

    SolidWorks:SolidWorks是一個專業(yè)的三維設計CAD軟件,其最大的優(yōu)點是它的友好界面為用戶提供易懂易用的建模方式,為數值仿真提供了一個良好的數據接口,更加適用于規(guī)則種植體的設計[21]。

    3.5? 醫(yī)學有限元分析常用軟件

    Hypermesh:醫(yī)學有限元分析中Hypermesh經常用于網格劃分、裝配等前處理中。該軟件功能強大,支持眾多的有限元前處理器。在骨肌系統(tǒng)建模與仿真中,應用較多。

    AnyBody軟件:AnyBody是一款用來模擬人體在環(huán)境中工作時內部骨肌系統(tǒng)生物力學響應的一款軟件。通過AnyBody軟件仿真可獲得人體的肌肉力、關節(jié)力、力矩、新陳代謝等豐富的生物力學數據。

    ANSYS Workbench軟件:美國ANSYS公司開發(fā)的工程仿真技術集成平臺,可進行線性靜態(tài)結構分析、模態(tài)分析、諧響應分析、響應譜分析、隨機振動分析、瞬態(tài)動力學分析、顯式動力學分析、結構非線性分析、流體動力學分析等。

    ABAQUS軟件:ABAQUS被廣泛地認為是功能最強的有限元軟件,可以分析復雜的固體力學結構力學系統(tǒng),特別是能夠駕馭非常龐大復雜的問題和模擬高度非線性問題。ABAQUS不但可以做單一零件的力學和多物理場的分析,同時還可以做系統(tǒng)級的分析和研究。

    COMSOL Multiphysics軟件:COMSOL Multiphysics是一款大型的高級數值仿真軟件。廣泛應用于各個領域的科學研究以及工程計算,模擬科學和工程領域的各種物理過程。COMSOL當前有一個基本模塊和八個專業(yè)模塊:結構力學模塊(Structural MechanicsModule)、以及反應工程實驗室(COMSOL Reaction Engineering LAB)、信號與系統(tǒng)實驗室(Signal&System LAB)、最優(yōu)化實驗室(Optimization LAB)、CAD導入模塊(CADImportModule)、二次開發(fā)模塊(COMSOL ScriptTM)。

    ADAMS軟件:ADMAS,即機械系統(tǒng)動力學自動分析軟件,是美國MDI公司開發(fā)的虛擬樣機分析軟件。ADAMS軟件中提供的多種多樣的關節(jié)和力的輸入方式,完全可以勝任頭頸部模型中的集總關節(jié)、韌帶和肌肉的建模與分析。

    4? 結束語

    建立醫(yī)學有限元實驗有兩個關鍵的問題:(1)醫(yī)用有限元模型快速準確的建立。模型的快速準確建立可以減少仿真實驗所需時間、降低費用、增加仿真的準確性和可信性。(2)建立通用的有限元模型庫,為進一步的實驗教學和科研打下堅實的基礎。因此需要在具體實驗實踐中逐步探索和積累。將工程有限元分析同醫(yī)學結合開設實驗課,屬于多學科之間的交叉領域,不僅可以提高學生對所學專業(yè)知識的綜合運用能力,增強學生就業(yè)與學習深造的競爭力,而且可以加強多學科教師的教學和科研合作,提高教師的教學科研水平。同時提高相關實驗室的利用率,為學生自主開展創(chuàng)新實驗提供平臺,加強學校和附屬醫(yī)院的教學科研合作,為醫(yī)學院校提供更為廣闊的教學和科學研究領域。

    參考文獻

    [1]Oftadeh R, Perez-Viloria M, Villa-Camacho J C, et al. Biomechanics and mechanobiology of trabecular bone: a review[J]. Journal of Biomechanical Engineering, 2015, 137(1): 010802.

    [2]Zhao X, Liu Y, Ding J L, et al. Hemodynamic effects of the anastomoses in the modified Blalock -Tassing shunt: a numerical study using A 0D/3D coupling method[J]. Journal of Mechanics in Medicine & Biology, 2015, 15(01): 2076-2086.

    [3]Seo, J., et al. A Highly Automated Computational Method for Modeling of Intracranial Aneurysm Hemodynamics. Frontiers in Physiology, 2018(9): 681.

    [4]Bonneville F, Sourour N, Biondi A. Intracranial Aneurysms: an Overview. Neuroimaging Clinics of North America. 2006; 16(3): 371-382.

    [5]Keedy A. An overview of intracranial aneurysms. McGill Journal of Medicine. 2006; 9(2): 141-146.

    [6]Wermer MJH, Van Der Schaaf IC, Algra A, Rinkel GJE. Risk of rupture of unruptured intracranial aneurysms in relation to patient and aneurysm characteristics: An updated meta- analysis. Stroke. 2007; 38 (4): 1404-1410.

    [7]Lasheras JC. The biomechanics of arterial aneurysms. Annual Review of Fluid Mechanics. 2007; 39: 293-319.

    [8]Sforza DM, Putman CM, Cebral JR. Hemodynamics of cerebral aneurysms. Annual Review of Fluid Mechanics. 2009; 41: 91-107.

    [9]Miura Y, Ishida F, Umeda Y, Tanemura H, Suzuki H, Matsushima S, et al. Low wall shear stress is independently associated with the rupture status of middle cerebral artery aneurysms. Stroke. 2013; 44(2): 519-521.

    [10]Fukazawa K, Ishida F, Umeda Y, et al. Using computational fluid dynamics analysis to characterize local hemodynamic features of middle cerebral artery aneurysm rupture points. World Neurosurg. 2015; 83(1): 80-86.

    [11]Tsuji M, Ishikawa T, Ishida F, et al. Stagnation and complex flow in ruptured cerebral aneurysms: a possible association with hemostatic pattern. J Neurosurg. 2016; 3: 1-7.

    [12]Hua, Y., J. H. Oh and Y. B. Kim, Influence of Parent Artery Segmentation and Boundary Conditions on Hemodynamic Characteristics of Intracranial Aneurysms. Yonsei Medical Journal, 2015. 56(5): p. 1328.

    [13]Paliwal, N., et al., Association between hemodynamic modifications and clinical outcome of intracranial aneurysms treated using flow diverters. Proc SPIE Int Soc Opt Eng, 2017. 10135.

    [14]Lu G, Huang L, Zhang XL, Wang SZ, et al. Influence of hemodynamic factors on rupture of intracranial aneurysms: patient-specific 3D mirror aneurysms model computational fluid dynamics simulation. AJNR American journal of neuroradiology. 2011; 32(7): 1255-61.

    [15]Xiang J, Natarajan SK, Tremmel M, et al. Hemodynamic- morphologic discriminants for intracranial aneurysm rupture. Stroke; a journal of cerebral circulation. 2011; 42(1): 144- 52.

    [16]Liu J, Fan J, Xiang J, et al. Hemodynamic characteristics of large unruptured internal carotid artery aneurysms prior to rupture: a case control study. J Neurointerv Surg. 2016; 8: 367-72.

    [17]Pentimalli L, Modesti A, Vignati A, et al. Role of apoptosis in intracranial aneurysm rupture. J Neurosurg. 2004; 101: 1018-25.

    [18]Jou LD, Wong G, Dispensa B, et al. Correlation between luminal geometry changes and hemodynamics in fusiform intracranial aneurysms. AJNR Am J Neuroradiol 2005; 26: 2357-63.

    [19]Acevedo-Bolton G, Jou LD, Dispensa BP, et al. Estimating the hemodynamic impact of interventional treatments of aneurysms: numerical simulation with experimental validation: technical case report. Neurosurgery 2006; 59: E429-30.

    [20]Meng H, Tutino VM, Xiang J, et al. High WSS or low WSS? Complex interactions of hemodynamics with intracranial aneurysm initiation, growth, and rupture: toward a unifying hypothesis. AJNR. American journal of neuroradiology. 2014; 35: 1254-1262. [2]

    [21]Wang C., et al., Hemodynamic alterations after stent implantation in 15 cases of intracranial aneurysm. Acta Neurochirurgica, 2016. 158(4): p. 811-819.

    猜你喜歡
    生物力學有限元
    自由式輪滑鞋的研究
    不同樁核材料對下頜第一磨牙殘根及樁核冠的生物力學影響
    可注射性磷酸鈣骨水泥對骨質疏松椎體壓縮骨折的生物力學研究
    真空攪拌骨水泥技術應用于椎體后凸成形術的生物力學研究
    磨削淬硬殘余應力的有限元分析
    基于SolidWorks的吸嘴支撐臂有限元分析
    箱形孔軋制的有限元模擬
    上海金屬(2013年4期)2013-12-20 07:57:18
    巨型總段吊裝中的有限元方法應用
    船海工程(2013年6期)2013-03-11 18:57:27
    欧美日韩av久久| 久久99热这里只频精品6学生| 2018国产大陆天天弄谢| 亚洲欧洲精品一区二区精品久久久| 夫妻午夜视频| 高潮久久久久久久久久久不卡| 中文字幕高清在线视频| 亚洲成a人片在线一区二区| 亚洲精品国产色婷婷电影| www.熟女人妻精品国产| 一级毛片精品| 成年人午夜在线观看视频| 黄色成人免费大全| 狠狠精品人妻久久久久久综合| 精品少妇内射三级| 高清欧美精品videossex| 黄色成人免费大全| av线在线观看网站| 十八禁高潮呻吟视频| 亚洲男人天堂网一区| 亚洲熟妇熟女久久| 精品福利永久在线观看| 香蕉国产在线看| a级毛片在线看网站| 成年人黄色毛片网站| 国产欧美日韩综合在线一区二区| 91麻豆av在线| 亚洲第一欧美日韩一区二区三区 | 久久 成人 亚洲| 亚洲中文日韩欧美视频| 国产老妇伦熟女老妇高清| 好男人电影高清在线观看| 亚洲成人免费av在线播放| 男女无遮挡免费网站观看| 欧美变态另类bdsm刘玥| 99久久99久久久精品蜜桃| 三上悠亚av全集在线观看| 丝袜美腿诱惑在线| 欧美性长视频在线观看| 十八禁网站免费在线| 两人在一起打扑克的视频| 成年女人毛片免费观看观看9 | 成人精品一区二区免费| 菩萨蛮人人尽说江南好唐韦庄| 精品久久久久久久毛片微露脸| videos熟女内射| 99精品在免费线老司机午夜| 桃花免费在线播放| 少妇裸体淫交视频免费看高清 | 午夜福利视频精品| 国产日韩一区二区三区精品不卡| 窝窝影院91人妻| 国产成人av激情在线播放| 两个人看的免费小视频| 最新在线观看一区二区三区| 美女扒开内裤让男人捅视频| 色尼玛亚洲综合影院| 91精品国产国语对白视频| 一区二区三区乱码不卡18| 午夜日韩欧美国产| 精品高清国产在线一区| 在线十欧美十亚洲十日本专区| 亚洲午夜精品一区,二区,三区| 老司机亚洲免费影院| 国产无遮挡羞羞视频在线观看| 999久久久精品免费观看国产| 黑人巨大精品欧美一区二区mp4| 侵犯人妻中文字幕一二三四区| 一区二区三区激情视频| 悠悠久久av| 国产日韩欧美亚洲二区| 久久国产精品男人的天堂亚洲| 久久国产精品男人的天堂亚洲| 考比视频在线观看| 亚洲av美国av| 新久久久久国产一级毛片| 免费日韩欧美在线观看| 王馨瑶露胸无遮挡在线观看| 免费一级毛片在线播放高清视频 | 精品午夜福利视频在线观看一区 | 一进一出好大好爽视频| 国产不卡一卡二| 蜜桃国产av成人99| 岛国在线观看网站| 久久青草综合色| 婷婷丁香在线五月| 久久久精品94久久精品| 精品一区二区三区av网在线观看 | 亚洲中文日韩欧美视频| 少妇猛男粗大的猛烈进出视频| 中文字幕av电影在线播放| 国精品久久久久久国模美| 国产精品电影一区二区三区 | 国产av国产精品国产| 成人精品一区二区免费| 一本综合久久免费| 12—13女人毛片做爰片一| 男人舔女人的私密视频| 久久精品aⅴ一区二区三区四区| 国产精品久久久人人做人人爽| 国产精品自产拍在线观看55亚洲 | 午夜精品国产一区二区电影| 亚洲国产看品久久| 亚洲av美国av| 亚洲性夜色夜夜综合| 午夜91福利影院| 最新美女视频免费是黄的| 窝窝影院91人妻| 欧美人与性动交α欧美软件| 麻豆成人av在线观看| 老司机在亚洲福利影院| 超碰97精品在线观看| 老熟妇仑乱视频hdxx| 女性被躁到高潮视频| 飞空精品影院首页| 国产亚洲精品一区二区www | 超碰97精品在线观看| 精品一区二区三区av网在线观看 | 亚洲熟女精品中文字幕| 大香蕉久久成人网| 亚洲一卡2卡3卡4卡5卡精品中文| 精品国产一区二区久久| 少妇的丰满在线观看| 最新在线观看一区二区三区| 香蕉丝袜av| 久久久久久久大尺度免费视频| 国产激情久久老熟女| 极品人妻少妇av视频| 久久精品亚洲av国产电影网| 久久香蕉激情| 成年女人毛片免费观看观看9 | 大香蕉久久网| videosex国产| 人妻一区二区av| 久久精品熟女亚洲av麻豆精品| 国产精品国产高清国产av | tocl精华| 欧美日韩视频精品一区| 国产免费av片在线观看野外av| 建设人人有责人人尽责人人享有的| 国产精品一区二区精品视频观看| xxxhd国产人妻xxx| 青青草视频在线视频观看| 国产亚洲精品一区二区www | 1024视频免费在线观看| 国产三级黄色录像| 国产免费福利视频在线观看| 国产伦人伦偷精品视频| 亚洲性夜色夜夜综合| 亚洲国产精品一区二区三区在线| 国产精品麻豆人妻色哟哟久久| 91字幕亚洲| av天堂在线播放| av欧美777| a级片在线免费高清观看视频| 国产成人系列免费观看| 少妇粗大呻吟视频| 天天躁狠狠躁夜夜躁狠狠躁| av在线播放免费不卡| 国产精品国产av在线观看| kizo精华| 一边摸一边做爽爽视频免费| 免费av中文字幕在线| 高清在线国产一区| 777米奇影视久久| 精品欧美一区二区三区在线| 午夜免费鲁丝| 美女国产高潮福利片在线看| 国产亚洲一区二区精品| 下体分泌物呈黄色| 亚洲一区中文字幕在线| 久久久国产一区二区| 亚洲精品乱久久久久久| 欧美中文综合在线视频| 精品熟女少妇八av免费久了| 国产熟女午夜一区二区三区| 精品福利永久在线观看| 狠狠婷婷综合久久久久久88av| 狠狠精品人妻久久久久久综合| 最近最新中文字幕大全电影3 | 久久香蕉激情| 三级毛片av免费| 两个人看的免费小视频| av国产精品久久久久影院| 曰老女人黄片| 国产欧美日韩一区二区精品| 高清毛片免费观看视频网站 | 日韩欧美国产一区二区入口| 好男人电影高清在线观看| 大型黄色视频在线免费观看| 国产深夜福利视频在线观看| 老熟女久久久| 国产一区二区在线观看av| 国产伦理片在线播放av一区| 热99国产精品久久久久久7| 亚洲av成人一区二区三| 国产日韩一区二区三区精品不卡| 久久久国产成人免费| 99久久国产精品久久久| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久国产一级毛片高清牌| 日本av免费视频播放| 视频区图区小说| 欧美国产精品一级二级三级| 久久国产亚洲av麻豆专区| 中文字幕精品免费在线观看视频| 高清欧美精品videossex| 亚洲精品一卡2卡三卡4卡5卡| 欧美 亚洲 国产 日韩一| 91成人精品电影| 中文字幕人妻丝袜制服| 国产单亲对白刺激| www日本在线高清视频| 国产成人系列免费观看| 母亲3免费完整高清在线观看| 亚洲精品乱久久久久久| 一二三四在线观看免费中文在| tube8黄色片| 精品国产乱码久久久久久男人| 亚洲伊人久久精品综合| 伦理电影免费视频| 蜜桃国产av成人99| 一二三四社区在线视频社区8| 欧美成狂野欧美在线观看| 日韩免费高清中文字幕av| 日韩成人在线观看一区二区三区| 露出奶头的视频| 99久久国产精品久久久| 色综合婷婷激情| 久久精品国产a三级三级三级| 午夜福利一区二区在线看| 久久天躁狠狠躁夜夜2o2o| 制服人妻中文乱码| 中文字幕制服av| 一个人免费在线观看的高清视频| 欧美人与性动交α欧美精品济南到| 国产亚洲午夜精品一区二区久久| 国产精品国产高清国产av | 亚洲欧美日韩另类电影网站| 亚洲第一青青草原| 麻豆国产av国片精品| 一区二区三区国产精品乱码| 久久狼人影院| 丝袜人妻中文字幕| 美女扒开内裤让男人捅视频| 欧美日韩av久久| 国产欧美日韩一区二区三| 欧美乱妇无乱码| 一本—道久久a久久精品蜜桃钙片| 久久午夜综合久久蜜桃| 国产老妇伦熟女老妇高清| videos熟女内射| 久久久久国内视频| 黄色毛片三级朝国网站| 日本欧美视频一区| 精品福利永久在线观看| 丰满人妻熟妇乱又伦精品不卡| 女性被躁到高潮视频| avwww免费| 一边摸一边抽搐一进一小说 | 国产福利在线免费观看视频| 久久久精品94久久精品| 午夜日韩欧美国产| 精品一品国产午夜福利视频| 亚洲人成电影观看| 男男h啪啪无遮挡| 欧美精品高潮呻吟av久久| 一级片免费观看大全| 精品国产乱码久久久久久小说| 久久青草综合色| 91字幕亚洲| 天堂中文最新版在线下载| 亚洲第一欧美日韩一区二区三区 | 午夜视频精品福利| 老司机影院毛片| 成人三级做爰电影| 热re99久久精品国产66热6| 黄色视频不卡| 亚洲色图 男人天堂 中文字幕| √禁漫天堂资源中文www| 日本黄色视频三级网站网址 | 久久亚洲精品不卡| 亚洲伊人久久精品综合| 搡老乐熟女国产| 午夜福利视频精品| 三上悠亚av全集在线观看| 免费在线观看完整版高清| 丝袜美腿诱惑在线| 69av精品久久久久久 | 91老司机精品| www日本在线高清视频| 国产深夜福利视频在线观看| kizo精华| 交换朋友夫妻互换小说| 成人国产av品久久久| 岛国在线观看网站| videos熟女内射| 大香蕉久久网| 久久精品91无色码中文字幕| 啪啪无遮挡十八禁网站| 久热爱精品视频在线9| 妹子高潮喷水视频| 女人被躁到高潮嗷嗷叫费观| 欧美人与性动交α欧美软件| 在线观看www视频免费| 女同久久另类99精品国产91| 中文字幕制服av| 亚洲男人天堂网一区| 国产精品久久久久久人妻精品电影 | 免费在线观看视频国产中文字幕亚洲| 一二三四在线观看免费中文在| 午夜视频精品福利| 久久av网站| 夜夜爽天天搞| 日韩中文字幕视频在线看片| 成人精品一区二区免费| 久久人人97超碰香蕉20202| 欧美老熟妇乱子伦牲交| 免费人妻精品一区二区三区视频| 国产一区二区在线观看av| 精品福利观看| 久久久久精品人妻al黑| 国产亚洲欧美精品永久| 飞空精品影院首页| 男女下面插进去视频免费观看| 色尼玛亚洲综合影院| 热re99久久精品国产66热6| 黄色视频,在线免费观看| 法律面前人人平等表现在哪些方面| 中文字幕av电影在线播放| 成年动漫av网址| 国产欧美日韩一区二区精品| 精品一区二区三区av网在线观看 | 久久99热这里只频精品6学生| 真人做人爱边吃奶动态| 又黄又粗又硬又大视频| 999久久久精品免费观看国产| netflix在线观看网站| 热re99久久国产66热| 午夜福利视频在线观看免费| 欧美日韩亚洲综合一区二区三区_| 国产黄色免费在线视频| 男女高潮啪啪啪动态图| 欧美激情 高清一区二区三区| 在线观看人妻少妇| 欧美 日韩 精品 国产| 婷婷成人精品国产| 久久久久久久大尺度免费视频| 欧美成人午夜精品| 日本五十路高清| 日韩欧美三级三区| 精品免费久久久久久久清纯 | 国产97色在线日韩免费| 一级毛片女人18水好多| 不卡一级毛片| 黄片小视频在线播放| 国产伦人伦偷精品视频| 大陆偷拍与自拍| 久久国产精品大桥未久av| 国产精品一区二区在线观看99| 天堂中文最新版在线下载| 日韩大码丰满熟妇| 色综合欧美亚洲国产小说| 欧美人与性动交α欧美软件| 老熟妇仑乱视频hdxx| 精品一品国产午夜福利视频| 大型av网站在线播放| 午夜激情久久久久久久| 亚洲午夜精品一区,二区,三区| 亚洲专区中文字幕在线| 十八禁网站网址无遮挡| 一区在线观看完整版| 国产精品久久电影中文字幕 | 色婷婷久久久亚洲欧美| 午夜福利欧美成人| 国产精品九九99| 午夜老司机福利片| 精品久久久精品久久久| 亚洲国产毛片av蜜桃av| 亚洲av第一区精品v没综合| 一本色道久久久久久精品综合| 50天的宝宝边吃奶边哭怎么回事| 五月天丁香电影| 亚洲色图av天堂| 国产不卡一卡二| 亚洲伊人久久精品综合| 在线观看人妻少妇| 自拍欧美九色日韩亚洲蝌蚪91| 伊人久久大香线蕉亚洲五| tocl精华| 在线亚洲精品国产二区图片欧美| 欧美日韩视频精品一区| 少妇粗大呻吟视频| 18禁黄网站禁片午夜丰满| 午夜精品国产一区二区电影| 老司机午夜十八禁免费视频| 亚洲欧美色中文字幕在线| 欧美黄色淫秽网站| 午夜福利欧美成人| 国产高清videossex| 国产xxxxx性猛交| 色综合欧美亚洲国产小说| 国产不卡一卡二| av网站在线播放免费| 少妇精品久久久久久久| 高清av免费在线| 午夜福利乱码中文字幕| 丰满人妻熟妇乱又伦精品不卡| 国产又色又爽无遮挡免费看| 免费av中文字幕在线| 日韩欧美一区视频在线观看| 欧美精品av麻豆av| 亚洲精品一二三| 久久国产精品人妻蜜桃| 国产一区二区三区综合在线观看| 久久人妻av系列| 国产欧美日韩精品亚洲av| 亚洲 欧美一区二区三区| 欧美激情极品国产一区二区三区| 91老司机精品| 国产亚洲欧美在线一区二区| 香蕉久久夜色| 纵有疾风起免费观看全集完整版| 午夜福利影视在线免费观看| 51午夜福利影视在线观看| 成年版毛片免费区| 高清毛片免费观看视频网站 | 欧美精品av麻豆av| 国产在线视频一区二区| 美女福利国产在线| 中文字幕另类日韩欧美亚洲嫩草| 又黄又粗又硬又大视频| 亚洲国产欧美一区二区综合| 日韩中文字幕欧美一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 三上悠亚av全集在线观看| 69精品国产乱码久久久| 人人妻,人人澡人人爽秒播| 亚洲成a人片在线一区二区| 免费黄频网站在线观看国产| 国产伦人伦偷精品视频| 中文字幕色久视频| 在线观看www视频免费| 少妇裸体淫交视频免费看高清 | www.熟女人妻精品国产| 黑人巨大精品欧美一区二区蜜桃| 国产有黄有色有爽视频| 国产高清videossex| 欧美+亚洲+日韩+国产| 亚洲黑人精品在线| 黄色怎么调成土黄色| 午夜视频精品福利| 亚洲欧美一区二区三区黑人| 多毛熟女@视频| 一区二区三区激情视频| 人人妻人人澡人人爽人人夜夜| 久久亚洲真实| 窝窝影院91人妻| 久久精品aⅴ一区二区三区四区| 国产成人影院久久av| 天天躁日日躁夜夜躁夜夜| 中文字幕人妻丝袜一区二区| 国产97色在线日韩免费| 超碰97精品在线观看| 18禁美女被吸乳视频| 日本av手机在线免费观看| 精品一区二区三卡| 欧美乱妇无乱码| 欧美激情久久久久久爽电影 | 91av网站免费观看| 一级毛片女人18水好多| 国产高清激情床上av| 欧美精品一区二区免费开放| 老司机靠b影院| 国产一区二区在线观看av| 国产精品久久久av美女十八| 国产高清videossex| 啦啦啦 在线观看视频| 亚洲专区字幕在线| 亚洲精品自拍成人| 老熟妇仑乱视频hdxx| 久久精品国产99精品国产亚洲性色 | 国产一区二区三区视频了| 国产精品美女特级片免费视频播放器 | 亚洲av日韩在线播放| 美女国产高潮福利片在线看| 中文欧美无线码| 十八禁高潮呻吟视频| 2018国产大陆天天弄谢| 亚洲成av片中文字幕在线观看| 国产99久久九九免费精品| 日韩欧美免费精品| 国产精品一区二区精品视频观看| a在线观看视频网站| 午夜精品国产一区二区电影| 亚洲综合色网址| 深夜精品福利| 在线观看舔阴道视频| 人成视频在线观看免费观看| 欧美日韩黄片免| 亚洲免费av在线视频| 久热这里只有精品99| 一级毛片女人18水好多| 亚洲欧美一区二区三区黑人| 超碰97精品在线观看| 十分钟在线观看高清视频www| 精品国产一区二区三区四区第35| 国产精品久久电影中文字幕 | videosex国产| 午夜激情av网站| 免费观看a级毛片全部| 丁香六月天网| 国产精品一区二区在线不卡| 视频区欧美日本亚洲| 欧美日韩av久久| 亚洲人成电影免费在线| 男人舔女人的私密视频| 少妇的丰满在线观看| 91麻豆精品激情在线观看国产 | 在线播放国产精品三级| 成人亚洲精品一区在线观看| 手机成人av网站| 麻豆国产av国片精品| 国产日韩欧美在线精品| 青青草视频在线视频观看| 香蕉丝袜av| 一边摸一边做爽爽视频免费| 精品国产一区二区久久| av视频免费观看在线观看| 黄色毛片三级朝国网站| 夜夜爽天天搞| 日韩免费av在线播放| 午夜久久久在线观看| 亚洲熟妇熟女久久| 下体分泌物呈黄色| 欧美日韩一级在线毛片| 国产主播在线观看一区二区| 久久 成人 亚洲| cao死你这个sao货| 男人舔女人的私密视频| 色尼玛亚洲综合影院| 国产精品一区二区在线观看99| 国产亚洲精品一区二区www | 欧美另类亚洲清纯唯美| 亚洲中文av在线| 97人妻天天添夜夜摸| av福利片在线| 久久久精品免费免费高清| 国产免费视频播放在线视频| 欧美在线一区亚洲| 三上悠亚av全集在线观看| 午夜久久久在线观看| 视频区欧美日本亚洲| 一级黄色大片毛片| 精品一区二区三区四区五区乱码| 久久精品aⅴ一区二区三区四区| 日本av手机在线免费观看| 满18在线观看网站| 亚洲专区中文字幕在线| 黄色丝袜av网址大全| 午夜福利在线免费观看网站| 侵犯人妻中文字幕一二三四区| 午夜精品久久久久久毛片777| 精品人妻1区二区| 搡老乐熟女国产| 久久国产精品大桥未久av| 一区二区三区激情视频| 久久久久国产一级毛片高清牌| 黄片小视频在线播放| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜福利欧美成人| 高清在线国产一区| 高清视频免费观看一区二区| 中文字幕人妻熟女乱码| 亚洲精品av麻豆狂野| av有码第一页| 黑人操中国人逼视频| 久久中文字幕人妻熟女| 啦啦啦中文免费视频观看日本| 另类亚洲欧美激情| 久久久久久久久免费视频了| 国产欧美亚洲国产| 国产成人精品久久二区二区免费| 国产精品.久久久| 精品国产国语对白av| 满18在线观看网站| 在线观看免费视频日本深夜| 桃红色精品国产亚洲av| 免费在线观看影片大全网站| 99久久精品国产亚洲精品| 99re在线观看精品视频| 搡老乐熟女国产| 亚洲精品美女久久av网站| 99国产精品99久久久久| 亚洲精品中文字幕一二三四区 | 精品少妇黑人巨大在线播放| 真人做人爱边吃奶动态| 国产免费福利视频在线观看| 变态另类成人亚洲欧美熟女 | av一本久久久久| 欧美老熟妇乱子伦牲交| 国产av国产精品国产| 91av网站免费观看| 亚洲专区国产一区二区| 亚洲伊人久久精品综合| 久久av网站| 在线 av 中文字幕| 亚洲伊人久久精品综合| 日本黄色日本黄色录像| 777久久人妻少妇嫩草av网站| 69精品国产乱码久久久| kizo精华| 热99国产精品久久久久久7| 如日韩欧美国产精品一区二区三区| 丁香六月天网|