基于跑姿控制训练模式下的跟腱生物力学研究

张希妮，王俊清，杨 洋，傅维杰

（上海体育学院 运动健身科技省部共建教育部重点实验室，上海200438）

从人类进化角度出发，现代人所具备的诸多形态（如骨盆狭窄、足弓出现等）是从 200万年前为了适应双足奔跑进化而来的（Davis et al.,2017）。在这些衍生的功能结构中，跟腱是连接足跟和足跖屈肌的关键结构，其主要承担运动时小腿三头肌肌力的传递，是下肢在跑、跳时有效完成力/能量储存和释放的关键（Bramble et al.,2004）。另一方面，随着跑步特别是马拉松运动的发展，有数据显示，我国参与马拉松的人次已增至 498万（林翠娟 等,2018），但相对应的却是高达 83.4%的下肢损伤率（肖梅, 2018）。而这其中，跟腱损伤是最常见且难以根治的跑步损伤，其发生率可达 8%～15%，且跑步者个体的终身患病风险高达52%（De Vos et al.,2010）。虽然现阶段学术界普遍认为，着地过程中跟腱的非同源性载荷（non-homologous loading）以及长时间非正常受力后，跟腱组织及结构的无法完全恢复是导致跟腱损伤的首要原因（Khan et al.,2002），但其真正的发病机制仍不清楚。

目前，针对跟腱损伤的保守治疗方法主要有小腿三头肌离心训练、平衡训练等（Alfredson et al.,1998; Fletcher et al.,2010）。但这些训练本身存在耗时长、可能会引起延迟性肌肉酸痛、以及25%～45%的患者最终仍需要借助手术治疗的不足（De Vos et al.,2010）。且单一的负重力量训练考虑更多的是人体内力问题，而跑步着地受力的整个过程则是地面反作用力等外力和下肢肌肉力等内力的交互综合效果，因此，仅依靠上述训练提高小腿三头肌肌力并希望使跟腱适应跑步时快速且重复的离/向心负荷，可能存在缺陷。

有研究发现长期采用前掌着地跑者相比于后跟着地跑者，其跟腱增加的横截面积和刚度都更利于提高长距离跑的跑步经济性，降低损伤风险（Histen et al.,2017）。McCarthy等（2014）的研究发现，12周跑姿转化训练，使后跟着地跑者更倾向采用前掌着地跑姿，同时表现出与裸足跑相似的运动学模式。这一结果得到了 Khowailed等（2015）的证实，其研究结果发现，6周的跑姿转化训练有可能改变习惯后跟着地跑者的运动模式。但这种跑姿转化训练对跟腱产生怎样的力学影响亟需探究。

综上所述，由于跑步着地时跟腱的内、外力学因素及其所受载荷的生物力学机理仍不清晰，我们应从跟腱力学本身出发，通过研究整个跑步着地中跟腱的载荷，包括肌力、跟腱力、应力、应变等，更深层次地理解跟腱受力在跑步着地时和跟腱损伤中所扮演的角色。因此，本研究旨在通过建立12周跑姿运动干预模型，并配合小腿三头肌离心训练，探究下肢姿态控制对跟腱生物力学的影响，以期增强跟腱力学性质，提高跟腱承受载荷的能力，从而达到预防和减小运动损伤风险的目的。

1 研究对象与方法

1.1 研究对象

选取30名男性跑步爱好者作为研究对象，根据随机表分为前掌跑训练组（SBR组，年龄：31.4±6.5岁、身高：175.7±4.9 cm、体重：70.6±8.2 kg）和极简跑鞋组（MIN组，年龄：28.4±5.8岁、身高：173.7±5.6 cm、体重：72.4±9.9 kg）各 15名。其中，剔除由于在训练干预过程中因伤（损伤原因与训练无关）、因病等不可抗拒的个人原因退出或中断1周以上人员数据，最后共17人符合标准，纳入结果统计（表 1）。所有受试者均习惯穿缓冲功能跑鞋进行后跟跑。近 4周每周跑步距离大于 20 km，且从未尝试过极简跑鞋，近 3个月内无明显下肢损伤。实验前告知受试者实验流程并签署知情同意书。

表1 本研究受试者基本情况Table 1 Basic Information of Subjects (M±SD)

1.2 实验仪器

