趙永玲，侯之超，黃友劍，劉建勛，

(1.清華大學(xué)汽車安全與節(jié)能國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京100084;2.株洲時(shí)代新材料科技股份有限公司，湖南株洲412007)

橡膠材料的一種5參數(shù)分?jǐn)?shù)導(dǎo)數(shù)模型

趙永玲1，侯之超1，黃友劍2，劉建勛1，2

(1.清華大學(xué)汽車安全與節(jié)能國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京100084;2.株洲時(shí)代新材料科技股份有限公司，湖南株洲412007)

通過(guò)在分?jǐn)?shù)導(dǎo)數(shù)Zener模型中增加Scott－Blair黏壺，建立了一種新的5參數(shù)分?jǐn)?shù)導(dǎo)數(shù)模型。分析了分?jǐn)?shù)導(dǎo)數(shù)算子參數(shù)對(duì)Cole－Cole曲線形狀的影響;基于某隔振橡膠材料的DMA實(shí)驗(yàn)，應(yīng)用溫—頻等效原理獲得了該材料的動(dòng)態(tài)模量和損耗因子的主曲線，并識(shí)別了新建模型的全部參數(shù)。結(jié)果表明，該模型能夠比較準(zhǔn)確地反映具有對(duì)稱或非對(duì)稱Cole－Cole曲線的橡膠材料的動(dòng)態(tài)特性。

Cole－Cole曲線;黏彈性;分?jǐn)?shù)導(dǎo)數(shù)模型;動(dòng)態(tài)力學(xué)實(shí)驗(yàn);參數(shù)識(shí)別

減振橡膠，作為一種黏彈性材料，具有緩沖減振等作用，在汽車運(yùn)輸、軌道交通、國(guó)防工業(yè)、民用建筑、醫(yī)療器械等領(lǐng)域均具有廣泛的應(yīng)用。隨著各行業(yè)對(duì)結(jié)構(gòu)可靠性和舒適性要求的提高，對(duì)橡膠元件進(jìn)行精細(xì)的結(jié)構(gòu)設(shè)計(jì)日益普及。為此，迫切需要對(duì)橡膠元件動(dòng)態(tài)特性進(jìn)行準(zhǔn)確表征與仿真分析。

橡膠材料的動(dòng)態(tài)特性一般采用儲(chǔ)能模量E'、損耗模量E″和損耗因子tanδ(=E'/E″)描述。其中，E'和E″二者之間的關(guān)系，可用Cole－Cole曲線［1］表示。作為一種描述橡膠材料動(dòng)態(tài)特性的一種有效方式，與其它參數(shù)一樣，黏彈性材料的Cole－Cole曲線與材料的分子結(jié)構(gòu)、加工成型條件、溫度和頻率等密切相關(guān)，其中溫度和頻率等因素至關(guān)重要。隨著溫度的改變，若材料內(nèi)部分子鏈只存在主轉(zhuǎn)變(玻璃化轉(zhuǎn)變)時(shí)，可以根據(jù)不同溫度下的E'和E″繪制一條光滑的Cole－Cole曲線:當(dāng)分子結(jié)構(gòu)為均一體系時(shí)，曲線具有對(duì)稱特性;當(dāng)分子結(jié)構(gòu)為共混體系時(shí)，曲線為非對(duì)稱［2］。

