下肢微重力模擬外骨骼系統(tǒng)的研究

吳岳橋, 梁國(guó)星, 李志利 , 戴小倩, 李文炯, 劉書(shū)娟, 張劍峰

(1.太原理工大學(xué)機(jī)械與運(yùn)載工程學(xué)院, 太原 030024; 2.中國(guó)航天員科研訓(xùn)練中心航天醫(yī)學(xué)基礎(chǔ)與應(yīng)用國(guó)家重點(diǎn)實(shí)驗(yàn)室, 北京 100094)

1 引言

長(zhǎng)期載人航天飛行過(guò)程中,太空微重力環(huán)境會(huì)對(duì)人體的心血管系統(tǒng)、骨骼肌肉系統(tǒng)等產(chǎn)生多種不良影響,影響航天員在軌健康和安全[1]。為確保航天員在太空的安全健康和任務(wù)的順利完成,國(guó)內(nèi)外航天機(jī)構(gòu)研究建立了多種地面微重力環(huán)境模擬方法與系統(tǒng),用于開(kāi)展科學(xué)研究和航天員訓(xùn)練。例如:拋物線飛行法、落塔法、氣浮法、水浮法、懸吊法以及外骨骼法,其中懸吊法出現(xiàn)的時(shí)間較早,是一種成熟且應(yīng)用廣泛的模擬方法。外骨骼法是一種新型的微重力模擬方法,目前仍在探索研究,尚未成熟應(yīng)用。

懸吊法可分為斜面懸吊法和垂直懸吊法。垂直懸吊法是通過(guò)豎直懸吊人體軀干降低足底壓力,從而達(dá)到模擬微重力的足底力學(xué)效果。該類(lèi)系統(tǒng)具有操作簡(jiǎn)單、易于實(shí)施和維護(hù)的優(yōu)點(diǎn)[2-3],但是該類(lèi)型系統(tǒng)無(wú)法補(bǔ)償人體四肢的重力載荷,受試者運(yùn)動(dòng)時(shí)肩、肘、髖、膝等大關(guān)節(jié)仍承受著1 G 的重力,系統(tǒng)無(wú)法模擬四肢的微重力環(huán)境受力狀態(tài)。

外骨骼法利用被動(dòng)重力平衡技術(shù)平衡人體的重力,其突出優(yōu)點(diǎn)是能夠模擬人體四肢的微重力環(huán)境受力狀態(tài)[4-5]。被動(dòng)重力平衡技術(shù)是指利用機(jī)械裝置平衡目標(biāo)系統(tǒng)的重力,而不使用任何的動(dòng)力組件和控制系統(tǒng)。在外骨骼法系統(tǒng)中,用于平衡人體四肢重力的外骨骼屬于串聯(lián)機(jī)器人,其常用的被動(dòng)重力平衡方法有配重法、純彈簧法和彈簧機(jī)構(gòu)法[6-8]。配重法是通過(guò)添加等量的配重,使機(jī)械臂的質(zhì)心位置保持不變,實(shí)現(xiàn)機(jī)器人的重力平衡。采用配重法的系統(tǒng),運(yùn)動(dòng)時(shí)會(huì)產(chǎn)生較大的慣性力,影響微重力環(huán)境模擬效果,因此該方法不適合應(yīng)用于微重力環(huán)境模擬系統(tǒng)。純彈簧法是通過(guò)在機(jī)械臂的適當(dāng)位置添加彈簧,依靠彈簧力耦合實(shí)現(xiàn)機(jī)器人的重力平衡。純彈簧法具有慣性小、結(jié)構(gòu)簡(jiǎn)單的優(yōu)點(diǎn),新墨西哥州立大學(xué)和南京航空航天大學(xué)研究建立的微重力環(huán)境模擬系統(tǒng),都采用了純彈簧法補(bǔ)償人體下肢重力。系統(tǒng)通過(guò)協(xié)同調(diào)整彈簧的安裝位置和彈性系數(shù),實(shí)現(xiàn)對(duì)重力補(bǔ)償水平的調(diào)節(jié)控制[9-11]。彈簧機(jī)構(gòu)法可分為凸輪機(jī)構(gòu)法、齒輪機(jī)構(gòu)法和輔助桿件法,三者的共同原理是通過(guò)彈簧機(jī)構(gòu)運(yùn)動(dòng)調(diào)節(jié)彈簧的變形量,從而精確平衡機(jī)器人的重力。彈簧機(jī)構(gòu)法的設(shè)計(jì)難度較大且安裝精度要求高,但同時(shí)也具有調(diào)節(jié)精度高、易于調(diào)控等優(yōu)點(diǎn)[12-13]。

