李旭鴻，郝衛(wèi)亞，肖曉飛，吳成亮

(1.浙江體育科學(xué)研究所，浙江 杭州 310004；2.國家體育總局 體育科學(xué)研究所，北京 100061)

·運動生物力學(xué)·

基于ADAMS.LifeMoD人體落地動作的動力學(xué)仿真與驗證

李旭鴻1，郝衛(wèi)亞2，肖曉飛2，吳成亮2

(1.浙江體育科學(xué)研究所，浙江 杭州 310004；2.國家體育總局 體育科學(xué)研究所，北京 100061)

摘要：基于人體運動仿真軟件LifeMoD建立人體個性化的三維多剛體動力學(xué)模型，該模型包括19個環(huán)節(jié)和50個自由度。對受試者的落地過程進行動作捕捉、運動學(xué)解析、動力學(xué)仿真，并完成仿真結(jié)果與三維測力臺的驗證。結(jié)果表明該模型具有很好的可信度，可以用來研究人體落地沖擊過程的動力學(xué)響應(yīng)特性。

關(guān)鍵詞：LifeMoD；地面反作用力；運動捕捉；仿真與驗證

收稿日期：2014-07-09

作者簡介：李旭鴻(1976-)，男，副研究員，博士，主要從事人體運動的生物力學(xué)仿真.

文章編號：1004-3624(2015)01-0114-03

中圖分類號:G804.66

文獻標(biāo)識碼：獻標(biāo)識碼：A

Abstract:A three dimensional multi-body dynamic model of the subject-specific human was developed using the software ADAMS and LifeMoD, and this model include 19 segments and 50 degree of freedom. Jump landing of the subject was captured, kinematic parsed and dynamics simulation. The results was tested and proved to valid and feasible with 3D force-plate. It is shown that the simulation model has a good reliability, and it can be used to study the kinematics and dynamics of the human landing under impact loads.

基金項目：浙江省體育局備戰(zhàn)2012奧運和2013全運會體育科研攻關(guān)課題項目

Simulation and Validation of Human Model in Landings

Based on ADAMS.LifeMoD

LI Xu-hong1, HAO Wei-ya2，XIAO Xiao-fei2,WU Chen-liang2

(1.Zhejiang Research Institute of Sport Science， Hangzhou 310004， China；

2.China Institute of Sport Science, Beijing 100061， China)

Key words:LifeMoD; ground reaction force; motion capture; simulation and validation

0前言

在眾多的體育活動中，不管是與人們?nèi)粘Ｉ钕⑾⑾嚓P(guān)的周期性落地(如走、跑、跳)，還是競技比賽場上的制動性落地(如體操)都面臨著一定的沖擊負(fù)荷。通常情況下，人體利用下肢關(guān)節(jié)屈曲過程中肌肉的主動收縮和骨骼的變形來吸收和衰減承受的沖擊負(fù)荷，進而避免運動損傷[1]。其實人體在體育鍛煉過程中每一次與地面接觸都需要承受1~3倍自身體重的地面反作用力[2]，而在高速跑跳過程中，這一沖擊力有時會高達自身體重的8倍以上。殊不知較大的沖擊力和較高的負(fù)載率(loading rate)有時會超出人體肌肉骨骼系統(tǒng)承受的生理極限，進而增加了下肢損傷的發(fā)生率[3]。所以，過去幾十年間體育科學(xué)和運動醫(yī)學(xué)領(lǐng)域的人們一直認(rèn)為跑跳過程中較高的沖擊載荷可能是導(dǎo)致下肢損傷的主要原因之一[4-6]。然而，近期的研究發(fā)現(xiàn)雖然較大的被動沖擊對人體的肌肉骨骼系統(tǒng)造成潛在的影響，但并不足以成為運動損傷的根本原因[7-8]。 因此現(xiàn)階段，一些專家學(xué)者針對人體落地沖擊方面的生物力學(xué)研究開始不僅僅局限于沖擊力本身，而是進一步考慮沖擊力的作用效果以及人體肌肉骨骼系統(tǒng)對此的綜合反應(yīng)[9]。但目前對人體肌骨系統(tǒng)在沖擊過程中承受內(nèi)在負(fù)荷的研究，通常都是建構(gòu)在人體模型的基礎(chǔ)之上。所以，首先要確定人體模型的合理性和有效性，方能進一步考慮獲取的內(nèi)在負(fù)荷是否真實、可信?；蛟S將理論模型與實驗研究整合在一起，能更好地理解和探索落地沖擊引起人體下肢損傷的生物力學(xué)機制[10]。