建立合適的黏彈性材料動(dòng)力學(xué)本構(gòu)模型是準(zhǔn)確預(yù)測(cè)橡膠元件動(dòng)態(tài)特性的前提。由于分?jǐn)?shù)導(dǎo)數(shù)模型在較寬的溫度和頻率范圍內(nèi)能夠很好地描述黏彈性材料的動(dòng)態(tài)特性，因此廣泛應(yīng)用于建立黏彈性材料的本構(gòu)模型。Pritz等［3－4］采用分?jǐn)?shù)導(dǎo)數(shù)Zener模型討論了黏彈性材料在較寬頻率范圍內(nèi)的力學(xué)特性。詹小麗等［5］采用分?jǐn)?shù)導(dǎo)數(shù)Maxwell模型研究了改性瀝青材料的動(dòng)態(tài)黏彈特性。趙永玲等［6］采用高階分?jǐn)?shù)導(dǎo)數(shù)模型對(duì)橡膠材料動(dòng)態(tài)力學(xué)性能進(jìn)行評(píng)估。其中分?jǐn)?shù)導(dǎo)數(shù)Zener模型被證明比較適合預(yù)測(cè)黏彈性材料在寬頻范圍內(nèi)的動(dòng)態(tài)特性，并且模型中的參數(shù)都具有明確的物理意義。但是這種模型只能描述具有對(duì)稱Cole－Cole曲線的黏彈性材料的動(dòng)態(tài)特性。

為描述具有非對(duì)稱Cole－Cole曲線的黏彈性材料的動(dòng)態(tài)特性，需要增加模型參數(shù)的數(shù)量。唐振寰等［7］采用應(yīng)力和應(yīng)變的分?jǐn)?shù)算子階次不相等的方法建立了一種五參數(shù)分?jǐn)?shù)導(dǎo)數(shù)模型。但是該模型不滿足熱力學(xué)定律。Pritz［8］、侯宏等［9］分別通過(guò)增加關(guān)于應(yīng)變的分?jǐn)?shù)導(dǎo)數(shù)項(xiàng)而構(gòu)建了一種5參數(shù)分?jǐn)?shù)導(dǎo)數(shù)模型，其中部分參數(shù)的物理意義還有待明確。

鑒于上述現(xiàn)狀，本文以分?jǐn)?shù)導(dǎo)數(shù)Zener模型為基礎(chǔ)，增加一個(gè)Scott－Blair(SB)黏壺［10］而建立了一種五參數(shù)分?jǐn)?shù)導(dǎo)數(shù)模型，并分析了分?jǐn)?shù)導(dǎo)數(shù)算子常數(shù)對(duì)Cole－Cole曲線形狀的影響。針對(duì)只存在主轉(zhuǎn)變(玻璃化轉(zhuǎn)變)的某橡膠材料，在不同溫度下進(jìn)行DMA頻率掃描實(shí)驗(yàn)，獲得了材料的儲(chǔ)能模量、損耗模量和損耗因子等實(shí)驗(yàn)數(shù)據(jù)，并利用溫—頻等效原理，得到寬頻范圍內(nèi)的主曲線［11］。進(jìn)而應(yīng)用試驗(yàn)數(shù)據(jù)，識(shí)別得到了所建立的5參數(shù)分?jǐn)?shù)導(dǎo)數(shù)模型的全部模型參數(shù)。通過(guò)對(duì)比分析，證明所提出模型可以描述具有對(duì)稱或非對(duì)稱Cole－Cole曲線的橡膠材料的動(dòng)態(tài)特性。

1 5參數(shù)分?jǐn)?shù)導(dǎo)數(shù)模型

分?jǐn)?shù)導(dǎo)數(shù)模型中，分?jǐn)?shù)導(dǎo)數(shù)算子的階次和項(xiàng)數(shù)反映材料黏彈性的強(qiáng)弱。依據(jù)Scott－Blair黏壺單元的應(yīng)力應(yīng)變關(guān)系，在頻域中導(dǎo)出對(duì)應(yīng)的存儲(chǔ)模量與損耗模量表達(dá)式。由此可以看出，當(dāng)頻率不變時(shí)，Scott－Blair黏壺單元的微分階算子階次決定其Cole－Cole曲線的斜率。顯然，對(duì)于非對(duì)稱的Cole－Cole曲線，材料模型中至少需要兩個(gè)Scott－Blair黏壺單元。因此，本文首先借鑒文獻(xiàn)［7－9］增加分?jǐn)?shù)導(dǎo)數(shù)項(xiàng)的作法，在分?jǐn)?shù)導(dǎo)數(shù)Zener模型中增加一個(gè)Scott－Blair黏壺，構(gòu)建出一個(gè)新的5參數(shù)分?jǐn)?shù)導(dǎo)數(shù)模型，以描述具有非對(duì)稱性Cole－Cole曲線的黏彈性特性。然后分析各分?jǐn)?shù)導(dǎo)數(shù)階次對(duì)動(dòng)態(tài)特性的影響，并通過(guò)模型退化闡述模型的適應(yīng)性。

