在汽车安全、军事工程和医疗诊断中，都需要对人体受冲击状态进行分析，研究人体组织的冲击力学性能，并发展能描述其动态性能的本构关系，以得到相关的参数，对于安全设计等方面有着重要的意义。肌肉是人体组织中一个重要的组成部分，按性质的不同通常分为3类：骨骼肌、平滑肌和心肌。骨骼肌在人体中占有很大的比重，躯体的各种运动和呼吸均离不开骨骼肌。骨骼肌由许多平行排列的肌纤维及纤维间的结缔组织组成，其力学性能由肌纤维和结缔组织共同决定。与其他工程材料有很大不同的是，它可以产生2种不同形式的力，即主动力和被动力。

研究骨骼肌在外部冲击载荷作用下的被动承载能力，对人体安全设计有重要意义。但由于肌肉组织具有一定的活性，试样获取困难，且极其柔软，给实验特别是动态实验带来了很大的困难，因而相关研究很少等［1］、C.Van Sligtenhorst等［2］曾利用SHPB装置获得了肌肉的动态压缩应力应变曲线，但其动态拉伸性能还无相关报道。为弥补这一空白，王宝珍等［3－4］利用改进的SHPB装置和套筒式SHTB装置对猪后腿肌肉的动态压缩和拉伸性能进行了研究，获得了不同应变率下的应力应变曲线。然而，描述肌肉材料动态力学性能的本构模型还需要进一步讨论。

关于骨骼肌的本构模型，以往的研究往往集中于描述肌纤维的主动收缩性能，如著名的Hill模型。而为了描述骨骼肌的被动力学性能，通常结合其特殊的纤维组织结构，将骨骼肌看作横向同性体 Loocke等［5］曾借用横向同性弹性理论，并引入应变相关杨氏模量，从而描述骨骼肌静态压缩下的非线性横向同性性能等［6］则从另一种思路出发，将肌肉看作纤维增强复合材料，在参考构形中，引入一个纤维方向的单位矢量，由该方向矢量与右Cauchy－Green变形张量可表述出5个应变不变量，他们将应变能密度函数用这5个应变不变量来表示，从而发展了描述心肌被动力学性能的横向同性超弹性本构模型，后被许多学者用于研究骨骼肌［7－8］的横向同性超弹性。然而，由于一直以来缺乏肌肉动态性能研究的实验结果，因而未能得到一个描述肌肉动态性能的率相关本构模型。

本文中，将结合已有的猪后腿肌肉动态压缩和拉伸实验结果，发展一个既能描述肌肉动态压缩性能又能描述其动态拉伸性能的率相关本构模型。

1 纤维增强材料应变能函数表示方法

肌肉、韧带等生物软组织通常可看作是柔性基质中分布了一簇簇连续纤维的复合结构。对于纤维增强复合材料，［9］提出了由5个应变不变量表示横向同性超弹性应变能函数的方法，该方法被一些研究者应用于生物软组织材料的本构模型研究中，并发现行之有效［10］。

质点在参考构形中的位置位于X，t时刻在当前构形中位于x，则变形梯度，右Cauchy－Green变形张量C＝FTF。C的5个应变不变量可定义为

式中：N0＝n0?n0为对称的二阶张量，单位方向矢量n0（X）表示各向异性材料的对称轴方向。其中I1、I2和I3表示材料的各向同性，在不可压假定下，I3＝1。不变量I4和I5用来表示材料的各向异性。则不可压纤维增强材料的应变能函数可表示为

通常可将式（2）进行解耦，以便与材料的结构特征相联系，如

式中：W1代表各向同性基体组织的力学性能，W2用来代表具有纤维集合的增强作用，W3表示纤维和基体之间的相互作用。

2 粘弹性性能

粘弹性材料最大的特点就是力学性能受之前变形的影响，应力依赖应变和应变率历史，对于不可压各向同性材料的粘弹性本构模型通常用下式来表示

式中：σv、pv分别为表征粘弹性性能的Cauchy应力和静水压力；Ω是张量泛函，反映应变历史对应力的影响，前人已经提出了多种形式的泛函的近似表达式，如BKZ模型、ZWT模型，这些模型通常采用积分形式来表示。

对于软组织材料，等［11］曾采用类似的积分形式来表示横向同性粘弹性性能，即

式中：φ是一个包含横向同性不变量的张量函数。不可压条件下，）通常定义为一个随时间衰减的函数，即这里α和θ为材料参数，θ通常被称为松弛时间。则有

3 肌肉的本构模型

按照前文的本构理论，将描述准静态下的本构模型与描述粘弹性性能的本构模型进行简单的相加，以表示肌肉的动态本构模型，那么总的Cauchy应力可分为2部分，一部分表征材料的准静态超弹性，一部分表征材料的粘弹性，即