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

    力量訓(xùn)練負荷和肌肉疲勞對表面肌電非線性參數(shù)的雙重調(diào)節(jié)

    2022-08-27 10:33:36許全盛白茹雪李世明
    中國體育科技 2022年8期
    關(guān)鍵詞:肌肉疲勞肌電肌力

    許全盛,白茹雪,李世明

    力量訓(xùn)練負荷和肌肉疲勞對表面肌電非線性參數(shù)的雙重調(diào)節(jié)

    許全盛1,白茹雪1,李世明2*

    1.燕山大學(xué),河北秦皇島 066004;2.中國海洋大學(xué),山東青島 266100

    目的:研究力量訓(xùn)練負荷和肌肉疲勞對表面肌電信號非線性參數(shù)的作用規(guī)律,探索精準力量訓(xùn)練中對肌力、肌肉做功和疲勞因素敏感的非線性指標。方法: 6名受試者在非疲勞和疲勞狀態(tài)下分別完成徒手至90%最大負荷共7個負重水平的下蹲起動作,采集股四頭肌的表面肌電、足底壓力和膝關(guān)節(jié)角數(shù)據(jù),提取表面肌電非線性參數(shù)分形維、多尺度熵、科爾莫戈諾夫熵、LZ復(fù)雜度,由足底壓力和膝關(guān)節(jié)角計算膝關(guān)節(jié)力矩,對疲勞前后負重蹲起動作過程的肌電參數(shù)和膝關(guān)節(jié)力矩進行回歸分析。結(jié)果:無論肌肉是否處于疲勞狀態(tài),表面肌電的非線性參數(shù)均隨負重水平增加而增大,但同等負重水平下,肌肉疲勞導(dǎo)致參數(shù)降低。其中,科爾莫戈諾夫熵的線性增長性最顯著,而分形維的數(shù)值穩(wěn)定性最高,其峰值出現(xiàn)在60%~75%最大負重。結(jié)論: 1)骨骼肌系統(tǒng)的復(fù)雜度受到肌肉疲勞和肌肉收縮的共同作用,在運用表面肌電非線性參數(shù)評價肌肉疲勞時,應(yīng)消除肌力或關(guān)節(jié)力矩的影響;2)表面肌電分形維能較好地反映肌肉做功能力,更全面地評估肌肉工作狀態(tài),可以作為精準化快速力量訓(xùn)練及肌肉傷病康復(fù)治療的評價指標。

    表面肌電信號;力量訓(xùn)練;肌肉疲勞;復(fù)雜度;遞增負重

    在健康中國背景下,人們的健身意愿日益強烈,參與健步走、廣場舞、慢跑等有氧運動越來越普遍。相較有氧運動,我國群眾無氧運動參與程度較低,尤其是力量練習(xí),盡管力量素質(zhì)對各年齡段人群的健身活動都很重要,但對多數(shù)人而言力量練習(xí)的量和強度較難把握。因此,科學(xué)的健身指導(dǎo)非常必要,一方面應(yīng)實時監(jiān)測肌肉活動狀態(tài),如果肌肉疲勞了仍繼續(xù)練習(xí)容易造成肌肉和關(guān)節(jié)的損傷;另一方面,應(yīng)控制好練習(xí)量和強度,量和強度過大容易受傷,過小則無法達到最佳的健身效果。檢驗力量練習(xí)對個體的肌肉力量的鍛煉效果,可通過表面肌電(surface electromyography, sEMG)信號分析對肌肉活動狀態(tài)進行評定,預(yù)防過度訓(xùn)練導(dǎo)致的肌肉損傷,并使力量練習(xí)的量和強度達到必要的水平,實現(xiàn)力量練習(xí)的科學(xué)化和精準化。

    有研究證實,sEMG時域和頻域指標在靜態(tài)收縮和中低強度動態(tài)收縮時對肌肉活動的解釋效果很好,比如,sEMG平均頻率、中值頻率隨肌肉疲勞的發(fā)展呈下降趨勢,而均方根振幅、積分肌電等時域指標則呈上升趨勢(王健,2001)。此外,肌肉疲勞與肌力、運動負荷、肌電振幅和頻譜間存在內(nèi)在聯(lián)系:運動負荷改變迫使肌力變化,負荷增大會加速疲勞發(fā)展,而肌電振幅和頻譜對肌力和疲勞狀態(tài)具有雙重依賴性(李世明,2014)。在中等負荷范圍內(nèi),肌電振幅隨肌力增大和疲勞的產(chǎn)生而增長;頻域指標隨肌力增大而增大,但隨疲勞的發(fā)生而降低。因此,當時域指標增大同時頻域指標下降可判定肌肉已經(jīng)疲勞。但當肌肉處于中高強度的動態(tài)收縮時,由電極位置、動作電位傳導(dǎo)速度、肌力以及肌纖維長度等因素帶來的變化使sEMG具有高度的非平穩(wěn)性、非線性和混沌背景(曲峰,2008;González et al., 2010)。此時,通過肌電時域和頻域分析來評估肌肉疲勞具有較大的局限性,在力量訓(xùn)練和健身活動指導(dǎo)方面無法取得滿意效果。

    有2種方法能夠克服生理信號的非平穩(wěn)性和非線性:1)時頻分析法,根據(jù)信號能量在時域和頻域的聯(lián)合分布計算出時變的頻域參數(shù),尤其適合檢測和評估動態(tài)肌肉疲勞(李世明等,2016; González et al., 2010),但時頻分析的計算量大,數(shù)值穩(wěn)定性也較差;2)非線性動力學(xué)分析,直接在時域提取信號的近似熵、樣本熵(sample entropy,SampEn)、多尺度熵(multi-scale entropy,MSE)、科爾莫戈諾夫熵(Kolmogorov entropy,KE)、LZ復(fù)雜度(LZ Complexity,LZC)、分形維數(shù)(fractal dimension,F(xiàn)D)、赫斯特指數(shù)以及李雅普諾夫指數(shù)等非線性特征(劉秉正等,2004)。這些特征從不同角度反映了系統(tǒng)的復(fù)雜度,由于生理系統(tǒng)往往表現(xiàn)出復(fù)雜的混沌模式和不規(guī)則行為,參考這些非線性特征參數(shù)有助于研究系統(tǒng)本身。SampEn是最早用于衡量生理系統(tǒng)的復(fù)雜性和無序性的參數(shù)之一;MSE是在多個離散時間尺度上的SampEn,能夠更完整的反映系統(tǒng)信息和時域演化進程(Costa et al., 2002; Silva et al., 2015);Kolmogorov把信息熵的概念精確化,用KE量度系統(tǒng)在運動過程中的混亂或無序程度;Lempel和Ziv提出了一種更一般性的復(fù)雜度定義LZC,以反映一個時間序列隨其長度的增長出現(xiàn)新模式的速率(劉秉正等,2004);FD是混沌系統(tǒng)中吸引子幾何結(jié)構(gòu)復(fù)雜性程度的完全表征;非線性參數(shù)與時頻參數(shù)也無須信號的平穩(wěn)性假設(shè),但計算量更小。骨骼肌運動單位、動作電位和神經(jīng)支配具有復(fù)雜性,因此,sEMG非線性分析成為研究肌肉系統(tǒng)的有力手段(Liu et al., 2005),尤其是肌肉在中高強度的動態(tài)收縮時,非線性分析相比傳統(tǒng)時域和頻域分析具有明顯優(yōu)勢。

    根據(jù)生理系統(tǒng)的復(fù)雜性理論,健康系統(tǒng)對復(fù)雜和不良環(huán)境的調(diào)節(jié)適應(yīng)能力更強,復(fù)雜度更高,相反,病理和衰老系統(tǒng)的信息量由于生理控制能力退化而下降,復(fù)雜度降低(Costa et al., 2002; Goldberger et al., 2002; Silva et al., 2015)。對肌肉系統(tǒng)而言,肌肉疲勞時力量輸出調(diào)節(jié)受損,系統(tǒng)處于“病態(tài)”,無論是靜態(tài)還是動態(tài)疲勞,sEMG非線性參數(shù)大體呈下降趨勢(Ancillao et al., 2014; Kathick et al., 2018)。但許全盛等(2017)發(fā)現(xiàn),MSE對肌肉疲勞和收縮強度同時敏感,提示在評估肌肉系統(tǒng)復(fù)雜性時,必須同時考慮這2個因素。

    目前對sEMG肌肉疲勞的研究已較為透徹,結(jié)論也比較明朗,但針對肌肉收縮力量(以下簡稱“肌力”)或收縮強度、運動負荷與sEMG非線性特征關(guān)系的研究仍然不足。特別是哪些或哪一種非線性參數(shù)能同時對肌肉疲勞、運動負荷以及肌肉做功能力的變化更為敏感尚不清楚,這在一定程度上限制了sEMG分析對力量練習(xí)和健身活動的指導(dǎo)作用。因此,本研究通過遞增負重下蹲實驗,研究力量訓(xùn)練中運動負荷、疲勞因素與sEMG 4類典型非線性參數(shù)(MSE、KE、LZC和FD)間的協(xié)同變化規(guī)律,探索力量訓(xùn)練中sEMG敏感指標。

    1 研究對象與方法

    1.1 實驗對象

    以體育專業(yè)6名男生為受試者,年齡(20.7±1.3)歲,身高(1.73±0.03)m,體質(zhì)量(65.90±2.28)kg。所有受試者均為右利手,身體狀況良好,無任何運動損傷,實驗前24 h未熬夜,未飲酒,未參與任何形式的劇烈運動,并在實驗前簽署了知情同意書。

    1.2 實驗數(shù)據(jù)采集

    1.2.1最大力量測試

    采用McNeely推薦的間接測試法測試受試者蹲起的最大負重,具體步驟為:1)采用受試者能輕松重復(fù)5~10次的重量進行熱身;2)休息2 min;3)重量增加10%~20%,盡量多的重復(fù),使受試者在2~20次達到疲勞;4)利用公式計算最大負重。計算公式為:最大力量(RM)=[(0.033×重復(fù)次數(shù))×重量]+重量,其中,重復(fù)次數(shù)指蹲起的次數(shù),重量指蹲起的重量,如果受試者的重復(fù)次數(shù)超過10次,則休息10 min后,將重量再增加10%~20%,再次測試。

    1.2.2非疲勞狀態(tài)下遞增負重實驗

    選擇負重下蹲(包括下蹲和起立2個階段)動作完成數(shù)據(jù)采集,便于控制實驗條件使肌力和關(guān)節(jié)力矩在較大范圍內(nèi)變化,便于大量重復(fù)操作。準備姿勢為:將杠鈴置于頸后,雙手抓握杠鈴,身體挺直,兩眼平視前方,兩腳與肩同寬,足尖外斜約15°,負重下蹲后接著負重起立還原。實驗員指導(dǎo)受試者依次在三維測力臺上完成從0% 1RM(徒手)、15% 1RM、30% 1RM、45% 1RM、60% 1RM、75% 1RM和90% 1RM共7個負重水平的下蹲起動作,其中30% 1RM、45% 1RM對應(yīng)力量耐力練習(xí)模式,60% 1RM、75% 1RM對應(yīng)爆發(fā)力量練習(xí)模式,90% 1RM對應(yīng)最大力量練習(xí)模式。每次動作重心都下蹲到最低位(深蹲)后,再全力快速蹬起,每個負重水平重復(fù)做3次,前2次為熟悉性操作,在實驗員的指導(dǎo)下感受動作速度和幅度,最后1次時采集實驗數(shù)據(jù),每個負重水平間休息2 min。實驗過程中完成如下測量:1)采用三維測力臺測試雙腳受力及壓心變化,采樣頻率1 000 Hz;2)采用Simi-Motion運動錄像解析系統(tǒng)在矢狀面內(nèi)對動作過程進行拍攝,拍攝頻率100 Hz,采集右側(cè)膝關(guān)節(jié)角和角速度;3)采用JE-TB0810八通道肌電系統(tǒng)采集右側(cè)腹直肌、股內(nèi)側(cè)肌、股外側(cè)肌的sEMG信號,采樣頻率1 000 Hz。實驗時將測力臺的觸發(fā)信號輸入肌電儀和運動錄像解析系統(tǒng),實現(xiàn)硬件同步。

    1.2.3疲勞實驗

    采用部義峰等(2013)研究中的蹬車方案,利用功率自行車給受試者施加遞增負荷至疲勞狀態(tài)。受試者手戴心率遙測表,首先在零負荷下蹬踏功率自行車(MONARK824型)1 min,適應(yīng)50轉(zhuǎn)/min的轉(zhuǎn)速后進入正式測試。起始負荷50 W,每3 min遞增50 W,轉(zhuǎn)速保持50轉(zhuǎn)/min,負荷遞增到200 W后不再增加。每一負荷結(jié)束前5 s記錄自覺用力評分指數(shù)(rating of perceived exertion,RPE)等級及心率。如果受試者經(jīng)鼓勵督促仍不能按規(guī)定強度運動持續(xù)10 s以上,并同時滿足達到預(yù)期最大心率(220-年齡)和RPE>18,便認為受試者股四頭肌已進入疲勞狀態(tài),停止運動。

    1.2.4非疲勞狀態(tài)下遞增負重實驗

    疲勞實驗后立即進行負重下蹲起實驗。疲勞后受試者很難連續(xù)力量耐力和爆發(fā)力量練習(xí),隨著實驗的進行,疲勞程度逐漸降低,而且不同受試者的恢復(fù)能力存在差異,因此,在疲勞狀態(tài)下只進行1~2個負重水平的測試,要求與非疲勞狀態(tài)下測試一致。完成后休息足夠長的時間,使機體完全恢復(fù)后,再次通過功率自行車達到疲勞狀態(tài),進行下一個負重水平的測試。

    1.3 肌電信號采集與預(yù)處理

    采用八通道肌電采集系統(tǒng)(JE-TB0810,中國)采集sEMG信號,肌電儀佩戴于受試者腰部,電極片采用Ag/AgCI一次性使用圓形心電電極(上海鈞康醫(yī)用設(shè)備有限公司),直徑2.5 cm;3個電極成等邊三角形貼于待測肌肉肌腹沿肌纖維走向的方向上。放置電極貼前分別用磨砂紙和醫(yī)用酒精擦拭皮膚以減小電阻。肌電儀配備截止頻率為10 Hz的2階巴特沃茲高通濾波器和截止頻率為500 Hz的低通濾波器,以消除環(huán)境噪聲和運動偽跡的影響,信號采樣頻率1 000 Hz,增益119 dB,共模抑制比>100 dB,輸入阻抗>1012Ω,輸入范圍±3 mV。將采集的sEMG信號導(dǎo)入計算機,利用MATLAB 2014a對股直肌、股內(nèi)側(cè)肌和股外側(cè)肌的疲勞前后7個負重級別的sEMG信號進行預(yù)處理,設(shè)計2階橢圓數(shù)字陷波器濾除50 Hz工頻干擾后,將每位受試者共42段原始sEMG信號(3塊肌肉在非疲勞和疲勞后的7級負荷)轉(zhuǎn)換為待處理的時間序列數(shù)據(jù),備后續(xù)分析使用。

    1.4 參數(shù)計算

    選擇MSE、KE、LZC和FD 4種具有代表性的非線性特征參數(shù),利用MATLAB直接在sEMG信號時域提取參數(shù)值,從不同角度反映肌肉系統(tǒng)復(fù)雜度。利用逆向動力學(xué)算法由采集的足底測力數(shù)據(jù)計算負重下蹲動作的膝關(guān)節(jié)力矩,以此反映肌肉收縮強度或肌力水平。

    1.4.1MSE

    本研究將SampEn在所有尺度上的平均值(mean multi-scale entropy,MMSE)作為MSE復(fù)雜度測量。研究表明,MMSE相比單一尺度的SampEn能更有效地識別肌肉動態(tài)疲勞(許全盛等,2017)。

    1.4.2KE

    1.4.3LZC

    3)當趨于∞時幾乎所有的()都會趨于一常數(shù),即:

    其中,()為序列的漸進行為,為字符串中不同字符的個數(shù)。用(對(進行歸一化得:

    1.4.4FD

    分形物體的基本特征是其長度、面積或體積取決于測量尺度。將sEMG信號視為一維曲線,F(xiàn)D就是曲線復(fù)雜性的特征,F(xiàn)D最常用的是盒型算法(box counting)(劉秉正等,2004)和Katz算法(Katz, 1988)。Liu等(2005)研究認為,Katz算法對信號的隨機性最敏感,對運動相關(guān)信號的量化能力最強。因此,本研究選擇Katz算法量化sEMG復(fù)雜性變化,其原理為:

    對一段連續(xù)平面曲線,F(xiàn)D的一般定義為:

    其中,是曲線的總長度,是曲線的直徑(平面范圍)。對一段波形,它們是由一系列點(x,y)的有序集,總長度是連續(xù)點之間距離的之和:

    其中,為序列長度,即xy的個數(shù)。由于波形有一個自然的起始點,平面范圍可以看作是起始點(1,1)與波形的任何其他點(點)間的最遠距離,即:

    對一維離散波形,F(xiàn)D計算需歸一化為平均步長/,代入公式(8)可得Katz的FD:

    其中,=/是曲線中的步數(shù)。對單調(diào)進行“向前”的平面曲線,F(xiàn)D介于1.0(直線)和1.5(高尖峰波)之間。真實波形永遠不會變得足夠復(fù)雜來填充一個平面,因此波形永遠不會有接近平面2.0維數(shù)FD。

    1.5 膝關(guān)節(jié)力矩與非線性參數(shù)和疲勞因素的回歸分析

    本研究通過膝關(guān)節(jié)角運動學(xué)數(shù)據(jù)、足底測力臺數(shù)據(jù),利用逆向動力學(xué)算法計算得到膝關(guān)節(jié)合力矩,膝力矩反映膝關(guān)節(jié)肌群的肌力貢獻,大小取決于肌肉收縮強度。膝力矩計算的核心是牛頓-歐拉方程,其中膝關(guān)節(jié)角度通過運動錄像解析系統(tǒng)根據(jù)髖膝踝骨性標志點的坐標計算得到,下肢關(guān)節(jié)轉(zhuǎn)動慣量根據(jù)受試者身高、體質(zhì)量,通過運動解析系統(tǒng)中人體環(huán)節(jié)參數(shù)模型獲得。測力臺數(shù)據(jù)包括足底壓力和轉(zhuǎn)矩,將測力臺關(guān)于矢狀面受力和冠狀軸轉(zhuǎn)矩數(shù)據(jù),和轉(zhuǎn)動慣量、膝關(guān)節(jié)角一起代入牛頓方程得到膝關(guān)節(jié)力矩,其數(shù)值求解過程通過MATLAB多剛體力學(xué)工具箱實現(xiàn)。

    利用MATLAB對所有受試者在非疲勞和疲勞態(tài)下分別完成負重下蹲動作過程膝關(guān)節(jié)力矩的平均值與股直肌、股內(nèi)側(cè)肌和股外側(cè)肌sEMG非線性參數(shù)值間的相關(guān)性進行統(tǒng)計分析。通過繪制非線性參數(shù)隨膝關(guān)節(jié)力矩的變化的散點圖,運用最小二乘法對散點圖進行線性回歸并計算其回歸系數(shù)、可決系數(shù)2和統(tǒng)計量,探究對力量練習(xí)負荷和疲勞程度敏感的sEMG指標,回歸方程如下:

    2 研究結(jié)果

    2.1 時域波形

    股直肌、股內(nèi)側(cè)肌和股外側(cè)肌是對負重蹲起動做貢獻最大的3塊肌肉。本研究中每位受試者在非疲勞狀態(tài)和疲勞狀態(tài)完成負重下蹲動作各產(chǎn)生21段肌電信號(7種負重水平,3塊肌肉),以受試者馬XX股內(nèi)側(cè)肌疲勞前后的原始肌電波形為例,無論是非疲勞狀態(tài)還是疲勞狀態(tài)時,肌電振幅和信號時長均隨負重水平從0% 1RM(徒手)至90% 1RM而遞增;與非疲勞狀態(tài)相比,負重較輕時振幅有所增大,但負重較大時并不顯著(圖1)。

    圖1 受試者馬XX在非疲勞狀態(tài)(左)和疲勞狀態(tài)(右)下不同負重時股內(nèi)側(cè)肌的原始肌電波形

    Figure 1. Raw sEMG of Vastus Medialis of Subject Ma XX in Fatigue-Free(Left) and Fatigue State(Right) at Different Loads

    2.2 sEMG參數(shù)

    6名受試者股直肌、股內(nèi)側(cè)肌和股外側(cè)肌sEMG 4個非線性參數(shù)在疲勞和非疲勞狀態(tài)下的計算結(jié)果(平均值±標準差)中MMSE計算時容許誤差取0.2,嵌入維取2,時間延遲取1,考慮到受試者完成下蹲起動作記錄的長度大多在1 500~2 000 ms,MSE最大時間尺度選擇50(表1)。

    2.3 膝關(guān)節(jié)力矩

    根據(jù)逆向動力學(xué)算法得到6名受試者的膝關(guān)節(jié)力矩在每個負重水平的所有最大值的平均值和標準差(表2)。

    以受試者馬XX疲勞狀態(tài)和非疲勞狀態(tài)的膝關(guān)節(jié)力矩隨時間變化為例,膝力矩最大值和平均值均隨負重增加而增加,由于完成動作的難度隨負重逐漸增大,膝關(guān)節(jié)角速度逐漸變慢,膝力矩發(fā)揮作用的時長逐漸增加(圖2)。

    2.4 膝力矩和疲勞對肌電參數(shù)的影響

    2.4.1sEMG非線性參數(shù)隨遞增負重水平的變化

    以股內(nèi)側(cè)肌為例,將疲勞狀態(tài)和非疲勞狀態(tài)時各參數(shù)的平均值和標準差數(shù)據(jù)繪制成按負重水平遞增的柱狀圖(圖3)。從徒手到最大負重90% 1RM全部7個負重水平,疲勞后的sEMG非線性參數(shù)均顯著低于非疲勞狀態(tài)下的數(shù)值。無論是疲勞還是非疲勞狀態(tài),4種參數(shù)變化趨勢基本類似,總體上均隨負重遞增而遞增。其中,KE的線性增長性最顯著,MMSE和LZC在30% 1RM負重時稍有下降,但LZC在75% 1RM時稍大于90% 1RM,而FD在60% 1RM和75% 1RM時數(shù)值接近,但在90% 1RM時有明顯下降;從數(shù)值穩(wěn)定性看,MMSE在不同受試者間個體差異最大,LZC次之,而FD的方差最小,表現(xiàn)最為穩(wěn)定。

    為觀察肌電參數(shù)在不同肌肉間的變化情況,將股直肌、股內(nèi)側(cè)肌和股外側(cè)肌的MMSE隨負重水平的變化繪制為柱狀圖(圖4),3塊肌肉的肌電參數(shù)在不同負重水平下的排序并不完全相同,但疲勞前后的總體變化趨勢基本一致,其他參數(shù)在不同肌肉間變化情況與MMSE類似。

    由圖5可知,肌電參數(shù)與膝關(guān)節(jié)力矩總體呈正相關(guān),同時,疲勞狀態(tài)下膝關(guān)節(jié)力矩和肌電參數(shù)均降低(數(shù)據(jù)點在非疲勞時的左下方)。由圖可見,MMSE和LZC的線性性較差,30% 1RM和45% 1RM負重下疲勞和非疲勞狀態(tài)MMSE和LZC數(shù)值有所交疊,而KE和FD增長的線性性更明顯。進一步觀察可見LZC和FD的峰值都出現(xiàn)在75% 1RM,但FD在60% 1RM時的數(shù)值比LZC更加接近75% 1RM。

    2.4.2sEMG非線性參數(shù)與膝關(guān)節(jié)力矩的回歸分析

    將6名受試者從徒手至90% 1RM共7個負重水平的sEMG參數(shù)和膝關(guān)節(jié)力矩的對應(yīng)關(guān)系繪制成散點圖(圖6)。其中疲勞狀態(tài)和非疲勞狀態(tài)的數(shù)據(jù)點各由42個散點組成,根據(jù)最小二乘法對散點數(shù)據(jù)進行線性回歸。結(jié)果發(fā)現(xiàn),疲勞狀態(tài)數(shù)據(jù)點的回歸線均位于非疲勞狀態(tài)的下方;從散點的分布范圍看,MMSE最大,F(xiàn)D最為集中,而KE和LZC介于兩者之間;從擬合的直觀效果看KE的線性性最顯著。

    表1 受試者肌電非線性參數(shù)匯總表

    表2 受試者的膝關(guān)節(jié)力矩最大值統(tǒng)計

    總結(jié)4類非線性參數(shù)與膝關(guān)節(jié)力矩的回歸方程的可決系數(shù)(2統(tǒng)計量)、統(tǒng)計量以及回歸系數(shù),結(jié)果1全部為正,而2全部為負,說明sEMG的非線性參數(shù)與肌力或力矩呈正相關(guān)而與疲勞程度的呈負相關(guān)。KE的2統(tǒng)計量和統(tǒng)計量顯著高于其他3類參數(shù),說明其線性擬合的效果最佳(表3)。

    圖2 受試者馬XX在非疲勞狀態(tài)(a)和疲勞狀態(tài)(b)遞增負重練習(xí)的膝關(guān)節(jié)力矩

    Figure 2. Knee Joint Torque of Subject Ma XX in Fatigue-Free (a) and Fatigue State (b) at Different Loads

    3 分析與討論

    力量素質(zhì)對各年齡段人群都很重要。力量訓(xùn)練過程中,練習(xí)者直觀感覺是運動負荷增大迫使肌肉加強收縮,以產(chǎn)生更大肌力或力矩來完成動作,同時,負荷增大也會加速疲勞發(fā)展。但無論是健身愛好還是專業(yè)訓(xùn)練,掌握力量練習(xí)的量和強度僅憑感覺是不夠的,因為人的主觀感受是整體性的,尤其在專項力量訓(xùn)練中很難形成某塊肌肉或某個肌肉群的精確感覺,所以需要有科學(xué)的量化指導(dǎo)。通過sEMG信號對肌肉狀態(tài)進行監(jiān)測評估可以具體到某一塊肌肉,由力量訓(xùn)練負荷和練習(xí)模式不同所帶來的肌肉收縮力、收縮速度、關(guān)節(jié)力矩、輸出功率、疲勞程度的變化信息都蘊含在sEMG信號中,通過時域、頻域、時頻域以及非線性參數(shù)定量得以反映。利用sEMG對肌肉活動狀態(tài)進行評定,能有效預(yù)防過度訓(xùn)練導(dǎo)致的肌肉損傷,制定最佳的練習(xí)量和強度,使肌肉力量獲得最大限度的超量恢復(fù)。

    圖3 股內(nèi)側(cè)肌sEMG非線性參數(shù)隨遞增負重水平的變化

    Figure 3. The Variation of sEMG Nonlinear Parameters in Medial Femoral with Increasing Load Level

    圖4 股直肌、股內(nèi)側(cè)肌和股外側(cè)肌在非疲勞狀態(tài)(a)和疲勞狀態(tài)(b)sEMG平均多尺度熵隨負重水平的變化

    Figure 4. The Variation of Mean Multiscale Entropy of sEMG with Load Level in Rectus Femoris, Medial Femoris and Lateral Femoris in Fatigue-Free (a) and Fatigue State (b)

    圖5  sEMG非線性參數(shù)隨膝關(guān)節(jié)力矩最大值的變化

    Figure 5. Variation of Nonlinear Parameters of sEMG with the Maximum Knee Joint Torque

    肌力與骨骼肌運動單元的募集程度呈正相關(guān),運動單元募集和發(fā)放率最大時形成最大肌力,sEMG信號是骨骼肌運動單元異步發(fā)放產(chǎn)生的多個脈沖串疊加的結(jié)果,運動單元越多,產(chǎn)生的sEMG信號越復(fù)雜。當肌肉被激活并施加負荷時,由此產(chǎn)生的sEMG信號的復(fù)雜性隨肌肉收縮加劇而增加。因此,可以推測肌肉系統(tǒng)為適應(yīng)收縮強度增大、肌力增加的變化使sEMG的復(fù)雜度也隨之增大。本研究選擇的sEMG 4種非線性參數(shù)從不同側(cè)面反映了肌肉系統(tǒng)的復(fù)雜度。股四頭肌是膝關(guān)節(jié)控制肌群中最重要的大肌肉,在屈膝和伸膝過程中分別完成退讓式收縮和克制性收縮。圖2表明,隨著負重等級增大,以股四頭肌為主的伸肌群收縮強度增大,以產(chǎn)生更大的肌力和膝力矩。表3、圖3、圖4、圖5中系數(shù)1都說明了sEMG 4種非線性參數(shù)值與負重水平或膝力矩正相關(guān),也證明了肌肉收縮強度和肌力水平的確與復(fù)雜度呈正相關(guān),但在圖5中FD和LZC的峰值并沒有出現(xiàn)在最大負重上。

    圖6  sEMG非線性參數(shù)與膝關(guān)節(jié)力矩的回歸結(jié)果

    Figure 6. Regression Results of sEMG Nonlinear Parameters and Knee Joint Torque

    表3  sEMG非線性參數(shù)對膝關(guān)節(jié)力矩回歸統(tǒng)計量

    健康系統(tǒng)對復(fù)雜和不良環(huán)境的調(diào)節(jié)適應(yīng)能力更強,而病態(tài)和衰老系統(tǒng)的信息量由于生理控制能力退化而下降,表現(xiàn)為復(fù)雜度降低(Costa et al., 2002; Goldberger et al., 2002; Silva et al., 2015)。對肌肉系統(tǒng)而言,肌肉疲勞時力量輸出調(diào)節(jié)受損,系統(tǒng)可視為處于“病態(tài)”。本實驗結(jié)果證實,疲勞時肌肉活動的復(fù)雜度的確降低。其中,圖3、圖4說明在同等負荷下,疲勞后的sEMG 4類非線性參數(shù)值均顯著低于非疲勞狀態(tài);疲勞抵消了一部分非線性參數(shù)隨膝力矩增長的趨勢(圖5),同時,疲勞引起了骨骼肌最大隨意收縮力量和輸出功率下降,同等負重下的膝關(guān)節(jié)力矩均有不同程度的降低;6名受試者的平均值的變化趨勢(圖3、圖4、圖5)和42個樣本點的線性回歸,疲勞狀態(tài)數(shù)據(jù)點的回歸線均位于非疲勞狀態(tài)的左下方(圖6),同樣說明疲勞抵消了非線性參數(shù)隨負重的一部分增長;系數(shù)2則說明,sEMG非線性參數(shù)與疲勞程度呈負相關(guān)(表3)。綜上,本研究認為sEMG非線性參數(shù)對力量訓(xùn)練負荷和疲勞狀態(tài)具有雙重依賴性,或者說肌肉系統(tǒng)復(fù)雜度對肌力和疲勞因素同時敏感。

    肌力的產(chǎn)生和肌肉收縮速度的結(jié)合將輸出肌肉功率。對sEMG非線性參數(shù)與肌肉做功之間聯(lián)系的研究較少。Liu等(2005)對受試者手持測力裝置時腦電信號分析,腦電FD與所測握力呈線性相關(guān);Gupta等(1997)研究表明,sEMG FD對肌力的大小和收縮速度都很敏感,與肌肉的運動潛力密切相關(guān);而Ancillao等(2014)發(fā)現(xiàn),縱跳高度與FD正相關(guān)。人體縱跳高度與下肢肌肉爆發(fā)力也就是肌肉功率密切相關(guān),因此,可推斷FD是一個對肌力和肌肉功率同時敏感的指標,本實驗結(jié)果也證實了這一點。從圖5可見,F(xiàn)D峰值都出現(xiàn)在60% 1RM附近,而LZC出現(xiàn)在75% 1RM,考慮到60%~75% 1RM的力量訓(xùn)練接近爆發(fā)力量練習(xí)模式,即肌力和肌肉收縮速度的乘積最大,肌肉輸出功率最大的模式,而中低負荷時(30% 1RM和45% 1RM)肌肉快速收縮和高負荷時(90% 1RM)肌肉慢速收縮分別屬于力量耐力和最大力量練習(xí)模式,這2種情況下肌肉輸出功率都無法達到最大值接近力量耐力模式,因此,F(xiàn)D和LZC屬于與肌力和肌肉功率同時敏感的指標,但FD更能反映爆發(fā)力量素質(zhì)。

    進一步分析可知,由于FD和LZC的峰值不在最大負重水平上,破壞了公式(12)的線性假設(shè),所以表3中FD和LZC的2和統(tǒng)計量均顯著小于線性擬合效果最好的KE,而圖5也更直觀地反映了FD對肌肉功率的敏感性。

    4 結(jié)論與建議

    骨骼肌系統(tǒng)的復(fù)雜度受到肌肉疲勞和收縮強度的雙重調(diào)制,在運用非線性參數(shù)評價肌肉動態(tài)疲勞時,應(yīng)同時考慮肌肉收縮強度的影響。

    FD能更全面地反映骨骼肌的工作狀態(tài)。在實際工作中,F(xiàn)D可作為快速力量訓(xùn)練效果檢驗、運動員選材以及肌肉傷病康復(fù)治療的評價指標。

    sEMG信號非線性參數(shù)提取在時域直接進行,計算復(fù)雜度較小,且無須信號的平穩(wěn)性假設(shè),有利于對力量練習(xí)者的肌肉活動狀態(tài)的實時監(jiān)測和評估,實現(xiàn)力量訓(xùn)練的科學(xué)化、定量化和精準化。

    部義峰,李世明,2013. 遞增負荷誘發(fā)肌肉疲勞過程中肌肉sEMG與負荷的非線性關(guān)系研究[J]. 山東體育學(xué)院學(xué)報,29(6): 47-52.

    李世明,2014.運動技術(shù)診斷概論[M].北京:科學(xué)出版社.

    李世明,許全盛,翟佳麗,等,2016.基于表面肌電時頻分析的遞增負荷誘導(dǎo)肌肉周期動態(tài)疲勞估計[J].中國體育科技, 52(3):48-56.

    劉秉正,彭建華, 2004.非線性動力學(xué)[M].北京:高等教育出版社.

    曲峰,2008.運動員表面肌電信號與分形[M].北京:北京體育大學(xué)出版社.

    王健,2001.靜態(tài)負荷肌肉疲勞過程中肌肉功率譜轉(zhuǎn)移特征[J].中國運動醫(yī)學(xué)雜志,20(2): 199-201.

    許全盛,李世明,季淑梅,2017.基于sEMG多尺度熵分析的遞增負荷誘導(dǎo)肌肉動態(tài)疲勞評估[J].天津體育學(xué)院學(xué)報,32(4): 344-349.

    ANCILLAO A, GALLI M, RIGOLDI C, et al., 2014. Linear correlation between fractal dimension of surface EMG signal from rectus femoris and height of vertical jump[J]. Chaos Soliton Fract, 66: 120-126.

    COSTA M, GOLDBERGER A L, PENG C K, 2002. Multiscale entropy analysis of complex physiologic time series[J]. Phys Rev Lett, 89(6): 1-4.

    GOLDBERGER A L, PENG C K, LIPSITZ L A, 2002. What is physiologic complexity and how does it change with ageing and disease?[J]. Neurobiol Aging, 23(1): 23-26.

    GONZáLEZ M, MALANDA A, NAVARRO A, et al., 2010. EMG spectral indices and muscle power fatigue during dynamic contractions[J]. J Electromyogr Kinesiol, 20(2): 233-240.

    GUPTA V, SURYANARAYANAN S, REDDY N P, 1997. Fractal analysis of surface EMG signals from the biceps[J]. Int J Med Inform, 45(3): 185-192.

    KATHICK P A, GHOSH D M, RAMAKRISHNAN S, 2018. Surface electromyography based muscle fatigue detection using high-resolution time-frequency methods and machine learning algorithms[J]. Comput Methods Programs Biomed, 154: 45-56.

    KATZ M J, 1988. Fractals and the analysis of waveforms[J]. Comput Biol Med, 18(3):145-156.

    LIU J Z, YANG Q, YAO B, et al., 2005. Linear correlation between fractal dimension of EEG signal and handgrip force[J]. Biol Cybern, 93(2):131-140.

    SILVA L E, CABELLA B C, NEVES U P, et al., 2015. Multiscale entropy-based methods for heart rate variability complexity analysis[J]. Physical A, 422(1): 143-152.

    Strength Training Load and Muscle Fatigue Cooperatively Modulates the Nonlinear Parameters of Surface Electromyography

    XU Quansheng1,BAI Ruxue1,LI Shiming2*

    1.Yanshan University, Qinhuangdao 066004, China;2.Ocean University of China, Qingdao 266100, China.

    Objective: This paper is aimed to study the effect of strength training load and muscle fatigue on nonlinear parameters of sEMG signal, and to explore the nonlinear indexes sensitive to muscle strength, muscle work and fatigue factors in precision strength training. Methods: Six volunteers were enrolled to perform weight-bearing crouching experiment under seven loading level from zero to 90% maximum load in fatigue and fatigue-free condition, respectively. The sEMG, plantar pressure and knee joint angle data of quadriceps femoris were collected; the nonlinear parameters of sEMG, including fractal dimension, multi-scale entropy, Kolmogorov entropy and LZ complexity, were extracted, and the knee joint torque was calculated from plantar pressure and knee joint angle by inverse dynamics; the nonlinear parameters of sEMG and knee joint torque in weight-bearing squatting were analyzed by regression analysis. Results: Whether the muscles are non-fatigued or fatigued, the nonlinear parameters increased with the increase of weighting load level, but at the same load level, muscle fatigue lead to lower parameters. The linear increase property of Kolmogonov entropy was the most significant, and the numerical stability of fractal dimension was the best, and the peak value appeared between 60% and 75% of the maximum load. Conclusions: 1) the complexity of skeletal muscle system is affected by muscle fatigue and muscle contraction, hence the influence of muscle strength when using nonlinear parameters to evaluate muscle fatigue should be eliminated; 2) fractal dimension can reflect the ability of muscle power to some extent, and evaluate the working state of skeletal muscle more comprehensively, therefore it could be used as an evaluation index of accurate explosive strength training and rehabilitation of muscle injuries.

    1002-9826(2022)08-0083-08

    10.16470/j.csst.2020081

    國家體育總局科技服務(wù)項目(2017B047)

    許全盛(1973-),男,博士,碩士研究生導(dǎo)師,主要研究方向為生物醫(yī)學(xué)信號處理、運動生物力學(xué),E-mail:123122122@qq.com。

    通信作者簡介:李世明(1969-),男,教授,博士,碩士研究生導(dǎo)師,主要研究方向為運動生物力學(xué),E-mail:leesm0503@ 163.com。

    G808.1

    A

    (2019-12-19;修訂日期:2022-03-23; 編輯:丁合)

    猜你喜歡
    肌肉疲勞肌電肌力
    等速肌力測試技術(shù)在踝關(guān)節(jié)中的應(yīng)用進展
    盆底肌電刺激聯(lián)合盆底肌訓(xùn)練治療自然分娩后壓力性尿失禁的臨床觀察
    BMI對拉力作業(yè)肌肉疲勞的影響研究
    高溫高濕環(huán)境長時間運動后神經(jīng)肌肉疲勞類型與PAP的關(guān)系
    28天肌力恢復(fù),告別夏日肌膚后遺癥
    好日子(2018年9期)2018-10-12 09:57:20
    散打訓(xùn)練對大學(xué)生肌力與肌耐力的影響
    造就美肌力 一日之際在于晨
    Coco薇(2017年12期)2018-01-03 21:18:37
    經(jīng)皮神經(jīng)肌電刺激治療上肢周圍神經(jīng)損傷的療效
    女生穿運動鞋與高跟鞋步行腰背肌電比較
    老年功能性便秘患者盆底表面肌電與臨床癥狀的相關(guān)性研究
    av在线app专区| av一本久久久久| 美女大奶头黄色视频| 国产成人freesex在线| 成人免费观看视频高清| 国产在线免费精品| 欧美老熟妇乱子伦牲交| 久久精品熟女亚洲av麻豆精品| 国产午夜精品一二区理论片| 青春草视频在线免费观看| 蜜桃久久精品国产亚洲av| 夜夜爽夜夜爽视频| 亚洲精品美女久久av网站| 波野结衣二区三区在线| 大话2 男鬼变身卡| 欧美xxⅹ黑人| 一级毛片黄色毛片免费观看视频| 日本黄色日本黄色录像| 麻豆成人av视频| 人妻一区二区av| 日本免费在线观看一区| 91国产中文字幕| 精品久久久久久久久亚洲| 亚洲国产av影院在线观看| 亚洲av.av天堂| 我的女老师完整版在线观看| 亚洲精品色激情综合| 十分钟在线观看高清视频www| 女人精品久久久久毛片| 欧美少妇被猛烈插入视频| 精品一区在线观看国产| 亚洲情色 制服丝袜| av国产精品久久久久影院| 黄色毛片三级朝国网站| 国产精品99久久99久久久不卡 | 黑人高潮一二区| 国产成人91sexporn| 亚洲av福利一区| 一级二级三级毛片免费看| 精品少妇久久久久久888优播| 国产av国产精品国产| 久久青草综合色| 国产精品熟女久久久久浪| 国产精品99久久99久久久不卡 | 黄色视频在线播放观看不卡| 久久人人爽人人片av| 欧美成人精品欧美一级黄| 国产69精品久久久久777片| 51国产日韩欧美| 飞空精品影院首页| 欧美成人午夜免费资源| 国产日韩欧美在线精品| 精品一品国产午夜福利视频| 亚洲精品久久成人aⅴ小说 | 少妇人妻久久综合中文| 大码成人一级视频| av一本久久久久| 99久久综合免费| 国产精品久久久久成人av| 亚洲欧美日韩卡通动漫| 亚洲欧美一区二区三区国产| 极品少妇高潮喷水抽搐| 丝袜喷水一区| 十八禁高潮呻吟视频| 亚洲人成网站在线观看播放| 欧美激情 高清一区二区三区| 国产av精品麻豆| 精品人妻偷拍中文字幕| 18在线观看网站| 欧美bdsm另类| 我的女老师完整版在线观看| 少妇的逼好多水| 午夜福利视频精品| 嫩草影院入口| 免费观看的影片在线观看| 人妻系列 视频| 内地一区二区视频在线| 国产亚洲av片在线观看秒播厂| 午夜视频国产福利| 亚洲欧美成人精品一区二区| 欧美激情 高清一区二区三区| 夫妻性生交免费视频一级片| 在线看a的网站| 亚洲精品aⅴ在线观看| 亚洲精品乱码久久久久久按摩| 天天躁夜夜躁狠狠久久av| 我的老师免费观看完整版| 日韩一本色道免费dvd| av网站免费在线观看视频| 美女大奶头黄色视频| 国产亚洲欧美精品永久| 久久ye,这里只有精品| 久久久久国产精品人妻一区二区| 人妻一区二区av| 欧美日韩在线观看h| 黑人欧美特级aaaaaa片| 中文字幕免费在线视频6| 五月开心婷婷网| 少妇丰满av| 久久精品久久精品一区二区三区| 国产在线一区二区三区精| 美女脱内裤让男人舔精品视频| www.色视频.com| 插阴视频在线观看视频| 精品亚洲成国产av| 肉色欧美久久久久久久蜜桃| 爱豆传媒免费全集在线观看| 简卡轻食公司| 国产精品99久久久久久久久| 久久免费观看电影| 校园人妻丝袜中文字幕| 久久久久久久亚洲中文字幕| 久久热精品热| 亚洲av电影在线观看一区二区三区| 少妇熟女欧美另类| 黄色一级大片看看| 国产成人91sexporn| 日韩免费高清中文字幕av| 日本91视频免费播放| 日本免费在线观看一区| 精品一区在线观看国产| 天堂俺去俺来也www色官网| 久久久a久久爽久久v久久| 乱人伦中国视频| 亚洲国产欧美日韩在线播放| 高清在线视频一区二区三区| 精品少妇黑人巨大在线播放| 国产成人精品久久久久久| 精品国产一区二区久久| 人人澡人人妻人| 熟女人妻精品中文字幕| 国产日韩欧美亚洲二区| 少妇高潮的动态图| 2022亚洲国产成人精品| 日本猛色少妇xxxxx猛交久久| 亚洲高清免费不卡视频| 性色av一级| h视频一区二区三区| 九九在线视频观看精品| 哪个播放器可以免费观看大片| 国产男女内射视频| 国产免费一区二区三区四区乱码| 久久99热6这里只有精品| 国产淫语在线视频| 免费av不卡在线播放| 日本色播在线视频| 制服人妻中文乱码| 久久久久国产精品人妻一区二区| 免费少妇av软件| 久久av网站| 91午夜精品亚洲一区二区三区| 亚洲综合色网址| av一本久久久久| 国产精品女同一区二区软件| 免费日韩欧美在线观看| 丰满饥渴人妻一区二区三| 欧美3d第一页| 国产高清三级在线| 久久精品人人爽人人爽视色| 欧美日韩精品成人综合77777| 国产视频首页在线观看| 亚洲综合精品二区| 自线自在国产av| 亚洲无线观看免费| a级毛片在线看网站| 777米奇影视久久| 一级毛片我不卡| 丝瓜视频免费看黄片| 亚洲欧美成人精品一区二区| 精品人妻偷拍中文字幕| 欧美精品国产亚洲| av不卡在线播放| 午夜影院在线不卡| 中文字幕人妻丝袜制服| 精品国产一区二区久久| 国精品久久久久久国模美| 一级黄片播放器| 日韩电影二区| 伊人久久精品亚洲午夜| 大香蕉久久网| 久久人妻熟女aⅴ| 26uuu在线亚洲综合色| 国产精品一区二区在线不卡| 嫩草影院入口| 王馨瑶露胸无遮挡在线观看| 少妇人妻精品综合一区二区| 亚洲欧美中文字幕日韩二区| 国产日韩欧美亚洲二区| 久久久久国产网址| 汤姆久久久久久久影院中文字幕| 精品久久蜜臀av无| 亚洲精品国产av成人精品| 亚洲av不卡在线观看| 大陆偷拍与自拍| 一边亲一边摸免费视频| 成年av动漫网址| 久久人妻熟女aⅴ| 久久综合国产亚洲精品| 新久久久久国产一级毛片| 男女啪啪激烈高潮av片| 成人无遮挡网站| 国产淫语在线视频| 观看av在线不卡| 免费观看的影片在线观看| 亚洲性久久影院| 精品久久久久久电影网| 热re99久久精品国产66热6| 国产成人精品在线电影| 多毛熟女@视频| 国产一级毛片在线| 日韩一本色道免费dvd| 91精品国产国语对白视频| 中文字幕久久专区| 亚洲国产av新网站| 国产精品国产三级国产专区5o| 亚洲av二区三区四区| 男女国产视频网站| 国产亚洲精品第一综合不卡 | 免费不卡的大黄色大毛片视频在线观看| 日韩 亚洲 欧美在线| 不卡视频在线观看欧美| 国产国拍精品亚洲av在线观看| 久久久久人妻精品一区果冻| 亚洲美女搞黄在线观看| 国产永久视频网站| 亚洲欧美清纯卡通| 成人亚洲精品一区在线观看| 搡老乐熟女国产| 九九爱精品视频在线观看| 高清黄色对白视频在线免费看| 色视频在线一区二区三区| 久久久国产欧美日韩av| 国产极品粉嫩免费观看在线 | 国产熟女欧美一区二区| 国产精品一区二区在线不卡| 精品一区二区三卡| 精品国产乱码久久久久久小说| 精品一区二区三区视频在线| 国产高清不卡午夜福利| av国产精品久久久久影院| 亚洲三级黄色毛片| 九九久久精品国产亚洲av麻豆| 亚洲精品第二区| 亚洲性久久影院| 亚州av有码| 18禁动态无遮挡网站| 人成视频在线观看免费观看| 久久99热这里只频精品6学生| 国产亚洲av片在线观看秒播厂| 国产精品成人在线| 狂野欧美白嫩少妇大欣赏| 亚洲av二区三区四区| 久久ye,这里只有精品| tube8黄色片| 成人手机av| 免费日韩欧美在线观看| 国产在线免费精品| 在线观看美女被高潮喷水网站| 夜夜爽夜夜爽视频| 搡女人真爽免费视频火全软件| 在线免费观看不下载黄p国产| 最近中文字幕2019免费版| 18禁在线播放成人免费| 午夜视频国产福利| 精品久久久久久久久av| 丁香六月天网| 日本vs欧美在线观看视频| 久久久精品区二区三区| 婷婷色麻豆天堂久久| 成人无遮挡网站| 伊人久久精品亚洲午夜| 青春草亚洲视频在线观看| av在线app专区| 日韩中字成人| 99热这里只有是精品在线观看| 国产av码专区亚洲av| 亚洲欧美日韩卡通动漫| 亚洲久久久国产精品| 久久久国产一区二区| 丝瓜视频免费看黄片| 日产精品乱码卡一卡2卡三| 色婷婷久久久亚洲欧美| 色网站视频免费| 黄色怎么调成土黄色| 欧美变态另类bdsm刘玥| 亚洲国产欧美在线一区| 国产一区有黄有色的免费视频| av福利片在线| 亚洲成人一二三区av| 91久久精品电影网| av女优亚洲男人天堂| 久久国产精品大桥未久av| 亚洲精品久久午夜乱码| 久久av网站| 丰满乱子伦码专区| 欧美97在线视频| 九九久久精品国产亚洲av麻豆| 妹子高潮喷水视频| 亚洲精品日韩在线中文字幕| 亚洲国产色片| 永久免费av网站大全| 男男h啪啪无遮挡| 91精品一卡2卡3卡4卡| 啦啦啦在线观看免费高清www| 国产深夜福利视频在线观看| 26uuu在线亚洲综合色| 国产精品一二三区在线看| 91久久精品电影网| 老司机影院成人| 国产免费现黄频在线看| 亚洲成人手机| 看免费成人av毛片| 国产视频首页在线观看| 国产免费又黄又爽又色| 女人精品久久久久毛片| 纯流量卡能插随身wifi吗| 18在线观看网站| 亚洲精品国产av成人精品| 亚洲中文av在线| 制服人妻中文乱码| 日韩在线高清观看一区二区三区| √禁漫天堂资源中文www| 欧美精品高潮呻吟av久久| 日本午夜av视频| 最近的中文字幕免费完整| 制服丝袜香蕉在线| 好男人视频免费观看在线| 99国产精品免费福利视频| √禁漫天堂资源中文www| 一级爰片在线观看| 色婷婷av一区二区三区视频| 丰满饥渴人妻一区二区三| 日韩一本色道免费dvd| 美女国产高潮福利片在线看| 肉色欧美久久久久久久蜜桃| videossex国产| 少妇高潮的动态图| 亚洲婷婷狠狠爱综合网| 成人免费观看视频高清| 王馨瑶露胸无遮挡在线观看| 免费高清在线观看日韩| 精品国产国语对白av| 不卡视频在线观看欧美| 最近中文字幕高清免费大全6| a级片在线免费高清观看视频| 又粗又硬又长又爽又黄的视频| 乱人伦中国视频| 观看av在线不卡| 妹子高潮喷水视频| 一区在线观看完整版| 精品卡一卡二卡四卡免费| 午夜视频国产福利| av一本久久久久| 内地一区二区视频在线| 欧美日韩成人在线一区二区| 少妇熟女欧美另类| 2022亚洲国产成人精品| 国产精品久久久久成人av| 伦理电影大哥的女人| 精品少妇久久久久久888优播| 成人午夜精彩视频在线观看| 激情五月婷婷亚洲| 午夜福利,免费看| 在线观看国产h片| av专区在线播放| 男女边吃奶边做爰视频| 在线观看免费视频网站a站| 五月玫瑰六月丁香| 国产乱人偷精品视频| 九九久久精品国产亚洲av麻豆| 黑人巨大精品欧美一区二区蜜桃 | 日日摸夜夜添夜夜爱| 伦精品一区二区三区| 国产综合精华液| 麻豆成人av视频| 日韩成人av中文字幕在线观看| 人人妻人人爽人人添夜夜欢视频| 久久97久久精品| 亚洲久久久国产精品| 久久99热这里只频精品6学生| 久久久久久久久久成人| 一二三四中文在线观看免费高清| 精品酒店卫生间| av卡一久久| 亚洲av二区三区四区| 精品少妇内射三级| 人体艺术视频欧美日本| 亚洲精品日本国产第一区| 欧美日韩一区二区视频在线观看视频在线| 欧美 亚洲 国产 日韩一| 欧美三级亚洲精品| 曰老女人黄片| 久久热精品热| 啦啦啦啦在线视频资源| 美女主播在线视频| 精品99又大又爽又粗少妇毛片| 最新中文字幕久久久久| 午夜精品国产一区二区电影| 人人妻人人澡人人看| 久久久久久久久久久丰满| 亚洲av成人精品一二三区| 免费看av在线观看网站| 汤姆久久久久久久影院中文字幕| 亚洲精品中文字幕在线视频| 亚洲美女视频黄频| 免费黄频网站在线观看国产| 一级a做视频免费观看| 久热这里只有精品99| 午夜久久久在线观看| 午夜免费观看性视频| 一级毛片电影观看| 在线播放无遮挡| 青青草视频在线视频观看| 亚洲国产日韩一区二区| 18禁在线无遮挡免费观看视频| 黑丝袜美女国产一区| av电影中文网址| 在线看a的网站| 老司机影院成人| 久久99一区二区三区| 欧美日韩av久久| 免费高清在线观看视频在线观看| 亚洲精品456在线播放app| 日本av免费视频播放| 欧美亚洲日本最大视频资源| 美女主播在线视频| 国产亚洲av片在线观看秒播厂| 日韩三级伦理在线观看| 精品视频人人做人人爽| a级毛片免费高清观看在线播放| 成人午夜精彩视频在线观看| 亚洲精品国产色婷婷电影| 少妇的逼水好多| 成人影院久久| av有码第一页| 丝袜美足系列| 亚洲欧美精品自产自拍| 97精品久久久久久久久久精品| 欧美精品一区二区免费开放| 如何舔出高潮| 97超视频在线观看视频| 人妻 亚洲 视频| 国产亚洲最大av| 高清视频免费观看一区二区| 毛片一级片免费看久久久久| 日本黄大片高清| 欧美精品亚洲一区二区| 欧美精品一区二区免费开放| 黑人巨大精品欧美一区二区蜜桃 | 一本大道久久a久久精品| 免费观看av网站的网址| 久久精品国产亚洲网站| 91aial.com中文字幕在线观看| 中文字幕制服av| 搡女人真爽免费视频火全软件| 久久久国产欧美日韩av| 97在线视频观看| 国产伦理片在线播放av一区| 啦啦啦中文免费视频观看日本| 欧美日韩综合久久久久久| 热99国产精品久久久久久7| 91午夜精品亚洲一区二区三区| 久久精品久久久久久久性| 黄片无遮挡物在线观看| 蜜桃久久精品国产亚洲av| 一本色道久久久久久精品综合| 一级毛片aaaaaa免费看小| 国产精品免费大片| 一个人看视频在线观看www免费| 色视频在线一区二区三区| 亚洲不卡免费看| 韩国高清视频一区二区三区| 国产精品久久久久久久电影| 一区二区三区免费毛片| 免费av不卡在线播放| 精品国产国语对白av| 免费av不卡在线播放| 日日摸夜夜添夜夜爱| 看十八女毛片水多多多| 黄色怎么调成土黄色| 少妇被粗大的猛进出69影院 | 免费黄网站久久成人精品| 亚洲精品乱码久久久久久按摩| 一本大道久久a久久精品| 亚洲av二区三区四区| 2021少妇久久久久久久久久久| 久久久国产精品麻豆| 亚洲欧美清纯卡通| 国产精品国产三级专区第一集| 女的被弄到高潮叫床怎么办| 七月丁香在线播放| 国产不卡av网站在线观看| 亚洲色图综合在线观看| 亚洲av日韩在线播放| 三级国产精品片| 最新中文字幕久久久久| 欧美97在线视频| 国产日韩一区二区三区精品不卡 | 欧美国产精品一级二级三级| 欧美激情国产日韩精品一区| 毛片一级片免费看久久久久| 午夜福利网站1000一区二区三区| 青春草国产在线视频| 纵有疾风起免费观看全集完整版| 老司机亚洲免费影院| 老司机影院成人| 欧美人与善性xxx| 亚洲精品色激情综合| 如何舔出高潮| 亚洲精品日韩av片在线观看| 久久精品国产亚洲网站| 久久久亚洲精品成人影院| 高清欧美精品videossex| 国产成人精品福利久久| 日韩成人av中文字幕在线观看| 免费看不卡的av| 激情五月婷婷亚洲| 国产不卡av网站在线观看| 国产视频首页在线观看| av视频免费观看在线观看| 国产成人精品婷婷| 日韩视频在线欧美| 日本猛色少妇xxxxx猛交久久| 青春草视频在线免费观看| 国产成人精品婷婷| 欧美一级a爱片免费观看看| 精品久久蜜臀av无| 国产成人精品一,二区| 黑人高潮一二区| 久久国产精品男人的天堂亚洲 | 亚洲av国产av综合av卡| 亚洲成人手机| 九草在线视频观看| av在线播放精品| 精品国产露脸久久av麻豆| 美女国产高潮福利片在线看| 午夜av观看不卡| videos熟女内射| 免费av中文字幕在线| 日韩,欧美,国产一区二区三区| 亚洲欧美清纯卡通| 2022亚洲国产成人精品| 蜜桃久久精品国产亚洲av| 国产精品免费大片| 精品酒店卫生间| 亚洲欧美中文字幕日韩二区| 哪个播放器可以免费观看大片| 丁香六月天网| 一级毛片电影观看| 欧美最新免费一区二区三区| 久久久久国产精品人妻一区二区| 又大又黄又爽视频免费| 免费人成在线观看视频色| 天美传媒精品一区二区| 狠狠婷婷综合久久久久久88av| 国产精品一二三区在线看| 色网站视频免费| 麻豆乱淫一区二区| 少妇人妻久久综合中文| 蜜桃在线观看..| 亚洲婷婷狠狠爱综合网| 久久免费观看电影| 久久久久久久国产电影| 国产精品欧美亚洲77777| 久久精品人人爽人人爽视色| 成年人午夜在线观看视频| 亚洲熟女精品中文字幕| 中文字幕久久专区| 卡戴珊不雅视频在线播放| 插逼视频在线观看| 99热全是精品| 日韩大片免费观看网站| 欧美日本中文国产一区发布| 五月天丁香电影| 亚洲天堂av无毛| 国产精品麻豆人妻色哟哟久久| 国产日韩一区二区三区精品不卡 | 亚洲成人手机| 亚洲av中文av极速乱| 内地一区二区视频在线| 亚洲精品国产色婷婷电影| 欧美日韩一区二区视频在线观看视频在线| 51国产日韩欧美| 免费播放大片免费观看视频在线观看| 成人午夜精彩视频在线观看| 久久午夜福利片| 久久99精品国语久久久| 国产一区有黄有色的免费视频| 久久午夜综合久久蜜桃| 久久人妻熟女aⅴ| 三上悠亚av全集在线观看| 97在线视频观看| 久久久国产精品麻豆| 一级黄片播放器| av又黄又爽大尺度在线免费看| 欧美人与善性xxx| 日本欧美国产在线视频| 26uuu在线亚洲综合色| 国产又色又爽无遮挡免| 久久午夜福利片| 超色免费av| 久久久a久久爽久久v久久| 国精品久久久久久国模美| 91午夜精品亚洲一区二区三区| 久久精品国产a三级三级三级| 精品一区在线观看国产| 亚洲欧美精品自产自拍| 18禁观看日本| 国产精品免费大片| 人人妻人人澡人人爽人人夜夜| 欧美3d第一页| 亚洲av福利一区| 大香蕉久久网| 日本黄色片子视频| 亚洲国产欧美在线一区|