1.2 5參數(shù)本構(gòu)方程與模量

在分?jǐn)?shù)導(dǎo)數(shù)Zener模型中增加一個(gè)Scott－Blair黏壺，即構(gòu)成一個(gè)新的模型(見圖1)。模型中含有5個(gè)參數(shù)α、β、E0、E∞和τ。其中α和β分別為兩個(gè)黏壺的分?jǐn)?shù)導(dǎo)數(shù)階次，E0為靜態(tài)模量，E∞為高頻動(dòng)態(tài)模量的極限值，τ為松弛時(shí)間。兩個(gè)黏壺對(duì)應(yīng)的應(yīng)變分別為ε3和ε2，兩個(gè)彈簧對(duì)應(yīng)的應(yīng)變分別為ε1和ε4。

5參數(shù)分?jǐn)?shù)導(dǎo)數(shù)模型(見圖1)，依據(jù)應(yīng)力應(yīng)變關(guān)系，得到本構(gòu)方程:

圖1 5參數(shù)分?jǐn)?shù)導(dǎo)數(shù)模型Fig.1 A fractional derivativemodel with five parameters

為滿足熱力學(xué)定律，應(yīng)該滿足0＜β≤α＜1。由于方程中可能包含導(dǎo)數(shù)階次大于1的項(xiàng)，因此該模型可稱為高階分?jǐn)?shù)導(dǎo)數(shù)模型。Dα為第α階分?jǐn)?shù)導(dǎo)數(shù)，可以定義如下:

式中:Γ為Gamma函數(shù)。

對(duì)式(1)作Fourier變換，可得到該模型所描述材料的復(fù)模量表達(dá)式:

式中:ωn=ωτ，ω=2πf。由此可以求得材料的儲(chǔ)能模量、損耗模量和損耗因子:

上述方程反映了材料各性能參數(shù)隨頻率f的變化關(guān)系。由式(4)可知，E'(ω→0)=E0，E'(ω→∞)= E∞。由于溫度和頻率對(duì)于黏彈性材料特性的影響可以通過(guò)溫度－頻率等效原理來(lái)轉(zhuǎn)化，因此式(4)～式(6)也可用于揭示E'、E″和tanδ隨溫度的變化規(guī)律。

1.3 分?jǐn)?shù)導(dǎo)數(shù)階次與Cole-Cole曲線

根據(jù)式(4)和式(5)可以計(jì)算出不同頻率時(shí)的儲(chǔ)能模量E'和損耗模量E″。以E'為橫坐標(biāo)E″為縱坐標(biāo)繪制曲線，即得到Cole－Cole曲線(見圖2)。由圖2可知，α≠β時(shí)，Cole－Cole曲線均為非對(duì)稱曲線。圖2(a)為α變化β不變時(shí)的Cole－Cole曲線。隨著α的變化，低頻時(shí)Cole－Cole曲線的斜率變化較大，高頻時(shí)其斜率均趨于一定值。這意味著參數(shù)α與低頻時(shí)的動(dòng)態(tài)特性有關(guān)，亦即α反映低頻時(shí)Cole－Cole曲線的變化率。圖2(b)為α不變?chǔ)伦兓瘯r(shí)的Cole－Cole曲線。隨著β的變化，低頻時(shí)Cole－Cole曲線的斜率變化不大，高頻時(shí)儲(chǔ)能模量和損耗因子的斜率隨β值的變化比較明顯，即β反映高頻時(shí)Cole－Cole曲線的變化率。這樣，5參數(shù)分?jǐn)?shù)導(dǎo)數(shù)模型中各參數(shù)的物理意義比較明確。

