鄒龍慶，曹義威，付海龍，王 玥

（東北石油大學(xué) 機(jī)械科學(xué)與工程學(xué)院，黑龍江 大慶市163000）

金屬橡膠阻尼性能好，可塑性強(qiáng)，可在惡劣的環(huán)境中保持良好性能，在航天航空、交通運(yùn)輸及工業(yè)生產(chǎn)等方面有著廣泛的應(yīng)用前景。金屬橡膠材料在加卸載過程中表現(xiàn)出較強(qiáng)的非線性遲滯特性，國(guó)內(nèi)外學(xué)者針對(duì)該特性進(jìn)行深入研究，并建立了多種力學(xué)模型，最具代表性的有跡法等效阻尼模型、雙折線恢復(fù)力模型和混合阻尼模型[1-2]等。其中將恢復(fù)力分解為彈性恢復(fù)力、黏性恢復(fù)力和干摩擦恢復(fù)力的模型能夠很好描述系統(tǒng)出現(xiàn)干摩擦?xí)r的情景，該類模型在實(shí)際工況中得到廣泛應(yīng)用。在眾多模型中，雙線性恢復(fù)力模型形式簡(jiǎn)單，參數(shù)識(shí)別精度較高，實(shí)際生產(chǎn)中應(yīng)用最為廣泛。該模型將黏性恢復(fù)力和干摩擦恢復(fù)力單獨(dú)表示，雖能夠很好地描述恢復(fù)機(jī)理，但模型中用線性方程表達(dá)黏性恢復(fù)力和干摩擦恢復(fù)力，很多情況下，不能準(zhǔn)確體現(xiàn)其非線性特征，特別在材料參數(shù)識(shí)別應(yīng)用時(shí)，受材料制備參數(shù)影響較大。有記憶型恢復(fù)力模型基于雙折線恢復(fù)力模型，將彈性恢復(fù)力改為5 次非線性方程，由非線性黏性阻尼力和干摩擦阻尼力構(gòu)造而成。有記憶恢復(fù)模型較之雙折線模型在描述高階非線性材料遲滯特性方面獲得較大進(jìn)步，可用于分析黏性恢復(fù)力和干摩擦恢復(fù)力所占比例及變化規(guī)律[3-4]。

有記憶型恢復(fù)力模型將恢復(fù)力分解為無記憶彈性力部分和有記憶阻尼力部分，考慮了變形對(duì)恢復(fù)力的影響，彌補(bǔ)了雙線性的不足。但因其干摩擦恢復(fù)力部分采用微分形式表達(dá)，導(dǎo)致該模型表達(dá)式復(fù)雜，求解困難。針對(duì)此問題本文旨在充分考慮干摩擦恢復(fù)力非線性的基礎(chǔ)上，利用微元結(jié)構(gòu)理論進(jìn)一步優(yōu)化有記憶型恢復(fù)力模型，并利用諧振動(dòng)試驗(yàn)驗(yàn)證。

1 有記憶型恢復(fù)力力學(xué)模型的優(yōu)化

1.1 金屬橡膠試塊恢復(fù)力分析

金屬橡膠有記憶型恢復(fù)力模型中，與整個(gè)變形有關(guān)的阻尼部分有著復(fù)雜的阻尼耗能機(jī)理，主要包括兩部分，一是部分金屬絲之間發(fā)生滑移和變形，原有形態(tài)及相對(duì)位置發(fā)生改變，并產(chǎn)生永久變形，該部分表現(xiàn)為黏性阻尼；另一部分為金屬絲接觸面在受力時(shí)發(fā)生相對(duì)滑動(dòng)，在摩擦力的作用下所引起的阻尼為干摩擦結(jié)構(gòu)阻尼[5]。以上兩種恢復(fù)力有著復(fù)雜的非線性特性，是恢復(fù)力模型中主要研究部分，文中著重對(duì)干摩擦恢復(fù)力部分進(jìn)行優(yōu)化研究。