1.2.1 实验用鞋

NIKE AIR ZOOM PEGASUS 34慢跑鞋（泡棉+气垫中底，前后落差12 mm，平均质量285 g）用于跑步实验（图1A）。INOV-8 BARE-XF 210极简跑鞋（3 mm橡胶外底，无中底，前后落差0 mm，平均质量227 g）用于干预训练（图1B）。

图1 本研究实验用鞋Figure 1.Shod Conditions in the Study

1.2.2 Vicon运动捕捉系统

16台英国 Vicon三维红外摄像头及运动捕捉系统（型号：T40）配套直径 14 mm的红外反光球，采集下肢髋、膝、踝3关节矢状面运动学参数，采样频率为100 Hz。

1.2.3 Kistler三维测力台

瑞士生产的Kistler三维测力台2块（90×60×10 cm，型号：9287C），外置信号放大器，采样频率为1 000 Hz。

1.2.4 Delsys肌电-加速度一体化测试仪

使用美国Delsys公司生产的64导EMG-Amplifier系统对下肢踝关节主要肌群（胫骨前肌、腓肠肌内/外侧头、比目鱼肌）的肌电信号进行测量。采样频率设为2 000 Hz。

1.2.5 光栅计时系统

采用意大利Microgate公司生产Witty-Manual光栅计时系统，测量受试者通过跑道的时间，从而控制平均速度。

1.2.6 Mindray M7 Super超声影像

Mindray公司生产的 M7 Super超声影像仪，本实验中采用 L6-14线阵探头，主要用于测量跟腱横截面面积以及踝关节等长跖屈时跟腱的在体运动影像，包括跟腱长度变化等。

1.2.7 Con-trex等速肌力测试仪

德国Physiomed公司生产的Con-Trex人体肌力评估和训练系统（型号：Con-Trex-Mj）。本实验主要针对踝关节的等长跖屈运动进行肌力测试。

1.2.8 Podoon压感智能鞋垫

使用EVA材料制造，内嵌0.2 mm柔性薄膜压力传感器，通过手机 APP内置算法对本实验训练干预中的前掌触地占比进行监控。

1.3 具体流程

1.3.1 超声影像学

利用超声影像仪获取跟腱横截面面积的成像，要求受试者俯卧于在治疗台上，使踝关节处于中立位（90°），将ML6-15-D线阵探头垂直放置在跟骨上 10 cm处以获取跟腱的横断面影像（Kernozek et al.,2018）。

受试者坐于Con-Trex等速肌力测试仪上，臀部和后背紧贴座位，膝关节完全伸展，踝关节处于中立位，进行踝关节跖屈/背屈的最大等长收缩（MVIC）（图 2A），并保持5 s，重复3次，同步踝关节跖屈力矩以及在体跟腱长度变化图像和数据（Joseph et al.,2012）（图2B、C）。

图2 跟腱超声影像学影像Figure 2.The Ultrasound Imaging of the Achilles Tendon

1.3.2 生物力学

跑步实验开始前，要求每名受试者先更换实验背心、短裤、袜子和慢跑鞋，而后在跑台上以 12 km/h的速度进行 5 min热身。热身结束后，操作人员为其贴上反光球和肌电电极，采集受试者静态模型。反光球的位置为髂前上棘、髂脊上缘等共36个点（Almonroeder et al.,2013）。实验所测试的肌群为胫骨前肌（TA）、腓肠肌外侧（LG）、腓肠肌内侧（MG）和比目鱼肌（SOL）（Suydam et al.,2015）。安置肌电电极前，先进行备皮工作，剃除汗毛，酒精擦拭之后，进行肌电电极的固定工作（图 3）。而后要求受试者坐于 Con-Trex等速肌力测试仪上，先进行踝关节跖屈/背屈的 MVIC，同步采集 TA、LG、MG、SOL肌电信号，以用于肌电标准化处理。

图3 受试者红外反光球和表面肌电的放置位置Figure 3.The Set-up of Infrared Ball Marker and EMG Position

实验开始时，要求受试者以 12 km/h（±5%）的速度通过跑道（利用光栅计时系统控制速度）并要求受试者优势侧（根据踢球法判定踢球距离更远的一侧为优势侧）（Zhang et al.,2018）成功踏上测力台3次，采集反光球轨迹、表面肌电等信号。