本文通過人體運動仿真軟件ADAMS. LifeMoD建立人體模型，利用紅外高速捕捉系統(tǒng)和高速攝像完成不同高度落地動作的采集，隨后計算機仿真獲取落地過程的地面反作用力，且與三維測力臺的數(shù)據(jù)進行比較和相似度研究，驗證基于ADAMS. LifeMoD建立人體模型的有效性。

1研究方法

1.1運動捕捉

受試對象是1名體育院校學(xué)生(24歲，170cm，70kg)，如圖1所示。利用高速紅外運動捕捉系統(tǒng)(Motion Analysis)和2臺高速攝像機(SONY)對其完成的落地跳實驗進行同步采集，兩者的采樣頻率均為200Hz，標(biāo)定框架為PEAK，落地高度分別為40cm和80cm。受試者首先穿鞋進行準(zhǔn)備活動，熟練后赤腳(barefoot)完成實驗且無向上和向前的跳躍動作。2塊三維測力臺(Kistler)置于厚度為20cm的落地墊下面記錄人體落地過程承受的地面反作用力，其采樣頻率為1 000Hz。人體關(guān)節(jié)點的選取如圖2所示，通過解析獲取人體運動學(xué)參數(shù)。

1.2建立模型

圖1　受試者落地實驗示意圖

利用人體運動仿真軟件LifeMoD，創(chuàng)建19個環(huán)節(jié)50個自由度的人體模型[11,12]，通過基于Python腳本語言的接口插件，把三維運動學(xué)數(shù)據(jù)導(dǎo)入人體模型中，并對初始狀態(tài)進行平衡分析和空間姿態(tài)的匹配，最終完成模型的建立。

圖2　人體關(guān)節(jié)點的選取示意圖

1.3仿真流程與簡易優(yōu)化

從數(shù)據(jù)文件中讀取初始位置、載荷和邊界條件，設(shè)置仿真時間與步長，進行逆向動力學(xué)運算。隨后模型在關(guān)節(jié)力矩的驅(qū)動下與地面產(chǎn)生交互作用，獲得人體落地的動力學(xué)響應(yīng)。仿真結(jié)束后，將仿真結(jié)果(GRF)與實驗結(jié)果進行比較，且對兩者的垂直方向上的GRF曲線開展相似度(復(fù)相關(guān)系數(shù)CMC)評價[13]。其中復(fù)相關(guān)系數(shù)(公式1)定義為：

(1)

式中是曲線的條數(shù)，是每條曲線中含有數(shù)據(jù)的個數(shù)，是第條曲線的第個數(shù)據(jù)，是條曲線的第個數(shù)據(jù)的平均值，是條曲線個數(shù)據(jù)的總體均值。最后，如果仿真結(jié)果不合理則需要調(diào)整相關(guān)參數(shù)，進行簡易的優(yōu)化最終得到最優(yōu)的仿真效果[12]。

2結(jié)果與分析

首先，利用高速攝像機和高速紅外運動捕捉系統(tǒng)對人體40cm高度的落地過程進行同步采集，分別展開仿真后得到垂直方向上的GRF(高速攝像機和motion系統(tǒng))并與三維測力臺上記錄的數(shù)據(jù)(Kistler)進行比較和相似度描述，結(jié)果發(fā)現(xiàn)兩者的仿真結(jié)果與實驗測試之間的相似度CMC分別為0.906和0.964。同時，兩者的GRF峰值與測力臺最大差值為13.93N，到達GRF峰值的時間平均相差3ms(圖3所示)。表明兩種攝像系統(tǒng)采集到的運動學(xué)數(shù)據(jù)可以有效的應(yīng)用到人體運動仿真研究中，且能較好地呈現(xiàn)出人體實際的落地動作。

其次，僅使用高速攝像機對人體40cm和80cm高度的落地過程進行運動捕捉，并對仿真獲取垂直方向上的GRF與三維測力臺上記錄的數(shù)據(jù)進行比較與相似度描述。結(jié)果發(fā)現(xiàn)， 40cm和80cm高度上的仿真結(jié)果與測力臺之間的相似度分別為0.935和0.967。在40cm 的高度上仿真的GRF峰值與測力臺的數(shù)據(jù)非常接近，僅相差6.29N，到達峰值的時間比測力臺提前6ms。而在80cm的高度上兩者的GRF峰值相差144.87N，是實測GRF峰值的3.3%，到達峰值的時間卻滯后2ms(圖4所示)。