非線性彈性恢復(fù)力采用5 次高階多項(xiàng)式表示，即

式中：Fk為非線性彈性恢復(fù)力；x為金屬橡膠加載位移；k1、k3、k5為彈性系數(shù)。

對(duì)于金屬橡膠非線性恢復(fù)力，傳統(tǒng)上通常采用兩種形式表達(dá)，即采用等效阻尼力整體表示金屬橡膠非線性阻尼力，另一種形式是將非線性阻尼力細(xì)化，分為黏性阻尼力和干摩擦阻尼力[6]。為進(jìn)一步研究干摩擦阻尼力產(chǎn)生機(jī)理，本文采用后者，用黏性阻尼力和干摩擦阻尼力之和表示非線性阻尼力。金屬橡膠黏性阻尼特性與振幅和振頻相關(guān)，其中一次黏性阻尼項(xiàng)占主要成分，高次速度項(xiàng)可忽略不計(jì)。因此黏性阻尼力部分可表示為

其中：C為黏性阻尼系數(shù)，A為振動(dòng)幅值。

1.2 干摩擦阻尼力

金屬橡膠試塊是金屬絲螺旋卷經(jīng)過拉伸、纏繞、沖壓而成型。在加卸載過程中，隨著運(yùn)動(dòng)空間的減少，金屬絲螺旋卷相互接觸并發(fā)生相對(duì)滑動(dòng)，產(chǎn)生干摩擦阻尼力。在加卸載過程中，金屬絲螺旋卷接觸狀態(tài)不斷變化，通過對(duì)不同接觸狀態(tài)金屬絲螺旋卷進(jìn)行力學(xué)分析，研究干摩擦阻尼力產(chǎn)生機(jī)理[7]。結(jié)合不同接觸狀態(tài)螺旋卷隨位移變化的規(guī)律，推導(dǎo)宏觀下金屬橡膠構(gòu)件的干摩擦阻尼力。

金屬橡膠基本結(jié)構(gòu)單元為金屬絲螺旋卷，如圖1所示。圖中LJ為螺旋卷螺距；DJ為螺旋卷直徑；ds為金屬絲直徑；α為螺旋升角。

圖1 螺旋卷結(jié)構(gòu)圖

在加卸載過程中，微元體相互擠壓并產(chǎn)生相對(duì)滑動(dòng)，引起內(nèi)摩擦耗能，是產(chǎn)生干摩擦阻尼力的主要形式。為此干摩擦阻尼力問題可轉(zhuǎn)化為螺旋卷滑移問題。由力學(xué)分析可知，在加載過程中，載荷首先克服微元體剛度下移，隨著微元體變形不斷增大，接觸點(diǎn)切向力增大，當(dāng)切向力大于臨界摩擦力時(shí)，螺旋卷產(chǎn)生相對(duì)滑動(dòng)；卸載過程中，隨著載荷減小，在結(jié)構(gòu)剛度的作用下，微元體恢復(fù)原有形狀。

如圖2所示，支撐力FN由下部橫向螺旋卷結(jié)構(gòu)剛度KT提供，KT根據(jù)材料和圓柱螺旋壓縮彈簧基本理論[8]來確定。

圖2 滑移接觸加載過程示意圖

E為金屬絲材料的彈性模量；υ為金屬絲材料的泊松比。

加載過程中，隨著金屬橡膠試塊體積不斷減小，金屬絲螺旋卷接觸狀態(tài)也隨之從未接觸狀態(tài)到滑移狀態(tài)轉(zhuǎn)化，隨體積進(jìn)一步減小，滑移狀態(tài)微元體轉(zhuǎn)化為擠壓狀態(tài)。微元體狀態(tài)數(shù)目比例與金屬橡膠試塊相對(duì)密度有著復(fù)雜的對(duì)應(yīng)關(guān)系，假設(shè)二者成二次非線性關(guān)系，得到如圖3所示的變化趨勢(shì)。