圖2 Cole－Cole曲線Fig.2 Cole－Cole curve

1.4 退化模型

對(duì)于所建立的5參數(shù)分?jǐn)?shù)導(dǎo)數(shù)模型，如果分?jǐn)?shù)階算子滿足α=β，式(1)就簡(jiǎn)化為:

另一方面，分?jǐn)?shù)導(dǎo)數(shù)Zener模型的本構(gòu)方程可表示為:

式(7)和式(8)二者的差別僅在于σ與ε一次項(xiàng)的常系數(shù)不同。換言之，當(dāng)α=β時(shí)，本文構(gòu)建的5參數(shù)分?jǐn)?shù)導(dǎo)數(shù)模型可以退化為與分?jǐn)?shù)導(dǎo)數(shù)Zener類似的模型。這說(shuō)明新建立的模型也適合于Cole－Cole曲線為對(duì)稱時(shí)的黏彈性材料，如圖2(a)中α=β=0.4和圖2 (b)中α=β=0.7兩種情況。

因此，當(dāng)材料具有對(duì)稱Cole－Cole曲線時(shí)，可應(yīng)用式(6)計(jì)算材料的儲(chǔ)能模量、損耗模量和損耗因子:

2 橡膠材料DMA實(shí)驗(yàn)

為了驗(yàn)證前文所建立的5參數(shù)分?jǐn)?shù)導(dǎo)數(shù)模型，對(duì)某減振橡膠進(jìn)行動(dòng)態(tài)機(jī)械分析(Dynamic Mechanical Analytic，DMA)實(shí)驗(yàn)。首先進(jìn)行恒定頻率下的溫度掃描實(shí)驗(yàn)，確定材料的玻璃化轉(zhuǎn)變區(qū)域;再通過(guò)不同恒定溫度下的頻率掃描實(shí)驗(yàn)，獲取材料的儲(chǔ)能模量、損耗模量和損耗因子隨頻率的變化曲線;然后利用溫－頻等效原理將這些曲線進(jìn)行平移，繪制出參考溫度下的主曲線。

2.1 DMA實(shí)驗(yàn)

采用德國(guó)GABO公司制造的EPLEXOR 500N型動(dòng)態(tài)熱力學(xué)分析儀(見圖3)，測(cè)定橡膠材料在不同溫度的一組頻率譜。橡膠試樣尺寸為20 mm×6 mm× 2 mm，加載方式采用軸向拉伸模式。

圖3 EPLEXOR 500N型動(dòng)態(tài)熱力學(xué)分析儀Fig.3 EPLEXOR 500N dynamic thermal analyzer

下面分別從溫度掃描和頻率掃描兩個(gè)方面對(duì)測(cè)試數(shù)據(jù)進(jìn)行分析，以闡述材料的力學(xué)性能。

2.1.1 溫度掃描

溫度掃描實(shí)驗(yàn)中選取的溫度范圍越寬，則通過(guò)溫－頻等效原理得到的主曲線跨越的頻率范圍就越寬。溫度范圍的選取可以通過(guò)溫度掃描實(shí)驗(yàn)來(lái)確定。圖4為頻率為2 Hz時(shí)的溫度掃描曲線，由圖4可知，此橡膠材料在－80～60℃范圍內(nèi)只存在主轉(zhuǎn)變(也稱玻璃化轉(zhuǎn)變);而在－60～－40℃之間E'、E″和tanδ變化均較大，此溫度范圍稱為玻璃化轉(zhuǎn)變區(qū)域。玻璃化轉(zhuǎn)變溫度Tg可以采用儲(chǔ)能模量轉(zhuǎn)折點(diǎn)對(duì)應(yīng)的溫度來(lái)定義，也可以采用損耗模量和損耗因子峰值對(duì)應(yīng)的溫度作為玻璃化轉(zhuǎn)變溫度。這三種方法得到的玻璃化溫度稍有不同，本文研究中以損耗因子峰值對(duì)應(yīng)的溫度定義為所研究橡膠材料的玻璃化轉(zhuǎn)變溫度，由圖4可知為－50℃。

