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

    股骨有限元模型構(gòu)建及鈍力性損傷研究

    2017-09-05 19:23:01汪浩陳彩鳳仲梁維
    軟件導(dǎo)刊 2017年7期
    關(guān)鍵詞:CT掃描

    汪浩+陳彩鳳+仲梁維

    摘 要:對(duì)股骨進(jìn)行多層螺旋CT掃描,利用Mimics軟件將CT醫(yī)學(xué)影像數(shù)據(jù)轉(zhuǎn)化為有效的股骨有限元模型。將模型在Workbench顯示動(dòng)力學(xué)Ls-Dyna模塊中進(jìn)行網(wǎng)格劃分、材料賦值及其它前處理,用碰撞物以不同的速度對(duì)其正面不同位置進(jìn)行模擬打擊,得到股骨被撞擊結(jié)果,從而為不同股骨骨折受傷源頭提供依據(jù)。實(shí)驗(yàn)結(jié)果表明,碰撞物的速度不同、被撞擊的位置不同,骨折發(fā)生時(shí)間和位置也不一樣。

    關(guān)鍵詞:股骨有限元;CT掃描;顯示動(dòng)力學(xué)

    DOIDOI:10.11907/rjdk.171231

    中圖分類(lèi)號(hào):TP319

    文獻(xiàn)標(biāo)識(shí)碼:A 文章編號(hào):1672-7800(2017)007-0137-03

    0 引言

    鈍力性損傷主要包括運(yùn)動(dòng)的制傷物作用于相對(duì)靜止的人體和運(yùn)動(dòng)的人體撞擊靜止的物體,同類(lèi)型鈍力性損傷程度和形態(tài)是推斷暴力來(lái)源和致傷特征分析的主要依據(jù)。隨著計(jì)算機(jī)技術(shù)的發(fā)展,有限元方法和CT掃描技術(shù)結(jié)合應(yīng)用在生物工程學(xué)上,以解決一些人體損傷的問(wèn)題[1]。利用軟件建立和人體類(lèi)似的三維模型,用鈍性物體以不同的速度對(duì)人體不同位置進(jìn)行撞擊,分析應(yīng)力、應(yīng)變、位移、速度、加速度等應(yīng)力響應(yīng)參數(shù)變化,為鈍力性損傷鑒定提供依據(jù)[2]。本文擬用鋼管以不同的速度撞擊股骨的不同位置,得到股骨在不同條件下各時(shí)間點(diǎn)的應(yīng)力、應(yīng)變和位移情況,為相關(guān)問(wèn)題研究提供參考。

    1 有限元模型創(chuàng)建及研究

    1.1 有限元模型創(chuàng)建

    利用CT機(jī)及AW圖像后處理工作站對(duì)尸體股骨頭上方至膝關(guān)節(jié)處進(jìn)行掃描,將掃描得到的數(shù)據(jù)導(dǎo)入Mimics中,利用閾值分割工具將股骨和股骨中段的軟組織提取出來(lái),然后使用蒙罩編輯和區(qū)域增長(zhǎng)等工具分離出股骨和軟組織,之后用填充工具將空洞填充和修補(bǔ),將提取出來(lái)的蒙罩三維重建,得到股骨及其中段軟組織的三維模型。

    1.2 研究方法

    根據(jù)相關(guān)文獻(xiàn)和實(shí)驗(yàn)得到的數(shù)據(jù),將股骨簡(jiǎn)化為與應(yīng)變率無(wú)關(guān)的塑性材料模型,密度ρ=1 990kg/m3,彈性模量E=14 635MPa,泊松比μ=0.343 5,屈服應(yīng)力δs=133MPa,切線模量1 024MPa。失效準(zhǔn)則采用最大主應(yīng)變準(zhǔn)則[3],包括它的拉伸閾值(εTmax=0.073)和壓縮閾值(εCmax=0.010 4)。中段軟組織采用Mooney-Rivlin材料模型建立,材料相關(guān)屬性[4]如下:密度ρ=3 600kg/m3,兩個(gè)Mooney-Rivlin常數(shù)為C01=4.1、C10=0.41。碰撞的鋼管材料采用常見(jiàn)的結(jié)構(gòu)鋼,密度ρ=7 850kg/m3,彈性模量E=200 000MPa,泊松比μ=0.3,體積模量K=166 670MPa,剪切模量G=76 923MPa。

    因?yàn)槭悄M股骨被鈍物撞擊,涉及到瞬態(tài)、大應(yīng)變、大型變、材料的破壞、材料的完全失效或伴隨復(fù)雜接觸的結(jié)構(gòu)問(wèn)題,所以通過(guò)Ansys顯示動(dòng)力學(xué)求解。Ansys顯示動(dòng)力學(xué)分析模塊包括Explicit Dynamics、ANSYS AUTODYN及Workbench 16.0的Ls-Dyna。

    Ls-Dyna是世界上最著名的通用顯式非線性有限元分析程序,能夠模擬真實(shí)世界的各種復(fù)雜問(wèn)題,特別適合求解各種二維、三維非線性結(jié)構(gòu)碰撞、金屬成型等非線性動(dòng)力沖擊問(wèn)題,它以Lagrange算法為主,兼有ALE和Euler算法;以顯式求解為主,兼有隱式求解功能;以非線性動(dòng)力分析為主,兼有靜力分析功能(如動(dòng)力分析前的預(yù)應(yīng)力計(jì)算和薄板沖壓成型后的回彈計(jì)算),是通用的結(jié)構(gòu)分析非線性有限元程序[5],操作和結(jié)果查看都十分快捷,本文采用Ls-Dyna模塊進(jìn)行股骨模擬碰撞[6]。

    為了模擬人體股骨受打擊情況,給股骨模型上端施加一個(gè)垂直的力,考慮人體情況,擬添加一個(gè)250N的壓力,股骨的下端固定,用直徑為40mm、壁厚為5mm、長(zhǎng)度為1.5m的結(jié)構(gòu)鋼鋼管正面打擊,分別以20m/s、10m/s的速度去撞擊股骨模型中段,然后以10m/s的速度撞擊股骨模型靠近膝關(guān)節(jié)的位置,用Workbench計(jì)算股骨和鋼管的應(yīng)力變化云圖,得到股骨的von mises應(yīng)力變化表和股骨在不同時(shí)間點(diǎn)的應(yīng)變、形變情況[7],得到股骨中段被鋼管以不同速度撞擊和股骨中段、下段被相同速度撞擊的生物力學(xué)響應(yīng),從而判斷股骨骨折情況,為司法鑒定中已知股骨的受傷情況來(lái)反推受傷源頭提供一定的科學(xué)依據(jù)。

    2 分析過(guò)程及結(jié)果

    2.1 添加材料

    使用Worbench16.0進(jìn)行分析。打開(kāi)Workbench16.0,啟動(dòng)Ls-Dyna模塊。在Engineering Data中選擇新建一種材料,這種材料用于股骨的各項(xiàng)彈塑性材料模型。在Tool Box中依次選擇密度、各項(xiàng)同性彈塑性材料、雙線性各項(xiàng)同性硬化和主要的應(yīng)變失效,密度中填入1 990kg/m3,彈性模量和泊松比分別為14 635MPa和0.343 5,體積模量和剪切模量分別為15 586MPa和5 446.6MPa,屈服強(qiáng)度和切線模量分別為133MPa和1 024MPa,在主要應(yīng)變失效中將拉伸閾值和壓縮閾值分別設(shè)置為0.073和0.010 4。然后再新建一種Mooney-Rivlin材料模型用于股骨中段的軟組織[8],密度設(shè)置為3 600kg/m3,在超彈性的材料屬性中選擇兩個(gè)Mooney-Rivlin常數(shù)選項(xiàng),在C01和C10中填入4.1和0.41。

    2.2 導(dǎo)入模型,劃分網(wǎng)格

    導(dǎo)入已經(jīng)建好的股骨和鋼管模型。因?yàn)槭菑?fù)雜曲面,首先需要對(duì)模型表面的復(fù)雜曲面進(jìn)行化簡(jiǎn),將短小細(xì)線進(jìn)行拓?fù)浣Y(jié)構(gòu)上的簡(jiǎn)化;其次對(duì)各面進(jìn)行段數(shù)控制,控制網(wǎng)格大小,在完成面網(wǎng)格劃分之后再進(jìn)行體網(wǎng)格劃分。

    2.3 分析前處理

    添加股骨、股骨中段軟組織還有鋼管材料,股骨上端施加300N的力,下端固定,撞擊終止時(shí)間0.01s。endprint

    2.4 結(jié)果

    當(dāng)鋼管以速度20m/s撞擊股骨中段時(shí),得到撞擊股骨的von mises等效應(yīng)力、時(shí)間變化曲線如圖1所示。

    由von mises等效應(yīng)力變化曲線可知,在時(shí)間為0.5ms時(shí)的von mises等效應(yīng)力就已經(jīng)超過(guò)了股骨的屈服極限133MPa,這時(shí)股骨已經(jīng)進(jìn)入塑性階段,該時(shí)間的應(yīng)力云圖如圖2所示。

    由圖2可以發(fā)現(xiàn),在這個(gè)時(shí)刻股骨靠近膝關(guān)節(jié)部分受到的應(yīng)力最大,為137.51GPa,超過(guò)了股骨的屈服極限,此處最易發(fā)生骨折。求解此時(shí)股骨總位移量,得到股骨位移云圖如圖3所示。

    在0.5ms時(shí),股骨被撞擊的中間段總位移量最大,達(dá)到了5.4445mm,伸長(zhǎng)比為1.09%,還未達(dá)到股骨的最大伸長(zhǎng)比1.41%[9],所以這時(shí)該位置還未拉斷。

    現(xiàn)在將撞擊速度改為10m/s,得到撞擊最后股骨的von mises等效應(yīng)力、時(shí)間變化曲線如圖4所示。

    當(dāng)撞擊時(shí)間為1.5ms時(shí),得到股骨von mises等效應(yīng)力云圖和位移云圖如圖5、圖6所示。

    由圖5、圖6可知,在1.5ms時(shí)股骨的von mises最大等效應(yīng)力在股骨兩端,為137.09MPa,超過(guò)了股骨的屈服極限,而此刻的最大位移發(fā)生在股骨中段,為7.4217mm,伸長(zhǎng)比為1.48%,已經(jīng)超過(guò)了股骨的最大伸長(zhǎng)比1.41%,所以這時(shí)股骨的中段和兩端都易發(fā)生骨折現(xiàn)象。

    撞擊速度依然為10m/s,現(xiàn)在用鋼管去撞擊股骨靠近膝關(guān)節(jié)的位置,得到撞擊最后股骨的von mises等效應(yīng)力、時(shí)間變化曲線如圖7所示。

    當(dāng)撞擊時(shí)間為0.5ms時(shí),股骨von mises的等效應(yīng)力云圖和位移云圖如圖8、圖9所示。

    由圖8、圖9可知,在0.5ms時(shí),股骨膝關(guān)節(jié)位置的von mises等效應(yīng)力最大,為156.93MPa,超過(guò)了股骨的屈服極限,此刻的最大位移發(fā)生在鋼管和股骨的碰撞位置,為2.3522mm,伸長(zhǎng)比為0.47%,還未達(dá)到股骨的最大伸長(zhǎng)比1.41%。

    當(dāng)股骨中段受到鋼管20m/s的高速正面撞擊時(shí),在0.5ms時(shí)股骨兩端的von mises等效應(yīng)力就已經(jīng)超過(guò)股骨屈服極限,而股骨被撞擊位置的伸長(zhǎng)比還未達(dá)到股骨的極限伸長(zhǎng)比,所以只有股骨兩端發(fā)生了骨折;當(dāng)股骨中段受到鋼管10m/s的低速正面撞擊時(shí),在1.5ms時(shí)股骨兩端的von mises等效應(yīng)力超過(guò)了股骨的屈服極限,而被撞擊位置的伸長(zhǎng)比也超過(guò)了股骨的極限伸長(zhǎng)比,所以股骨兩端和被撞擊位置都發(fā)生了骨折;當(dāng)股骨靠近膝關(guān)節(jié)的位置受到鋼管10m/s的低速正面撞擊時(shí),在0.5ms時(shí)膝關(guān)節(jié)位置的von mises等效應(yīng)力超過(guò)了股骨的屈服極限,而股骨被撞擊位置的伸長(zhǎng)比還未達(dá)到股骨的極限伸長(zhǎng)比,所以只有股骨的膝關(guān)節(jié)位置發(fā)生骨折。

    股骨被撞擊的骨折情況和撞擊物的速度、撞擊位置都有關(guān)系。當(dāng)撞擊物以高速撞擊股骨中段時(shí),骨折發(fā)生得快,而且骨折位置主要集中在股骨兩端;當(dāng)撞擊物以低速撞擊股骨中段時(shí),骨折發(fā)生得慢,而且在股骨兩端和中間會(huì)發(fā)生骨折; 當(dāng)撞擊物以低速撞擊股骨靠近膝關(guān)節(jié)位置時(shí),骨折發(fā)生得快,而且骨折位置集中在膝關(guān)節(jié)位置。

    3 結(jié)語(yǔ)

    本研究對(duì)股骨和其它關(guān)節(jié)部位只做了簡(jiǎn)單固定,而人體骨骼是非線性、粘彈性以及非連續(xù)的材料[10],在有限元建模中,受到軟件材料庫(kù)的限制,只能簡(jiǎn)化為各項(xiàng)同性彈塑性材料和Mooney-Rivlin材料。

    本文通過(guò)CT掃描圖像和Mimics軟件制作出股骨的有限元模型,并進(jìn)行了損傷生物力學(xué)分析,可從股骨的骨折情況大概推算出撞擊物的速度、撞擊位置等,為司法鑒定提供依據(jù)。

    參考文獻(xiàn):

    [1]陳憶九,鄒冬華,邵煜.有限元方法在法醫(yī)損傷鑒定中的應(yīng)用[J].中國(guó)司法鑒定,2010(5):28-32.

    [2]劉子建,張建華,楊濟(jì)匡.碰撞生物力學(xué)基礎(chǔ)及其應(yīng)用[J].中華創(chuàng)傷雜志,2015,17(5):261-263.

    [3]BAYRAKTAR HH,MORGAN EF,NIEBUR GL.Comparison of the elastic and yield properties of human femoral trabecular and cortical bone tissue[J].J Biomech,2004,37(1): 27-35.

    [4]ANDERSON AE,PETERS CL,TUTTLE BD.Subjectspecific finite element model of the pelvis: development,validation and sensitivity studies[J]. Biomech Eng,2005,127(3):364-373.

    [5]唐長(zhǎng)剛.Ls-dyna有限元分析及仿真[M].北京:電子工業(yè)出版社,2013.

    [6]張冠軍,魏搖嵬,曹立波,等.行人膝關(guān)節(jié)有限元模型的建立與驗(yàn)證[J]. 汽車(chē)工程,2014,36(9):1069-1073.

    [7]李正東,劉寧國(guó),黃平,等.下肢有限元模型的建立及損傷機(jī)制重建[J]. 中國(guó)司法鑒定,2012(6):37-42.

    [8]周家澤.工程材料彈塑性應(yīng)力應(yīng)變模型分析[J]. 武漢科技學(xué)院學(xué)報(bào),2003,16(6):42-45.

    [9]馮元楨.生物力學(xué)[M].戴克剛,鞠烽熾,譯.長(zhǎng)沙:湖南科學(xué)技術(shù)出版社,1986.

    [10]CHEN Y B,MIAO Y Y,XU C A,et al.Wound ballistics of the pig mandibular angle:a preliminary finite element anal-ysis and experimental study[J].J Biomech,2010,43(6):1131-1137.endprint

    猜你喜歡
    CT掃描
    創(chuàng)傷骨科中CT掃描分析研究
    商情(2018年24期)2018-07-06 13:38:28
    CT及MRI在診斷三叉神經(jīng)瘤中的臨床價(jià)值分析
    MRI、CT單獨(dú)或聯(lián)合應(yīng)用在肝細(xì)胞癌中的診斷價(jià)值
    腦梗死及腦血腫患者復(fù)查中采用CT低劑量掃描的可行性
    探究CT影像診斷在臨床醫(yī)學(xué)中的應(yīng)用價(jià)值
    科技資訊(2017年4期)2017-04-23 12:27:25
    X線平片及CT掃描在脊柱爆裂骨折診斷中的應(yīng)用價(jià)值
    今日健康(2016年7期)2017-04-12 01:54:45
    淺談CT在先天性心血管疾病診斷中的應(yīng)用
    支氣管肺癌的影像診斷研究
    卷宗(2016年4期)2016-05-30 11:01:28
    CT掃描診斷腦血管病的療效分析
    中老年人心源性肺水腫影像診斷
    69av精品久久久久久 | 久久综合国产亚洲精品| 无限看片的www在线观看| 一级黄色大片毛片| 国产亚洲av片在线观看秒播厂| 日本黄色日本黄色录像| 久久九九热精品免费| 亚洲欧洲日产国产| 欧美日韩精品网址| 国产一区有黄有色的免费视频| 国产免费现黄频在线看| 精品久久久久久电影网| 天天躁日日躁夜夜躁夜夜| 中文字幕制服av| 亚洲av美国av| 亚洲国产欧美一区二区综合| 欧美另类一区| 丝袜脚勾引网站| 天天躁狠狠躁夜夜躁狠狠躁| 日本a在线网址| 视频区图区小说| 亚洲精品久久午夜乱码| 在线亚洲精品国产二区图片欧美| 亚洲av美国av| 欧美日韩一级在线毛片| 亚洲精品久久午夜乱码| 黄色 视频免费看| 亚洲,欧美精品.| 亚洲精华国产精华精| 18禁裸乳无遮挡动漫免费视频| 在线观看免费高清a一片| 国产一级毛片在线| 99国产精品免费福利视频| 汤姆久久久久久久影院中文字幕| 法律面前人人平等表现在哪些方面 | 精品福利观看| 少妇被粗大的猛进出69影院| 免费人妻精品一区二区三区视频| 亚洲 欧美一区二区三区| 性色av一级| 久热这里只有精品99| 男女下面插进去视频免费观看| 亚洲成国产人片在线观看| netflix在线观看网站| 狠狠精品人妻久久久久久综合| 亚洲第一青青草原| 欧美日韩福利视频一区二区| 国产免费一区二区三区四区乱码| 亚洲少妇的诱惑av| 色94色欧美一区二区| 日韩有码中文字幕| 9191精品国产免费久久| 午夜免费成人在线视频| 中文字幕人妻丝袜制服| 免费久久久久久久精品成人欧美视频| 欧美日韩精品网址| 欧美黑人欧美精品刺激| 黄色片一级片一级黄色片| 中文字幕av电影在线播放| 飞空精品影院首页| 国产在线免费精品| av福利片在线| 无限看片的www在线观看| 丝袜喷水一区| 亚洲国产欧美一区二区综合| 国产精品 欧美亚洲| 亚洲精品一卡2卡三卡4卡5卡 | 女人精品久久久久毛片| 岛国在线观看网站| 亚洲第一青青草原| 波多野结衣一区麻豆| 久久99热这里只频精品6学生| 亚洲精品国产区一区二| 国产欧美亚洲国产| 国产精品亚洲av一区麻豆| 精品卡一卡二卡四卡免费| 老熟女久久久| 51午夜福利影视在线观看| 久久久久久亚洲精品国产蜜桃av| 在线 av 中文字幕| 亚洲,欧美精品.| 国产一区二区三区av在线| 黄色a级毛片大全视频| 悠悠久久av| 19禁男女啪啪无遮挡网站| 黄色片一级片一级黄色片| 国产免费一区二区三区四区乱码| 欧美精品高潮呻吟av久久| 成年美女黄网站色视频大全免费| 大码成人一级视频| 久久久久国产精品人妻一区二区| 他把我摸到了高潮在线观看 | 国产精品一区二区在线不卡| 黄网站色视频无遮挡免费观看| 美女高潮喷水抽搐中文字幕| 男人操女人黄网站| 午夜福利影视在线免费观看| 精品人妻在线不人妻| 一级,二级,三级黄色视频| 搡老乐熟女国产| 黄网站色视频无遮挡免费观看| 99re6热这里在线精品视频| 久9热在线精品视频| av国产精品久久久久影院| 久久久久久久国产电影| 国产精品自产拍在线观看55亚洲 | 波多野结衣一区麻豆| 99国产精品99久久久久| 亚洲国产看品久久| 一区在线观看完整版| 中文字幕av电影在线播放| 国产精品影院久久| 一级片免费观看大全| 国产精品 欧美亚洲| 亚洲成人免费电影在线观看| 国产成人精品在线电影| 日本猛色少妇xxxxx猛交久久| 水蜜桃什么品种好| 高清视频免费观看一区二区| 夜夜夜夜夜久久久久| 超碰成人久久| 欧美午夜高清在线| 国产麻豆69| 9191精品国产免费久久| 在线永久观看黄色视频| 国产欧美日韩一区二区精品| 99国产精品一区二区三区| 日韩欧美一区视频在线观看| 波多野结衣一区麻豆| 亚洲成人手机| 亚洲第一av免费看| 一区二区三区精品91| 亚洲精品久久午夜乱码| 在线观看一区二区三区激情| 桃红色精品国产亚洲av| 精品高清国产在线一区| 精品一品国产午夜福利视频| av福利片在线| 亚洲午夜精品一区,二区,三区| 在线永久观看黄色视频| 日韩制服骚丝袜av| 97精品久久久久久久久久精品| www.999成人在线观看| bbb黄色大片| 老司机深夜福利视频在线观看 | av天堂久久9| 亚洲一区中文字幕在线| 50天的宝宝边吃奶边哭怎么回事| 美女主播在线视频| 亚洲色图 男人天堂 中文字幕| 国产成人欧美在线观看 | 巨乳人妻的诱惑在线观看| 欧美日本中文国产一区发布| av网站免费在线观看视频| 国产精品秋霞免费鲁丝片| 久久人人爽av亚洲精品天堂| 曰老女人黄片| 国产亚洲精品第一综合不卡| 日韩三级视频一区二区三区| 一级,二级,三级黄色视频| 欧美黄色淫秽网站| 黄色毛片三级朝国网站| 天堂俺去俺来也www色官网| 国产av一区二区精品久久| 亚洲九九香蕉| 欧美久久黑人一区二区| 欧美97在线视频| 日韩大码丰满熟妇| 一本久久精品| 男女床上黄色一级片免费看| 大码成人一级视频| 男人舔女人的私密视频| 永久免费av网站大全| 久久人人爽av亚洲精品天堂| 日韩视频在线欧美| 别揉我奶头~嗯~啊~动态视频 | 日韩人妻精品一区2区三区| 亚洲情色 制服丝袜| 老熟女久久久| 一区二区三区精品91| 91精品伊人久久大香线蕉| 美女午夜性视频免费| 母亲3免费完整高清在线观看| 久久青草综合色| 新久久久久国产一级毛片| 亚洲精品自拍成人| 亚洲精品国产av成人精品| 国产欧美日韩一区二区三 | 又紧又爽又黄一区二区| 免费高清在线观看日韩| 亚洲情色 制服丝袜| 国产精品影院久久| 青草久久国产| 一本久久精品| 大香蕉久久网| 国产高清国产精品国产三级| 欧美精品高潮呻吟av久久| 在线观看免费日韩欧美大片| 女性被躁到高潮视频| 精品少妇黑人巨大在线播放| 国产激情久久老熟女| 老熟妇仑乱视频hdxx| 亚洲一码二码三码区别大吗| 亚洲av片天天在线观看| www.999成人在线观看| 日韩中文字幕视频在线看片| 亚洲专区中文字幕在线| 国产黄频视频在线观看| 亚洲国产精品成人久久小说| 久久久水蜜桃国产精品网| 老司机午夜福利在线观看视频 | videos熟女内射| 在线观看免费日韩欧美大片| 老司机在亚洲福利影院| 日韩视频一区二区在线观看| 黄色怎么调成土黄色| 亚洲欧美色中文字幕在线| 亚洲美女黄色视频免费看| 亚洲自偷自拍图片 自拍| 99国产精品免费福利视频| 精品欧美一区二区三区在线| 免费在线观看影片大全网站| 99国产精品一区二区蜜桃av | 国产人伦9x9x在线观看| 搡老岳熟女国产| 国产一区二区三区av在线| 国产精品熟女久久久久浪| 国产一区二区三区综合在线观看| 亚洲欧洲精品一区二区精品久久久| 精品亚洲乱码少妇综合久久| 久9热在线精品视频| 国产精品久久久久久精品电影小说| 黄色视频不卡| 两个人看的免费小视频| 精品一品国产午夜福利视频| 中文字幕色久视频| av线在线观看网站| 男人操女人黄网站| 国产精品国产三级国产专区5o| 叶爱在线成人免费视频播放| 午夜福利视频精品| 老汉色av国产亚洲站长工具| 亚洲欧洲精品一区二区精品久久久| 国产亚洲精品第一综合不卡| 三级毛片av免费| 亚洲三区欧美一区| 伊人亚洲综合成人网| 色播在线永久视频| 国产三级黄色录像| 精品国产乱子伦一区二区三区 | 精品国产乱子伦一区二区三区 | 精品国内亚洲2022精品成人 | 97精品久久久久久久久久精品| 成在线人永久免费视频| 人人妻人人添人人爽欧美一区卜| 欧美午夜高清在线| 午夜成年电影在线免费观看| 欧美97在线视频| 欧美另类一区| 国产xxxxx性猛交| 三级毛片av免费| 久久久精品区二区三区| 中国美女看黄片| 日韩免费高清中文字幕av| 十八禁网站网址无遮挡| 丰满人妻熟妇乱又伦精品不卡| 日本五十路高清| 欧美日韩视频精品一区| 国产视频一区二区在线看| 男女之事视频高清在线观看| 国产99久久九九免费精品| av视频免费观看在线观看| 十分钟在线观看高清视频www| 免费在线观看影片大全网站| 欧美精品一区二区免费开放| 国产精品国产三级国产专区5o| 最近最新中文字幕大全免费视频| 亚洲中文日韩欧美视频| 美女高潮到喷水免费观看| 国产高清videossex| 国产在线观看jvid| 精品一区二区三区四区五区乱码| 成年av动漫网址| 青草久久国产| 天堂俺去俺来也www色官网| 欧美在线一区亚洲| 色老头精品视频在线观看| 成在线人永久免费视频| 王馨瑶露胸无遮挡在线观看| av有码第一页| 亚洲国产精品一区三区| 18禁国产床啪视频网站| av一本久久久久| 巨乳人妻的诱惑在线观看| 菩萨蛮人人尽说江南好唐韦庄| 不卡一级毛片| 日本av免费视频播放| 亚洲精品一卡2卡三卡4卡5卡 | 精品少妇内射三级| 亚洲男人天堂网一区| 免费久久久久久久精品成人欧美视频| 欧美97在线视频| 十分钟在线观看高清视频www| 亚洲精品国产av蜜桃| 免费高清在线观看视频在线观看| 国产视频一区二区在线看| 纵有疾风起免费观看全集完整版| 久久久久网色| 纵有疾风起免费观看全集完整版| 欧美激情久久久久久爽电影 | 亚洲av男天堂| 99久久人妻综合| 后天国语完整版免费观看| 久热爱精品视频在线9| 亚洲七黄色美女视频| 精品久久久精品久久久| 黑丝袜美女国产一区| 国产免费视频播放在线视频| 一二三四社区在线视频社区8| 欧美另类一区| 国产又色又爽无遮挡免| 国产高清videossex| 久久精品熟女亚洲av麻豆精品| 欧美黄色淫秽网站| 欧美av亚洲av综合av国产av| 在线av久久热| 国产精品久久久av美女十八| 老司机午夜福利在线观看视频 | 国产欧美日韩一区二区三 | 日本91视频免费播放| 亚洲国产av影院在线观看| 飞空精品影院首页| 在线观看免费午夜福利视频| 热re99久久精品国产66热6| 日韩人妻精品一区2区三区| 老司机亚洲免费影院| 亚洲性夜色夜夜综合| 一级毛片精品| 性高湖久久久久久久久免费观看| 在线观看免费高清a一片| 欧美日本中文国产一区发布| 中文字幕av电影在线播放| www.999成人在线观看| 久久久久精品国产欧美久久久 | 十分钟在线观看高清视频www| av线在线观看网站| 亚洲av男天堂| 亚洲欧美色中文字幕在线| 十分钟在线观看高清视频www| av天堂在线播放| 99热国产这里只有精品6| 丝袜在线中文字幕| 女警被强在线播放| 好男人电影高清在线观看| 免费观看a级毛片全部| 久久久精品区二区三区| 一区在线观看完整版| 亚洲精品中文字幕在线视频| av福利片在线| 国产深夜福利视频在线观看| 色婷婷av一区二区三区视频| 午夜福利视频在线观看免费| 亚洲国产欧美在线一区| 久久狼人影院| 高清在线国产一区| 国产精品久久久av美女十八| 一区在线观看完整版| 午夜激情av网站| 欧美老熟妇乱子伦牲交| 99精品久久久久人妻精品| 免费一级毛片在线播放高清视频 | 亚洲情色 制服丝袜| 18禁观看日本| 久久免费观看电影| 国产1区2区3区精品| 蜜桃国产av成人99| 亚洲av日韩在线播放| 脱女人内裤的视频| 日韩精品免费视频一区二区三区| 午夜精品久久久久久毛片777| 精品久久久精品久久久| 成年动漫av网址| 老司机午夜十八禁免费视频| 老司机在亚洲福利影院| 精品人妻熟女毛片av久久网站| 爱豆传媒免费全集在线观看| 欧美 日韩 精品 国产| 国产欧美日韩一区二区精品| 自线自在国产av| 日韩有码中文字幕| 极品少妇高潮喷水抽搐| 亚洲精品乱久久久久久| 在线观看免费视频网站a站| www.精华液| 国产精品香港三级国产av潘金莲| 男女国产视频网站| 亚洲久久久国产精品| av国产精品久久久久影院| 中文字幕人妻丝袜一区二区| 黄色a级毛片大全视频| 高清视频免费观看一区二区| 亚洲色图综合在线观看| 精品久久久久久久毛片微露脸 | 99国产精品免费福利视频| 看免费av毛片| 天天添夜夜摸| 欧美性长视频在线观看| 精品福利永久在线观看| 久久人人爽人人片av| 国产黄频视频在线观看| 男人舔女人的私密视频| 18禁裸乳无遮挡动漫免费视频| 亚洲第一av免费看| av一本久久久久| 欧美黄色片欧美黄色片| 免费一级毛片在线播放高清视频 | 成人国产av品久久久| 国产一卡二卡三卡精品| 久久青草综合色| videosex国产| 国产在视频线精品| 一级片免费观看大全| 国产精品国产三级国产专区5o| 如日韩欧美国产精品一区二区三区| 美女视频免费永久观看网站| 少妇人妻久久综合中文| 国产日韩欧美在线精品| 国产高清videossex| 欧美黑人欧美精品刺激| 大片免费播放器 马上看| 69av精品久久久久久 | 在线永久观看黄色视频| 免费在线观看日本一区| 精品久久久久久久毛片微露脸 | 亚洲中文av在线| 9191精品国产免费久久| 69av精品久久久久久 | 国产真人三级小视频在线观看| 国产精品免费大片| 一本一本久久a久久精品综合妖精| 国产高清videossex| 午夜免费观看性视频| 青春草亚洲视频在线观看| 91麻豆精品激情在线观看国产 | 久久国产精品大桥未久av| 亚洲专区字幕在线| 亚洲精品自拍成人| 蜜桃国产av成人99| 国产区一区二久久| 两个人免费观看高清视频| 久久午夜综合久久蜜桃| 国产一区二区激情短视频 | 麻豆av在线久日| 999精品在线视频| 国产精品秋霞免费鲁丝片| 国产精品麻豆人妻色哟哟久久| 国产精品 欧美亚洲| 亚洲久久久国产精品| 日日夜夜操网爽| 在线天堂中文资源库| 精品一区二区三区av网在线观看 | 性色av一级| 黄频高清免费视频| 中文字幕高清在线视频| 久久久欧美国产精品| 自拍欧美九色日韩亚洲蝌蚪91| 丝瓜视频免费看黄片| 99精品久久久久人妻精品| 亚洲中文字幕日韩| a级毛片在线看网站| 国产精品久久久久成人av| 黄色视频,在线免费观看| 国产成人a∨麻豆精品| 久久久水蜜桃国产精品网| 亚洲情色 制服丝袜| 国产一卡二卡三卡精品| 一级毛片电影观看| 欧美黑人精品巨大| 狂野欧美激情性xxxx| 精品亚洲成a人片在线观看| 欧美日韩福利视频一区二区| 超碰97精品在线观看| 18在线观看网站| 日韩 亚洲 欧美在线| 纵有疾风起免费观看全集完整版| 我要看黄色一级片免费的| 国产有黄有色有爽视频| 欧美乱码精品一区二区三区| 90打野战视频偷拍视频| 在线观看免费高清a一片| 久久久水蜜桃国产精品网| 黄色怎么调成土黄色| 亚洲成av片中文字幕在线观看| 国产精品一区二区免费欧美 | 正在播放国产对白刺激| 日韩电影二区| 国产激情久久老熟女| 美女中出高潮动态图| 久久ye,这里只有精品| 国产91精品成人一区二区三区 | 日韩免费高清中文字幕av| 免费看十八禁软件| 老熟女久久久| 一级毛片电影观看| 色94色欧美一区二区| 老司机在亚洲福利影院| 国产欧美日韩综合在线一区二区| 美女中出高潮动态图| 久久久久久久精品精品| 狠狠婷婷综合久久久久久88av| 国产精品免费视频内射| 一区二区av电影网| 国产真人三级小视频在线观看| 欧美国产精品一级二级三级| 十八禁网站免费在线| 国产精品久久久久成人av| 大香蕉久久成人网| 久久久国产欧美日韩av| 亚洲第一欧美日韩一区二区三区 | 91av网站免费观看| 女人爽到高潮嗷嗷叫在线视频| 叶爱在线成人免费视频播放| 亚洲专区字幕在线| 狠狠狠狠99中文字幕| 久久精品亚洲熟妇少妇任你| 亚洲欧美精品自产自拍| 王馨瑶露胸无遮挡在线观看| 久久国产精品大桥未久av| 一区二区三区乱码不卡18| 成人亚洲精品一区在线观看| 中文字幕高清在线视频| 精品久久久久久久毛片微露脸 | 黄色 视频免费看| 69精品国产乱码久久久| 每晚都被弄得嗷嗷叫到高潮| 菩萨蛮人人尽说江南好唐韦庄| 国产xxxxx性猛交| 成人18禁高潮啪啪吃奶动态图| 国产一区二区三区综合在线观看| 9191精品国产免费久久| 中文字幕人妻丝袜制服| 手机成人av网站| 视频在线观看一区二区三区| 国产日韩欧美视频二区| 在线天堂中文资源库| 中文字幕人妻丝袜一区二区| 国产成人精品久久二区二区免费| 狠狠狠狠99中文字幕| 国产在视频线精品| 香蕉丝袜av| 日韩免费高清中文字幕av| 欧美在线一区亚洲| 国产精品一二三区在线看| 美女主播在线视频| 午夜免费成人在线视频| 可以免费在线观看a视频的电影网站| 国产一区有黄有色的免费视频| 可以免费在线观看a视频的电影网站| 日韩 欧美 亚洲 中文字幕| 91麻豆av在线| 男人添女人高潮全过程视频| 男女高潮啪啪啪动态图| 高清视频免费观看一区二区| 黄色毛片三级朝国网站| 免费高清在线观看日韩| 我要看黄色一级片免费的| 亚洲黑人精品在线| 国产精品国产三级国产专区5o| 国产av又大| av线在线观看网站| 亚洲第一欧美日韩一区二区三区 | 两性午夜刺激爽爽歪歪视频在线观看 | 妹子高潮喷水视频| 久久久久视频综合| 日韩三级视频一区二区三区| 老司机在亚洲福利影院| 老司机午夜十八禁免费视频| 人妻一区二区av| 久久精品亚洲熟妇少妇任你| 99热国产这里只有精品6| 国产精品国产av在线观看| 亚洲精品国产av成人精品| 久久久国产成人免费| 国产淫语在线视频| 王馨瑶露胸无遮挡在线观看| 国产精品成人在线| 日韩电影二区| 黑人操中国人逼视频| 青草久久国产| 成年动漫av网址| 久久精品国产综合久久久| 国产一区二区三区av在线| 精品一品国产午夜福利视频| 亚洲欧美一区二区三区久久| 中文欧美无线码| 成年动漫av网址| 咕卡用的链子| 国产一区二区三区av在线| 欧美黑人精品巨大| 一本久久精品| 欧美日韩一级在线毛片| kizo精华| 国产老妇伦熟女老妇高清| 国产一区二区在线观看av| 日本精品一区二区三区蜜桃| 老熟妇乱子伦视频在线观看 | 欧美在线黄色| 91精品伊人久久大香线蕉| 在线观看免费日韩欧美大片| 日韩电影二区| 搡老熟女国产l中国老女人| 69精品国产乱码久久久| 国产欧美日韩一区二区三 | 国产成人精品在线电影| 国产无遮挡羞羞视频在线观看|