圖3 不同接觸狀態(tài)微元體數(shù)目比例圖

橫坐標(biāo)為金屬橡膠試塊相對(duì)密度，縱坐標(biāo)n為不同接觸狀態(tài)微元體數(shù)目比例。

圖中為隨位移加載量不斷變化的金屬橡膠試塊相對(duì)密度。其中曲線1為未接觸狀態(tài)微元體數(shù)目和擠壓狀態(tài)微元體數(shù)目比例，在之前為未接觸狀態(tài)微元體比例之后為擠壓狀態(tài)微元體比例。曲線2為隨加載位移而不斷變化的滑移狀態(tài)微元體數(shù)目的比例。兩曲線在相同相對(duì)密度下比例系數(shù)和為1。其中為采用拉伸編織工藝后及毛坯狀態(tài)的相對(duì)密度，毛坯狀態(tài)的微元體的接觸狀態(tài)為未接觸。根據(jù)毛坯結(jié)構(gòu)可得毛坯狀態(tài)密度

為金屬橡膠構(gòu)件在彈性范圍內(nèi)能達(dá)到的最大相對(duì)密度，通常情況下其值在0.3～0.4 變化。金屬橡膠構(gòu)件微元體總數(shù)目N可根據(jù)其成型參數(shù)確定

式中：Ls為金屬絲長(zhǎng)度；ρ為金屬橡膠試塊密度；V為金屬橡膠試塊體積；ρs為金屬絲密度。

根據(jù)相對(duì)密度變化確定3種接觸狀態(tài)微元體數(shù)目，由于干摩擦耗能只與滑移狀態(tài)微元體貢獻(xiàn)數(shù)目相關(guān)［9］，得滑移狀態(tài)的微軟體數(shù)目比例ns

結(jié)合以上公式可得到隨位移不斷變化的干摩擦力Ff

其中：μ為干摩擦系數(shù)。

結(jié)合式（1）、式（2）和式（7）得到金屬橡膠構(gòu)件的本構(gòu)關(guān)系

上文干摩擦阻尼力分析中只考慮了一組螺旋卷的相互作用，實(shí)際狀態(tài)下周圍螺旋卷對(duì)該組螺旋卷產(chǎn)生擠壓、摩擦等影響，隨著振幅的增大，影響程度隨之增加。因此在干摩擦阻尼力部分引入與振動(dòng)幅值A(chǔ)相關(guān)的修正系數(shù)B，得到以下數(shù)學(xué)模型

式中干摩擦阻尼修正系數(shù)B可由試驗(yàn)獲得。

2 金屬橡膠準(zhǔn)靜態(tài)試驗(yàn)研究

金屬橡膠恢復(fù)力在宏觀上表現(xiàn)為具有記憶特性的遲滯回線[10]。圖4中1為加載曲線，2為卸載曲線。加載時(shí)，黏性阻尼力和干摩擦阻尼力阻止內(nèi)部金屬絲的變形，表現(xiàn)為金屬橡膠塊剛度變大；卸載時(shí)，黏性阻尼力和干摩擦阻尼力方向改變，阻止變形恢復(fù)，變?yōu)榻饘傧鹉z實(shí)際剛度減小。

圖4 金屬橡膠遲滯回線

加載曲線和卸載曲線所圍面積代表黏性阻尼力和干摩擦阻尼力在單個(gè)振動(dòng)周期內(nèi)耗能，體現(xiàn)金屬橡膠減振能力。金屬橡膠遲滯曲線可由加卸載時(shí)力-位移曲線得到。

金屬橡膠試塊結(jié)構(gòu)參數(shù)中，相對(duì)密度對(duì)金屬橡膠遲滯特性有較大影響，為驗(yàn)證優(yōu)化后有記憶型恢復(fù)力模型準(zhǔn)確性及適用性，通過試驗(yàn)研究不同相對(duì)密度金屬橡膠試塊的遲滯回線。