1.4 12周跑姿控制训练方案

1.4.1 MIN组

每位受试者发放一双极简跑鞋。要求受试者穿着极简跑鞋，并以现有的跑姿习惯进行跑步训练，时长同 SBR组。方案外的时间可以穿着原本的跑步鞋进行跑步。

1.4.2 SBR组

每位受试者发放一双极简跑鞋和一双Podoon压感智能鞋垫并通过手机app用以配合检测为期12周的跑姿再训练干预方案（Tam et al.,2015）（表 2）。训练干预时要求受试者穿着内置智能鞋垫极简跑鞋，采用中等强度的自选速度进行跑步，要求其在着地时采用前掌触地方式，即需要受试者用前掌跖球部（the ball of the foot）触地，并允许后跟的随后着地，同时尽量使足落于髋关节下方（Rice et al.,2017）。受试者通过智能鞋垫所提供的数据反馈获知其触地方式是否为前掌。干预方案只是替代实验组部分的训练量，每周总距离保持不变。实验期间禁止受试者参加任何相关赛事。除此之外，指导受试者加强足部和下肢肌群的功能及力量练习，从而使下肢肌-骨和软组织系统适应长期不同的跑姿控制及着地负荷（Warne et al.,2015）。

表2 12周跑姿控制干预方案Table 2 12-week Running Retraining Protocol

所有受试者每完成一次训练需要记录一次训练日志。其中，MIN组需要记录训练开始时间、结束时间、地点、距离以及损伤情况。SBR组除需要记录上述内容外，还需另外记录前掌跑开始时间、结束时间以及距离。所有受试者每次开始跑步前都需要告知实验人员具体地点和开始时间，以方便实验人员进行不定期检查，并通过智能鞋垫的云端数据与记录的训练日志进行比对。

1.5 数据处理和实验参数

本实验研究的落地阶段为触地瞬间至脚尖离地的阶段。采用美国 C-Motion公司开发的 V3D三维步态分析软件（版本：3.21.0）对采集到的运动学数据进行滤波和逆向动力学分析，其中，运动学的截止频率分别为7 Hz（Zhang et al.,2018）。采用美国 NI公司开发的 DASYLab信号采集分析软件（版本：v.8.0）分析表面肌电（sEMG）数据。滤波范围为10～400 Hz （ 燊张 等, 2016）。具体参数如下：

1）跟腱横截面面积（CSA）：采用超声影像仪获取跟腱横截面面积的成像，并使用 Image J软件（NIH，USA）确定跟腱图像中的跟腱横截面面积（Histen et al.,2017）。

2）跟腱力（FAT）≈ 肌肉Fmus= TQ / MA，其中，TQ为跖屈力矩，由肌力测试仪获得；MA为跟腱力臂，由从踝至跟腱前边界的平均垂直距离得出，默认值为 0.05 m（Komi, 1990）。

3）跟腱应力 σ = FAT/ CSA （Joseph et al.,2012）。

4）跟腱应变 ε =ΔL / L0，其中，ΔL为是跟腱近端束的位移，L0为跟腱静息长度（Joseph et al.,2012）。

5）跟腱力峰值（FMAX）= Mankle/ LAT，其中，Mankle为跑步时的踝关节力矩峰值，由逆向动力学计算获得；LAT为跟腱力臂长度，本团队基于 Lyght等（2016）用于计算小腿三头肌肌肉力臂的多项式算法及 Rugg等（Rugg et al.,1990） 在体跟腱影像学数据所得（图4）。

图4 踝关节角度与跟腱力臂长度的关系拟合曲线Figure 4.The Fitted Curve of the Relationship between Ankle Angle and Achilles Tendon Force Arm Length

6）均方根振幅值（RMS）：计算着地冲击前 50 ms肌肉预激活阶段（pre-activation）、着地后50 ms后激活阶段（post-activation）、离心阶段和向心阶段的 RMS，其中，根据膝关节屈曲到最大值的时刻划分离心阶段和向心阶段。与MVIC下的RMS进行比值以获得标准化肌电数据。计算方法如下（ 燊张 等, 2016）：

其中，t为 EMG信号开始时间，t＋T为信号结束时间。

1.6 统计学