圖4 溫度掃描曲線Fig.4 Curves from temperature scanning

2.1.2頻率掃描

由圖4可知，＜－60℃和＞20℃時(shí)E'、E″和tanδ曲線均趨于穩(wěn)定。因此選?。?0～20℃之間的溫度分別作為實(shí)驗(yàn)溫度，進(jìn)行頻率掃描，頻率掃描范圍為0.5～15 Hz。圖5為在恒定溫度分別為－60℃、－55℃、－50℃、－45℃、－40℃、－30℃、－20℃、－10℃、0℃、10℃、20℃時(shí)頻率掃描實(shí)驗(yàn)得到的E'、E″和tanδ曲線。從圖5(a)可知，儲(chǔ)能模量隨著頻率的升高而增大、隨著溫度的升高而減小。由圖5(b)可知，－60℃時(shí)損耗模量隨著頻率的升高而減小，T＞－55℃時(shí)，損耗模量隨著頻率的升高而增大，說(shuō)明在－60～－55℃時(shí)損耗模量具有最大值。從圖5(c)可知，當(dāng)T＜－55℃時(shí)，損耗因子隨著頻率的升高而減小;當(dāng)T＞－45℃時(shí)，損耗因子隨著頻率的升高而增大，說(shuō)明在－55～－45℃時(shí)損耗因子具有最大值，與溫度掃描得到的玻璃化轉(zhuǎn)變溫度是一致的。

圖5 頻率掃描實(shí)驗(yàn)結(jié)果Fig.5 Curves from frequency scanning

2.2 繪制主曲線

為準(zhǔn)確描述橡膠材料在寬頻率范圍內(nèi)的動(dòng)態(tài)特性，將經(jīng)由DMA實(shí)驗(yàn)得到的不同溫度下的頻率掃描數(shù)據(jù)轉(zhuǎn)換為參考溫度下的黏彈性函數(shù)主曲線。為此，需要應(yīng)用合適的溫－頻等效模型。目前最常用的溫－頻等效模型是WLF模型，其具體形式為:

式中:αT為溫－頻轉(zhuǎn)換因子，T0為參考溫度，C1和C2是經(jīng)驗(yàn)常數(shù)。選擇不同溫度作為參考溫度，該式的形式不變，只是參數(shù)C1和C2不同。

選玻璃化轉(zhuǎn)變溫度－50℃作為參考溫度，利用非線性最小二乘法擬合獲得相對(duì)該參考溫度的溫頻轉(zhuǎn)換因子和經(jīng)驗(yàn)常數(shù)。其中，材料常數(shù)C1=7.2，C2=37.5，而不同溫度下的溫－頻轉(zhuǎn)換因子見圖6。

根據(jù)溫－頻轉(zhuǎn)換因子繪制橡膠材料主曲線(見圖7)。由圖7可知，不同溫度下的儲(chǔ)能模量、損耗模量和損耗因子經(jīng)過(guò)平移后，得到參考溫度下的光滑主曲線，且頻率跨越接近10個(gè)數(shù)量級(jí)。

圖6 不同溫度下的溫－頻轉(zhuǎn)換因子Fig.6 Temperature－frequency conversion factor under different temperatures

圖7 主曲線Fig.7 Themaster curve

圖8 Cole－Cole曲線Fig.8 The Cole－Cole curve

3 模型參數(shù)識(shí)別

根據(jù)圖7所示主曲線中的儲(chǔ)能模量E'和損耗模量E″的數(shù)值，可繪制出Cole－Cole曲線(見圖8)。