2.1 金屬橡膠試塊的制備

采用直徑為0.2 mm、密度為7.93 kg/cm3的奧氏體不銹鋼303 金屬絲為原材料，利用金屬橡膠卷簧機(jī)床機(jī)制備金屬絲螺旋卷，并經(jīng)等距拉伸至恒定螺距2 mm，通過模具繞制與沖壓成型，制備得到高為20 mm、內(nèi)外徑分別為16 mm和36 mm的柱型構(gòu)件，如圖5所示。

圖5 金屬橡膠試塊

保持形狀因子β為0.32，制備出2種不同相對(duì)密度的試塊，詳細(xì)參數(shù)如表1所示。

表1 試塊基本參數(shù)

2.2 理論與試驗(yàn)對(duì)比分析

為驗(yàn)證理論模型可靠性和在靜態(tài)激勵(lì)下的適用性，采用準(zhǔn)靜態(tài)試驗(yàn)對(duì)模型進(jìn)行驗(yàn)證。試驗(yàn)采用CMT5105型全數(shù)字化金屬材料電子萬能試驗(yàn)機(jī)，通過可編程控制器實(shí)施加載和卸載，如圖6所示。

圖6 試驗(yàn)裝置及金屬橡膠試塊安裝

分別對(duì)2組金屬橡膠試塊進(jìn)行準(zhǔn)靜態(tài)加卸載試驗(yàn)，為保證初始載荷與位移的準(zhǔn)確性，試驗(yàn)開始前，對(duì)試塊預(yù)加2 mm的壓縮位移，試驗(yàn)中控制加卸載速率為5 mm/min，最大加載位移為6 mm。對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行處理得到如圖7所示結(jié)果。

為驗(yàn)證數(shù)學(xué)模型的可靠性，利用試驗(yàn)所得數(shù)據(jù)對(duì)式（9）進(jìn)行參數(shù)識(shí)別。優(yōu)化后的模型無復(fù)雜的微分形式，結(jié)構(gòu)較簡(jiǎn)單，可采用普通的非線性最小二乘法對(duì)式中參數(shù)進(jìn)行一次性識(shí)別[11]。

圖7 不同相對(duì)密度試塊加卸載試驗(yàn)

利用相對(duì)密度為0.2 和相對(duì)密度為0.25 的試塊進(jìn)行理論驗(yàn)證，得到的試驗(yàn)值和理論值如圖8所示。

可以看出在不同的相對(duì)密度下，該理論模型與試驗(yàn)值都能很好擬合，尤其當(dāng)相對(duì)密度較大時(shí)擬合精度更高。雖有誤差，但模型一定程度上可以準(zhǔn)確描述金屬橡膠試塊遲滯恢復(fù)力特性。誤差產(chǎn)生的原因?yàn)樵谟美碚撃Ｐ瓦M(jìn)行微元體結(jié)構(gòu)力學(xué)分析時(shí)，根據(jù)材料和結(jié)構(gòu)性形式加載剛度設(shè)定為定剛度，而在實(shí)際的金屬橡膠試塊內(nèi)部，螺旋卷剛度不僅隨結(jié)構(gòu)變形而改變，并且受周圍微元體的影響，剛度表現(xiàn)為較強(qiáng)的非線性。

3 數(shù)學(xué)模型物理意義辨析

數(shù)學(xué)模型中參數(shù)有明確的物理意義，對(duì)于指導(dǎo)工程應(yīng)用有著重要的作用。為辨析本文建立的金屬橡膠遲滯特性本構(gòu)模型中參數(shù)物理意義，選取相對(duì)密度為0.2金屬橡膠試塊進(jìn)行試驗(yàn)，參數(shù)識(shí)別中采用加載速度為5 mm/min、加載幅值分別為2 mm、4 mm、6 mm 和8 mm 的試驗(yàn)值。識(shí)別結(jié)果如表2所示。