所得参数值均用平均值±标准差（M±SD）表示。本研究采用重复测量双因素方差（two-way repeated measures ANOVA，SPSS 21.0）观察自变量（组别×训练前后）对于各因变参数（跟腱力、应力、应变、肌电信号等）的影响，显著性水平α设为0.05。

2 研究结果

2.1 跟腱横截面积、跟腱力、应力、应变

12周跑姿控制训练后，两组的跟腱力、应力和应变均显著增加（P＜0.05），而跟腱横截面积没有显著变化（图5）。

图5 12周跑姿控制训练对跟腱力、横截面积、应力、应变的影响Figure 5.The Effect of 12-week Running Retraining Protocol on Achilles Tendon Force, Cross Section Area, Stress and Strain

2.2 跟腱力峰值

12周跑姿控制训练后，两组在 12 km/h（±5%）速度下跑步时的跟腱力峰值均表现出显著增加（P＜0.05），且 SBR组的跟腱力峰值显著高于 MIN组（P＜0.05）（图6）。

2.3 肌电图学

12周跑姿控制训练后，两组胫骨前肌的 RMS在预激活阶段均表现出显著减小（P＜0.05），SBR组的比目鱼肌RMS显著高于 MIN组（P＜0.05）；两组腓肠肌外侧头的RMS在后激活和离心收缩阶段均表现出显著减小（P＜0.05）（图 7）。

图6 12周跑姿控制训练对跑步时跟腱力峰值的影响Figure 6.The Effect of 12-week Running Retraining Protocol on Achilles Tendon Peak Force during Running

图7 12周跑姿控制训练对触地前预激活阶段、触地后激活阶段、离心阶段和向心阶段RMS的影响Figure 7.The Effect of 12-week Running Retraining Protocol on RMS of Pre-50ms, Post-50 ms,Eccentric and Concentric Phase during Running

3 讨论

本研究从跟腱力学本身出发，通过建立12周跑姿控制训练的运动干预模型，并配合小腿三头肌离心训练，以期达到改变跟腱力学，提高跟腱承受载荷能力的目的，从而为预防和减小跟腱损伤提供新的理论依据。目前研究普遍认为，由于循环应变所引起的力传导机制会影响结缔组织（如跟腱等）的内稳态，以调节其适应过程，从而通过力学刺激的调节影响跟腱对运动的适应性反馈（Wang et al.,2006）。Arampatzis等（2007）的研究发现，与非负重项目运动员相比，负重项目运动员具有更强的跟腱力学性质。类似地，Waugh等（2012）在比较儿童与成人跟腱力学性能差异时发现，成人跟腱刚度增加的主要原因是由于在负重运动中跟腱负荷刺激的增加，而与年龄并无关系。因此，跟腱的训练需要通过一定的力学刺激，促进肌腱胶原蛋白的排列，以提高其力学性质，从而适应应力的需要。本研究通过极简鞋诱导受试者进行前掌跑训练，刺激其每次触地时的跟腱受力，以提高其力学性质。

本研究结果支持前人的研究（Joseph et al.,2017），12周跑姿控制训练干预使相应的跟腱力学性质发生改变，即跟腱力、应力和应变均显著增加。与后跟触地相比，刚触地时，前掌触地模式在踝关节周围会产生更大的跖屈内力矩以对抗背屈外力矩（Lieberman et al.,2010）。此时，小腿三头肌收缩，跟腱被拉长。这种重复性的“拉长-收缩”循环，可以刺激跟腱增强力学性质，以提高其承受载荷的能力。另一方面，前人研究表明，穿着极简鞋跑步时，下肢主要采取软着陆（soft landing）的触地方式，同时可以增加跟腱应力、应变以及足部肌肉力量（Almonroeder et al.,2013;Lyght et al.,2016）。然而，由于本研究受试者均为习惯穿着传统跑鞋的跑者，穿着极简鞋跑步被作为一种训练跟腱的方法，本研究不仅仅希望其在穿着极简鞋跑步时能提高跟腱力学性质，更希望其在回归传统跑鞋跑步后能继续维持较高的跟腱力学性质。本研究发现12周跑姿控制训练后，两组在穿着慢跑鞋跑步时的跟腱力峰值也表现出显著增加，且SBR组的跟腱力峰值显著高于MIN组。这与抗阻训练（Bonacci et al.,2009）和最大自主等长训练（Arampatzis et al.,2010） 的干预结果类似，这类训练大多每周进行2～3次，并在 8～12周后观察到跖屈力矩、跟腱力、弹性模量等显著增加。跟腱传递小腿三头肌收缩产生的力量。通过我们的研究结果，可以假设，前掌跑训练后，力量的增加可能发生在小腿三头肌，且与MIN组相比，SBR组的力量增加更加显著。我们推测，系统的前掌跑训练可以使受试者的跟腱能够在跑步期间进行更有效地储存和释放弹性能。这表明，12周跑姿控制训练是一种可以有效提高跟腱力学性质的训练方式，而更优的跟腱力学特性，很可能是预防跟腱损伤的根本（Henriksen et al.,2009） 。