圖3　兩種采集系統(tǒng)運動學(xué)數(shù)據(jù)仿真與測力臺上的力值

圖4　三維人工解析數(shù)據(jù)仿真與測力臺的力值

最后在運動學(xué)方面的驗證上，課題組還利用Motion計算80cm高度落地過程中膝關(guān)節(jié)(以左側(cè)為例)角度變化與其仿真結(jié)果進行對比，結(jié)果發(fā)現(xiàn)兩者CMC值為0.98。同時，兩者均在0.24s達到最小膝角，且屈曲角度相差8°左右[15]。

綜上所述，無論利用紅外高速運動捕捉系統(tǒng)還是高速攝像機都可以準(zhǔn)確地解析出人體關(guān)節(jié)點的運動軌跡，并應(yīng)用到人體運動仿真研究中，且在運動學(xué)和動力學(xué)上都高度相似。另外，從不同的落地高度上也能較好地反映人體運動過程，從而進一步表明建立的人體模型是合理、可信的。

3結(jié)論與建議

3.1利用高速攝像或高速紅外運動捕捉系統(tǒng)采集到的運動學(xué)參數(shù)均可以較好地應(yīng)用到人體運動仿真研究中，而基于LifeMoD建立19個環(huán)節(jié)的人體模型是可行和合理的。

3.2通過簡易優(yōu)化，結(jié)合外部環(huán)境與人體形成系統(tǒng)模型能較好地反映人體落地動作的實際運動，并對高水平運動員在實際訓(xùn)練或比賽中的GRF進行有效估算。

參考文獻

[1]Nigg BM. Biomechanics, load analysis and sports injuries in the lower extremites[J].Sports Med,1985,2:367-379.

[2]Challis JH and Pain MTG. Soft tissue motion influences skeletal loads during impacts[J].Exerc Sport Sci Rev,2008,36:71-75.

[3]O’Leary K, Vorpahl KA, Heiderscheit B. Effect of cushioned insoles on impact forces during running[J].J Am Podiatr Med Assoc,2008,98(1):36-41.

[4]Clarke TE, Frederick EC, Cooper LB. Effects of shoe cushioning upon ground reaction forces in running[J].Int J Sports Med,1983a,4:247-251.

[5]Clarke TE, Frederick EC, Hamill CL. The effects of shoe design parameters on rearfoot control in running[J].Med Sci Sports Exerc,1983b,15:376-381.

[6]Nigg BM. Impact forces in running [J].Current Opinion in Orthopaedics,1997,8:43-47.

[7]Wakeling JM, Von Tscharner V, Nigg BM, Stergiou P. Muscle activity in the leg is tuned in response to ground reaction forces[J].J Appl Physiol,2001,91:1307-1317.

[8]Boyer KA and Nigg BM. Muscle activity in the leg is tuned in response to impact force charactertics[J].J Biomech,2004,37:1583-1588.

[9]傅維杰.運動鞋和緊身裝備對沖擊力、軟組織振動及肌肉活動影響的研究[D].上海體育學(xué)院，2011.

[10]李旭鴻，郝衛(wèi)亞.落地沖擊引起下肢損傷的生物力學(xué)研究進展[J].中國運動醫(yī)學(xué)雜志，2013,32(2)：186-191.

[11]李旭鴻.基于LifeMoD對跳馬過程中人-器械動力學(xué)關(guān)系的計算機仿真[D].上海體育學(xué)院，2013.

[12]李旭鴻，郝衛(wèi)亞，于佳彬，等.基于LifeMoD對跳馬過程中體操運動員-落地墊動力學(xué)關(guān)系的計算機仿真[J].體育科學(xué)，2013,33(3)：81-87.

[13]李旭鴻，等.基于LifeMoD對跳馬運動員推手動作動力學(xué)的仿真研究[J].中國體育科技，2014,50(6):53-58.

[14]Kadaba MP, Ramakrishnan HK, Wootten ME, et al. Repeatability of kinematic, kinetic, and electromyographic data in normal adult gait[J].J Orthop Res,1989,7(6):849-860.

[15]于佳彬，郝衛(wèi)亞，周興龍.基于LifeMoD系統(tǒng)對縱跳落地動作仿真方法的研究[C].第十五屆全國運動生物力學(xué)學(xué)術(shù)交流大會論文摘要匯編，2012:44.