從圖8可知，Cole－Cole曲線具有明顯的非對(duì)稱性。采用前文構(gòu)建的5參數(shù)分?jǐn)?shù)導(dǎo)數(shù)模型進(jìn)行參數(shù)識(shí)別:將11組不同溫度下的頻率掃描試驗(yàn)數(shù)據(jù)代入式(4)～式(6)，采用非線性最小二乘法擬合得到全部模型參數(shù)(見表1)。

表1 分?jǐn)?shù)導(dǎo)數(shù)模型擬合參數(shù)Tab.1 The fitting param etersof fractional derivativemodel

應(yīng)用識(shí)別獲得的模型參數(shù)，根據(jù)式(4)～式(6)計(jì)算所研究橡膠材料的儲(chǔ)能模量、損耗模量和損耗因子。預(yù)測(cè)精度應(yīng)用下式定義的相對(duì)誤差進(jìn)行評(píng)價(jià):

式中:EiE和EiP分別為第i個(gè)實(shí)驗(yàn)值和擬合值，n為擬合點(diǎn)數(shù)。計(jì)算顯示，儲(chǔ)能模量、損耗模量和損耗因子的擬合誤差分別為2.27%、8.72%、7.51%，均不超過(guò)10%，滿足工程精度要求。

繪制E'－lg f、E″－lg f、tanδ－lg f曲線和Cole－Cole曲線，并與試驗(yàn)曲線對(duì)比(見圖9)。由圖9可知，5參數(shù)分?jǐn)?shù)導(dǎo)數(shù)模型計(jì)算得到的儲(chǔ)能模量、損耗模量和損耗因子，與試驗(yàn)數(shù)據(jù)吻合良好。這驗(yàn)證了模型的合理性。

圖9 模型預(yù)測(cè)結(jié)果與試驗(yàn)數(shù)據(jù)的對(duì)比Fig.9 Comparison between model prediction and experimental data

4 結(jié)論

本文提出了橡膠材料的一種分?jǐn)?shù)導(dǎo)數(shù)黏彈性本構(gòu)模型。針對(duì)某種橡膠材料，設(shè)計(jì)并進(jìn)行了不同溫度下的DMA頻率掃描實(shí)驗(yàn)，應(yīng)用實(shí)測(cè)數(shù)據(jù)，繪制了材料的儲(chǔ)能模量、損耗模量及損耗因子的主曲線;進(jìn)而通過(guò)非線性最小二乘擬合，識(shí)別得到了全部模型參數(shù)，模型預(yù)測(cè)結(jié)果與實(shí)驗(yàn)數(shù)據(jù)吻合良好。分析與計(jì)算顯示，所提出的分?jǐn)?shù)導(dǎo)數(shù)模型既適合具有對(duì)稱Cole－Cole曲線的材料，也適合具有非對(duì)稱Cole－Cole曲線的材料。利用該模型可以預(yù)測(cè)不同頻率(或溫度)下的黏彈性材料動(dòng)力學(xué)特性，具有重要的應(yīng)用價(jià)值。

［1］范仰才，李景德，符德勝.廣義導(dǎo)納圓［J］.中山大學(xué)學(xué)報(bào):自然科學(xué)版，1995，34(3):31－35.

FAN Yang－cai，LI Jing－de，F(xiàn)U De－sheng.Generalized admittance circles［J］.Acta Scientiarum Naturalium Universitatis Sunyatseni:Natural Sciences，1995，34(3):31－35.

［2］于逢源.聚烯烴流動(dòng)誘導(dǎo)結(jié)晶的流變學(xué)研究［D］.上海:上海交通大學(xué)，2009.

［3］Pritz T.Analysis of four－parameter fractionalderivativemodel for polymeric damping materials［J］.Journal of Sound and Vibration，2003，265:935－953.

［4］Charles Lu.Characterize the high－frequency dynamic properties of elastomers using fractional calculus for FEM［C］// SAE2007 Noise and Vibration Conference and Exhibition，SAE Technical Paper，2007－01－2417.