跟腱力与跟腱横截面积是影响跟腱应力的主要原因（Joseph et al.,2012）。本研究结果发现，跟腱横截面积在干预训练前后变化并不显著，表明干预训练后，跟腱受其载荷影响而产生适应性增生，横截面积有增加的趋势，但不会不断增粗。有研究认为跟腱横截面积对训练后跟腱力学性质的增强并不起主要作用（沈勇伟 等,2010）。另一种解释认为与肌肉相比，跟腱内血管稀疏，供血能力差，其代谢速度较慢，重塑可能需要更长的时间（顾文奇 等,2009）。因此，跟腱应力的变化可能更多地与跟腱力学性质的增强有关，表明跟腱承受载荷的能力提高。另一方面，训练过程中跟腱需要承受较大的外力，经历较大的长度变化，因此具有较高的依从性（compliance）。本研究发现，经过 12周跑姿控制训练干预后，跟腱应变显著增加。跟腱在触地早期通过拉长储存能量，在支撑相后期通过释放能量产生推进力。较高的依从性有利于小腿肌纤维通过改变长度以更有效地完成力/能量储存和释放（Kubo et al.,2015）。

跑步时，下肢肌肉触地前被预先激活，而后迅速进入离心收缩，紧接着迅速进入向心收缩。本研究发现12周跑姿控制训练后，两组胫骨前肌的RMS在预激活阶段均表现出显著减小，表明胫骨前肌在这一阶段被较少激活。这一结果支持前人研究。郑义等（2018）对习惯后跟触地的跑步爱好者进行30 min的裸足跑训练后发现，裸足跑训练后着鞋跑时预激活阶段的胫骨前肌肌电平均振幅显著减小。类似地，Khowailed等人（2015） 研究表明，经过 6周仿裸足跑训练后，在预激活阶段受试者胫骨前肌的激活肌电活动显著降低。人体在穿着极简跑鞋时更倾向于前掌触地，这一跑步姿态使跑步时踝关节处于较跖屈的状态，因此胫骨前肌的激活水平减弱。两组腓肠肌外侧头的RMS在后激活和离心收缩阶段均表现出显著减小。Albracht等（2006）利用希尔肌肉模型模拟肌肉收缩时发现，小腿三头肌具有较大收缩力量和较高刚度的特性，导致其在高激活水平下发力比较困难，而在较低的激活水平下发力更具有优势。另外，Suydam等（2015）的研究也证明，跟腱断裂患者累及侧的腓肠肌内、外侧头激活程度均高于未累及侧，以补偿术后跟腱松弛所引起的肌肉无力问题。本研究显示，作为影响跟腱力的主要肌群之一，虽然腓肠肌外侧头在后激活和离心收缩阶段的激活水平表现出显著降低，但跑步时跟腱力峰值表现出显著增加，表明12周跑姿训练干预后，腓肠肌外侧头在较低的激活水平下表现出了更高的激活效率。

