王寶珍，胡時勝

（中國科學(xué)技術(shù)大學(xué)中國科學(xué)院材料力學(xué)行為和設(shè)計重點(diǎn)實驗室，安徽 合肥 230026）

在汽車安全、軍事工程和醫(yī)療診斷中，都需要對人體受沖擊狀態(tài)進(jìn)行分析，研究人體組織的沖擊力學(xué)性能，并發(fā)展能描述其動態(tài)性能的本構(gòu)關(guān)系，以得到相關(guān)的參數(shù)，對于安全設(shè)計等方面有著重要的意義。肌肉是人體組織中一個重要的組成部分，按性質(zhì)的不同通常分為3類：骨骼肌、平滑肌和心肌。骨骼肌在人體中占有很大的比重，軀體的各種運(yùn)動和呼吸均離不開骨骼肌。骨骼肌由許多平行排列的肌纖維及纖維間的結(jié)締組織組成，其力學(xué)性能由肌纖維和結(jié)締組織共同決定。與其他工程材料有很大不同的是，它可以產(chǎn)生2種不同形式的力，即主動力和被動力。

研究骨骼肌在外部沖擊載荷作用下的被動承載能力，對人體安全設(shè)計有重要意義。但由于肌肉組織具有一定的活性，試樣獲取困難，且極其柔軟，給實驗特別是動態(tài)實驗帶來了很大的困難，因而相關(guān)研究很少。B.Song等［1］、C.Van Sligtenhorst等［2］曾利用SHPB裝置獲得了肌肉的動態(tài)壓縮應(yīng)力應(yīng)變曲線，但其動態(tài)拉伸性能還無相關(guān)報道。為彌補(bǔ)這一空白，王寶珍等［3－4］利用改進(jìn)的SHPB裝置和套筒式SHTB裝置對豬后腿肌肉的動態(tài)壓縮和拉伸性能進(jìn)行了研究，獲得了不同應(yīng)變率下的應(yīng)力應(yīng)變曲線。然而，描述肌肉材料動態(tài)力學(xué)性能的本構(gòu)模型還需要進(jìn)一步討論。

關(guān)于骨骼肌的本構(gòu)模型，以往的研究往往集中于描述肌纖維的主動收縮性能，如著名的Hill模型。而為了描述骨骼肌的被動力學(xué)性能，通常結(jié)合其特殊的纖維組織結(jié)構(gòu)，將骨骼肌看作橫向同性體。M.Van Loocke等［5］曾借用橫向同性彈性理論，并引入應(yīng)變相關(guān)楊氏模量，從而描述骨骼肌靜態(tài)壓縮下的非線性橫向同性性能。J.D.Humphrey等［6］則從另一種思路出發(fā)，將肌肉看作纖維增強(qiáng)復(fù)合材料，在參考構(gòu)形中，引入一個纖維方向的單位矢量，由該方向矢量與右Cauchy－Green變形張量可表述出5個應(yīng)變不變量，他們將應(yīng)變能密度函數(shù)用這5個應(yīng)變不變量來表示，從而發(fā)展了描述心肌被動力學(xué)性能的橫向同性超彈性本構(gòu)模型，后被許多學(xué)者用于研究骨骼?。?－8］的橫向同性超彈性。然而，由于一直以來缺乏肌肉動態(tài)性能研究的實驗結(jié)果，因而未能得到一個描述肌肉動態(tài)性能的率相關(guān)本構(gòu)模型。

本文中，將結(jié)合已有的豬后腿肌肉動態(tài)壓縮和拉伸實驗結(jié)果，發(fā)展一個既能描述肌肉動態(tài)壓縮性能又能描述其動態(tài)拉伸性能的率相關(guān)本構(gòu)模型。

1 纖維增強(qiáng)材料應(yīng)變能函數(shù)表示方法

肌肉、韌帶等生物軟組織通?？煽醋魇侨嵝曰|(zhì)中分布了一簇簇連續(xù)纖維的復(fù)合結(jié)構(gòu)。對于纖維增強(qiáng)復(fù)合材料，A.T.M.Spencer［9］提出了由5個應(yīng)變不變量表示橫向同性超彈性應(yīng)變能函數(shù)的方法，該方法被一些研究者應(yīng)用于生物軟組織材料的本構(gòu)模型研究中，并發(fā)現(xiàn)行之有效［10］。

式中：N0＝n0?n0為對稱的二階張量，單位方向矢量n0（X）表示各向異性材料的對稱軸方向。其中I1、I2和I3表示材料的各向同性，在不可壓假定下，I3＝1。不變量I4和I5用來表示材料的各向異性。則不可壓纖維增強(qiáng)材料的應(yīng)變能函數(shù)可表示為

通?？蓪⑹剑?）進(jìn)行解耦，以便與材料的結(jié)構(gòu)特征相聯(lián)系，如