本文選擇采用彈簧機(jī)構(gòu)法補(bǔ)償人體下肢的重力,設(shè)計(jì)了一種帶有輔助桿件的下肢外骨骼系統(tǒng),用于解決垂直懸吊式系統(tǒng)無(wú)法補(bǔ)償人體四肢重力載荷的問(wèn)題,提高垂直懸吊式系統(tǒng)的微重力環(huán)境模擬真實(shí)度。新型外骨骼通過(guò)協(xié)同調(diào)節(jié)輔助桿件的安裝位置、彈簧的彈性系數(shù)和位置參數(shù),實(shí)現(xiàn)對(duì)重力補(bǔ)償水平的調(diào)節(jié)控制,具有調(diào)控范圍大、調(diào)控方便、補(bǔ)償效果好的優(yōu)點(diǎn)。本文對(duì)所設(shè)計(jì)的系統(tǒng)進(jìn)行了重力補(bǔ)償推導(dǎo)計(jì)算,并從動(dòng)力學(xué)角度研究分析了重力平衡對(duì)人體下肢系統(tǒng)動(dòng)力學(xué)性能的影響。

2 外骨骼系統(tǒng)概念設(shè)計(jì)

2.1 外骨骼概念設(shè)計(jì)

圖1 為改進(jìn)后安裝有下肢外骨骼系統(tǒng)的垂直懸吊式微重力模擬系統(tǒng),主要由恒拉力懸吊系統(tǒng)、上身束縛系統(tǒng)、下肢外骨骼系統(tǒng)和跑步機(jī)等組成。圖2 為帶有輔助桿件的新型下肢微重力模擬外骨骼系統(tǒng)。下肢外骨骼由腰部連桿、大腿桿、小腿桿、輔助連桿、髖關(guān)節(jié)、膝關(guān)節(jié)和彈簧組成,所有連桿均為伸縮桿,可根據(jù)佩戴者腿部的幾何形狀和慣性參數(shù)進(jìn)行調(diào)節(jié)。輔助連桿將2 根彈簧的張力轉(zhuǎn)換為髖關(guān)節(jié)和膝關(guān)節(jié)的平衡扭矩,以克服重力的影響。通過(guò)調(diào)節(jié)彈簧輔助桿件機(jī)構(gòu),實(shí)現(xiàn)任意重力補(bǔ)償水平的精確調(diào)節(jié)控制,滿足不同體重的受試者,模擬不同微重力環(huán)境狀態(tài)。下肢外骨骼在踝關(guān)節(jié)之前結(jié)束,對(duì)踝關(guān)節(jié)沒(méi)有限制,髖關(guān)節(jié)有2 個(gè)自由度、膝關(guān)節(jié)有1 個(gè)自由度,能夠?qū)崿F(xiàn)髖關(guān)節(jié)的屈伸、外展、內(nèi)收和膝關(guān)節(jié)的屈伸[12-14]。

圖1 微重力模擬系統(tǒng)示意圖Fig.1 Schematics of low gravity simulation system

圖2 被動(dòng)外骨骼概念設(shè)計(jì)圖Fig.2 Conceptual design of passive exoskeleton

2.2 重力補(bǔ)償原理

被動(dòng)重力平衡的設(shè)計(jì)主要有力矩對(duì)等和能量守恒2 種方式[4-7],本文選擇能量守恒的方法以過(guò)髖關(guān)節(jié)的水平面為零勢(shì)能面。由外骨骼和人體下肢組成的人機(jī)系統(tǒng),其重力勢(shì)能會(huì)隨著髖關(guān)節(jié)和膝關(guān)節(jié)角度變化而變化,彈簧的作用是使系統(tǒng)的重力勢(shì)能和彈簧的彈性勢(shì)能相互轉(zhuǎn)換,系統(tǒng)的總勢(shì)能在任意位姿保持恒定。當(dāng)人機(jī)系統(tǒng)總勢(shì)能保持不變時(shí),系統(tǒng)處于重力平衡狀態(tài),人體下肢重力得到補(bǔ)償,其數(shù)學(xué)表達(dá)式見(jiàn)式(1)。