虽然本研究仅要求 SBR组在穿着极简跑鞋时采用前掌触地模式，而并未对MIN组触地模式作出要求，但结果发现，12周训练后两组的跟腱生物力学特性均发生显著性变化。这一结果支持前人研究。McCarthy等（2014）要求习惯后跟触地的跑步爱好者穿着极简跑鞋进行12周训练，但并未对其跑姿作出要求，结果发现训练后，100%的受试者在裸足跑时采用非后跟触地模式，56%的受试者在穿鞋跑时依然采用非后跟触地模式。我们推测，12周重复性的周期运动使 MIN组有足够的时间进行自我运动控制学习（Bonacci et al.,2009）。穿着极简跑鞋时，MIN组跑者可能通过主动改变跑姿减小触地冲击以产生自我保护性反馈（Lieberman et al.,2010）。因此，有理由推测MIN组在穿着极简跑鞋跑步时可能倾向于非后跟触地模式从而导致跟腱力学性质及相关肌群的变化，但由于本研究仅关注 SBR组在训练中的触地方式，因此，在未来的研究中应该同时关注训练中两组的触地方式变化。本研究结果表明，没有受试者因为训练造成受伤而退出的情况，因此，12周跑姿控制训练可以作为一种有效增强跟腱力学性质，提高跟腱承受载荷能力，从而为预防和降低跟腱损伤提供了可能的干预方式。

4 结论

本研究结果表明，12周跑姿控制训练会对跟腱生物力学性质产生显著性影响，且 12周训练后 SBR组跑步时的跟腱力显著高于MIN组，这使跟腱在跑步期间可以更有效地完成能量储存和释放。而增加的跟腱力学性质和降低的肌电均方根振幅表明，相同跑速下肌肉激活程度降低但效率更具有优势，提示，12周跑姿控制训练可以积极改善跟腱承受载荷的能力，从而预防和降低损伤风险。此外，对于习惯后跟触地的跑者而言，进行跑姿转化训练时建议采用循序渐进的方式并尽量配合下肢肌群的强化练习。

顾文奇, 施忠民, 柴益民, 2009.慢性跟腱断裂的治疗进展 [J].中国骨与关节损伤杂志, 24 (10): 953-955.

林翠娟, 杨海, 2018.我国城市马拉松赛事问题的社会学分析 [J].辽宁体育科技, 40 (5): 34-37.

沈勇伟, 张林, 刘劲松, 2010.短跑训练对跟腱横截面积和弹性模量的影响 [J].17 (1): 92-95.

肖梅, 2018.马拉松运动损伤的预防康复措施研究 [J].南方职业教育学刊, 43 (1): 107-113.

张燊 , 傅维杰, 刘宇, 2016.不同着地冲击模式的下肢生物力学研究[J].体育科学, 36 (1): 59-66.

郑义, 王宇, 2018.裸足跑时人体下肢运动学和肌电特征的变化 [J].科学技术与工程, 18 (5): 173-179.

ALBRACHT K, ARAMPATZIS A, 2006.Influence of the mechanical properties of the muscle-tendon unit on force generation in runners with different running economy [J].Biol Cybern, 95 (1): 87-96.

ALFREDSON H, PIETILA T, JONSSON P, et al.,1998.Heavy-load eccentric calf muscle training for the treatment of chronic achilles tendinosis [J].Am J Sports Med, 26 (3): 360-366.

ALMONROEDER T, WILLSON J D, KERNOZEK T W, 2013.The effect of foot strike pattern on achilles tendon load during running [J].Ann Biomed Eng, 41 (8): 1758-1766.

ARAMPATZIS A, KARAMANIDIS K, ALBRACHT K, 2007.Adaptational responses of the human achilles tendon by modulation of the applied cyclic strain magnitude [J].J Exp Biol, 210 (Pt 15):2743-2753.

ARAMPATZIS A, PEPER A, BIERBAUM S, et al.,2010.Plasticity of human achilles tendon mechanical and morphological properties in response to cyclic strain [J].J Biomech, 43 (16): 3073-3079.

BONACCI J, CHAPMAN A, BLANCH P, et al.,2009.Neuromuscular adaptations to training, injury and passive interventions: Implications for running economy [J].Sports Med, 39 (11): 903-921.

BRAMBLE D M, LIEBERMAN D E, 2004.Endurance running and the evolution of homo [J].Nature, 432 (7015): 345-352.

DAVIS I S, RICE H M, WEARING S C, 2017.Why forefoot striking in minimal shoes might positively change the course of running injuries [J].J Sport Health Sci, 6 (2): 154-161.

DE VOS R J, WEIR A, VAN SCHIE H T, et al.,2010.Platelet-rich plasma injection for chronic achilles tendinopathy: A randomized controlled trial [J].JAMA, 303 (2): 144-149.

