头顶高速运动对肌肉和关节力的依赖性：肌肉骨骼模型及比较分析

i本文的最新情况见最后医学信息学解锁20（2020）100415运动员在高速头顶任务中的比例和技能水平对肌肉和关节力的依赖性：肌肉骨骼模拟研究Hamidreza Barnameheia，*，Farhad Tabatabai Ghomsheha，b，Afsaneh Safar Cherati c，Majid Pouladianaa伊朗德黑兰伊斯兰阿扎德大学生物医学工程系科学和研究处b伊朗德黑兰社会福利和康复科学大学人类工程学系c伊朗医学科学大学，Rasoul Akram医院，伊朗A R T I C L EI N FO关键词：头顶肌肉骨骼建模仿真肩关节羽毛球A B S T R A C T羽毛球过顶扣球是一项高速过顶动作，由不同的肌肉-肌腱单元产生运动力。直接检查运动过程中所有肌肉和关节的力是不可能的。本研究的目的是：（1）估计OFS过程中的肌肉和关节负荷，（2）比较通用和缩放通用模型之间的结果。20名羽毛球运动员分为业余和专业两组。使用运动分析系统记录标记轨迹。在OpenSim中建立了5节段、10自由度、26个肌腱单元的肩关节肌肉骨骼模型。运动学逆解用于评价OFS过程中关节运动学。静态优化和利用关节反作用力工具分别通过OFS估计关节和肌肉力。结果通过OFS为业余和职业球员的关节和肌肉的力量模式。在大多数肌肉和关节力中，通用和缩放通用肌肉骨骼模型之间检测到显着差异。此外，业余运动员和专业运动员在OFS期间存在显著差异。肌肉和关节力的结果是有用的运动员，教练员，教练员，和运动医学。识别哪些肌肉以及如何发挥有效作用以产生正确方向的力对于教练员和训练员来说至关重要，因为设计有效的训练以提高运动员的表现。另一方面，关节载荷的模式可以帮助评估哪些关节和哪些序列是高风险的。本研究是第一次全面的研究，研究了两个不同技术水平的头顶高速运动。1. 介绍羽毛球是世界上许多运动员参加的奥运会项目[1]。羽毛球是一项球拍运动，由不同技术水平的人执行，作为单打（每方一名球员）和双打（每方两名球员）[2]。它有时被非精英人士作为休闲户外娱乐[3]。它包括下肢和上肢的一些缓慢和快速运动[4]。不可预测性是羽毛球运动的主要特点之一，这些突然的运动可能在很短的时间内发生在每个运动方向上。因此，羽毛球运动员不仅需要优秀的上肢体能和心理能力，还需要优秀的下肢体能和心理能力。其中最受欢迎和关键的技术之一，羽毛球比赛是头顶扣球（OFS）[4]。OFS是头顶上的，是肩部和牛鞭运动的近端到远端序列的样本，包括上肢和下肢高速运动。对网球发球、投掷和踢球运动学的各种评价显示了端点的角速度或线速度[5，6]。远中近中运动是指手臂和前臂绕纵轴旋转发生在运动终末部位的复杂运动。一些先前的文献表明，学习可以减少运动范围和运动持续时间[7OFS中风是一种高风险损伤运动，可能涉及肘关节和肩关节复合体，因为这些在OFS的执行中起关键作用[11]。根据乌梅托大学医院的研究结果，在所有需要急救的运动损伤中，* 通讯作者。电子邮件地址：h. wsu.edu（H.Barnamehei）。https://doi.org/10.1016/j.imu.2020.100415接收日期：2020年4月21日;接收日期：2020年8月3日;接受日期：2020年8月11日在线预订2020年2352-9148/©2020的 自行发表通过Elsevier 公司这是一个开放接入文章下的CCBY-NC-ND许可证（http://creativecommons.org/licenses/by-nc-nd/4.0/）。可在ScienceDirect上获得目录列表医学信息学期刊主页：http://www.elsevier.com/locate/imuH. Barnamehei等人医学信息学解锁20（2020）1004152是打羽毛球引起的。在90.7%的情况下，受伤的受试者是娱乐运动员或初学者[10]。51.3%的损伤为轻度（AIS 1），48.7%的损伤为中度（AIS2）。32%的病例累及上肢[12]。最常见的肩关节损伤是肩关节不稳定、肩袖损伤和肌肉失衡。肩关节复合体由盂肱关节（肩关节）、胸锁关节、肩锁关节和肩胛胸关节组成。盂肱关节是一种三维球窝关节，在人类生活中便于上肢三维运动，用于执行体育活动或生理活动[13]。肩关节复合体包括一些肌腱、韧带和肌肉，在那里它是一个冗余系统，并包括更多的26个肌肉-肌腱单位（MTU）[14]。这些软组织在维持肩关节复合体的稳定和平衡方面发挥着重要作用。盂肱关节周围产生的肌力和力矩的平衡对盂肱关节稳定性有影响[15]。为了限制盂肱关节的平移运动，肩袖和其他肩部肌肉的协调激活产生了从肱骨头到关节盂的压缩力[16]。Rathi等人和Wattanaprakornkul等人表明，为了肩关节的伸展-屈曲，肩关节的前部和后部肌肉同时产生的力量，以压缩头部的幽默在垂直运动平面上的关节盂窝[17，18]。Lucas在1973年以及Ackland和Pandy在2009年指出，在肩关节外展过程中，三角肌产生向上的激动力，下旋袖肌产生向下的力，将肱骨头压缩到关节盂窝中[19，20]。先前的研究揭示了肌肉在盂肱关节稳定性中的作用。Yousif和Bicos提到，肌肉活动不平衡的影响是疼痛、部分脱位、过度剪切载荷和早发性骨关节炎[21]。一些先前的研究评估了运动或生理活动中肩部的肌肉和关节负荷[22，23]。肌肉和关节负荷的一些方面，如运动学变量之间的关系，肌-拮抗剂关系，每个肌肉产生关节力矩和关节负荷的贡献是模糊的。了解上肢肌肉激活、肌力和关节反作用负荷之间的关系，对于预防运动损伤、认识运动损伤机理、设计康复处方、提高运动成绩和技术水平具有重要意义。尽管如此，非侵入性地量化肩部复合体的神经肌肉变量是困难的。此外，如前所述，肩关节复合体包括引起复杂运动的不同关节[24]。此外，由于肩部的肌肉数量众多，并且其中一些肌肉是多关节的，因此不可能评估每个肌肉在三维运动中的作用[25]。肌肉骨骼建模是一种估计肌肉的方法[26，27]和关节载荷[28该方法不仅可以用来模拟人体运动，也可以用来模拟动物运动。此外，它经常用于评估神经肌肉手术对人体或动物运动的影响。肌肉骨骼模型使研究人员能够估计特定的肌肉力量和每块肌肉产生力矩的作用，然后是特定的运动。肌肉骨骼建模有助于我们理解肌肉骨骼行为，并用于运动应用[32，33]、康复应用[34]、手术效果[32，35]、肩关节植入物效果[36]和人体工程学应用[37]。常见的肌肉骨骼模型通常包括骨骼模型（由几何形状、骨骼文件和属性组成）、关节模型（包括分段、关节、广义坐标、运动学函数、材料和包裹对象）和肌肉（Hill型）模型（由最大等长肌力（FM0）和最大纤维缩短速度最佳肌纤维长度（LM0）、张角（α）和肌腱松弛长度（LTS）。这些特性不仅随受试者的体型、性别、年龄而变化，还随神经肌肉状况而变化[38，39]。肌肉骨骼模型可以用作通用模型或受试者特异性模型。由通用数据先前研究（尸体数量）信息的测量值[40，41]; A.A.[29，42在大多数情况下，模型测量是困难的，在某些情况下是不可能的，由于大量的主题，而与泛型模型的工作更容易和更快。尽管在某些情况下，观察到受试者之间存在很大差异[46，47]。以前的一些文献研究了通用模型和特定对象模型之间的差异[26，47受试者特定的准确性非常重要，可能会影响模拟结果[52受试者特异性肌肉-肌腱缩放通过缩放特定受试者的肌肉-肌腱测量来进行，以达到准确的结果。缩放通常分三部分进行，包括几何、惯性和肌肉[56]。 精确的缩放方法提供精确的几何参数，以确认对反向运动学结果的非负面影响。先前的研究提出了不同的方法，以建立受试者特异性模型，例如标记物基础缩放和成像信息缩放[57，58]。目前，医学图像缩放是从医学图像提供准确几何数据的黄金标准[56，59在该缩放方法中，通常来自于以下的三维几何信息： ct扫描 [30，63 MRI [50，70 EOS [74]或超声波[75]手动[62]或半自动[61]重建。Bartels等人和Lenaerts、Bartels等人分别指出了受试者特定模型对关节力矩、肌肉力和接触力的影响[30，65]。CT图像通常用于提取骨几何结构[63]，而MRI图像通常用于提取肌肉路径[47，52，76，77]。准确的肌肉路径数据直接影响肌肉和关节力[52]，因为与通用模型相比，准确的受试者特定信息会影响肌肉力矩臂[47]。其他方法之一是与光电数据缩放[32，78，79]。这种方法是常见的，因为不需要高层次的知识，时间和金钱，用于肌肉肌腱缩放是改变肌腱松弛长度和最佳纤维长度。由于它们之间的比率是恒定的，因此可以估计给定关节角位置处的总肌肉-肌腱单位（MTU）长度[60，80]。肌肉骨骼模型的主要应用之一是估计-使用数值优化匹配肌肉和关节力。有各种数值优化方法，例如最小化肌肉激活[80-82 ]、能量消耗[ 83 ]、疲劳[ 84 ]和肌肉应力[ 85-87 ]的Brand等人和Pedersen等人指出了成本函数公式对关节载荷的影响[88，89]。通过最小化肌肉活动平方和的静态优化是估计个体肌肉力的常用方法。羽毛球过顶扣球是一项高速过顶动作，是由不同的肌肉-肌腱单元产生力量。肌肉和关节负荷模式不仅用于提高技能和表现，而且还用于通过检查哪些关节和哪些序列是高风险的来防止受伤。另一方面，它是不可能检查所有的肌肉和关节的力量，在这个运动直接而肌肉骨骼建模使可能检查所有的肌肉。此外，不同技术水平的羽毛球运动员使用特定的策略来执行OFS。因此，本研究的目的是：（1）估计OFS期间的肌肉和关节力(2)比较通用模型和比例通用模型的结果。识别哪些肌肉以及如何发挥有效作用以产生正确方向的力对于教练员和训练员来说至关重要，因为设计有效的训练以提高运动员的表现。另一方面，关节载荷的模式可能有助于评估哪些关节和哪些序列是高风险的。本研究的主要创新之处是通过肌肉骨骼模拟比较两种不同技能水平之间的肌肉和关节负荷，而以前的工作仅通过有限肌肉的EMG数据比较不同技能水平。H. Barnamehei等人医学信息学解锁20（2020）10041532. 材料和方法2.1. 科目20名男子职业和业余羽毛球运动员参加了这项研究。研究对象为羽毛球专业运动员，其羽毛球运动经验均在十年以上。以业余羽毛球运动员为研究对象，以业余羽毛球运动员的身份参与本实验。受试者在实验前五年处于健康状态，没有任何上肢或下肢受伤或疼痛。此外，在10年的实验前，他们没有做过任何外科手术.所有业余和专业运动员都阅读并签署了谢里夫理工大学道德委员会的道德表格和知情同意书。业余志愿者和专业志愿者的信息见表1。2.2. 文书在当前实验中使用十台Vicon高速运动分析相机（Vicon MX，OXford Mexico，United Kingdom，采样频率200 Hz）和两个测力板（Kistler，Winterthur，Switzerland，采样频率2000 Hz），以分别记录运动学变量和地面反作用力表面标记安装在Wu等人[90]使用的骨平台标记处。那辆车2.4. 肌肉骨骼建模选择Wu肩关节模型作为基础肌肉骨骼模型[92]。该模型包括五个段和十个自由度（DOF）。使用质量和人体测量数据制作受试者特定模型，并通过每个受试者的表面标记信号完成缩放程序[80]。手腕和前臂段定义为一个刚体段。肘关节定义为连接前臂和臂（肱骨）节段的2自由度万向关节。胸锁关节定义为2自由度万向关节，而肩锁关节和盂肱关节定义为提供旋转运动的3自由度（球窝）关节。该肌肉模型包括由26个肌肉-肌腱单位（muscle- tendon units，MTU）驱动的轴肩胛大肌群、轴肱肌群和肩胛肱肌群。通常，我们使用的肌肉骨骼模型包括肩胛提肌、上斜方肌、中斜方肌（上）、中斜方肌（下）、下斜方肌、菱形小肌、菱形大肌（上）、菱形大肌（下）、上前锯肌、中前锯肌、下前锯肌、锁骨下肌、胸小肌、前三角肌、中三角肌、后三角肌、棘上肌、冈下肌、肩胛下肌、小圆肌、大圆肌、上胸大肌、中胸大肌、下胸大肌、背阔肌和喙肱肌。2.5. 缩放喷射物 的 SIXTY 表面 反射 标记 是 过滤 经由四阶，低通，巴特沃兹，截止频率10 Hz [90]。由于肩关节的复杂性和深部肌肉的存在，采用瑞士产的无线Myon EMG系统，采样频率为2000 Hz，对肱二头肌、棘下肌、冈上肌、中三角肌和前锯肌5块肌肉的肌电活动进行在放置电极前，先刮除皮肤并用酒精清洗，以减少实验误差。EMG信号通过二阶Butterworth高通滤波，截止频率为25Hz ， 并 进 行 整 流 [90] 。 通 过 Vicon Nexus 软 件 （ Vicon ， OxfordMetric，UK）同步记录肌电图、ki-肌电（运动捕获）和动力学（测力板）数据2.3. 实验程序热身15 min后，根据Wu肩模型将表面标记物贴附于骨骼标志物上，然后在选择的肌肉上安装EMG电极。我们遵循了之前关于EMG电极放置的研究[91]。在安装电极和标记后，每个玩家分别对每块肌肉进行MVC测试。最后，每个运动员执行5次OFS，每次试验之间休息2分钟。计算每名受试者五次试验的平均值以获得真实数据。定义了羽毛球后备后卫在当前研究中，为了开发按比例缩放的通用模型，通过静态中性姿势下的表面骨标志（表面标记）位置对肌肉-肌腱单位长度、力臂、节段质量、肌张力和惯性特性进行线性缩放[80，93]。OpenSim缩放工具通过比较虚拟标记之间的距离和实验标记之间的距离来缩放模型中每个片段的尺寸。通过比较模型上真实和虚拟标记中实验标记之间的距离来计算比例因子。此外，模型中依赖于距离或长度的属性和特征，如韧带和肌肉致动器，也是个性化的。2.6. 逆动力学利用逆运动学方法计算了OFS过程中关节的角位置。IK工具通过优化肌肉骨骼模型上的虚拟标记与对象上的真实表面标记之间的差异来计算关节角度。虚拟标记物安装在模型的不同节段上，而真实标记物附着在骨标志皮肤上。IK计算加权最小二乘方程，以最小化每个时间间隔中模型标记数据（Xi）和表面反射标记数据（Xi）之间的长度[94]：作为起始位置，当手下降并垂直于躯干时，通过MOKKA软件中的目视检查定义结束相位。 图 1显示了OFS的不同阶段，最小q（∑N我wi xi-xi（q）2）羽毛球运动员在运动分析实验室中进行的当前研究（见图1）。 2）的情况。表1每组受试者人数、羽毛球经验和人体测量学信息。专业（精英）业余（非精英）受试者人数年龄（岁）24± 2.5 26± 3.2身高（cm）173± 6.1 175± 4.4重量（kg）62± 5.5 68± 7.2经验（年）13±2 2±1所有值均为平均值±SD。其中q表示模型的广义坐标的向量，wi表示第i个标记的权重逆动力学（ID）工具用于计算广义关节力，包括关节力和力矩[95]。ID计算动力学运动方程：M（q）q-+C（q，qstec）+G（q）=τ其中q、qstec和q？是广义坐标的位置、速度和加速度，例如线性或角位置。M、C、G分别是质量（或惯性矩）、离心力、科里奥利力和重力矩阵，而τ是未知广义力（力或力矩，[23]）的矢量H. Barnamehei等人医学信息学解锁20（2020）1004154∑Fig. 1. 利用Mokka软件和实验表面标记数据，在运动分析实验室中对羽毛球运动员完成的头顶扣球的各个阶段进行了可视化。图二. 前视图（左）、侧视图在当前研究中使用的肌肉骨骼模型的后视图侧（中视图）和后视图侧包括5个节段、10个自由度、26个Hill型肌肉。（包括肩胛提 肌 、 上 斜 方 肌 、 中 斜 方 肌 （ 上 ） 、 中 斜 方 肌（下）、下斜方肌、小菱形肌、大菱形肌（上）、大菱形肌（下），上前锯肌、中前锯肌、下前 锯肌、锁骨下肌、胸小肌、前三角肌、中三角肌、后三角肌、棘上肌、冈下肌、肩胛下肌、小圆肌、大圆肌、上胸大肌、中胸大肌、下胸大肌、背阔肌和喙肱肌）和表面标记。2.7. 静态优化静态优化（SO）用于估计OFS运动期间的肌肉激活和力[96]。SO工具通过最小化平方肌肉激活的总和来估计所有26个肌肉激活和力。SO工具使用Hill型肌肉模型和力-长度-速度属性。肌肉激活的结果，将SO方法与实验肌肉活动进行比较，以验证静态优化结果[97][42]。nJ=（am）pm=1其中，a表示肌肉m的激活水平，n表示肌肉骨骼模型中的致动器数量，并且p表示幂H. Barnamehei等人医学信息学解锁20（2020）1004155++由运算符定义的等式。在本研究中，功效值定义为2。2.8. 关节反利用OpenSim分析工具中的JointReaction工具估计OFS期间上肢的关节反应负荷关节反作用载荷包括每五个节段估计的力和力矩（在父节段参考系中施加力将运行联合反作用力分析，因此联合反作用力估计如下：Rj=mjaj-（∑Fext+∑Fcons+Rj+1）其中，Rj表示关节“j”处的反作用力通过比较模拟的肌肉活动结果和EMG的实验肌肉活动，验证了肌肉骨骼模拟的结果[98，99]。2.9. 数据分析同步运动分析（标记轨迹）、测力板和EMG数据，并从NEXUS软件（Vicon，OX ford Metric，UK）中导出为C3D文件。使用Mokka软件从C3D文件中区分开集、接触点和偏移。在Mokka可视化之后，裁剪的数据被合并到MATLAB（版本2019b）软件中。将OFS的持续时间标准化为百分比，以便于组间的时间比较。 此外，将EMG数据标准化为最大自愿等长收缩（MVIC），以促进组和肌肉之间的比较。将肌肉和关节反作用力标准化为体重百分比。此外，为了验证模拟结果[1022.10. 统计平均值和标准差用于呈现静态优化和联合反应分析的结果。Shapiro-Wilk检验用于检查正态性和方差等同性。双向方差分析用于区分业余组和专业组之间的显着差异，以及通用和标度通用模型的显着差异。显著性水平设定为p值低于0.05。所有统计计算均采用MATLAB 2019b软件进行。3. 结果图3显示了羽毛球专业组和业余组在实施头顶扣球时，肩关节和肘关节在屈伸、外展、内收和旋转轴上的角位置和角速度的关系。右轴表示角速度，左轴表示角位置。图中的画笔部分显示了接触时间的平均范围，这在专业组和业余组中是不同的。外展-内收运动仅定义为肩关节，因为我们的模型限制了肘关节的外展-内收。 显著 差异 是 发现 之间业余运动员与职业运动员肘关节和肩关节各方向角速度差异有统计学意义（P0.05）。<专业组的OFS角位置和角速度在不同时间有一定的变化，而业余组则服从稳态运动模式。的图三. 羽毛球专业组与业余组扣球时肩肘关节屈伸、外展内收、旋转轴角度位置与角速度关系之研究。右轴表示角速度，左轴表示角位置。图中的画笔部分显示了接触时间的平均范围，这在专业组和业余组中是不同H. Barnamehei等人医学信息学解锁20（2020）1004156=====-==-=-=-==-=业余组和专业组在接触时间之前观察到外展的峰值角位置。两组手臂的旋转角速度峰值均出现在接触时。肩关节最大外展角速度在业余组和职业组之间存在差异，职业组肩关节最大外展角速度出现在接触后，业余组肩关节最大外展角速度出现在接触前图4显示了在由专业和业余羽毛球运动员执行的OFS期间，来自通用和缩放通用肌肉骨骼模型的平均选择性肌肉力。这些子图形分别代表三角肌前肌、三角肌中肌、三角肌后肌、冈上肌、冈下肌、肩胛下肌、小圆肌、胸大肌上肌、胸大肌下肌、背阔肌和喙肱肌。统计学意义 差异之间 通用 和 标度类属 结果 为所有业余组和职业组的肌肉均观察到（p<0.05）。此外，观察到业余组和专业组之间存在显著差异。表2列出了这些数据之间的比较统计报告（专业通用、专业通用分级、业余通用分级和业余通用分级）。职业组和业余组在波幅和时间上存在差异。这些差异在各种肌肉力量上是不同的。例如，肩胛下肌、小圆肌和大圆肌的通用和缩放通用结果之间的差异在业余和专业组中都很小。业余组和专业组之间的大部分肌肉，除了冈上肌，小圆肌，喙肱肌和冈上肌的差异。业余组和专业组之间在胸大肌上下、大圆肌、三角肌前后两表2列出了平均值、标准差和峰值肌肉在OFS期间，业余和专业团体的力量，规模通用模型。最后三列显示两组之间差异的p值。业余组的肌力结果显示，除大圆肌外，通用模型和比例通用模型之间存在显著差异（p0.229）。此外，除菱形大（下）肌外，专业组的肌力结果在通用和比例通用肌肉骨骼模型之间存在显著差异（p0.701）、前中锯肌（p0.287）、胸小肌（p 0.927）、三角肌中部（p 0.239）、冈上肌（p0.114）、冈下肌（p 0.103）、肩胛下肌（p 0.512）、小圆肌（p0.085）、胸大肌上部（p 0.134）、胸大肌下部（p 0.203）、背阔肌（p0.241）和喙肱肌（p 0.095）。此外，业余和专业团体之间的显着差异，观察到所有肌肉（p 0.05），三角肌前部除外（p 0.558），<肩胛肌（p0.772）、小圆肌（p0.965）和喙肱肌（p0.097）。图5显示了业余和专业运动员在OFS期间，与肩关节伸展-屈曲、外展-内收和内外旋角度相关的平均上、压缩、前和合成关节反作用力，来自通用和缩放通用肌肉骨骼模型。仅在专业组中观察到通用模型和比例通用模型之间的显著差异。业余和专业组之间的显着差异，所有联合反应部队进行了观察。最大的联合反作用力的变化被发现在外展-内收变量。表3列出了业余和专业组以及每组通用和缩放通用模型的上、压缩、前和合成表示的平均值、标准差和峰值关节力。对于业余选手来说，通用模型和比例通用模型估计的联合反应力是相似的见图4。在职业和业余羽毛球运动员执行OFS期间，来自通用和缩放通用肌肉骨骼模型的平均选择性标准化肌肉力量与体重的关系。子图表示前三角肌、中三角肌、后三角肌、冈上肌、冈下肌、肩胛下肌、小圆肌、大圆肌、上胸大肌、下胸大肌、背阔肌和喙肱肌。表2平均值，标准差，和峰值肌肉标准化力量体重在业余和专业团体在OFS的通用和缩放通用模型。最后三列显示两组之间差异的p值。专业业余专业业余分级通用分级通用分级通用分级通用平均STD峰值平均STD峰值平均STD峰值平均STD峰值一般VS标度-一般一般VS标度-一般业余VS专业肩胛提肌0.028 0.004 0.038 0.026 0.006 0.042 0.020 0.005 0.037 0.012 0.004 0.018 p 0.05 p 0.05 p 0.05<<<上斜方肌0.059 0.008 0.073 0.051 0.013 0.067 0.043 0.012 0.087 0.025 0.010 0.045 p 0.05 p 0.05 p 0.05<<<中斜方肌（上）0.033 0.007 0.038 0.031 0.005 0.041 0.025 0.006 0.051 0.015 0.007 0.028 p 0.05 p 0.05 p 0.05<<<中斜方肌（下）0.032 0.004 0.045 0.031 0.005 0.038 0.024 0.006 0.049 0.014 0.006 0.025 p 0.05 p 0.05 p 0.05<<<下斜方肌0.195 0.043 0.342 0.177 0.058 0.265 0.144 0.050 0.298 0.084 0.043 0.165 p 0.05 p 0.05 p 0.05<<<小菱形0.044 0.005 0.049 0.042 0.009 0.061 0.032 0.010 0.061 0.020 0.010 0.040 p 0.05 p 0.05 p 0.05<<<大菱形（上）0.021 0.004 0.026 0.022 0.006 0.041 0.017 0.005 0.035 0.010 0.004 0.016 p 0.05 p 0.05 p 0.05<<<大菱形（下）0.012 0.002 0.015 0.012 0.004 0.023 0.009 0.003 0.017 0.005 0.002 0.009 p= 0.701 p 0.05 p 0.05<<前上锯肌0.061 0.012 0.096 0.068 0.014 0.093 0.058 0.023 0.115 0.031 0.015 0.056 p 0.05 p 0.05<<<0.0100.0350.0080.0060.0200.1800.0670.0450.0490.1050.1170.0350.1000.1050.0610.1600.025p 0.05 p 0.05<