式中,UToal、ULG、UMG、US分別表示系統(tǒng)總勢(shì)能、人體下肢重力勢(shì)能、機(jī)構(gòu)重力勢(shì)能、彈簧彈性勢(shì)能。式(1)表示完全補(bǔ)償重力的失重狀態(tài),如果需要模擬月球、火星等微重力狀態(tài),則如式(2)所示:

式中,ρ為比例系數(shù),表示被補(bǔ)償?shù)南轮亓εc原始下肢重力的比值。例如,月球的重力是地球重力的1/6,因此模擬月球的微重力環(huán)境時(shí),需要平衡5/6 的重力,此時(shí)ρ為5/6。

3 重力平衡計(jì)算

系統(tǒng)重力平衡設(shè)計(jì)的關(guān)鍵在于如何確定彈簧的參數(shù)和安裝位置。本文基于能量守恒使用剛度矩陣法[15-16]推導(dǎo)重力補(bǔ)償被動(dòng)外骨骼的設(shè)計(jì)公式,確定彈簧輔助桿件機(jī)構(gòu)的幾何參數(shù)。下肢系統(tǒng)數(shù)學(xué)模型如圖3 所示。

圖3 下肢人機(jī)系統(tǒng)數(shù)學(xué)模型Fig.3 Mathematical model of lower limb man-machine system

圖中,m1、m2、m3、m4分別表示大腿段AB(大腿和大腿連桿)、小腿段BC(小腿和小腿連桿)、輔助連桿DE、輔助連桿EF的質(zhì)量。mp1、mp2分別表示髖關(guān)節(jié)A、膝關(guān)節(jié)B的質(zhì)量,mp3表示踝關(guān)節(jié)和腳的質(zhì)量。l1、l2分別表示大腿段AB、小腿段BC的長(zhǎng)度。s1表示髖關(guān)節(jié)A到鉸接點(diǎn)D的距離,s2表示膝關(guān)節(jié)B到鉸接點(diǎn)F的距離。r1表示大腿段的質(zhì)心到髖關(guān)節(jié)A的距離,r2表示小腿段質(zhì)心到膝關(guān)節(jié)B的距離,r3表示輔助連桿EF質(zhì)心到鉸接點(diǎn)E的距離,r4表示輔助連桿DE的質(zhì)心到鉸接點(diǎn)D的距離。k1、k2表示2 個(gè)彈簧的彈性系數(shù),u1表示彈簧k1連接點(diǎn)G到髖關(guān)節(jié)A的距離,u2表示表彈簧k2連接點(diǎn)H到鉸接點(diǎn)D的距離。q1表示沿連桿AG的單位方向向量,q2表示沿連桿AB的單位方向向量,q3表示沿連桿BC的單位方向向量。

首先計(jì)算系統(tǒng)重力勢(shì)能的剛度矩陣,整個(gè)系統(tǒng)的質(zhì)心向量如式(3)所示。

式中,M表示系統(tǒng)總質(zhì)量,C表系統(tǒng)的質(zhì)心向量,C可表示為式(4)。

式中,K2∈R6×6為分塊矩陣,I為2× 2 的單位矩陣,K2即為系統(tǒng)彈性勢(shì)能的剛度矩陣。最后計(jì)算系統(tǒng)總勢(shì)能的剛度矩陣,根據(jù)式(5)、(7)系統(tǒng)的總勢(shì)能可表示為式(9)。

K即為系統(tǒng)總勢(shì)能的剛度矩陣,當(dāng)剛度矩陣的非對(duì)角線分量全為零時(shí),系統(tǒng)總勢(shì)能守恒系統(tǒng)處于重力平衡狀態(tài)[6,15-16],如式(10)所示。

由式(10)可知,設(shè)計(jì)方程中的變量有輔助桿件的安裝位置u1、u2,彈簧的彈性系數(shù)k1、k2以及位置參數(shù)s1、s2。由于在實(shí)際應(yīng)用中很難獲取特定彈性系數(shù)的彈簧,因此可以通過(guò)協(xié)同調(diào)節(jié)輔助連桿的安裝位置、彈簧的彈性系數(shù)和位置參數(shù),滿足不同體重的受試者和模擬不同微重力環(huán)境狀態(tài)的需求。

4 動(dòng)力學(xué)分析