FLETCHER J R, ESAU S P, MACINTOSH B R, 2010.Changes in tendon stiffness and running economy in highly trained distance runners [J].Eur J Appl Physiol, 110 (5): 1037-1046.

HENRIKSEN M, AABOE J, BLIDDAL H, et al.,2009.Biomechanical characteristics of the eccentric achilles tendon exercise [J].J Biomech, 42 (16): 2702-2707.

HISTEN K, ARNTSEN J, L'HEREUX L, et al.,2017.Achilles tendon properties of minimalist and traditionally shod runners [J].J Sport Rehabil, 26 (2): 159-164.

JOSEPH M F, HISTEN K, ARNTSEN J, et al.,2017.Achilles tendon adaptation during transition to a minimalist running style [J].J Sport Rehabil, 26 (2): 165-170.

JOSEPH M F, LILLIE K R, BERGERON D J, et al.,2012.Measuring achilles tendon mechanical properties: A reliable, noninvasive method [J].J Strength Cond Res, 26 (8): 2017-2020.

KERNOZEK T W, KNAUS A, RADEMAKER T, et al.,2018.The effects of habitual foot strike patterns on achilles tendon loading in female runners [J].Gait Posture, 66 (9): 283-287.

KHAN K M, COOK J L, KANNUS P, et al.,2002.Time to abandon the“tendinitis” myth [J].BMJ, 324 (7338): 626-627.

KHOWAILED I A, PETROFSKY J, LOHMAN E, et al.,2015.Six weeks habituation of simulated barefoot running induces neuromuscular adaptations and changes in foot strike patterns in female runners [J].Med Sci Monit, 21: 2021-2030.

KOMI P V, 1990.Relevance of in vivo force measurements to human biomechanics [J].J Biomech, 23 (1): 23-34.

KUBO K, MIYAZAKI D, TANAKA S, et al.,2015.Relationship between achilles tendon properties and foot strike patterns in longdistance runners [J].J Sport Sci, 33 (7): 665-669.

LIEBERMAN D E, VENKADESAN M, WERBEL W A, et al.,2010.Foot strike patterns and collision forces in habitually barefoot versus shod runners [J].Nature, 463 (7280): 531-535.

LYGHT M, NOCKERTS M, KERNOZEK T W, et al.,2016.Effects of foot strike and step frequency on achilles tendon stress during running [J].J Appl Biomech, 32 (4): 365-372.

MCCARTHY C, FLEMING N, DONNE B, et al.,2014.12 weeks of simulated barefoot running changes foot-strike patterns in female runners [J].Int J Sports Med, 35 (5): 443-450.

RICE H, PATEL M, 2017.Manipulation of foot strike and footwear increases achilles tendon loading during running [J].Am J Sports Med, 45 (10): 2411-2417.

RUGG S G, GREGOR R J, MANDELBAUM B R, et al.,1990.In vivo moment arm calculations at the ankle using magnetic resonance imaging (mri) [J].J Biomech, 23 (5): 495,499-497,501.

SUYDAM S M, BUCHANAN T S, MANAL K, et al.,2015.Compensatory muscle activation caused by tendon lengthening postachilles tendon rupture [J].Knee Surg Sport Traumatol Arthrosc, 23(3): 868-874.

TAM N, TUCKER R, ASTEPHEN WILSON J L, et al.,2015.Effect on oxygen cost of transport from 8-weeks of progressive training with barefoot running [J].Int J Sports Med, 36 (13): 1100-1105.

WANG J H, THAMPATTY B P, 2006.An introductory review of cell mechanobiology [J].Biomech Model Mechanobiol, 5 (1): 1-16.

WARNE J P, KILDUFF S M, GREGAN B C, et al.,2015.A 4-week instructed minimalist running transition and gait-retraining changes plantar pressure and force [J].Scand J Med Sci Sports, 24 (6): 964-973.

WAUGH C M, BLAZEVICH A J, FATH F,.,et al.,2012.Age-related changes in mechanical properties of the achilles tendon [J].Journal of Anatomy, 220 (2): 144-155.

ZHANG X, XIA R, DAI B, et al.,2018.Effects of exercise-induced fatigue on lower extremity joint mechanics, stiffness, and energy absorption during landings [J].J Sports Sci Med, 17 (4): 640-649.