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

    脊柱手法的力學(xué)研究概況

    2014-11-29 00:49:04郭偉李藝韓磊田飛鵬閔亞青周衛(wèi)趙平
    頸腰痛雜志 2014年5期
    關(guān)鍵詞:有限元模型研究

    郭偉,李藝,韓磊,田飛鵬,閔亞青,周衛(wèi),趙平

    (1.空軍總醫(yī)院正骨科,100142;2.望京醫(yī)院,北京 100102)

    脊柱手法治療脊柱損傷退變性疾病在我國已有數(shù)千年的歷史。對其作用機制曾有著豐富的論述。時至今日,脊柱手法仍然是眾多國家衛(wèi)生防治指南中最常建議使用的治療方法之一。目前手法主要有脊柱關(guān)節(jié)手法和椎旁軟組織手法,其中脊柱關(guān)節(jié)手法主流手法分為關(guān)節(jié)調(diào)整或沖擊手法和關(guān)節(jié)松動手法。關(guān)節(jié)調(diào)整手法(Manipulation)或沖擊手法(Thrusting)[1]是利用作用于特殊解剖位置的長或短杠桿手段達(dá)到的關(guān)節(jié)手法。其特征是在一定的速率、振幅和方向下振動推力,既國際手法學(xué)術(shù)界慣常稱謂的(High Velocity Low Amplitude)HVLA特征。這種手法經(jīng)常伴有彈響聲(空腔效應(yīng)),比如傳統(tǒng)中醫(yī)的斜扳手法和美式整脊手法。關(guān)節(jié)松動術(shù)(Mobilization)是一種被動運動療法,可以提升關(guān)節(jié)活動度,但不超過其解剖范圍,一個非推力的關(guān)節(jié)運動通常不會超過關(guān)節(jié)的彈性界限,并不常見明顯的關(guān)節(jié)彈響聲。脊柱松動術(shù)是一個在關(guān)節(jié)活動范圍內(nèi)或在關(guān)節(jié)極度運動時更溫和且更易于接受的可操作外力,比如馮氏手法和美式整脊的關(guān)節(jié)松動手法。

    隨著計算機仿真技術(shù)的發(fā)展,脊柱手法的模擬力學(xué)檢測及脊柱力學(xué)模型仿真技術(shù)越來越趨于成熟,并得到了很大應(yīng)用和發(fā)展。很多研究者運用力傳感測量技術(shù)、光學(xué)定位捕捉技術(shù),通過腰椎CT建立腰椎多體動力學(xué)模型,通過力學(xué)測量系統(tǒng)測量腰部手法操作的力學(xué)、運動學(xué)變化過程等,為揭示手法治療脊柱損傷性疾病的機理,分析脊柱手法技術(shù)操作特點及影響因素提供了科學(xué)量化的方法,并對臨床治療提供了寶貴的基礎(chǔ)研究數(shù)據(jù),本文總結(jié)了國內(nèi)外關(guān)于脊柱手法的力學(xué)機制研究和計算機仿真的研究進(jìn)展,綜述如下。

    1 脊柱手法活體力學(xué)研究

    1.1關(guān)節(jié)突關(guān)節(jié)壓力的測量

    以往相關(guān)學(xué)者研究,測量關(guān)節(jié)突關(guān)節(jié)面壓力的方法很多。分為直接測量法和間接測量法。早期多以間接測量法為主。

    1.1.1間接測量法 通過關(guān)節(jié)突關(guān)節(jié)的受力情況靠測量有關(guān)數(shù)據(jù)來推算。沒有破壞關(guān)節(jié)突關(guān)節(jié)的完整性是間接法的優(yōu)點,但其結(jié)果的可靠程度有待商榷。

    1.1.2直接測量法 常見如壓力傳感器法。近年來隨著相關(guān)技術(shù)的高速發(fā)展,壓敏片[2]和其它測量方法相比,方法相對簡單,能同時獲得較完善的相關(guān)信息,能反映出應(yīng)力分布的整體趨勢,這是其他壓力傳感器法所無法比擬的。

    1.1.3有限元分析法 通過脊柱三維數(shù)學(xué)模型的建立來模擬承載的實際情況和關(guān)節(jié)突關(guān)節(jié)的運動(具體總結(jié)見5)。Triano等[3]作用時間在200-300 ms情況下,使用反動力模型估算旋轉(zhuǎn)手法的旋扳力是111N和123N。

    1.2軟組織按摩手法的力學(xué)研究

    很多研究者做了有益的嘗試,鄧玫等[4]應(yīng)用Ergocheck系統(tǒng)壓力傳感器進(jìn)行測試。從而得出由于著力面積越小,壓強值最大。在力學(xué)指標(biāo)方面國外研究者進(jìn)行了深入定量的研究。Hessell[5]發(fā)現(xiàn)預(yù)加載力與較大的作用力密切相關(guān)。而Herzog等[6]則實測施術(shù)者旋轉(zhuǎn)手法時的各種力學(xué)變化。手法的力學(xué)研究應(yīng)多多參考以上研究方法。

    1.3頸椎旋轉(zhuǎn)手法活體測量

    目前國內(nèi)外研究比較多的是頸椎手法的生物力學(xué)分析。Jimmw等[7]對頸椎扳轉(zhuǎn)時的平均加速度、最大加速度與持續(xù)時間均進(jìn)行了量化測量,為手法的力學(xué)測量提供了很多參考指標(biāo)。Van zoest[8]等直接測量手法操作時醫(yī)生肢體與患者間三維作用力,研究顯示三維的力學(xué)參數(shù)對單方向的力學(xué)參數(shù)存在明顯優(yōu)勢,且證明療效的可靠性可通過手法操作的動態(tài)標(biāo)準(zhǔn)化保證。李義凱[9]等應(yīng)用壓力傳感器,檢測旋轉(zhuǎn)手法作用過程中術(shù)者拇指頂推患者頸椎棘突的最大推力。Triano[10]提出腰椎“高速低幅沖擊手法(High Velocity Low Amplitude thrusting manipulation)”的力學(xué)參數(shù),分別對手法操作時峰值、預(yù)加載力、扳動速度、谷值、扳動時間等力學(xué)參數(shù)進(jìn)行了具體的闡釋,并在作力學(xué)曲線上標(biāo)示(圖1)。他的研究為很多研究者所公認(rèn)和應(yīng)用。

    圖1 Triano模擬旋轉(zhuǎn)手法力學(xué)曲線圖:1、預(yù)加載力(Preload):手法起始階段作用于軟組織的準(zhǔn)靜態(tài)載荷,2、谷值(Downward incisural point)沖擊發(fā)力前預(yù)加載力下降最低值,3、扳動速度(Thrust speed)沖擊相的斜率(△y/△x),4、峰值(Peak force)沖擊作用力升高最大值,5、扳動時間(Thrust duration)沖擊相開始至峰值的時間段。

    國內(nèi)的較深入的研究者有朱立國等[11]選擇神經(jīng)根型頸椎病患者施行旋轉(zhuǎn)手法,活體實測手法操作者力學(xué)參數(shù),測量左右手施行旋轉(zhuǎn)手法時預(yù)加載力、最大作用力、扳動力都是相近的,旋轉(zhuǎn)手法從緩慢上牽到扳動操作結(jié)束整個過程具有一定的規(guī)律性。王立恒[12]等對51名頸椎病患者施行旋轉(zhuǎn)手法,測量不同痛閾值下力學(xué)參數(shù)以及手法治療后的軟組織張力。Wherzog[13]運用固定于患者頸部治療區(qū)域的彈性壓力薄板,測量過程中作用力的變化。其最大作用力范圍從99 N-140 N。操作的作用時間為90-120 ms。房敏等[14]用電測技術(shù)觀察頸部拔伸法的應(yīng)力,測試不同施力點的拔伸力度,發(fā)現(xiàn)小角度前屈時手法安全,效果最佳。

    1.4腰椎的扳動手法活體測量

    力學(xué)效應(yīng)是手法行為的本質(zhì),所以很多人對腰椎的扳動手法進(jìn)行生物力學(xué)測量。Triano[15]指出手法給脊柱和它周圍組織施加外部負(fù)荷以影響脊柱及其周圍組織,但目前還不能預(yù)測哪種手法的優(yōu)劣,應(yīng)該將測量脊柱推拿生物力學(xué)參數(shù),應(yīng)用于手法療效定量研究。姜宏[16]認(rèn)為手法運動力學(xué)的動態(tài)變化對手法效應(yīng)途徑與作用環(huán)節(jié)起作用。陳守吉[17]指出手法的作用由手法刺激性、局部組織生物學(xué)特性、人體生理病理狀態(tài)三個因素構(gòu)成。徐海濤等[18]采用壓力傳感器方法主要研究力、力矩、時間以及力-時間曲線,通過壓力檢測系統(tǒng),分析手法影響因素,測量發(fā)現(xiàn)左右肩部扳力分別為12.552±1.715 kg和12.748±1.645 kg,為該手法的臨床操作提供了定量的、具體的力學(xué)參數(shù)。張軍[19]等采集和分析了腰椎手法時,常規(guī)推拿手法,掌根或者肘部的壓力信息。潘良春等[20]歸納正骨手法的力學(xué)特點為“軸向用力”、“旋轉(zhuǎn)用力”、“側(cè)方用力”,將旋轉(zhuǎn)用力與側(cè)方用力進(jìn)行力學(xué)分析,這項研究非常有助于對不同手法力學(xué)的差異進(jìn)行量化比較。Marcotte等[21]通過將微型壓力傳感器安裝在施術(shù)者手上,記錄推扳的壓力在仰臥定點旋轉(zhuǎn)手法過程中的變化,發(fā)現(xiàn)推扳力波動在4.0-41.0 N。

    脊柱手法主要效應(yīng)點在于腰椎關(guān)節(jié)突關(guān)節(jié),很多研究者非常關(guān)注角度和姿勢變化對關(guān)節(jié)突關(guān)節(jié)的影響,手法施加在關(guān)節(jié)突關(guān)節(jié)的載荷也與手法的療效和副作用密切相關(guān)。這些研究者主要運用攝影技術(shù)方法研究手法或位移傳感器的生物力學(xué)分析,觀察力學(xué)指標(biāo)主要是速度、時間、角度以及位移。侯筱魁[22]認(rèn)為關(guān)節(jié)突關(guān)節(jié)錯縫可以通過腰椎旋轉(zhuǎn)手法得到調(diào)整。但有人認(rèn)為手法過重傷及關(guān)節(jié)突關(guān)節(jié)和關(guān)節(jié)囊使一些患者經(jīng)腰椎手法后,早期出現(xiàn)腰痛癥狀加重等現(xiàn)象。

    2 數(shù)字化人體模型研究

    隨著計算機數(shù)字時代的到來,虛擬現(xiàn)實與仿真技術(shù)也突顯出其重要的作用,人體建模仿真也日漸成為人機工程、計算機圖形學(xué)、生物力學(xué)、交通事故分析、康復(fù)工程等多領(lǐng)域的研究熱點。用于損傷與防護(hù)研究的數(shù)字化人體模型的發(fā)展經(jīng)歷了集總參數(shù)模型、多剛體動力學(xué)模型、多體動力學(xué)模型、有限元模型,目前正在探討在虛擬人的基礎(chǔ)上建立動力學(xué)生理模型。

    2.1腰椎三維有限元分析

    2.1.1腰椎三維有限元模型(Finite element models)關(guān)節(jié)突關(guān)節(jié)力學(xué)特征分析

    腰椎三維有限元模型對于分析脊柱損傷及其功能紊亂的潛在機制,具有重大意義。葉淦湖[23]對腰椎關(guān)節(jié)突關(guān)節(jié)的模型進(jìn)行各種手法的力學(xué)分解模擬,發(fā)現(xiàn)腰椎關(guān)節(jié)突關(guān)節(jié)在旋轉(zhuǎn)推拿中具有主要的抵抗作用。關(guān)節(jié)在抗旋轉(zhuǎn)運動時,主要通過上、下關(guān)節(jié)突關(guān)節(jié)面的直接接觸,而關(guān)節(jié)囊可防止兩個關(guān)節(jié)突關(guān)節(jié)面分開過大。他把關(guān)節(jié)突關(guān)節(jié)處理為接觸模型,正確分析了關(guān)節(jié)突關(guān)節(jié)的力學(xué)功能。由于腰椎關(guān)節(jié)突關(guān)節(jié)面在下腰椎與矢狀面呈45°的角,這樣使其旋轉(zhuǎn)的余地很小,對于正常腰椎來說,在有壓縮載荷的條件下,椎間盤旋轉(zhuǎn)時關(guān)節(jié)突關(guān)節(jié)在旋轉(zhuǎn)1-2°時主要抵抗力為旋轉(zhuǎn)力。

    2.1.2三維有限元模型對手法的療效及安全性評價

    UenoK實驗發(fā)現(xiàn)[24]關(guān)節(jié)軟骨表層骨骨折的應(yīng)力范圍為14-59 MPa,平均為35.7 MPa。而尸體實驗和有限元模型證明[25]前屈側(cè)彎旋轉(zhuǎn)手法使腰椎關(guān)節(jié)突關(guān)節(jié)的活動幅度最大,直立旋轉(zhuǎn)手法次之,關(guān)節(jié)突關(guān)節(jié)最大應(yīng)力均小于14 MPa,腰椎關(guān)節(jié)突關(guān)節(jié)的應(yīng)力為21.98 N,不會造成關(guān)節(jié)突關(guān)節(jié)骨折。向左旋轉(zhuǎn)時,右側(cè)關(guān)節(jié)突關(guān)節(jié)間隙大,向右側(cè)旋轉(zhuǎn)時,左側(cè)關(guān)節(jié)突關(guān)節(jié)間隙增大。當(dāng)前屈最大時,關(guān)節(jié)突關(guān)節(jié)的上、下關(guān)節(jié)突可相對滑動5-7 mm,由于關(guān)節(jié)突關(guān)節(jié)突的反復(fù)活動,關(guān)節(jié)囊的伸縮,可松動神經(jīng)根周圍的粘連[26]。

    從整個脊椎單元來看,除了手法使旋轉(zhuǎn)側(cè)椎間盤突出程度加重外,旋轉(zhuǎn)力亦不大可能使椎間盤纖維環(huán)破裂,從這一點考慮,手法是安全的,除去椎間盤突出的位移,兩側(cè)椎間盤后外側(cè)角位移基本一致,而旋轉(zhuǎn)對側(cè)的椎間孔在手法過程中,由于上位關(guān)節(jié)突關(guān)節(jié)的上移而擴大了椎間孔。畢勝[26]用有限元模型模擬手法作用時,發(fā)現(xiàn)椎間盤的髓核內(nèi)壓力牽扳手法最小且為負(fù)值,坐位旋轉(zhuǎn)手法時最大。坐位旋轉(zhuǎn)手法和斜扳手法外層纖維環(huán)應(yīng)力最大,內(nèi)層纖維環(huán)應(yīng)力最大的是牽扳手法,手法均可造成椎間盤的后外側(cè)與神經(jīng)根之間移動,并且有限元方法與生物力學(xué)研究結(jié)果基本一致。他還對腰椎關(guān)節(jié)突關(guān)節(jié)的模型進(jìn)行各種手法的分解模擬,發(fā)現(xiàn)腰椎關(guān)節(jié)突關(guān)節(jié)在旋轉(zhuǎn)推拿中具有主要的抵抗作用。

    2.1.3三維有限元模型的局限性

    目前有限元模型分為兩類:靜力學(xué)和動力學(xué)研究模型。用于靜力學(xué)分析的模型更加注重椎體的幾何形態(tài),模型往往只包含一兩個運動節(jié)段。三維有限元的研究主要集中在骨折和椎體變形分析方面。在分析肌肉、韌帶等軟組織對脊柱系統(tǒng)的整體穩(wěn)定、載荷分布、退變等因素的影響時,存在一定偏差。亦即,單純采用有限元法分析脊柱局部的線性或非線性的變形尚不足以獲得所分析單元的整體變形、載荷分布和局部位移特征。

    2.2人體多剛體模型

    多剛體模型能夠避免人體內(nèi)部的復(fù)雜調(diào)控機制,簡單而較為真實的反應(yīng)出人體運動的基本機理,具有廣泛的應(yīng)用領(lǐng)域。著名的Hanavan模型[27]將人體分為15個獨立的剛體,各個剛體之間通過球形鉸鏈連接在一起,建立力學(xué)方程組,通過求解得到運動模型參數(shù)。美國國家空軍實驗室[28]早在1975年就建立了三維多剛體假人ABT(Articulated Total Body)得到了廣泛的應(yīng)用。荷蘭應(yīng)用科學(xué)研究院(TNO)[29]采用多剛體系統(tǒng)動力學(xué)的方法構(gòu)建了人體運動系統(tǒng),使用橢球面模擬人體幾何外形,建立了Hannavan多剛體假人,被廣泛應(yīng)用于包括汽車碰撞在內(nèi)的眾多領(lǐng)域[30]。

    2.3神經(jīng)-肌肉-骨骼綜合模型

    人體運動系統(tǒng)是一個非常復(fù)雜的綜合非線性系統(tǒng),需要神經(jīng)、肌肉和骨骼的共同調(diào)控和相互協(xié)調(diào)。南非的Hatze[31]構(gòu)建了一個2段5塊肌肉3自由度的腿平面模型解決了踢腳動作的動力學(xué)最優(yōu)化問題。Naomichi[32]等提出了更符合人體解剖學(xué)的肌肉-骨骼-神經(jīng)系統(tǒng)的集成模型,考慮了中樞神經(jīng)系統(tǒng)運動神經(jīng)元的控制,然而這種建模方法十分復(fù)雜,同時人們對人體神經(jīng)系統(tǒng)對肌肉的具體機理認(rèn)識還不十分清楚,模型建立的難度很大,模型的計算也需要適當(dāng)?shù)暮喕?/p>

    2.4頸椎多柔體動力學(xué)模型

    國外科研工作者開展的頸椎多柔體動力學(xué)模型主要針對汽車碰撞問題。Tien[33]等開發(fā)了人體頭頸部系統(tǒng)數(shù)值模擬的算法和程序,然而無法模擬肌肉等軟組織的力學(xué)行為。Merrill[34]等建立了包括2個胸椎、7個頸椎與頭部的三維動力學(xué)模型,并使用無質(zhì)量的彈簧和阻尼連接這些剛體,但該模型是二維的,肌肉模型也存在不足。Deng等[35]在Merrill基礎(chǔ)上,將15組主要頭頸部肌肉從參數(shù)模型里分離出來,使用3點肌肉模型模擬繞過頸椎時形成的彎曲,使肌肉模型更仿真。Jager[36,37]等建立了1、不單獨考慮肌肉的全局模型2、單獨考慮肌肉作用的細(xì)化模型,從而使仿真度有了部分改善。

    國內(nèi)多體動力學(xué)仿真方面的研究相對較少,張竟等[38]使用多體動力學(xué)軟件建立了簡單的人體頭頸部多剛體動力學(xué)模型,但該模型僅考慮了4組主要肌肉,對肌肉力學(xué)行為的模擬不足。柳松楊等[39]對艦載機著陸過程中飛行員頭頸部的動力學(xué)響應(yīng)進(jìn)行了仿真試驗,但由于肌肉模型的限制,模型與實驗結(jié)果還存在一定差距。魯廷等[40]建立包括頭部、7個頸椎和2個胸椎共10個剛體的頭頸部多體動力學(xué)模型;采用集總參數(shù)法描述韌帶、椎間盤等軟組織的力學(xué)特性;采用最新的非線性應(yīng)力-應(yīng)變關(guān)系方法來描述人體頭頸部15組肌肉的力學(xué)特性,并用多組汽車碰撞過程動力學(xué)響應(yīng)實驗數(shù)據(jù)對模型進(jìn)行驗證。

    2.5腰椎多體動力學(xué)模型的建立

    數(shù)字模擬仿真使得我們能夠?qū)θ梭w的復(fù)雜結(jié)構(gòu)和功能進(jìn)行簡化,從而避免了對人體內(nèi)部多樣運動的關(guān)注,而對主要影響因素和主要問題進(jìn)行分析。其次,仿真模型突破了人體自身的局限性和測量儀器的局限性,使我們能夠獲得那些無法直接測量或者很難有效測量的參數(shù)、指標(biāo),便于我們對人體損傷機理及其防護(hù)的進(jìn)一步研究。目前多剛體動力學(xué)仿真軟件在車輛碰撞分析和沖擊生物力學(xué)研究領(lǐng)域得到了廣泛的使用。多體系統(tǒng)模型適合于人體動力學(xué)響應(yīng)的模擬。專門針對軟組織研究的多柔體動力學(xué)方法不僅具備三維有限元方法的變形和應(yīng)力分析功能,而且可以更加方便地處理骨骼體與肌肉、韌帶、間盤等軟組織的相互作用、變形約束與協(xié)調(diào)以及相對運動。多體動力學(xué)模型可以被應(yīng)用于人體動力學(xué)的各個方面,可以通過動力學(xué)仿真的方法獲得人體的各個部分在運動中的狀態(tài)以及相互作用或載荷分布的狀態(tài),完全可以勝任脊柱力學(xué)分析任務(wù)。能夠獲得探究各種載荷下脊柱及其周圍軟組織在各種力學(xué)條件下的受力狀態(tài)、失衡狀態(tài)等??傊?,建立脊柱系統(tǒng)的多體動力學(xué)模型,就可以根據(jù)動力學(xué)仿真獲其受力和失衡狀態(tài)。

    脊柱的力學(xué)模型在生物力學(xué)研究領(lǐng)域中已經(jīng)有一段發(fā)展過程,學(xué)者利用脊柱力學(xué)模型[41,42]研究了脊柱載荷、脊柱的病變、動力學(xué)損傷、脊柱穩(wěn)定性等等。有限元方法和多剛體模型被廣泛用來建立脊柱力學(xué)模型。有限元方法可以建立較為細(xì)致的分析變形-應(yīng)力關(guān)系的力學(xué)模型,但不能很好地描述脊柱運動學(xué)和復(fù)雜連接關(guān)系。隨著計算力學(xué)和計算機運算能力的發(fā)展,多體動力學(xué)已經(jīng)由初期的多剛體動力學(xué)模型發(fā)展到了目前的多剛體-柔性體耦合的動力學(xué)模型。因此可以說多柔性體動力學(xué)融合了有限元法和多剛體動力學(xué)的優(yōu)勢。在生物多體動力學(xué)方面,國外較為著名的軟件有OpenSim[43]和Anybody[44]兩款軟件。UC伯克利大學(xué)的Miguel Christophy[45]在OpenSim軟件中建立了包含238條肌肉的脊柱模型,但是沒有考慮關(guān)節(jié)突、韌帶等組織的影響。國內(nèi)上海交通大學(xué)的聶文忠[46]利用有限元建立了一個脊柱胸腰部生物力學(xué)模型,但是其肌肉是以力邊界的形式給出。目前國內(nèi)外研究中,鮮有包含較為完整軟組織的脊柱模型。

    郭偉[47]使用清華大學(xué)航天航空學(xué)院開發(fā)的多柔體動力學(xué)分析軟件,能分析非線性的多個運動體間相互關(guān)系,而且在彈性體的模型構(gòu)建上非常有優(yōu)勢。目前已經(jīng)研制一套以絕對節(jié)點坐標(biāo)建模方法為理論基礎(chǔ)的柔性體動力學(xué)求解器-THUSolver。該方法能夠進(jìn)行各種條件下力學(xué)分析(包括失重和過載狀態(tài)),已經(jīng)用于飛機各種飛行狀態(tài)下動力學(xué)分析。

    2.5.1骨骼模型建立

    建立骨骼模型的動力學(xué)模型,需要其骨骼相應(yīng)的三維幾何、質(zhì)量信息。隨著計算機技術(shù)的發(fā)展,可以通過人體兩維斷層影像來重構(gòu)人體骨骼的三維數(shù)字模型,如Minics軟件。OpenSim的源代碼[48]是開源的,用戶可以使用其建立相關(guān)的人體肌骨模型,也可以使用其相關(guān)的文件、信息在其他軟件中建立模型。

    2.5.2關(guān)節(jié)突接觸

    椎骨和椎骨之間,通過上下關(guān)節(jié)突的接觸約束了腰部的扭轉(zhuǎn)自由度。關(guān)節(jié)突之間的接觸應(yīng)力是關(guān)心的內(nèi)容,其接觸力可以使用赫茲接觸模型[49]。

    2.5.3椎間盤模型

    椎間盤是一個被動結(jié)構(gòu),在維持脊柱的穩(wěn)定性方面起到重要作用。椎間盤的力學(xué)模型是脊柱力學(xué)分析中的一個核心內(nèi)容。一般認(rèn)為,椎間盤對椎骨的反作用力與椎骨間的相對位移有關(guān)系[50][51]。Panjabi的線性剛度陣模型表現(xiàn)出一種類似于梁的性質(zhì)。Gardner-Morse[52]提出了用剪切梁模型來代替K矩陣。取剪切梁的椎間盤模型的彈性模量為E=1E7,泊松比=0.49。擬合了在軸壓為150 N的幾何數(shù)據(jù)[11]。

    2.5.4韌帶建模

    不同的韌帶,其力學(xué)性質(zhì)有著明顯的不同。韌帶的應(yīng)變-應(yīng)力曲線通常是非線性的。下圖是黃韌帶的力學(xué)性質(zhì)曲線。

    圖2黃韌帶應(yīng)變-力關(guān)系

    2.5.5肌肉建模

    人體肌肉結(jié)構(gòu)及其力學(xué)特性十分復(fù)雜。目前大部分生物力學(xué)研究中,骨骼肌模型通常采用直線或曲線來構(gòu)建,其中直線或曲線稱為肌肉路徑[55-57]。肌腱的力學(xué)行為則等效為非線性的彈簧阻尼器。

    3 目前研究中的問題

    當(dāng)前有關(guān)脊柱生物力學(xué)實驗的研究主要是圍繞椎體附件進(jìn)行研究,尤其是骨折、植入物的垂直剪切力關(guān)注較多,對肌肉、韌帶、關(guān)節(jié)囊等椎旁軟組織關(guān)注少。局部(如腰椎4-5、腰5-骶骨,椎弓根的應(yīng)力應(yīng)變)關(guān)注多,而對整體(脊柱-骨盆-肌肉-韌帶之間的關(guān)系)關(guān)注較少。更未發(fā)現(xiàn)脊柱失衡條件下不同手法對于脊柱軟組織、椎間盤和關(guān)節(jié)突關(guān)節(jié)影響的實驗研究。脊柱椎旁肌和韌帶組織的實驗研究之所以比較少,除了基礎(chǔ)實驗的研究者關(guān)注點不同,還可能緣于應(yīng)力指標(biāo)的選擇比較局限。但這一點恰好是臨床上脊柱損傷病理過程及脊柱手法安全實施的關(guān)鍵環(huán)節(jié)。

    4 展望

    總之,當(dāng)前常見的脊柱力學(xué)模型的建立雖然為脊柱創(chuàng)傷帶來了一些研究成果,甚至在如何規(guī)范、量化手法操作方面也提出了一些原則性看法,但很難具備指導(dǎo)臨床應(yīng)用的作用。主要原因就是既往脊柱建模的目標(biāo)仍然主要局限于椎體和間盤的創(chuàng)傷,相關(guān)肌肉和韌帶組織的損傷應(yīng)力指標(biāo)研究并不精確。而這一點恰好是脊柱損傷及脊柱手法安全實施的關(guān)鍵環(huán)節(jié)。我們知道,關(guān)節(jié)穩(wěn)定結(jié)構(gòu)——韌帶與肌肉等軟組織的力學(xué)失衡是導(dǎo)致最終脊柱結(jié)構(gòu)力學(xué)失衡(骨折或間盤突出)的關(guān)鍵。所以,我們不僅要重視骨折、植入物的固定等剛體的變化,還應(yīng)該特別重視軟組織在脊柱運動應(yīng)力損傷中的臨床意義。多柔體動力學(xué)分析法與計算機可視化技術(shù)的出現(xiàn)正是此領(lǐng)域所迫切需要的技術(shù)。所以,有關(guān)脊柱手法的生物力學(xué)研究重點應(yīng)該是1、各種脊柱力學(xué)失衡狀態(tài)下和實施手法時的椎旁肌及韌帶組織的應(yīng)力測定,2、基于以上研究的基礎(chǔ)結(jié)合數(shù)字化人體仿真技術(shù)構(gòu)建包括腰椎主要肌肉和韌帶的多柔體動力學(xué)模型的建立。以上兩方面的研究不僅豐富和完善關(guān)于脊柱運動應(yīng)力損傷的機制研究,還對確立脊柱手法實施的安全原則具有重要的臨床意義,是脊柱生物力學(xué)研究未來發(fā)展的方向。

    [1]Tomas F.Chiropractic Technique:Principles and Procedures[M].Elservier(Singapore),2013.87-89.

    [2]MaekoVW,MatthewsLW,ZwirkoskiP,et al.Thejoint ontactarea of the ankle:the contribution of the posterior malleolus[J].Bone Joint Surg,1997,3:347-351.

    [3]John J Triano,Marion McGregor,Dennis R Skogsbergh.Use of chiropractic manipulation in lumbar rehabilitation[J].Journal of Rehabilitation Research and Development,1997,34(4):394-404.

    [4]鄧玫,劉紅萍,張曉剛,等.五種按法最大壓強的比較研究[J].按摩與導(dǎo)引,2005,21(7):14-16.

    [5]HessellBw,Herzogw,eonwayPJ,et al.Experiment measurement of the fore during spinal manipulation using the Thompson technique[J].Journal of manipulative and physiological therapeutics,1990,13(8):448-453.

    [6]HerzogP.J,Conway,G.N,Kawehuw,et al.Forces exerted during spinal elative therapy[J].Spine,1993,18(9):1207-1208.

    [7]JimMW,Daniel Chow,Andrew D.The kinematics and intra-and inter-the rapisteonsisteneies of lower cervical rotational manipulation[J].Medical Engineeringffamp;Physies,2005,27:396-400.

    [8]Van zoest GG,GosselinG.Three-dimension analyze forces in chiropractic spinal manipulative therapy[J].Journal of Manipulative and Physiological Therapy,2003,26(9):549-556.

    [9]李義凱,徐海濤,王國林,等.頸椎定點旋提手法所致咔噠聲響與最大扳動力的量效關(guān)系研究[J].中國康復(fù)醫(yī)學(xué)雜志,2004,19(9):644-646.

    [10]TrianoJJ,Carol M.Rogers,SarahCombs,et al.Quantitative feedback versus standard training for cervical and thoracic manipulation[J].Manipulative Physical Ther,2003,26(3):131.

    [11]朱立國,馮敏山,畢方杉,等.頸椎旋轉(zhuǎn)(提)手法的活體力學(xué)測量[J].中國康復(fù)醫(yī)學(xué)雜志,2007,22(8):673-674.

    [12]王立恒,孫樹椿,陳朝暉,等.手法治療頸椎病的在體力學(xué)測量[J].中國中醫(yī)骨傷科雜志,2009,17(3):22-24.

    [13]Wherzog.Force secreted during Spinalmanipulative therapy[J].spine,1993,18(9):1206-1212.

    [14]房敏,嚴(yán)雋陶,沈國權(quán),等.頸部推拿拔伸手法的 在體 研究[J].頸腰痛雜志,2000,21(3):202-204.

    [15]TrianoJ.Biomechanics of spinal manipulative[J].The spinal journal,2001,l:121-131.

    [16]姜宏,施祀.頸椎手法的生物力學(xué)研究與探索[J].中國中醫(yī)骨傷科雜志,2014(2):51-53.

    [17]陳守吉,許世雄,史一蓬,等.中醫(yī)推拿擺動類手法的動力學(xué)研究(工)生物力學(xué)模型及方程[J].醫(yī)用生物力學(xué),1996,2(2):112.

    [18]徐海濤,徐達(dá)傳,李義凱,等.腰椎斜扳手法所致“咔嗒”聲時推扳力的研究[J].中國康復(fù)醫(yī)學(xué)雜志,2005,23(3):202-204.

    [19]張軍,韓磊,胡磊,等.推拿手法培訓(xùn)系統(tǒng)的研制及應(yīng)用[J].中國中醫(yī)骨傷科雜志,2009,17(2):60.

    [20]潘良春,周太安,周奉皋,等.中醫(yī)正骨手法中的力學(xué)機制[J].海軍醫(yī)學(xué)雜志2009,28(l):694-696.

    [21]Marcotte,Normand.Measurement of the Pressure applied during motion palpation and reliability for cervical spine rotation[J].Manipulative Physiol ther,2005,28(8):591-600.

    [22]侯筱魁,董凡,趙文成.斜扳時完整腰椎三維立體運動的研究[J].中國骨傷,1996,9(4):5-7.

    [23]葉淦湖,張美超,李義凱.模擬推拿時腰椎關(guān)節(jié)突關(guān)節(jié)有限元模型的生物力學(xué)分析[J].廣州中醫(yī)藥大學(xué)學(xué)報,2003,20(3):195-200.

    [24]UenoK,LiuYK.A3-D nonlinear finite element model of lumbar intervertebral joint intorsion[J].Journal of Biomechanical Engineering,1987,109:200-209.

    [25]高允海.有限元分析法研究脊柱生物力學(xué)的新進(jìn)展[J].國外醫(yī)學(xué)-生物醫(yī)學(xué)工程分冊,2003,26(6):281-284.

    [26]畢勝,李義凱,趙衛(wèi)東,等.腰部推拿手法生物力學(xué)和有限元比較研究[J].中華物理醫(yī)學(xué)與康復(fù)雜志,2002,24(9):525-527.

    [27]Wojciech Blajera,Adam Czaplickib.Modeling and inverse simulation of somersaults on the trampoline[J].Journal of Biomechanics,2001,3(4):68-69.

    [28]Zhou Xiaodong,Louis F Draganich,F(xiàn)arid Amirouche.A dynamic model for simulating a trip and fall during gait[J].Medical Engineering and Physics,2002,24(4):29-30.

    [29]陳健、姚頌平.虛擬現(xiàn)實技術(shù)活體育運動技術(shù)仿真中的應(yīng)用[J].體育科學(xué),2006,26(9):34-39.

    [30]張曉云,金先龍,申杰.基于人體損傷特征的車—人碰撞事故再現(xiàn)模擬與分析[J].機械工程學(xué)報,2010,46(12):78.

    [31]Hatze,Gentaro Taga.A model of neuro-muscule-skeletal system for anticipatory adjustment of human locomotion during obstacle avoidance[J].Biological Cybernetics,1988,78:89-90.

    [32]Naomichi Ogihara Naomichi Ogihara,Nobutoshi Yamazaki.Generation ofhuman bipedal locomotion by a bio-mimetic neuro-musculo-skeletal system[J].Biological Cybernetics,2001,84:90-91.

    [33]Tien CS,Huston RL.Biodynamic modeling of the head/neck system.In:Field accidents data collections,analysis,methodologiesand crash injury reconstructions[J].Warrendale,USA:Society of Automotive Engineers,1995,359-364.

    [34]Merrill T,Goldsmith W,Deng YC.Three dimensional response of a lumped parameter head neck model due to impact and impulsive loading[J].JBiomech,1984,17(2):81-85.

    [35]Deng YC,Goldsmith W.Response of a human head/neck/upper-torso replica to dynamic loading-Ⅱanalytical/numerical model[J].J Biomech,1987,20(5):487-497.

    [36]Jager M,Sauren A,Thunnissen J,et al.A 3D head-neck model:Validation for frontal and lateral impacts[C].Proceedings of the 38th Stapp Car Crash Conference.Fort Lauderdale,USA,1994.9,46-47.

    [37]Jager M,Sauren A,Thunnissen J.A global and a detailed mathematicalmodel for head-deck dynamics[C].Proceedingsof the40th Stapp CarCrash Conference.Albuquerque,USA,1996.9,89-92.

    [38]張竟,韓旭,文桂林.基于ADAMS LifeMOD的人體頭頸部動力學(xué)仿真與驗證[J].系統(tǒng)仿真學(xué)報,2008,20(10):2718-2721.

    [39]柳松楊,叢紅,王鶴,等.軍機飛行員的頸部損傷研究[J].醫(yī)用生物力學(xué),2010,25(4):262-265.

    [40]魯廷,王亞偉,柳松楊.攔阻著艦過程中飛行員頭頸部的動力學(xué)響應(yīng)[J].醫(yī)用生物力學(xué),2012,27(6):642-646.

    [41]N.Bogduk,J.E.Macintosh,and M.J.Pearcy.A universal model of the lumbar back muscles in the upright position[J].Spine,1992,17(8):897-913.

    [42]Aubin,C.E,Petit,Y,StokesI.A.F,Poulin F,Gardner-MorseMG,LabelleH.Biomechanical modelingof posterior instrumentation of the scoliosisspine[J].ComputerMethods in Biomechanics and BiomedicalEngineering,2003,6:27-32.

    [43]Delp S.L,et al.OpenSim:Open-Source software to create and analyze dynamic simulations of movement[J].Biomedical Engineering,2012,54(11):1940-1950.

    [44]Michael Damsgaarda,John Rasmussena,Christensena,Egidijus Surmaa,Mark de Zee.Analysis of musculoskeletal systems in the anybody modeling system[J].Simulation Modelling Practice and Theory,2006,14(8):1100-1111.

    [45]Miguel Christophy,Nur Adila Faruk Senan,Jeffrey C.A musculoskeletal model for the lumbar spine[J].Biomech Model Mechanobiol,2012,11:19-34.

    [46]聶文忠.脊柱胸腰部的生物力學(xué)建模與應(yīng)用研究[D].北京:上海交通大學(xué)機械與動力學(xué)院,2009.38-39.

    [47]郭偉,陳立君,王飛.腰椎三維有限元模型的建立及高載荷(Gz)環(huán)境的受力分析[J].空軍醫(yī)學(xué)雜志,2013,29(2):74-75.

    [48]http://opensim.stanford.education.

    [49]虞磊.基于絕對節(jié)點坐標(biāo)法的柔性多體系統(tǒng)建模與計算方法研究[D].北京:清華大學(xué)航天航空學(xué)院,2010.39-40.

    [50]PanjabiM.M.,BrandR.A.,White IIIA.A.Three-dimensionalflexibility and stiffness properties of the human thoracic spine[J].Journal of Biomechanics,1976,9:185-192.

    [51]Miguel Christophy,Maurice Curtin et al.On the modeling of the intervertebral joint in multibody models for the spine[J].Multibody Syst Dyn,2013,30(4):413-432.

    [52]Mack G.Gardner-Morse,Ian A.F.Stokes.Structural behavior of human lumbar spinal motion segments[J].Journal of Biomechanics,2004,37:205-212.

    [53]Alf L.Nachemson.Some mechanical propert ies of the third human lumbar interlaminar ligament(ligamentum flavum)[J].Journal of Biomechanics,1968,1(3):211-214.

    [54]Eberlein,Holzapfel,F(xiàn)rohlich.Multi-segmentFEA of the human lumbar spine including the heterogeneity of the annulus fibroses[J].Computational Mechanics,2004,34:147-163.

    [55]Myer Kutz.Biomedical Engineering and Design Handbook Volume 1:Fundamentals(Second Edition)[M].McGraw-Hill Professional,2009.164-169.

    [56]TangG,WangC.A muscle-path-plane method for representing muscle contraction during joint movement[J].Computer Methods in Biomechanics and Biomedical Engineering,2010,13(6):723-729.

    [57]GarnerB,PandyM.The obstacle-set method for representing muscle paths in musculoskeletalmodels[J].Computermethods in biomechanics and biomedical engineering,2000,3(1):1-30.

    猜你喜歡
    有限元模型研究
    一半模型
    FMS與YBT相關(guān)性的實證研究
    遼代千人邑研究述論
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    視錯覺在平面設(shè)計中的應(yīng)用與研究
    科技傳播(2019年22期)2020-01-14 03:06:54
    EMA伺服控制系統(tǒng)研究
    3D打印中的模型分割與打包
    磨削淬硬殘余應(yīng)力的有限元分析
    基于SolidWorks的吸嘴支撐臂有限元分析
    22中文网久久字幕| 99热全是精品| av专区在线播放| 99视频精品全部免费 在线| 精品久久久久久久久亚洲| 久久久久久国产a免费观看| 男人的好看免费观看在线视频| 国产片特级美女逼逼视频| 成人漫画全彩无遮挡| 精品久久久久久久久久久久久| 欧洲精品卡2卡3卡4卡5卡区| 99在线视频只有这里精品首页| 日本撒尿小便嘘嘘汇集6| 欧美激情在线99| 亚洲四区av| 国产精品久久久久久精品电影| 久久精品国产鲁丝片午夜精品| 深夜a级毛片| 国产大屁股一区二区在线视频| 亚洲av第一区精品v没综合| 精品久久久久久久久av| 国产精品人妻久久久影院| 国产精品人妻久久久影院| 欧美一区二区国产精品久久精品| 99在线人妻在线中文字幕| 日日摸夜夜添夜夜爱| 麻豆乱淫一区二区| 搡老熟女国产l中国老女人| 嫩草影视91久久| 少妇的逼好多水| 国产精华一区二区三区| 国产在视频线在精品| 精品福利观看| 久久精品夜夜夜夜夜久久蜜豆| 最近中文字幕高清免费大全6| 夜夜夜夜夜久久久久| 中文字幕av在线有码专区| 在线天堂最新版资源| 欧美一区二区精品小视频在线| 黄色一级大片看看| 欧美在线一区亚洲| 寂寞人妻少妇视频99o| 校园春色视频在线观看| 午夜亚洲福利在线播放| 国产av在哪里看| 人妻夜夜爽99麻豆av| 日韩av不卡免费在线播放| 97人妻精品一区二区三区麻豆| 午夜a级毛片| 色吧在线观看| 最近视频中文字幕2019在线8| 午夜亚洲福利在线播放| 直男gayav资源| .国产精品久久| 国产伦精品一区二区三区视频9| а√天堂www在线а√下载| 男人狂女人下面高潮的视频| 国内精品一区二区在线观看| 男女下面进入的视频免费午夜| 两个人视频免费观看高清| 欧美成人精品欧美一级黄| 国产精品99久久久久久久久| 国内揄拍国产精品人妻在线| 亚洲精品国产av成人精品 | www日本黄色视频网| 最近中文字幕高清免费大全6| 99久国产av精品| 搞女人的毛片| 美女大奶头视频| 搡老妇女老女人老熟妇| 在线播放国产精品三级| 精品一区二区三区人妻视频| 中文字幕熟女人妻在线| 最新在线观看一区二区三区| 大香蕉久久网| 国产成人91sexporn| 亚洲精品国产av成人精品 | 99九九线精品视频在线观看视频| 午夜福利成人在线免费观看| 在线天堂最新版资源| 亚洲精品色激情综合| 1024手机看黄色片| 在线播放无遮挡| 国产av在哪里看| 午夜亚洲福利在线播放| 亚洲人成网站在线播| 精品福利观看| 国产国拍精品亚洲av在线观看| 别揉我奶头~嗯~啊~动态视频| 亚洲丝袜综合中文字幕| 亚洲美女黄片视频| 免费搜索国产男女视频| 色尼玛亚洲综合影院| 在线观看免费视频日本深夜| 国产一区二区激情短视频| 干丝袜人妻中文字幕| 国产精品伦人一区二区| 黄色配什么色好看| 日日撸夜夜添| av女优亚洲男人天堂| 真人做人爱边吃奶动态| 18禁黄网站禁片免费观看直播| 亚洲av熟女| 成年av动漫网址| 在线a可以看的网站| 免费电影在线观看免费观看| 能在线免费观看的黄片| 伊人久久精品亚洲午夜| 色吧在线观看| 亚洲成人久久爱视频| 18+在线观看网站| 国产午夜精品久久久久久一区二区三区 | 国产精品久久视频播放| 啦啦啦韩国在线观看视频| 国产精品综合久久久久久久免费| 搡老妇女老女人老熟妇| 亚洲不卡免费看| 全区人妻精品视频| 97在线视频观看| 99热这里只有是精品在线观看| 日韩,欧美,国产一区二区三区 | 亚洲成人精品中文字幕电影| 午夜老司机福利剧场| 亚洲欧美日韩高清在线视频| 国产一区二区三区av在线 | 一卡2卡三卡四卡精品乱码亚洲| 黄色日韩在线| 久久久久久久午夜电影| 免费黄网站久久成人精品| 小说图片视频综合网站| 久久久精品94久久精品| 国产久久久一区二区三区| 永久网站在线| 日韩成人av中文字幕在线观看 | 看免费成人av毛片| 人妻夜夜爽99麻豆av| 国产精品免费一区二区三区在线| 成人av在线播放网站| 蜜桃久久精品国产亚洲av| 人人妻人人澡人人爽人人夜夜 | 亚洲精品乱码久久久v下载方式| 国产亚洲精品久久久久久毛片| 天堂√8在线中文| 免费在线观看成人毛片| 成人高潮视频无遮挡免费网站| 此物有八面人人有两片| 国产一区二区三区在线臀色熟女| 国产真实乱freesex| 久久久a久久爽久久v久久| 成人欧美大片| 成人亚洲欧美一区二区av| av中文乱码字幕在线| 亚洲av免费在线观看| 欧美成人一区二区免费高清观看| 日本欧美国产在线视频| 欧美成人精品欧美一级黄| 天堂av国产一区二区熟女人妻| 亚洲av五月六月丁香网| 色在线成人网| 亚洲内射少妇av| 日韩亚洲欧美综合| 久久精品国产亚洲av天美| 色综合色国产| a级毛片a级免费在线| 中出人妻视频一区二区| 两个人的视频大全免费| 亚洲aⅴ乱码一区二区在线播放| 老师上课跳d突然被开到最大视频| 人妻夜夜爽99麻豆av| 深爱激情五月婷婷| 乱码一卡2卡4卡精品| 欧美一区二区精品小视频在线| 天美传媒精品一区二区| 中文字幕久久专区| 日本色播在线视频| 18禁裸乳无遮挡免费网站照片| 国产精品野战在线观看| 男女做爰动态图高潮gif福利片| 免费观看的影片在线观看| 91在线观看av| 国产男人的电影天堂91| 精品日产1卡2卡| 寂寞人妻少妇视频99o| 久久综合国产亚洲精品| 乱码一卡2卡4卡精品| 日日啪夜夜撸| 国产高清有码在线观看视频| 午夜激情福利司机影院| 精品久久久久久久久久免费视频| 九九在线视频观看精品| 日本五十路高清| or卡值多少钱| 亚洲四区av| 色视频www国产| 精品一区二区三区av网在线观看| 久久久久国产网址| 久久久午夜欧美精品| 五月玫瑰六月丁香| 校园人妻丝袜中文字幕| 最近在线观看免费完整版| 国产精品一区二区免费欧美| 国产高清不卡午夜福利| 亚洲av.av天堂| 久久6这里有精品| 亚洲成a人片在线一区二区| 免费人成在线观看视频色| 日日啪夜夜撸| 一级毛片久久久久久久久女| 欧美另类亚洲清纯唯美| 国产爱豆传媒在线观看| 男女视频在线观看网站免费| 尾随美女入室| 床上黄色一级片| 人妻久久中文字幕网| 亚洲成人久久性| 亚洲av第一区精品v没综合| 午夜精品一区二区三区免费看| 久久人人精品亚洲av| 一级毛片我不卡| 国产伦一二天堂av在线观看| 久久精品综合一区二区三区| 欧美国产日韩亚洲一区| 日韩精品青青久久久久久| 午夜福利在线观看吧| 在线国产一区二区在线| 国产亚洲91精品色在线| 国产大屁股一区二区在线视频| 日产精品乱码卡一卡2卡三| 18禁裸乳无遮挡免费网站照片| av在线亚洲专区| 亚洲自偷自拍三级| 成人亚洲精品av一区二区| 久久久久久九九精品二区国产| 男女视频在线观看网站免费| 免费看光身美女| 有码 亚洲区| 色尼玛亚洲综合影院| 99国产极品粉嫩在线观看| 男插女下体视频免费在线播放| av在线亚洲专区| 波多野结衣高清作品| 如何舔出高潮| 精品人妻偷拍中文字幕| 禁无遮挡网站| 亚洲成人av在线免费| 激情 狠狠 欧美| 国产大屁股一区二区在线视频| 国产精品无大码| 久久久国产成人免费| 黄色一级大片看看| 深爱激情五月婷婷| 看黄色毛片网站| 久久久久久伊人网av| 成年免费大片在线观看| 国产乱人视频| 最近在线观看免费完整版| h日本视频在线播放| a级毛片a级免费在线| 男女那种视频在线观看| 久久久久久伊人网av| 一级av片app| 国产高清视频在线播放一区| 美女内射精品一级片tv| 国内精品一区二区在线观看| 综合色丁香网| 日韩欧美在线乱码| 国产男人的电影天堂91| 国产色爽女视频免费观看| 亚洲精品成人久久久久久| 简卡轻食公司| 三级国产精品欧美在线观看| 亚洲精品国产av成人精品 | 美女cb高潮喷水在线观看| 伊人久久精品亚洲午夜| 亚洲av一区综合| 国产伦精品一区二区三区四那| 男人舔女人下体高潮全视频| 蜜桃亚洲精品一区二区三区| 最好的美女福利视频网| 最近的中文字幕免费完整| 国产精品亚洲一级av第二区| 嫩草影院入口| 欧美不卡视频在线免费观看| 亚洲欧美日韩卡通动漫| 国产精品久久久久久久久免| 午夜影院日韩av| 国产片特级美女逼逼视频| 欧美日韩综合久久久久久| 我要搜黄色片| 亚洲无线观看免费| 欧美一区二区国产精品久久精品| av专区在线播放| 亚洲av美国av| 亚洲不卡免费看| 一级a爱片免费观看的视频| 秋霞在线观看毛片| 欧美又色又爽又黄视频| 午夜激情福利司机影院| 久久精品人妻少妇| 在线免费十八禁| 久久精品91蜜桃| 亚洲高清免费不卡视频| 亚洲精品一卡2卡三卡4卡5卡| 免费大片18禁| 亚洲性久久影院| 97在线视频观看| 亚洲欧美日韩高清在线视频| 亚洲国产欧洲综合997久久,| 欧美一级a爱片免费观看看| 一进一出好大好爽视频| 精品人妻一区二区三区麻豆 | 99热全是精品| 99国产极品粉嫩在线观看| 成人av在线播放网站| 久久精品91蜜桃| 午夜福利在线观看吧| 国产精品一区二区三区四区免费观看 | 又爽又黄a免费视频| 少妇丰满av| 欧美3d第一页| 欧美日韩在线观看h| 国产中年淑女户外野战色| 国产女主播在线喷水免费视频网站 | 国产午夜精品论理片| 丰满乱子伦码专区| 成人精品一区二区免费| 成人特级av手机在线观看| 成人综合一区亚洲| 久久精品影院6| 99热这里只有是精品50| 波多野结衣高清作品| 丝袜喷水一区| 亚洲av美国av| 国内精品久久久久精免费| 在线免费观看的www视频| 久久精品国产鲁丝片午夜精品| 少妇丰满av| 晚上一个人看的免费电影| 欧美日韩在线观看h| 97超级碰碰碰精品色视频在线观看| 中文字幕av成人在线电影| 最近在线观看免费完整版| 日本-黄色视频高清免费观看| 国模一区二区三区四区视频| 国产亚洲精品久久久久久毛片| 中国美女看黄片| 白带黄色成豆腐渣| 久久天躁狠狠躁夜夜2o2o| 国产黄片美女视频| 久久天躁狠狠躁夜夜2o2o| 看十八女毛片水多多多| 嫩草影院入口| 插逼视频在线观看| 最近手机中文字幕大全| 亚洲电影在线观看av| 亚洲精品国产av成人精品 | 在线免费观看不下载黄p国产| 麻豆久久精品国产亚洲av| 国产亚洲av嫩草精品影院| 亚洲成人av在线免费| 97热精品久久久久久| av天堂中文字幕网| 97热精品久久久久久| 精品不卡国产一区二区三区| 男插女下体视频免费在线播放| 国产熟女欧美一区二区| 日韩亚洲欧美综合| 国产乱人视频| 夜夜夜夜夜久久久久| 最好的美女福利视频网| 国产v大片淫在线免费观看| 精华霜和精华液先用哪个| 干丝袜人妻中文字幕| 国产真实伦视频高清在线观看| 亚洲成a人片在线一区二区| 人人妻人人澡欧美一区二区| 干丝袜人妻中文字幕| 国产真实伦视频高清在线观看| 一级av片app| 亚洲成人中文字幕在线播放| 色av中文字幕| 免费av毛片视频| 精品人妻一区二区三区麻豆 | 色av中文字幕| 丰满的人妻完整版| 男人狂女人下面高潮的视频| 国产黄片美女视频| 色综合色国产| 在线免费观看的www视频| 一进一出抽搐动态| 99视频精品全部免费 在线| 欧美一区二区亚洲| 亚洲在线观看片| 国产高清激情床上av| 啦啦啦啦在线视频资源| 麻豆久久精品国产亚洲av| 一个人观看的视频www高清免费观看| 午夜免费男女啪啪视频观看 | 美女高潮的动态| 国产女主播在线喷水免费视频网站 | 久久国内精品自在自线图片| 特级一级黄色大片| 色哟哟·www| 乱码一卡2卡4卡精品| 日韩中字成人| 免费在线观看成人毛片| 国产极品精品免费视频能看的| 99久久无色码亚洲精品果冻| 内射极品少妇av片p| 欧美性猛交黑人性爽| 精品久久久久久久久久免费视频| 午夜福利18| 久久九九热精品免费| 伦理电影大哥的女人| 免费不卡的大黄色大毛片视频在线观看 | 男女视频在线观看网站免费| 免费av观看视频| 亚洲国产欧洲综合997久久,| 欧洲精品卡2卡3卡4卡5卡区| 免费看a级黄色片| 久久久国产成人免费| av福利片在线观看| 俄罗斯特黄特色一大片| 免费人成在线观看视频色| 亚洲综合色惰| 欧美成人精品欧美一级黄| 国产三级中文精品| 亚洲av免费高清在线观看| 变态另类成人亚洲欧美熟女| 校园人妻丝袜中文字幕| 人人妻人人看人人澡| 3wmmmm亚洲av在线观看| 欧美高清性xxxxhd video| 国产精品一及| 99精品在免费线老司机午夜| 黄色视频,在线免费观看| 亚洲欧美日韩东京热| 丝袜美腿在线中文| 免费观看在线日韩| 99热6这里只有精品| 久久国产乱子免费精品| 尤物成人国产欧美一区二区三区| 国产国拍精品亚洲av在线观看| 看黄色毛片网站| 男女之事视频高清在线观看| 欧美日韩综合久久久久久| 国产欧美日韩一区二区精品| 国产男靠女视频免费网站| 精品日产1卡2卡| 91麻豆精品激情在线观看国产| 欧美最黄视频在线播放免费| 欧美性猛交╳xxx乱大交人| 国产精品乱码一区二三区的特点| 一级毛片我不卡| 成人毛片a级毛片在线播放| 国产精品美女特级片免费视频播放器| 亚洲美女黄片视频| 亚洲精品粉嫩美女一区| 成人永久免费在线观看视频| 天天躁日日操中文字幕| 1000部很黄的大片| 性欧美人与动物交配| 久久婷婷人人爽人人干人人爱| 国产精品一区二区免费欧美| 久久久久久久久中文| 久久综合国产亚洲精品| 变态另类成人亚洲欧美熟女| 国产精品一区www在线观看| 午夜影院日韩av| 蜜桃久久精品国产亚洲av| 久久久久久久久大av| 成人一区二区视频在线观看| 午夜亚洲福利在线播放| 在线天堂最新版资源| videossex国产| 精品福利观看| 男女视频在线观看网站免费| 亚洲国产精品合色在线| 成人漫画全彩无遮挡| 亚洲av中文字字幕乱码综合| 久久人人爽人人爽人人片va| 久久人妻av系列| 久久久久久久久大av| 别揉我奶头 嗯啊视频| 欧美bdsm另类| 亚洲av不卡在线观看| 少妇人妻一区二区三区视频| 狂野欧美白嫩少妇大欣赏| 免费在线观看影片大全网站| 久久久色成人| 欧美日韩一区二区视频在线观看视频在线 | av在线老鸭窝| 国产私拍福利视频在线观看| 尤物成人国产欧美一区二区三区| 国产精品三级大全| a级毛色黄片| 男人舔奶头视频| 国产精品一区www在线观看| 精品一区二区三区人妻视频| 午夜精品在线福利| 国产精品1区2区在线观看.| 观看美女的网站| 欧美人与善性xxx| 国产精品女同一区二区软件| 日日摸夜夜添夜夜添av毛片| a级毛片a级免费在线| 在线播放无遮挡| 国产精品精品国产色婷婷| 午夜日韩欧美国产| 国产精品爽爽va在线观看网站| 精品久久久久久久人妻蜜臀av| 国产精品三级大全| 亚洲av中文字字幕乱码综合| 国国产精品蜜臀av免费| 国产欧美日韩精品一区二区| 久久人人爽人人爽人人片va| 一级毛片aaaaaa免费看小| 九色成人免费人妻av| 精品国内亚洲2022精品成人| 国产一区二区在线av高清观看| 国产激情偷乱视频一区二区| 男女边吃奶边做爰视频| 看片在线看免费视频| 我要看日韩黄色一级片| 成人亚洲欧美一区二区av| 美女 人体艺术 gogo| 在线观看一区二区三区| 国产欧美日韩精品一区二区| 色播亚洲综合网| 国产精品综合久久久久久久免费| 成人亚洲精品av一区二区| 亚洲熟妇熟女久久| 国产成人freesex在线 | 婷婷精品国产亚洲av在线| 久久久久久久久久久丰满| 欧美一区二区国产精品久久精品| 婷婷精品国产亚洲av| 亚洲人与动物交配视频| 国产蜜桃级精品一区二区三区| 亚洲中文日韩欧美视频| 亚洲五月天丁香| 亚洲av中文av极速乱| 在线播放国产精品三级| 中国国产av一级| 亚洲精品一卡2卡三卡4卡5卡| 亚洲人与动物交配视频| 色综合亚洲欧美另类图片| 国产亚洲精品综合一区在线观看| 久久久久国内视频| 国产男靠女视频免费网站| 老女人水多毛片| 舔av片在线| 全区人妻精品视频| 久久精品国产99精品国产亚洲性色| 亚州av有码| 久久精品国产亚洲av涩爱 | 国产色婷婷99| 搡老妇女老女人老熟妇| 国产女主播在线喷水免费视频网站 | 国产乱人视频| 久久精品国产鲁丝片午夜精品| 亚洲最大成人av| 国产老妇女一区| 美女高潮的动态| 波多野结衣高清作品| 欧美激情久久久久久爽电影| 在线免费十八禁| 99久久精品热视频| 日日摸夜夜添夜夜添av毛片| 国产亚洲av嫩草精品影院| av国产免费在线观看| 日日摸夜夜添夜夜添小说| 少妇猛男粗大的猛烈进出视频 | 我要搜黄色片| 国产免费一级a男人的天堂| 老司机福利观看| 变态另类丝袜制服| 成年版毛片免费区| 麻豆乱淫一区二区| 国产精品人妻久久久影院| 久久久久国产精品人妻aⅴ院| 国产成人aa在线观看| 日韩制服骚丝袜av| 尾随美女入室| 1024手机看黄色片| 草草在线视频免费看| 中文资源天堂在线| 日本欧美国产在线视频| 国产不卡一卡二| 赤兔流量卡办理| 俺也久久电影网| 18禁黄网站禁片免费观看直播| 99久久精品一区二区三区| 六月丁香七月| 亚洲成人精品中文字幕电影| 国产精品电影一区二区三区| 六月丁香七月| 国产午夜精品久久久久久一区二区三区 | 三级国产精品欧美在线观看| 精品久久久久久成人av| 亚洲性夜色夜夜综合| 欧美三级亚洲精品| 99热全是精品| 男人舔女人下体高潮全视频| 精品一区二区三区av网在线观看| 色综合站精品国产| 久久鲁丝午夜福利片| 日韩欧美精品免费久久| 免费在线观看影片大全网站| 免费在线观看成人毛片| 看非洲黑人一级黄片| 丝袜喷水一区| 最近视频中文字幕2019在线8| 白带黄色成豆腐渣| 最近在线观看免费完整版|