下肢系統(tǒng)屬于多關(guān)節(jié)空間連桿機(jī)構(gòu),故將系統(tǒng)簡(jiǎn)化為多連桿模型。下肢的運(yùn)動(dòng)主要集中于矢狀面內(nèi),為了方便計(jì)算在建模時(shí)只考慮在矢狀面內(nèi)的運(yùn)動(dòng),忽略下肢在額狀面和水平面的運(yùn)動(dòng)。下肢運(yùn)動(dòng)的一個(gè)步態(tài)周期可分為支撐相和擺動(dòng)相2 個(gè)階段,由于本文主要目的是研究重力補(bǔ)償對(duì)系統(tǒng)動(dòng)力學(xué)性能的影響,因此本文只建立單腿支撐的擺動(dòng)相動(dòng)力學(xué)模型[17-18]。圖4 為矢狀面內(nèi)擺動(dòng)相階段的多連桿動(dòng)力學(xué)模型。

圖4 中,θ1表示大腿段AB與豎直方向的夾角,θ2表示大腿段AB與小腿段BF的夾角,Q1、Q2分別表示髖關(guān)節(jié)A、膝關(guān)節(jié)B的關(guān)節(jié)力矩,其余參數(shù)與前文數(shù)學(xué)模型參數(shù)定義相同。利用拉格朗日方程進(jìn)行動(dòng)力學(xué)分析,下肢系統(tǒng)連桿模型的拉格朗日方程為式(11)。

圖4 下肢人機(jī)系統(tǒng)動(dòng)力學(xué)模型Fig.4 Dynamic model of lower limb man-machine system

式中,EK表示人機(jī)系統(tǒng)總動(dòng)能,EP表示人機(jī)系統(tǒng)的總勢(shì)能,計(jì)算各連桿質(zhì)心的位置和速度。

桿AB的質(zhì)心坐標(biāo)(x1,y1) 為式(12)。

將式 (14)、(15)代入式(11)整理得系統(tǒng)的拉格朗日方程。根據(jù)拉格朗日方程,計(jì)算下肢系統(tǒng)的動(dòng)力學(xué)方程,系統(tǒng)動(dòng)力學(xué)方程為式(16)。

將式(17)、(18)代入式(16)中,整理可得人機(jī)系統(tǒng)髖關(guān)節(jié)和膝關(guān)節(jié)的關(guān)節(jié)力矩,為了方便研究分析將結(jié)果寫(xiě)成矩陣形式,如式(19)。

分析計(jì)算結(jié)果可知,關(guān)節(jié)力矩與重力項(xiàng)矩陣和彈力項(xiàng)矩陣有關(guān),選擇合適的彈簧參數(shù)能夠抵消或者減小重力項(xiàng)矩陣,從而減小髖關(guān)節(jié)和膝關(guān)節(jié)的關(guān)節(jié)力矩,實(shí)現(xiàn)人體下肢的微重力環(huán)境模擬。

5 仿真驗(yàn)證

根據(jù)國(guó)家標(biāo)準(zhǔn)GB/T17245-2004《成年人人體慣性參數(shù)》和GB 1000-88《中國(guó)成年人人體尺寸》確定人體下肢和外骨骼的幾何和慣性參數(shù)。根據(jù)上文推導(dǎo)的設(shè)計(jì)公式(10)確定彈簧輔助桿件機(jī)構(gòu)的幾何參數(shù)。在動(dòng)力學(xué)建模的基礎(chǔ)上利用三維建模軟件SOLIDWORKS 構(gòu)建人機(jī)系統(tǒng)下肢模型,并將建立好的模型導(dǎo)入ADAMS 中添加約束,不添加關(guān)節(jié)驅(qū)動(dòng),設(shè)置系統(tǒng)重力為9.8 N/kg,設(shè)置模型的參數(shù)見(jiàn)表1,設(shè)置仿真時(shí)間為10 s,step步數(shù)為1000 步。將仿真模型初始位姿設(shè)置為提膝狀態(tài),仿真驗(yàn)證外骨骼系統(tǒng)能否平衡下肢重力,仿真結(jié)果如圖5 所示。

表1 系統(tǒng)參數(shù)表Table 1 Values of system parameters