式中：W1代表各向同性基體組織的力學(xué)性能，W2用來代表具有纖維集合的增強(qiáng)作用，W3表示纖維和基體之間的相互作用。

2 粘彈性性能

粘彈性材料最大的特點(diǎn)就是力學(xué)性能受之前變形的影響，應(yīng)力依賴應(yīng)變和應(yīng)變率歷史，對于不可壓各向同性材料的粘彈性本構(gòu)模型通常用下式來表示

式中：σv、pv分別為表征粘彈性性能的Cauchy應(yīng)力和靜水壓力；Ω是張量泛函，反映應(yīng)變歷史對應(yīng)力的影響，前人已經(jīng)提出了多種形式的泛函的近似表達(dá)式，如BKZ模型、ZWT模型，這些模型通常采用積分形式來表示。

對于軟組織材料，Z.A.Taylor等［11］曾采用類似的積分形式來表示橫向同性粘彈性性能，即

3 肌肉的本構(gòu)模型

按照前文的本構(gòu)理論，將描述準(zhǔn)靜態(tài)下的本構(gòu)模型與描述粘彈性性能的本構(gòu)模型進(jìn)行簡單的相加，以表示肌肉的動態(tài)本構(gòu)模型，那么總的Cauchy應(yīng)力可分為2部分，一部分表征材料的準(zhǔn)靜態(tài)超彈性，一部分表征材料的粘彈性，即

3.1 肌肉準(zhǔn)靜態(tài)本構(gòu)模型

取n0＝ （1 ，0 ，0） ，對于纖維向加載，假定纖維向伸長比為λ，則3個主伸長比為λ1＝λ，λ2＝λ3＝λ－1／2。 變形梯度F和右Cauchy－Green變形張量C可表示為

5個基本不變量可分別表示為

選用Yeoh模型［12］來描述基體的各向同性性能，即

假定肌纖維不能承受壓縮載荷而僅能承受拉伸載荷，選用這樣一個描述纖維伸長的應(yīng)變能函數(shù)

式中：符號 〈x〉表示當(dāng)x≤0時，〈x〉＝0，當(dāng)x＞0時，〈x〉＝x。

對于纖維與基體的相互作用部分，參考J.C.Criscione等［13］提出的2個表示纖維向剪切和橫向剪切性能的新變形不變量，并將其改寫為

纖維與基體相互作用的那部分應(yīng)變能函數(shù)表示為

將式（9）代入式（12），有

因此

則表征準(zhǔn)靜態(tài)超彈性性能的Cauchy應(yīng)力分量可以表示為

圖1 準(zhǔn)靜態(tài)肌肉纖維方向壓縮和拉伸的應(yīng)力應(yīng)變曲線Fig.1Compressive and tensile stress－strain curves of muscle under quasi－static loading along the fiber direction

采用式（18）對纖維向單軸壓縮和單軸拉伸曲線進(jìn)行擬合，如圖1所示?？傻玫絽?shù)C10＝821Pa，C20＝1.81kPa，C30＝－72.4Pa，a＝7.09kPa，b＝13.94。發(fā)現(xiàn)實驗曲線和擬合曲線有很好的一致性。

3.2 動態(tài)本構(gòu)模型

對纖維向動態(tài)壓縮，仍假定肌纖維不承壓，式（6）表示為

這里取

由實驗得到的應(yīng)力應(yīng)變曲線，發(fā)現(xiàn)3 200s－1下的應(yīng)力應(yīng)變曲線與另2種應(yīng)變率（1 100、2 000s－1）下的應(yīng)力應(yīng)變曲線間距較大，即在高應(yīng)變率區(qū)，應(yīng)變率增加時應(yīng)變率敏感性增加，為描述這一現(xiàn)象，令

式中：為伸長速率，壓縮為負(fù)，拉伸為正；A1、A2和均為材料參數(shù)。由式（19）～（21）可得到

同樣，由σ22＝σ33＝0，可以得到

則動態(tài)單軸壓縮時沿纖維方向的工程應(yīng)力σ（＝σ11／λ）為

圖2 高應(yīng)變率下肌肉沿纖維方向壓縮的應(yīng)力應(yīng)變曲線Fig.2Compressive stress－strain curves at high strain rates along the fiber direction

式中：由式（18）給出，采用上式對高應(yīng)變率下纖維方向的動態(tài)壓縮工程應(yīng)力－伸長比曲線進(jìn)行擬合，如圖2，參數(shù)A1＝12.1kPa，A2＝15.4kPa＝1 237s－1，θ＝67.11μs。結(jié)果發(fā)現(xiàn)，擬合曲線與實驗結(jié)果一致性較好。