［5］詹小麗，張肖寧，王端宜，等.改性瀝青非線性黏彈性本構(gòu)關(guān)系研究及應(yīng)用［J］.工程力學(xué)，2009，26(4):187－191.

ZHAN Xiao－li，ZHANG Xiao－ning，WANG Duan－yi，et al.Study on nonlinear viscoelastic constitutive equation ofmodified asphalt and its applications［J］.Engineering Mechanics，2009，26(4):187－191.

［6］趙永玲，侯之超.基于分?jǐn)?shù)導(dǎo)數(shù)的橡膠材料兩種黏彈性本構(gòu)模型［J］.清華大學(xué)學(xué)報(bào):自然科學(xué)版，2013，53(3): 378－383.

ZHAO Yong－ling，HOU Zhi－chao.Two viscoelastic constitutive models of rubber materials using fractional derivation［J］.Jounal of Tsinghua University:Science and Technology，2013，53(3):378－383.

［7］唐振寰，羅貴火，陳偉，等.橡膠隔振器黏彈性5參數(shù)分?jǐn)?shù)導(dǎo)數(shù)并聯(lián)動(dòng)力學(xué)模型［J］.航空動(dòng)力學(xué)報(bào)，2013，28(2): 275－282.

TANG Zhen－yuan，LUO Gui－h(huán)uo，CHENWei，et al.Parallel dynamic model of rubber isolator about five－parameter fractional derivative［J］.Journal of Aerospace Power，2013，28 (2):275－282.

［8］Pritz T.Five－parameter fractional derivative model for polymeric damping materials［J］.Journal of Sound and Vibration，2003，265:935－953.

［9］侯宏，高星，孫亮，等.變頻變溫條件黏彈材料參數(shù)化數(shù)學(xué)模型［J］.振動(dòng)與沖擊，2013，32(21):209－213.

HOU Hong，GAO Xing，SUN Liang，et al.Parameterized mathematicalmodel of a viscoelasticmaterial under variations of temperature and frequency［J］.Journal of Vibration and Shock，2013，32(21):209－213.

［10］Scott Blair G W.Analytical and integrative aspects of the stress－strain－time problem［J］.Journal of ScientificInstruments，1994，21:80－84.

［11］楊挺青，羅文波，徐平，等.黏彈性理論與應(yīng)用［M］.北京:科學(xué)出版社，2004.

A fractional derivativemodelw ith five parameters for rubber materials

ZHAO Yong－ling1，HOU Zhi－chao1，HUANG You－jian2，LIU Jian－xun1，2
(1.State Key Laboratory of Automotive Safety and Energy Saving，Tsinghua University，Beijing 100084，China; 2.Zhuzhou Time New Material Technology Co.，Ltd.，Zhuzhou 412007，China)

A fractional derivative model with five parameters was constructed by adding a Scott－Blair pot to the fractional derivative Zenermodel.Numerical analysiswas conducted to analyze the influence of parameters of the fractional derivative operators on Cole－Cole curves.Dynamic analysis testswere then performed on a rubbermaterial usually used for vibration isolation.Using the measured data，the storage modulus，loss modulus and loss factor of the material were derived by adopting the temperature－frequency equivalence principle，and all parameters of the new fractional derivative model were identified.The results demonstrated that the newly established fractional derivativemodelwith five parameters can accurately describe the dynamic characteristics of rubber materials with a symmetric or non－symmetric Cole－Cole curve.

Cole－Cole curve;visco－elasticity;fractional derivativemodel;dynamic test;parameter identification

TB303

A

10.13465/j.cnki.jvs.2015.23.007

中國(guó)南車科技資助項(xiàng)目(2011NCK098)

2014－09－15修改稿收到日期:2014－12－26

趙永玲女，博士生，1974年生

侯之超男，博士，教授，博士生導(dǎo)師，1966年生