表2 加載幅值為2 mm至8 mm的參數(shù)識(shí)別結(jié)果

從表2中識(shí)別結(jié)果可以看出，非線性彈性恢復(fù)力系數(shù)隨振幅改變呈現(xiàn)出明顯規(guī)律性，其中低階剛度系數(shù)相對(duì)較大，表明在彈性恢復(fù)力中低階剛度系數(shù)占主要成分。系數(shù)隨振幅的增大而減小，表明彈性恢復(fù)力在整個(gè)恢復(fù)力中貢獻(xiàn)比例隨振幅增大而變小。

該模型中，非線性阻尼力系數(shù)也具有明顯的物理意義。從識(shí)別結(jié)果中看出黏性阻尼力系數(shù)隨著振幅增加而變大，表明黏性阻尼耗能能力隨振幅增加而變大，減振特性也會(huì)越明顯。從表2參數(shù)識(shí)別結(jié)果中可以看出干摩擦阻尼力修正系數(shù)B與振幅正相關(guān)。如圖9所示。

圖8 理論值與試驗(yàn)值結(jié)果對(duì)比

圖9 干摩擦阻尼修正系數(shù)B和黏性阻尼系數(shù)C隨加載位移變化曲線

圖9中對(duì)比黏性阻尼系數(shù)C和干摩擦阻尼力修正系數(shù)B隨振幅變化速率，可看出干摩擦修正系數(shù)變化速率較快，表明隨著振幅的增加，干摩擦阻尼力貢獻(xiàn)比例也在不斷增加。

4 結(jié)語

（1）從微觀角度研究金屬橡膠遲滯恢復(fù)力特性，通過對(duì)金屬橡膠螺旋卷結(jié)構(gòu)力學(xué)分析及滑移狀態(tài)微元體數(shù)目的變化得出干摩擦阻尼力產(chǎn)生機(jī)理：隨著位移的增加，金屬橡膠螺旋卷接觸并發(fā)生滑移產(chǎn)生干摩擦，同時(shí)滑移接觸微元體數(shù)目也隨位移的增加而增多，增大了干摩擦阻尼力在非線性力中的比重。

（2）提出一種金屬橡膠遲滯恢復(fù)力模型，該模型考慮了干摩擦阻尼力對(duì)遲滯特性的影響，代替原有模型中微分表達(dá)形式。優(yōu)化后的模型中干摩擦阻尼力隨位移變化特性清晰且表達(dá)形式簡(jiǎn)練，具有求解簡(jiǎn)單、適用范圍廣的優(yōu)點(diǎn)，彌補(bǔ)了原有微分模型表達(dá)式復(fù)雜、物理意義不清晰和參數(shù)識(shí)別過程中求解困難的不足。

（3）提出的金屬橡膠遲滯恢復(fù)力模型將非線性阻尼力進(jìn)一步用黏性阻尼力和干摩擦阻尼力準(zhǔn)確表示。從參數(shù)識(shí)別結(jié)果可知，非線性彈性力部分1 階剛度占主要成分，且隨著位移的增大剛度減小，說明進(jìn)入到軟特性階段；非線性阻尼力中，振幅對(duì)黏性阻尼力和干摩擦阻尼力都有較大影響，其中干摩擦阻尼力隨振幅增加貢獻(xiàn)系數(shù)不斷增大。

本文雖優(yōu)化了模型表達(dá)式，但在描述干摩擦阻尼力時(shí)將滑移支撐剛度假設(shè)為定剛度，與金屬橡膠中復(fù)雜的非線性剛度存在差異，與原有微分模型相比降低了計(jì)算精度。因此在后續(xù)工作中亟待對(duì)干摩擦阻尼力分量中的剛度做進(jìn)一步研究，以提高模型精度。