為了描述纖維向動態(tài)拉伸的情況，需對式（25）進(jìn)行修改，即在等號右側(cè)增加2項：一項考慮肌纖維在拉伸時對粘彈性性能的貢獻(xiàn)，另一項考慮動態(tài)拉伸應(yīng)力應(yīng)變曲線與動態(tài)壓縮曲線形狀的不同。動態(tài)拉伸曲線的初始階段不存在低應(yīng)力幅值區(qū)，而是有一類似屈服點(diǎn)的轉(zhuǎn)折點(diǎn)，且動態(tài)實驗結(jié)果在小應(yīng)變時也并不可靠，因此擬合曲線時并非從零點(diǎn)開始，而是僅考慮該點(diǎn)以后的情況。最終得到動態(tài)單軸拉伸沿纖維方向的工程應(yīng)力－伸長比本構(gòu)關(guān)系為

式中：f（I4）＝A3〈I4－1〉／2，那么

圖3 高應(yīng)變率下肌肉沿纖維方向拉伸的應(yīng)力應(yīng)變曲線Fig.3Tensile stress－strain curves at high strain rates along the fiber direction

采用上式對高應(yīng)變率下纖維方向的動態(tài)應(yīng)力應(yīng)變曲線進(jìn)行擬合，如圖3所示，參數(shù)A3＝300kPa，A4＝3kPa，A5＝3.09，發(fā)現(xiàn)理論模型與實驗結(jié)果有很好的一致性。

4 結(jié) 論

肌肉軟組織可看作是一種纖維增強(qiáng)復(fù)合材料。應(yīng)用相關(guān)本構(gòu)理論，提出了一個既能描述動態(tài)壓縮力學(xué)性能又能描述動態(tài)拉伸力學(xué)性能的橫向同性率相關(guān)本構(gòu)模型，該模型由描述準(zhǔn)靜態(tài)響應(yīng)的超彈性模型和描述高應(yīng)變率響應(yīng)的粘彈性模型組成，理論曲線與實驗得到的應(yīng)力應(yīng)變曲線有很好的一致性。說明提出的模型能夠很好地描述肌肉的力學(xué)性能，可為人體安全防護(hù)設(shè)計提供理論依據(jù)。

［1］Song B，Chen W，Ge Y，et al.Dynamic and quasi－static compressive response of porcine muscle［J］.Journal of Biomechanics，2007，40（13）：2999－3005.

［2］Sligtenhorst V C，Cronin D S，Brodland G W.High strain rate compressive properties of bovine muscle tissue determined using a split Hopkinson bar apparatus［J］.Journal of Biomechanics，2006，39（10）：1852－1858.

［3］王寶珍，胡時勝.豬后腿肌肉的沖擊壓縮特性實驗［J］.爆炸與沖擊，2010，30（1）：33－38.

WANG Bao－zhen，HU Shi－sheng.Dynamic compression experiments of porcine ham muscle［J］.Explosion and Shock Waves，2010，30（1）：33－38.

［4］王寶珍，鄭宇軒，胡時勝.豬后腿肌肉的動態(tài)拉伸性能研究［J］.爆炸與沖擊，2010，30（5）：449－455.

WANG Bao－zhen，ZHENG Yu－xuan，HU Shi－sheng.Dynamic tensile properties of porcine ham muscle［J］.Explosion and Shock Waves，2010，30（5）：449－455.

［5］Van Loocke M，Lyons C G，Simms C K.Viscoelastic properties of passive skeletal muscle in compression：Stressrelaxation behaviour and constitutive modelling［J］.Journal of Biomechanics，2008，41（7）：1555－1566.

［6］Humphrey J D，Strumpf R K，Yin F C P.Determination of a constitutive relation for passive myocardium［J］.Journal of Biomechanical Engineering，1990，112：333－346.

［7］Gielen A W，Oomens C W，Bovendeerd P H，et al.A finite element approach for skeletal muscle using a distributed moment model of contraction［J］.Computer Methods in Biomechanics ＆Biomedical Engineering，2000，3（3）：231－244.

［8］Martins J A C，Pires E B，Salvado R，et al.A numerical model of passive and active behavior of skeletal muscles［J］.Computer Methods in Applied Mechanics and Engineering，1998，151（3）：419－433.

［9］Spencer A J M.Continuum theory of the mechanics of fibre－reinforced composites［M］.New York：Springer－Verlag，1984：1－32.

［10］Limbert G，Middleton J.A constitutive model of the posterior cruciate ligament［J］.Medical Engineering ＆Phys－ics，2006，28（2）：99－113.

［11］Taylor Z A，Comas O，Cheng M，et al.On modelling of anisotropic viscoelasticity for soft tissue simulation：Numerical solution and GPU execution［J］.Medical Image Analysis，2009，13（2）：234－244.

［12］Yeoh O H.Characterization of elastic properties of carbon－black－filled rubber vulcanizates［J］.Rubber Chemistry and Technolog，1990，63（5）：792－805.

［13］Criscione J C，Douglas A S，Hunter W C.Physically based strain invariant set for materials exhibiting transversely isotropic behavior［J］.Journal of the Mechanics and Physics of Solids，2001，49（4）：871－897.