1. 引言
行走是一种日常生活中常见的人类运动,是人体瞬时动态平衡连续传递的运动过程,这个过程需要肌肉和骨骼的协同作用控制各关节的运动进行实时动态调整人体的质心。在正常行走过程中往往受到外界环境干扰而使步态失稳。据调查,目前导致人体步态失稳最常见的外界干扰是地面湿滑[1]。由于湿滑的外界干扰而引起的跌倒是职业环境中严重的安全问题之一[2]。有文献指出,滑倒是工业部门主要的意外事故,达到职业意外伤害的一半[3]。在社会生活中滑倒造成的损伤也占到家庭事故的20%,公共场所事故的18%。美国国家安全委员会的调查统计显示,2002年仅在美国由于跌倒造成至少12,900名成年人死亡,并且至少1,016,700名因滑倒而造成损伤的成年人接受治疗[4]。
为了制定有效的防滑策略,有必要了解滑动后成功恢复平衡和终止跌倒之间的人体内在本质属性的区别[5]。
本文将下肢骨骼简化为六杆的倒立摆物理模型,通过建立数学模型和计算,得出人体水平直线步道行走时单支撑相水平方向与垂直方向的惯性力表达式,然后再对单支撑相进行受力分析得出摩擦力与人体平衡的关系,最后对身体下肢运动状态的分析,对人体发生滑跌过程中下肢关节力矩和关节角度变化趋势和机理作了详细的研究和分析。试验结果和模型推论相一致。因此,验证了本文所建立模型的有效性。本文的研究结果可为下肢助行装置的研发、防滑跌损伤和康复训练提供一定的参考依据。
2. 步态模型构建
2.1. 步态运动学特性
本文假设人体各个环节都是刚体,人体质心运动轨迹可以看作以右足踝关节为圆心,身体质心m到踝关节的距离为半径l的圆弧,质心位置先升高再降低,其简化模型如图1所示。
图1中,G为竖直向下的重力,Fgt为切向惯性力,Fgn为法向惯性力,N为垂直向上的地面支反力,f为沿X轴负方向的摩擦力。
根据人体行走过程质心的运动轨迹,本文将人体平衡站立假设为倒立摆模型[6] [7]。下面运用数学模型对此进行讨论分析。在此运动过程中,质心在圆弧轨迹上切向受力为Gsinθ,θ为人体质心与足底压力中心的连线与铅垂线的夹角。质心的切向加速度为at,根据牛顿第二定律F = ma可以得出
,下面对此式进行计算。式中at为:
Figure 1. Simplified model of centroid motion in walking
图1. 行走过程质心运动简化模型图
(1)
所以得到人体质心运动的微分方程为:
(2)
取θ的最大值
,r为步长的一半,一般成年人的步长为(0.5~0.8 m),中国人的平均身高为1.67 m [8]。所以θ的最大值为
,当
时误差小于1%,所以这里可以认为
。所以式(2)可以写为:
(3)
计算此方程得其通解为:
(4)
两个边界条件分别看作为
时
,
时
,并且令
,
。代入上市可以得到常数项为:
,
(5)
再把常数项A1和A2带回式(4)中求出θ的表达式为:
(6)
对式(6)进行简化,由双曲正弦函数
,因此可简化为:
(7)
2.2. 步态动力学特性
人体步行运动是一个连续的动态平衡过程,可以根据达朗贝尔原理将动力学问题用静力学的方法进行求解。在运动的任一瞬间可以在物体上虚加一个惯性力,与主动力和约束力形成一个平衡力系。设质心质量为m,加速度为a,作用在质心上的主动力为F,约束力为FN,由牛顿第二定律有
(8)
上式中的−ma就是惯性力,惯性力的大小为质心质量与加速度的乘积,方向与加速度方向相反。因此人体在步行过程中受到的外力有重力、地面摩擦力、地面支反力、空气阻力及虚加的惯性力。由于试验环境中空气流速小,受试者步行速度慢,因此空气阻力的数值很小,在本文的分析中可以忽略不计。根据达朗贝尔原理在人体步行的过程中的静力学受力表达式为
(9)
式(9) G为重力,FN为地面支反力,f为地面摩擦力,Fg为惯性力。惯性力Fg可分解为切向惯性力和法向惯性力,切向惯性力为Fgt = mat,法向惯性力为Fgn = man。
切向加速度:
(10)
法向加速度:
(11)
求得切向加速度与法向加速度后便可求得切向惯性力与法向惯性力,其大小为:
切向惯性力:
(12)
法向惯性力:
(13)
本文要研究人体步行过程中矢状面的运动,要将惯性力分解为X方向与Z方向的分力,X方向的惯性力为Fgx,Z方向的惯性力为Fgz。
(14)
(15)
从(14)和(15)两式中可以得出Fgx和Fgz均与θ成正相关。
2.3. 滑跌步态动力学分析
本文要研究人体在行走时的滑跌步态,因此要先分析在滑跌时的受力情况。在本文研究所选的人体步态的单支撑相内,受力情况如图2所示。以质心O为坐标系原点建立坐标系XOZ,X轴的正方向为人体运动方向,Z轴的正方向为垂直于X轴向上。
Figure 2. Force analysis of single support phase
图2. 单支撑相受力分析图
人体所受到的外力如图2所示:竖直向下的重力G,垂直向上的惯性力Fgz,沿X轴负方向的惯性力Fgx,垂直向上的地面支反力FN,沿X轴正方向的摩擦力f。要保持身体的平衡,这一阶段人体所受的外力和力矩均要处于平衡,因此可得力学方程如下:
(16)
在人体步行过程中若不发生打滑,则要求式中的摩擦力f小于地面与受试者足部接触面所提供的最大静摩擦力fmax。因此平衡行走所需的摩擦力f越小,接触面提供的静摩擦力越大,则行走发生打滑情况的可能性就越小,否则发生打滑情况的可行性就会越大。由于
,因此在行走过程中惯性力越大或加速度越大,则需要的摩擦力越大,发生打滑情况的可能性越大。
(17)
式(17)中的RCOF是在人体行走过程中不出现打滑情况时,摩擦力f与正压力FN的比值,也就是安全行走所需要的摩擦系数。由RCOF = tanθ,可知此摩擦系数的大小取决于人体质心与足底压力中心的连线与铅垂线的夹角θ,其中θ角越大则所需的摩擦力越大,也就越容易出现打滑现象。
2.4. 滑跌步态模型构建
本文研究内容为人体下肢,人体下肢由髋骨、股骨、髌骨、腓骨、胫骨及足骨构成的框架,并且骨骼之间形成三大关节:髋关节Hip joint (H)、膝关节Knee joint (K)、踝关节Ankle joint (A)。各个骨骼与关节由肌肉连接,由肌肉的收缩驱动关节的转动使身体运动。对于简化的下肢多杆运动学模型,关节的内驱动力矩对整个步态模型来说为内力,其外在主动力为测力台或者步道的反作用力及自身重力。
Figure 3. Gait model
图3. 步态模型图
由于身体在进入一个步态周期时足跟触地时身体为制动状态,加速度向后惯性力向前,水平方向忽略其他阻力后只有摩擦力。所以摩擦力在触地时是向后的,这时若地面所提供的摩擦力小于平衡行走所需的摩擦力,足部便向前滑动身体相对向后倾倒[9]。所以以下的分析均假设身体发生打滑时为相对踝关节向后倾倒。如图3所示的步态模型,将运动方向为右设为正方向,分析各关节运动对身体质心位置的改变情况。设定RA点为定点相对于地面为静止不动,模型便可以看为绕定点转动的多杆机构,杆RK-RA为主动杆,通过其运动参数可以计算其他杆件的运动情况。设其有初始角速度ωra,角加速度为αra。然后分析H点的运动情况,H点为髋关节的简化点,髋关节与身体质心的相对运动较小,因此H点的运动情况可以近似看作身体质心的运动。
首先忽略RK点与H点的驱动力矩对右膝关节的运动进行分析,RK-RA为绕定点转动的刚体长度为l1,RK点的切向速度为vrkt,切向加速度为arkt,法向加速度为arkn。对这三个参数求解为:切向速度
,切向加速度
,法向加速度
。现在假设MRA的方向为顺时针方向,根据RK点的运动学参数求此时的H点的运动学参数,对H点的速度分析如图4所示。
Figure 4. Velocity analysis of hip joint
图4. 髋关节速度分析图
如图4所示,可以运用刚体在平面运动基点法将H点随RK的平移与绕RK转动的速度合成H点的运动速度vh的值。其大小为:
(18)
vh的大小与方向会出现图4中A、B、C三种情况,本文分开进行讨论。A情况时,H点相对于RK点做逆时针方向的转动,但是其产生的向左的速度分量不足以改变vh正负,H点的运动趋势仍然是向右方向,为正方向。B情况时,H点相对于RK做顺时针方向的转动,H点的运动趋势是向右方向,为正方向。C情况时,H点相对于RK点做逆时针方向的转动,但是其产生的向左的速度分量可以改变vh正负,H点的运动趋势是向左方向,为负方向。C情况髋关节H带动身体质心向左运动,会使人体质心与足底压力中心的连线与铅垂线的夹角θ变大,使RCOF增大,因此人体容易滑倒的可能性增大。如假设MRA的方向为逆时针方向,由于人体的膝关节的生理特征限制H-RK杆不能相对RK-RA杆向右弯曲,这就说明RK-RA杆向左转动时会使θ角变大,同样是发生打滑的不稳定现象。
下面对H点的加速度进行分析,分析由于身体的内力作用改变运动状态对加速度有何影响,从而推论出对人体行走时惯性力与所需摩擦力的影响。这里分析MRA的方向为顺时针方向时H点的加速度,如图5所示。
ah的大小与方向会出现图5中A、B、C三种情况,本文分开进行讨论。当出现A情况时,角加速度a为逆时针,但是其产生的向左的加速度分量不足以改变ah正负,但是这个速度分量会增大行走时的惯性力从而增大发生打滑的可能性。B情况角加速度a为顺时针,产生正方向的加速度可以减小惯性力。C情况时,ah的方向发生改变为负方向,相较B情况更大程度上的增加了行走时的惯性力。并且这个加速度的大小可以改变H点的加速度反向与图4中的C情况一致,是一种不稳定的情况。
由图4和图5可以分析出当与地面接触脚发生打滑时,为了保持身体平衡要使H点速度有相对向前运动的趋势,但是由于人体的膝关节的生理特征限制H-RK杆不能相对RK-RA杆向右弯曲。因此可以认为MRA为顺时针方向即踝关节屈曲运动,H点相对于RK顺时针转动即髋关节伸展运动,有助于身体保持平衡。
Figure 5. Acceleration analysis of hip joint
图5. 髋关节加速度分析图
3. 试验
3.1. 试验对象
本试验选取10名健康年轻男性作为试验对象,年龄(24.5 ± 0.5)岁,体重(72.78 ± 8.6 kg),身高(175.3 ± 5.4 cm)。受试者均为在校学生,身体健康,没有显著的肌肉骨骼病史及神经疾病、心血管疾病。在试验前无剧烈运动、肌肉拉伤和肌肉疲劳等情况,保证每名受试者自愿参加试验并且能够按照试验要求进行试验。
3.2. 试验方法
本试验拟采集人体在行走过程中发生滑跌时关节角度及关节力矩参数。试验平台由VICON三维光学运动捕捉系统(VICON Motion Systems, Oxford, UK)、AMTI三维测力台(AMTI, Watertown, US)、步态行走试验台以及随行保护装置四部分组成。步态行走试验台长8 m,宽3 m,将两块AMTI测力台嵌入平台内并进行固定,在测力台上放置光滑的大理石并保证其高度与试验平台一致。试验中要求受试者穿戴随行保护装置以便保证受试者在坠落冲击时的安全受到保护。
本实验采用对比试验法,在水平步道上涂油创造湿滑条件,将人体滑移后成功恢复平衡与滑移后终致跌倒两试验组相互对照。对比分析成功恢复平衡与跌倒之间人体下肢关节角度和关节力矩的变化情况及差异性,从而探讨人体下肢关节运动的力学效应与平衡恢复的联系。在进行数据采集前要给予受试者适当的准备时间,受试者可以在试验平台上进行试验过程演练以适应身上装配的试验设备,并调整自己的步伐使之双脚在步行过程中尽量分别踩踏在两块测力台的中心位置。受试者要多次行走在普通步道和涂油的湿滑步道,为了不让受试者准确辨认地面的变化,在实验室中使光线昏暗并且指示受试者在行走中目光向前看。这样的试验条件是合理的,因为大量的跌倒发生在光线昏暗的环境中。每位受试者要在无外界扰动与有外界扰动的两种情况下完成,两种试验随机进行以尽量减小受试者对滑跌事件的心理预防。目视前方,并且未告知其地面情况,使其在无意识的情况下完成试验,使得试验更接近实际情况。试验中由工作人员操作VICON系统软件采集两种试验环境下受试者的步态参数、运动学参数以及动力学参数。
3.3. 试验参数
本文是基于步态分析,对人体行走步态滑跌现象进行研究,对比分析滑移恢复平衡与跌倒之间人体下肢关节角度和关节力矩参数的变化情况及差异性,从而探究在步态滑跌过程中人体下肢关节运动的力学效应。
本文运用VICON采集油面水平步道行走环境下,髋、膝、踝关节各关节力矩和关节角度的变化规律。“关节力矩”定义为一个力和关节中心点到此力作用线垂直距离的矢量积。在步态研究中,关节力矩均指关节内力矩。举例说明:伸膝力矩是指当发生在胫骨、股骨间反作用力的作用线向后经过关节屈伸轴时,即当外部力矩驱使膝关节发生屈曲时,膝关节内部所产生的伸膝力矩。Liu等[10]调查了平衡恢复过程中年龄对下肢关节力矩的影响,结果表明在反应发生至平衡恢复过程中,老年人依靠额状面和矢状面的关节力矩较多,而年轻人则主要依靠矢状面的关节力矩。因此针对本试验中受试者的具体情况,主要考虑受试者矢状面内的关节力矩与关节角度作为研究指标。
4. 试验结果与分析
本文是基于步态分析,对人体行走滑跌现象进行研究,对比分析恢复平衡与跌倒之间人体下肢关节角度和关节力矩参数的变化情况及差异性,探究在步态滑跌过程中人体下肢关节运动的力学效应对恢复平衡有何作用与影响。所有受试者全部完成实验,对实验所得下肢髋关节角度(angle of hip joint)、膝关节角度(angle of knee joint)、踝关节角度(angle of ankle joint)和髋关节力矩(hip joint torque)、膝关节力矩(knee joint torque)、踝关节力矩(ankle joint torque)进行标准化和求解平均值,所得结果均采用(均数±标准差)的形式表示,所得结果采用Origin进行绘图。
4.1. 关节角度
将人体步态滑移恢复平衡与跌倒两试验组髋关节、膝关节和踝关节矢状面的关节角度数据进行时间标准化,求多组试验的数据进行均值处理,做成变化曲线图。
Figure 6. Comparison of joint angle of perturbed leg
图6. 滑动腿关节角度对比图
如图6所示,发生意外滑移时,在滑动腿足跟接触步道时刻到非滑动腿足跟接触步道时刻(第一双支撑相与单支撑相)期间,自主恢复平衡试验组(I组)和滑倒试验组(II组)进行比较相,可得结果如下:
(1) I组的髋关节屈曲角度大于II组,当这个步态滑跌的过程进行到所选区间大于40%时,I组髋关节的屈曲角度快速下降。可以得知此时为滑动腿单支撑期,为了使身体恢复平衡滑动腿会激活髋关节的伸展肌群使其收缩减小髋关节的屈曲角度,而跌倒试验组的屈曲角度变化较慢,说明髋关节的伸展肌群在这个过程中并未激活或者激活较小终至滑倒。
(2) 膝关节传递负荷,参与运动,辅助动量守恒,并为腿部活动提供力偶。膝关节是人体最大,也是最复杂的关节,它由胫骨关节和髌骨关节组成双关节机构。I组和II组的膝关节角度变化趋势基本一致,但是II组膝关节角度小于I组,对比结果表明滑跌过程中膝关节角度减小。
(3) 踝关节和足部的生物力学相当复杂且彼此关系紧密。足部是整个下肢运动系统的一部分,对步态平稳顺畅有重要作用。踝关节将承重由下肢传递到足部,决定着足部在地面上的定位。I组和II组踝关节角度变化趋势基本一致,在0%~20%时间内I组踝关节角度小于II组,在20%~80%时间内I组踝关节角度大于II组,然后II组大于I组。
4.2. 关节力矩
试验对所有关节力矩幅值进行标准化,标准化的方法为用采集到的试验数据除以受试者的体重和腿长,以便对10组受试者的关节力矩进行统计分析。
Figure 7. Comparison of joint torque of perturbed leg
图7. 滑动腿关节力矩对比图
如图7所示,发生意外滑移时,在滑动腿足跟接触步道时刻到非滑动腿足跟接触步道时刻(第一双支撑相与单支撑相)期间,自主恢复平衡试验组(I组)和滑倒试验组(II组)相比较,可得结果如下:
(1) I组髋关节关节力矩和II组变化趋势基本相同,但是I组需要更大的屈曲力矩。进入滑动腿单支撑时其髋关节的屈曲力矩以更快的速率降低,并且在50%时转变为伸展力矩并最终其值更大。
(2) I组膝关节关节力矩和II组变化趋势大致相同,但是I组膝关节力矩变化比II组稍慢。20%后的时刻膝关节的力矩为屈曲力矩,在40%后的时刻滑动腿单独支撑时,成功恢复平衡的试验组较跌倒试验组拥有更大膝关节屈曲力矩。这说明在步态滑跌过程中,为了保持身体的平衡需要更大的膝关节屈曲力矩,带动膝关节的屈曲运动趋势。
(3) I组踝关节关节力矩和II组踝关节力矩在步态周期的前20%期间变化趋势相同且大小基本相等,在20%之后I组明显大于II组。这说明在步态滑跌过程中踝关节的背屈运动更加有助于平衡的恢复。
5. 数据对比分析
现有的人体行走平衡理论认为,当人体的质心垂直投影COM在双脚的支撑区间BOS内便可维持身体平衡[9]。通过建立人体行走的运动模型,分析得出在发生滑跌时判断身体滑跌程度。之后进行试验检测人体滑跌时下肢关节运动的变化情况,对比成功恢复与跌倒试验之间数据的差异性。将检测分析结果与理论分析结果相比较是否相符。
对髋关节关节角度对比会发现:I组(自主恢复平衡)与II组(滑倒)相比,两者髋关节均呈屈曲状态,但II组屈曲角度更小。滑倒时,质心向下后移,导致髋关节屈曲角度减小,质心垂直投影超出BOS,引发滑倒。
对髋关节关节力矩对比会发现:髋关节屈曲力矩有小幅波动且呈下降趋势,伸展力矩增大,使髋关节从屈曲转为伸展,质心向上前移,维持身体直立。进入单支撑阶段,屈曲力矩快速降低,50%时转为伸展力矩并达最大。说明滑跌时,髋关节屈曲角度和力矩均减小。
对膝关节角度变化对比会发现:膝关节角度表现为屈曲角度,当滑动腿足跟着地时,膝关节处于伸膝状态。在滑移瞬间,人体质心向后下移动,膝关节开始屈曲。
对膝关节力矩变化对比会发现:滑移瞬间,膝关节力矩仍为伸膝力矩。但随着人体通过屈膝运动维持质心在BOS内,膝关节屈曲角度增加。20%时刻后,膝关节转为屈曲力矩,结合髋关节伸展,使人体恢复直立。自主恢复平衡组(I组)需更大屈曲力矩,表明步态滑跌时,保持平衡需更大膝关节屈曲角度和力矩。
对踝关节角度对比分析会发现:自主恢复平衡组(I组)与滑倒组(II组)对比,I组前20%背屈角度小,之后增大,并在17%时背屈角度、屈曲力矩开始上升。滑移时,踝关节由背屈转为跖屈,背屈角度减小,跖屈力矩增大。大约17%时,随着髋关节伸展和膝关节屈曲,人体质心前移,背屈角度增大,跖屈力矩减小,转为背屈力矩,增大足地接触面积和接触力,以减缓或停止滑动。
6. 结论
本文从步态模型分析入手,结合试验结果得出以下结论:
(1) 本文建立的人体下肢六杆刚体步态模型有效。滑跌时,惯性力随质心与足底压力中心连线与铅垂线夹角增大而增大,导致RCOF增大,引发滑跌。
(2) 滑移瞬间,髋关节屈曲角度和力矩减小。质心偏移出BOS,滑动腿髋关节伸展肌群激活,伸展力矩增大,维持身体平衡。
(3) 膝关节在滑移瞬间开始屈曲,屈膝肌群使伸膝向屈曲转变,结合髋关节伸展,维持身体平衡。
(4) 滑移瞬间,踝关节背屈,质心后移激活跖屈肌群,跖屈力矩增大,背屈角度减小。随髋关节伸展和膝关节屈曲,质心前移,踝关节背屈角度增大,跖屈力矩减小,直至背屈力矩增大,增大足地接触面积和力,减小滑动速度,维持身体平衡。
NOTES
*通讯作者。