• <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
    久久亚洲真实| 女警被强在线播放| 亚洲精品粉嫩美女一区| 亚洲 欧美一区二区三区| 亚洲熟妇中文字幕五十中出| 熟女少妇亚洲综合色aaa.| 亚洲九九香蕉| 97碰自拍视频| 嫩草影视91久久| 亚洲中文字幕日韩| 欧美国产精品va在线观看不卡| 黄色丝袜av网址大全| 操出白浆在线播放| 嫁个100分男人电影在线观看| 变态另类丝袜制服| 亚洲第一青青草原| 国产高清视频在线播放一区| 真人做人爱边吃奶动态| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美国产日韩亚洲一区| 亚洲一区高清亚洲精品| 午夜福利欧美成人| 99久久久亚洲精品蜜臀av| 无限看片的www在线观看| 在线观看免费日韩欧美大片| 男女床上黄色一级片免费看| 黄片小视频在线播放| 女人被躁到高潮嗷嗷叫费观| 我的亚洲天堂| 国产精品亚洲av一区麻豆| 国产熟女午夜一区二区三区| 19禁男女啪啪无遮挡网站| 少妇熟女aⅴ在线视频| 精品高清国产在线一区| 国产成人啪精品午夜网站| 久久久久国内视频| 久久午夜综合久久蜜桃| 丝袜人妻中文字幕| 精品久久久久久久毛片微露脸| 精品久久久久久,| 大型黄色视频在线免费观看| 亚洲午夜理论影院| 19禁男女啪啪无遮挡网站| 亚洲精品国产精品久久久不卡| 国产欧美日韩一区二区三| 国产精品一区二区三区四区久久 | 老汉色∧v一级毛片| 国产午夜福利久久久久久| 国产精品亚洲一级av第二区| 亚洲一区二区三区色噜噜| 日韩精品中文字幕看吧| 一级黄色大片毛片| 精品午夜福利视频在线观看一区| 国产精品一区二区精品视频观看| 亚洲欧美精品综合一区二区三区| 午夜成年电影在线免费观看| 久久精品国产99精品国产亚洲性色 | 99在线人妻在线中文字幕| 国产精品99久久99久久久不卡| 高清在线国产一区| 99香蕉大伊视频| 这个男人来自地球电影免费观看| 欧美激情高清一区二区三区| 欧美色欧美亚洲另类二区 | 美女 人体艺术 gogo| av超薄肉色丝袜交足视频| 黑人操中国人逼视频| 50天的宝宝边吃奶边哭怎么回事| 国产精品一区二区免费欧美| 国产成人一区二区三区免费视频网站| 中文字幕另类日韩欧美亚洲嫩草| 亚洲狠狠婷婷综合久久图片| 露出奶头的视频| 天天躁狠狠躁夜夜躁狠狠躁| 国产亚洲精品久久久久久毛片| 久久精品亚洲熟妇少妇任你| 久久久久精品国产欧美久久久| 熟女少妇亚洲综合色aaa.| www国产在线视频色| 国产麻豆69| 琪琪午夜伦伦电影理论片6080| 亚洲人成77777在线视频| 女人被躁到高潮嗷嗷叫费观| 中文字幕人成人乱码亚洲影| 国产91精品成人一区二区三区| 美女高潮到喷水免费观看| 亚洲人成77777在线视频| 视频在线观看一区二区三区| 99久久国产精品久久久| 欧美一级a爱片免费观看看 | 欧美乱妇无乱码| 欧美色视频一区免费| 国产成人精品无人区| 免费人成视频x8x8入口观看| 俄罗斯特黄特色一大片| 亚洲aⅴ乱码一区二区在线播放 | 久久久久久免费高清国产稀缺| 精品免费久久久久久久清纯| 午夜免费激情av| 午夜精品在线福利| 成年女人毛片免费观看观看9| 最新美女视频免费是黄的| 三级毛片av免费| 国产区一区二久久| 亚洲第一欧美日韩一区二区三区| 日韩一卡2卡3卡4卡2021年| 亚洲欧美精品综合久久99| 久久久精品国产亚洲av高清涩受| 麻豆国产av国片精品| 亚洲在线自拍视频| 两人在一起打扑克的视频| 制服人妻中文乱码| 成人国语在线视频| 欧美色视频一区免费| 久久久久久亚洲精品国产蜜桃av| 91麻豆av在线| 精品一区二区三区av网在线观看| 91九色精品人成在线观看| 精品一区二区三区视频在线观看免费| 免费少妇av软件| ponron亚洲| 麻豆久久精品国产亚洲av| 黄网站色视频无遮挡免费观看| www.自偷自拍.com| 一二三四社区在线视频社区8| 女生性感内裤真人,穿戴方法视频| 1024视频免费在线观看| 久久精品国产99精品国产亚洲性色 | 久久这里只有精品19| 国产精品一区二区免费欧美| 757午夜福利合集在线观看| 午夜福利高清视频| 国产精品久久久久久精品电影 | 伊人久久大香线蕉亚洲五| 亚洲国产高清在线一区二区三 | 黄色片一级片一级黄色片| 亚洲精品国产一区二区精华液| 久久中文看片网| 欧美黑人欧美精品刺激| 黄色视频不卡| 日韩高清综合在线| 性欧美人与动物交配| 一级毛片女人18水好多| 黄色毛片三级朝国网站| 亚洲国产高清在线一区二区三 | 亚洲午夜理论影院| 手机成人av网站| 麻豆久久精品国产亚洲av| 一区福利在线观看| 啦啦啦 在线观看视频| 一本久久中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 香蕉久久夜色| 女人被躁到高潮嗷嗷叫费观| 日韩成人在线观看一区二区三区| 视频区欧美日本亚洲| 中亚洲国语对白在线视频| 极品人妻少妇av视频| 精品第一国产精品| 看片在线看免费视频| 亚洲一区中文字幕在线| 亚洲色图 男人天堂 中文字幕| 两个人视频免费观看高清| 亚洲国产精品成人综合色| 一级片免费观看大全| 真人一进一出gif抽搐免费| 亚洲av电影不卡..在线观看| 热re99久久国产66热| 变态另类成人亚洲欧美熟女 | 国产精品九九99| 满18在线观看网站| 中文字幕人妻熟女乱码| 午夜福利成人在线免费观看| 成熟少妇高潮喷水视频| 最近最新免费中文字幕在线| 亚洲精品在线观看二区| 91麻豆精品激情在线观看国产| 国产一区二区三区在线臀色熟女| 一区二区三区精品91| 老司机午夜福利在线观看视频| 国内久久婷婷六月综合欲色啪| 欧美色欧美亚洲另类二区 | 激情视频va一区二区三区| 国产成人系列免费观看| 久久精品人人爽人人爽视色| 一个人观看的视频www高清免费观看 | 一本久久中文字幕| 日本免费a在线| 久久亚洲真实| 亚洲成a人片在线一区二区| 丁香六月欧美| 午夜日韩欧美国产| 欧美精品亚洲一区二区| 久久国产乱子伦精品免费另类| 大陆偷拍与自拍| 日韩大码丰满熟妇| 欧美色欧美亚洲另类二区 | 国产精品免费一区二区三区在线| 婷婷丁香在线五月| 国产精品,欧美在线| 长腿黑丝高跟| 国产91精品成人一区二区三区| 国产免费男女视频| 久久久久久久久免费视频了| 在线十欧美十亚洲十日本专区| 亚洲欧美激情在线| 免费少妇av软件| 国产成人精品在线电影| 久久精品人人爽人人爽视色| 黄片大片在线免费观看| 欧美另类亚洲清纯唯美| 中文字幕精品免费在线观看视频| 久久欧美精品欧美久久欧美| 国产一区二区激情短视频| 久久精品国产亚洲av高清一级| 亚洲人成伊人成综合网2020| 成人精品一区二区免费| 99久久99久久久精品蜜桃| 久久国产精品男人的天堂亚洲| 欧美乱妇无乱码| 欧美绝顶高潮抽搐喷水| 岛国视频午夜一区免费看| aaaaa片日本免费| 亚洲专区国产一区二区| 成在线人永久免费视频| 免费观看精品视频网站| 亚洲激情在线av| 免费在线观看亚洲国产| 老汉色∧v一级毛片| 天天添夜夜摸| 亚洲 欧美一区二区三区| 亚洲国产毛片av蜜桃av| 别揉我奶头~嗯~啊~动态视频| 757午夜福利合集在线观看| 长腿黑丝高跟| 高清在线国产一区| 久久精品国产综合久久久| 丁香欧美五月| 久久这里只有精品19| 精品不卡国产一区二区三区| 精品国产乱码久久久久久男人| 国产日韩一区二区三区精品不卡| 一a级毛片在线观看| 一进一出抽搐动态| 午夜两性在线视频| 亚洲aⅴ乱码一区二区在线播放 | 夜夜爽天天搞| 免费在线观看视频国产中文字幕亚洲| 国产欧美日韩一区二区三| 日本撒尿小便嘘嘘汇集6| 国产精品亚洲一级av第二区| 男女下面插进去视频免费观看| 国产激情欧美一区二区| 亚洲,欧美精品.| 欧美中文日本在线观看视频| 国产精品影院久久| 午夜福利视频1000在线观看 | 国产欧美日韩一区二区三| 久久久久久久午夜电影| or卡值多少钱| 操出白浆在线播放| 亚洲 欧美一区二区三区| 最好的美女福利视频网| 亚洲七黄色美女视频| 给我免费播放毛片高清在线观看| 校园春色视频在线观看| 97人妻天天添夜夜摸| 88av欧美| 久久久国产成人免费| 正在播放国产对白刺激| 国产精品一区二区精品视频观看| 一区二区日韩欧美中文字幕| 亚洲aⅴ乱码一区二区在线播放 | 免费少妇av软件| 最近最新中文字幕大全免费视频| 精品久久蜜臀av无| 免费在线观看黄色视频的| 俄罗斯特黄特色一大片| 性欧美人与动物交配| netflix在线观看网站| 丁香欧美五月| 黑丝袜美女国产一区| 色播亚洲综合网| 大码成人一级视频| 在线播放国产精品三级| 亚洲人成77777在线视频| 给我免费播放毛片高清在线观看| 国产精品电影一区二区三区| 一二三四在线观看免费中文在| 91麻豆精品激情在线观看国产| a级毛片在线看网站| 99精品久久久久人妻精品| 中文字幕色久视频| 国产成人一区二区三区免费视频网站| 精品一区二区三区视频在线观看免费| 欧美不卡视频在线免费观看 | 久久精品国产综合久久久| 国内精品久久久久精免费| 精品日产1卡2卡| 国产精品久久视频播放| 久99久视频精品免费| 久久精品成人免费网站| 在线观看免费日韩欧美大片| 黄色片一级片一级黄色片| 久久香蕉激情| 可以免费在线观看a视频的电影网站| 欧美日韩一级在线毛片| 久久久久久亚洲精品国产蜜桃av| 日日摸夜夜添夜夜添小说| 熟妇人妻久久中文字幕3abv| 午夜免费成人在线视频| 久久香蕉国产精品| 好男人在线观看高清免费视频 | 亚洲成人免费电影在线观看| 伊人久久大香线蕉亚洲五| 国产免费av片在线观看野外av| 最近最新中文字幕大全免费视频| www.熟女人妻精品国产| 在线观看66精品国产| 国产精品,欧美在线| 又大又爽又粗| 午夜久久久久精精品| 99久久久亚洲精品蜜臀av| 涩涩av久久男人的天堂| 亚洲三区欧美一区| 亚洲激情在线av| 欧美乱妇无乱码| 中文字幕精品免费在线观看视频| 免费av毛片视频| 精品免费久久久久久久清纯| 精品午夜福利视频在线观看一区| 变态另类成人亚洲欧美熟女 | 成人手机av| 99国产精品99久久久久| 欧美日本亚洲视频在线播放| 老熟妇乱子伦视频在线观看| 亚洲狠狠婷婷综合久久图片| 久久久久久亚洲精品国产蜜桃av| 欧美绝顶高潮抽搐喷水| 国产一区在线观看成人免费| 性色av乱码一区二区三区2| 女人被狂操c到高潮| 亚洲自偷自拍图片 自拍| 亚洲第一av免费看| 午夜福利视频1000在线观看 | 午夜福利一区二区在线看| 午夜视频精品福利| 国产精品,欧美在线| 两性午夜刺激爽爽歪歪视频在线观看 | 精品第一国产精品| 国产亚洲精品第一综合不卡| 两性夫妻黄色片| 亚洲情色 制服丝袜| 黄网站色视频无遮挡免费观看| 露出奶头的视频| 欧美日韩亚洲国产一区二区在线观看| 欧美国产精品va在线观看不卡| 精品人妻在线不人妻| 国产精品九九99| 精品日产1卡2卡| 亚洲中文字幕一区二区三区有码在线看 | 精品国产超薄肉色丝袜足j| 久久精品国产清高在天天线| 亚洲色图 男人天堂 中文字幕| 99久久99久久久精品蜜桃| 日韩 欧美 亚洲 中文字幕| 人成视频在线观看免费观看| 国产精品秋霞免费鲁丝片| 宅男免费午夜| 高清毛片免费观看视频网站| 久久人妻av系列| 亚洲一区二区三区色噜噜| 啦啦啦免费观看视频1| 久久久久久久久久久久大奶| 窝窝影院91人妻| 久久狼人影院| 欧美人与性动交α欧美精品济南到| 少妇熟女aⅴ在线视频| 精品国内亚洲2022精品成人| 精品国产超薄肉色丝袜足j| 色综合站精品国产| 亚洲欧美精品综合久久99| 亚洲色图av天堂| 叶爱在线成人免费视频播放| 我的亚洲天堂| avwww免费| bbb黄色大片| 91老司机精品| 亚洲精品av麻豆狂野| 97碰自拍视频| 国产成人精品无人区| 久久热在线av| 巨乳人妻的诱惑在线观看| 高清黄色对白视频在线免费看| 色综合婷婷激情| 长腿黑丝高跟| 我的亚洲天堂| 国产私拍福利视频在线观看| 欧美成人免费av一区二区三区| 757午夜福利合集在线观看| 久久久久久国产a免费观看| 成人18禁高潮啪啪吃奶动态图| 韩国精品一区二区三区| 满18在线观看网站| 啦啦啦韩国在线观看视频| 午夜福利欧美成人| 一进一出抽搐动态| 又紧又爽又黄一区二区| 亚洲欧美日韩无卡精品| 波多野结衣一区麻豆| 欧美日韩精品网址| 成在线人永久免费视频| 人人澡人人妻人| 国产xxxxx性猛交| 三级毛片av免费| 黄网站色视频无遮挡免费观看| 国产精品爽爽va在线观看网站 | 精品电影一区二区在线| 一边摸一边做爽爽视频免费| av在线天堂中文字幕| 亚洲欧美日韩高清在线视频| 国产av又大| 一二三四社区在线视频社区8| 久久精品国产99精品国产亚洲性色 | 天堂√8在线中文| 国产成人av教育| 午夜激情av网站| 久久中文字幕人妻熟女| 亚洲欧美日韩高清在线视频| 精品国产超薄肉色丝袜足j| 亚洲av成人av| 色综合欧美亚洲国产小说| 激情视频va一区二区三区| 最近最新中文字幕大全免费视频| 国产激情欧美一区二区| av视频免费观看在线观看| 操出白浆在线播放| 女人爽到高潮嗷嗷叫在线视频| 美女免费视频网站| 国产免费av片在线观看野外av| 99久久国产精品久久久| 精品少妇一区二区三区视频日本电影| 亚洲视频免费观看视频| 国产三级黄色录像| 黄片大片在线免费观看| 麻豆一二三区av精品| 88av欧美| 国产日韩一区二区三区精品不卡| 婷婷六月久久综合丁香| 欧洲精品卡2卡3卡4卡5卡区| 亚洲国产精品sss在线观看| 可以在线观看的亚洲视频| 中文亚洲av片在线观看爽| 日本a在线网址| 免费在线观看日本一区| av视频免费观看在线观看| 亚洲第一av免费看| 色精品久久人妻99蜜桃| 女人爽到高潮嗷嗷叫在线视频| 免费在线观看完整版高清| 少妇的丰满在线观看| 国产高清视频在线播放一区| 18禁裸乳无遮挡免费网站照片 | 美女 人体艺术 gogo| 国产一级毛片七仙女欲春2 | 后天国语完整版免费观看| 久久久久国产精品人妻aⅴ院| 91精品三级在线观看| 亚洲男人的天堂狠狠| 日韩成人在线观看一区二区三区| 精品熟女少妇八av免费久了| 一边摸一边做爽爽视频免费| 免费人成视频x8x8入口观看| 欧美 亚洲 国产 日韩一| 欧美av亚洲av综合av国产av| 日本vs欧美在线观看视频| 韩国精品一区二区三区| 日韩有码中文字幕| 色哟哟哟哟哟哟| 久久天躁狠狠躁夜夜2o2o| 亚洲国产精品合色在线| 男女之事视频高清在线观看| 中文字幕高清在线视频| 欧美黄色片欧美黄色片| 国产精品一区二区在线不卡| 男女做爰动态图高潮gif福利片 | 午夜a级毛片| 九色国产91popny在线| videosex国产| 成年版毛片免费区| 国产欧美日韩一区二区精品| 久久久国产欧美日韩av| 又大又爽又粗| 国产精品影院久久| 日本在线视频免费播放| 亚洲中文av在线| 国产精品乱码一区二三区的特点 | 亚洲 国产 在线| 久久草成人影院| 1024香蕉在线观看| 精品卡一卡二卡四卡免费| cao死你这个sao货| 精品国产亚洲在线| 久久九九热精品免费| 女人被躁到高潮嗷嗷叫费观| 最近最新免费中文字幕在线| 两性夫妻黄色片| 一本综合久久免费| 亚洲av电影在线进入| 黑人欧美特级aaaaaa片| 巨乳人妻的诱惑在线观看| 在线观看日韩欧美| 九色亚洲精品在线播放| 亚洲黑人精品在线| 熟女少妇亚洲综合色aaa.| 亚洲成人精品中文字幕电影| 99久久精品国产亚洲精品| 91字幕亚洲| 一进一出抽搐动态| 久久精品国产99精品国产亚洲性色 | 无限看片的www在线观看| 叶爱在线成人免费视频播放| 黑人操中国人逼视频| 一个人免费在线观看的高清视频| 午夜日韩欧美国产| 久久久国产精品麻豆| 丝袜美足系列| 天天添夜夜摸| 欧美乱码精品一区二区三区| 俄罗斯特黄特色一大片| 中文字幕精品免费在线观看视频| 国产亚洲精品久久久久5区| 免费在线观看视频国产中文字幕亚洲| 亚洲七黄色美女视频| 欧美不卡视频在线免费观看 | 欧美+亚洲+日韩+国产| 久久这里只有精品19| 免费看a级黄色片| 国产97色在线日韩免费| 美女大奶头视频| 美女扒开内裤让男人捅视频| 亚洲精品久久成人aⅴ小说| 午夜久久久在线观看| 亚洲精品久久成人aⅴ小说| 日韩有码中文字幕| 免费在线观看黄色视频的| 精品电影一区二区在线| 老鸭窝网址在线观看| 曰老女人黄片| 日本 欧美在线| 天天一区二区日本电影三级 | 怎么达到女性高潮| 精品国产美女av久久久久小说| 久久久久精品国产欧美久久久| 亚洲精品久久成人aⅴ小说| 国产激情久久老熟女| 国产精品久久久久久精品电影 | 两个人看的免费小视频| 99在线人妻在线中文字幕| 啦啦啦观看免费观看视频高清 | 国产99久久九九免费精品| 国产黄a三级三级三级人| 国产成人精品久久二区二区免费| 成在线人永久免费视频| 老熟妇乱子伦视频在线观看| 亚洲av日韩精品久久久久久密| 成年人黄色毛片网站| 亚洲精品国产区一区二| 欧美人与性动交α欧美精品济南到| 国产欧美日韩一区二区三区在线| 操美女的视频在线观看| 日本欧美视频一区| 国内久久婷婷六月综合欲色啪| 天堂动漫精品| 亚洲中文日韩欧美视频| 嫩草影院精品99| 不卡av一区二区三区| 麻豆av在线久日| 黄色 视频免费看| 国产麻豆成人av免费视频| 男女午夜视频在线观看| 中文字幕色久视频| 一区二区日韩欧美中文字幕| 国产在线精品亚洲第一网站| 变态另类丝袜制服| 亚洲国产欧美网| 免费人成视频x8x8入口观看| 精品不卡国产一区二区三区| 国产色视频综合| 午夜日韩欧美国产| 亚洲电影在线观看av| 亚洲成av人片免费观看| 老汉色∧v一级毛片| 老汉色av国产亚洲站长工具| 亚洲色图综合在线观看| 国产午夜精品久久久久久| 亚洲午夜理论影院| а√天堂www在线а√下载| 国语自产精品视频在线第100页| 亚洲五月天丁香| 久久久久久免费高清国产稀缺| 免费在线观看完整版高清| 国产成人av教育| 日韩成人在线观看一区二区三区| 欧美色视频一区免费| 女人被躁到高潮嗷嗷叫费观| 免费在线观看亚洲国产| 精品人妻1区二区| 性欧美人与动物交配| avwww免费| 精品第一国产精品| 国产亚洲精品第一综合不卡| 大陆偷拍与自拍| 午夜福利影视在线免费观看| 男女之事视频高清在线观看|