圖5 中,虛線代表未添加彈簧時(shí)的仿真結(jié)果,實(shí)線代表添加彈簧時(shí)的仿真結(jié)果。結(jié)果表明,在未添加彈簧力時(shí),受重力影響各部分的質(zhì)心位置會(huì)發(fā)生變化,模型無(wú)法在重力場(chǎng)中保持原始位姿狀態(tài),在添加彈簧力時(shí),各部分的質(zhì)心位置保持靜止,下肢所受重力得到精確補(bǔ)償,模型能在重力場(chǎng)中保持平衡狀態(tài),證明外骨骼系統(tǒng)系統(tǒng)能夠平衡下肢重力。

圖5 質(zhì)心位置變化曲線Fig.5 Change curve of centroid position

依據(jù)文獻(xiàn)[19],確定髖關(guān)節(jié)和膝關(guān)節(jié)在一個(gè)步態(tài)周期內(nèi)的角度變化曲線。利用ADAMS 的STEP 函數(shù)編寫(xiě)髖關(guān)節(jié)和膝關(guān)節(jié)的驅(qū)動(dòng)函數(shù),設(shè)置系統(tǒng)重力為9.8 N/kg,設(shè)置仿真時(shí)間為5 s,step步數(shù)為500 步。

圖6 為模擬不同重力環(huán)境下,下肢擺動(dòng)時(shí)髖關(guān)節(jié)和膝關(guān)節(jié)的關(guān)節(jié)力矩變換曲線,圖6(a)代表未平衡系統(tǒng)重力,模擬地球重力環(huán)境的仿真結(jié)果,圖6(b)代表平衡系統(tǒng)部分重力,模擬月球和火星低重環(huán)境的仿真結(jié)果,圖6(c)代表平衡系統(tǒng)全部重力,模擬太空失重環(huán)境的仿真結(jié)果。結(jié)果表明,隨著重力補(bǔ)償水平即彈簧力的不斷增加,下肢的關(guān)節(jié)力矩逐漸減小,證明彈簧彈力補(bǔ)償了下肢重力,減小或消除了重力對(duì)下肢運(yùn)動(dòng)的影響,與本文動(dòng)力學(xué)分析結(jié)果具有一致性。并且仿真結(jié)果證明在外骨骼模擬的低重力環(huán)境和實(shí)際低重力環(huán)境中下肢運(yùn)動(dòng)的生物力學(xué)特性具有相似性,外骨骼能夠模擬真實(shí)低重力環(huán)境[10,17-18]。

圖6 系統(tǒng)關(guān)節(jié)力矩變化曲線Fig.6 Curve of system joint torque

圖7 為失重狀態(tài)下,系統(tǒng)勢(shì)能的變化曲線,結(jié)果表明,在仿真運(yùn)動(dòng)過(guò)程中,系統(tǒng)的重力勢(shì)能不斷下降,彈性勢(shì)能在不斷上升,系統(tǒng)的總勢(shì)能保持恒定不變,證明能夠?qū)崿F(xiàn)重力勢(shì)能和彈性之間的相互轉(zhuǎn)換,維持系統(tǒng)的總勢(shì)能保持不變[20],結(jié)果與本文重力補(bǔ)償原理一致,驗(yàn)證了重力補(bǔ)償原理的可行性和設(shè)計(jì)公式的正確性。

圖7 系統(tǒng)勢(shì)能變化曲線Fig.7 Curve of system potential energy

6 結(jié)論

本文針對(duì)懸吊式微重力環(huán)境模擬系統(tǒng)無(wú)法模擬人體下肢低重力環(huán)境的問(wèn)題,提出利用外骨骼機(jī)器人模擬下肢低重力環(huán)境的解決方案。采用被動(dòng)重力平衡的方法設(shè)計(jì)了外骨骼,以實(shí)現(xiàn)模擬下肢的低重力環(huán)境受力狀態(tài),并利用剛度矩陣法對(duì)重力平衡條件進(jìn)行了推導(dǎo),得到了外骨骼機(jī)器人重力補(bǔ)償設(shè)計(jì)公式。通過(guò)數(shù)學(xué)分析和計(jì)算機(jī)仿真對(duì)外骨骼的重力補(bǔ)償設(shè)計(jì)方案和設(shè)計(jì)公式的有效性和正確性進(jìn)行驗(yàn)證。結(jié)果表明,采用彈簧輔助桿件法設(shè)計(jì)的外骨骼能夠模擬下肢低重力環(huán)境。