填充橡膠材料動(dòng)態(tài)特性與激振頻率及激振振幅的關(guān)系

填充橡膠材料動(dòng)態(tài)特性與激振頻率及激振振幅的關(guān)系

0本構(gòu)模型的建立橡膠分離器廣泛應(yīng)用于航空航天、汽車(chē)、船舶等減振器。試驗(yàn)表明,橡膠隔振器的動(dòng)態(tài)特性具有較強(qiáng)的非線性特征,并且與激振頻率及激振振幅相關(guān)。因此,對(duì)橡膠隔振器的動(dòng)態(tài)性能的預(yù)測(cè)難度較大,一般需要反復(fù)地設(shè)計(jì)和試驗(yàn)才能獲得滿意的橡膠隔振器動(dòng)態(tài)性能。近年來(lái),已有不同的模型,如Maxwell模型、Kelvin-Voigt模型以及分?jǐn)?shù)導(dǎo)數(shù)模型來(lái)表征橡膠隔振器動(dòng)態(tài)特性的頻率相關(guān)性。與Kelvin-Voigt模型相比,Maxwell模型提高了描述橡膠隔振器動(dòng)特性的精度。BERG引入分?jǐn)?shù)導(dǎo)數(shù)概念,可用較少的模型參數(shù)來(lái)表征橡膠材料的頻率相關(guān)性,但在時(shí)域上計(jì)算的低效率限制了分?jǐn)?shù)導(dǎo)數(shù)模型的工程應(yīng)用。將若干Maxwell單元并聯(lián)形成廣義Maxwell模型,可建立橡膠材料的本構(gòu)模型,這種模型能較好地預(yù)測(cè)橡膠隔振器在高頻、小振幅激勵(lì)下的動(dòng)態(tài)性能,被商業(yè)有限元軟件廣泛采用。在橡膠材料中加入碳黑填充劑,使得橡膠材料表現(xiàn)出類(lèi)似于摩擦現(xiàn)象的特征,其動(dòng)態(tài)特性表現(xiàn)出較強(qiáng)地激振振幅相關(guān)性。BERG建立了描述橡膠隔振器動(dòng)態(tài)特性的模型,用庫(kù)侖摩擦力單元描述其振幅相關(guān)性。CHRISTIAN等及LIN等從粘彈性和彈塑性理論出發(fā),給出了描述橡膠材料激振頻率和激振振幅相關(guān)性的本構(gòu)關(guān)系模型,但在實(shí)際工程應(yīng)用時(shí)還有諸多不便和困難。OLSSON將若干個(gè)理想彈塑性單元并聯(lián),建立了填充橡膠材料的本構(gòu)模型,該本構(gòu)模型可較好地描述填充橡膠材料動(dòng)態(tài)特性的振幅相關(guān)性。文中利用已有的超彈性模型、廣義Maxwell模型和理想彈塑性模型,采用模型疊加方法,得到了表征填充橡膠材料動(dòng)態(tài)特性的激振頻率和激振振幅相關(guān)性的超彈性—粘彈性—彈塑性本構(gòu)關(guān)系模型,模型中的參數(shù)都具有明確的物理意義。利用橡膠試件的簡(jiǎn)單剪切試驗(yàn)數(shù)據(jù)擬合得到了模型參數(shù)。計(jì)算了某橡膠懸置的動(dòng)態(tài)特性,并與試驗(yàn)數(shù)據(jù)進(jìn)行了對(duì)比分析。文中建立的填充橡膠材料動(dòng)態(tài)特性計(jì)算的本構(gòu)模型、模型參數(shù)的識(shí)別方法和計(jì)算方法,可用于橡膠隔振器動(dòng)態(tài)特性的計(jì)算分析。1橡膠材料彈性、粘彈性和彈塑性單元的非線性流變模型填充橡膠材料的動(dòng)態(tài)特性與激振頻率和激振振幅相關(guān)。對(duì)于填充碳黑橡膠材料,一般可假設(shè)其彈性特性、動(dòng)態(tài)特性的激振頻率相關(guān)性和激振振幅相關(guān)性三者相互獨(dú)立。因此可以建立如圖1所示的彈性、粘彈性和彈塑性單元并聯(lián)而形成的一維非線性流變模型。其中,彈性模型、粘彈性模型和彈塑性模型分別用于描述橡膠材料的彈性特性、動(dòng)態(tài)特性的頻率相關(guān)性和振幅相關(guān)性。圖1中,τ為應(yīng)力,ε為應(yīng)變,G∞為彈性單元的準(zhǔn)靜態(tài)彈性模量,Give和tr_i(i=1,2,,N)分別為粘彈性單元中的彈性模量和松弛時(shí)間,Gjep和κy_j(j=1,2,,M)分別為彈塑性單元中的彈性模量和屈服應(yīng)變。由圖1,模型總的應(yīng)力可以表示為式中τe——彈性單元應(yīng)力τep——M個(gè)彈塑性單元疊加后的總應(yīng)力τve——N個(gè)粘彈性單元疊加后的總應(yīng)力2試驗(yàn)結(jié)果和表征圖2、3分別為橡膠試片的簡(jiǎn)單剪切試驗(yàn)裝置圖和試片示意圖,剪切試片長(zhǎng)度a=25mm,寬度b=20mm,厚度δ0=5.8mm。橡膠試片的靜、動(dòng)態(tài)特性在MTS831液壓伺服振動(dòng)試驗(yàn)臺(tái)上測(cè)試。測(cè)試時(shí),在作動(dòng)端施加不同振幅和頻率的正弦位移激勵(lì),記錄固定端力傳感器的反力。橡膠試片受到簡(jiǎn)諧位移激勵(lì)時(shí),填充橡膠材料的阻尼引起滯后效應(yīng),其力—位移的關(guān)系形成遲滯環(huán),將力—位移關(guān)系轉(zhuǎn)換為切應(yīng)力—切應(yīng)變關(guān)系,如圖4所示。橡膠材料的動(dòng)態(tài)特性由動(dòng)態(tài)切變模量和阻尼表征。動(dòng)態(tài)切變模量和阻尼的定義為式中τ0——切應(yīng)力幅值ε0——切應(yīng)變幅值Uc——滯回環(huán)的損耗能3動(dòng)態(tài)切變模量和阻尼在時(shí)域中采用時(shí)間—步長(zhǎng)方法計(jì)算彈性單元、粘彈性單元和彈塑性單元的動(dòng)態(tài)切變模量和阻尼。將激振周期T進(jìn)行n等分,每段時(shí)間為?t=T/n,初始時(shí)刻時(shí),整個(gè)模型的應(yīng)力和應(yīng)變均為零值。以下計(jì)算彈性、粘彈性及彈塑性單元在各時(shí)間段?lt(l=1,2,,n)內(nèi)的應(yīng)力增量。3.1彈性單元在時(shí)間段?lt內(nèi),彈性單元的應(yīng)力增量式中,?εl為?lt內(nèi)切應(yīng)變的增量3.2共粘彈性單元在時(shí)間段?lt內(nèi),第i個(gè)Maxwell粘彈性單元的應(yīng)力為由式(4)可知,在時(shí)間段?lt內(nèi),第i個(gè)粘彈性單元的應(yīng)力增量為3.3模型參數(shù)的確定第j個(gè)彈塑性單元典型的應(yīng)力—應(yīng)變遲滯環(huán)如圖5所示。設(shè)第j個(gè)彈塑性單元的在第l個(gè)時(shí)刻的應(yīng)力為則其在第l+1個(gè)時(shí)刻的應(yīng)力如式(6)所示,當(dāng)γ>1時(shí),令γ=1。時(shí)間段?lt內(nèi),第j個(gè)彈塑性單元的切應(yīng)力增量在時(shí)間段?lt內(nèi),整個(gè)模型的切應(yīng)力增量為所有單元切應(yīng)力增量的總和。這樣,周期T中各時(shí)間點(diǎn)的切應(yīng)力和切應(yīng)變均已求出,可以得到如圖4所示的切應(yīng)力—切應(yīng)變遲滯環(huán)。由式(2)可以計(jì)算得到模型的動(dòng)態(tài)切變模量和阻尼。模型參數(shù)可以通過(guò)最小化動(dòng)態(tài)切變模量和阻尼的計(jì)算值和試驗(yàn)值的相對(duì)誤差來(lái)得到式中,dd,i、Gd,i為計(jì)算值,de,i、Ge,i為測(cè)試值,m為測(cè)量點(diǎn)總數(shù)。為縮小優(yōu)化范圍,假設(shè)tr_1>tr_2>tr_3>0,max(ε0)>εy_1>εy_2>εy_3>0。試驗(yàn)值和計(jì)算值的擬合效果見(jiàn)圖6,得到的模型參數(shù)見(jiàn)下表。由圖6可見(jiàn),采用三個(gè)粘彈性單元和三個(gè)彈塑性單元,可較好地?cái)M合橡膠試片的動(dòng)態(tài)特性。4本構(gòu)模型三維建模的實(shí)現(xiàn)在進(jìn)行橡膠隔振器動(dòng)態(tài)特性的有限元分析時(shí),不必構(gòu)造新的橡膠材料本構(gòu)模型,只需根據(jù)已識(shí)別出的一維超彈性、粘彈性和彈塑性本構(gòu)模型參數(shù),進(jìn)行簡(jiǎn)單地有限元模型疊加。本構(gòu)模型疊加的有限元實(shí)現(xiàn)的具體步驟如下所述。首先,分別建立對(duì)應(yīng)于超彈性、粘彈性和彈塑性本構(gòu)模型的有限元分析模型,三種模型具有相同的有限元網(wǎng)格劃分形式,如圖7所示。然后,將三種有限元模型的結(jié)點(diǎn)進(jìn)行“粘結(jié)”,疊加三種有限元模型的網(wǎng)格,如圖8所示。網(wǎng)格疊加的結(jié)果是,結(jié)點(diǎn)數(shù)沒(méi)有增加,但網(wǎng)格單元數(shù)為原單元數(shù)的3倍。最終,形成了類(lèi)似于圖1中一維流變模型的三維彈性—粘彈性—彈塑性有限元分析模型。需要注意的是,在結(jié)構(gòu)載荷的作用下,被“粘結(jié)”的各結(jié)點(diǎn)不會(huì)產(chǎn)生分離。5橡膠懸置動(dòng)態(tài)特性的頻率相關(guān)使用擬合得到的本構(gòu)模型參數(shù),由有限元分析方法計(jì)算一橡膠懸置的動(dòng)態(tài)特性,橡膠懸置的有限元分析模型如圖9所示。約束面為金屬外管,分析時(shí)對(duì)其進(jìn)行位移固定約束,即金屬外管上結(jié)點(diǎn)的各方向位移為零。加載面為金屬內(nèi)管,用剛性連接單元關(guān)聯(lián)到中間某一結(jié)點(diǎn),以該結(jié)點(diǎn)作為加載點(diǎn)。在加載點(diǎn)施加沿z方向的簡(jiǎn)諧位移激勵(lì),經(jīng)有限元分析可獲得該點(diǎn)的反力,從而可以計(jì)算該懸置在z方向的動(dòng)剛度和阻尼。橡膠懸置動(dòng)態(tài)特性的頻率相關(guān)性如圖10所示。由圖10a,在某一激振振幅下,動(dòng)剛度隨著頻率的增大而增大,動(dòng)剛度的計(jì)算值與試驗(yàn)值吻合較好,經(jīng)計(jì)算可知相對(duì)誤差小于8%。由圖10b,滯后角隨著頻率的增大而增大,滯后角的計(jì)算值與試驗(yàn)值也比較吻合,基本滿足工程精度要求。在記錄試驗(yàn)數(shù)據(jù)時(shí),不可避免的會(huì)存在偏差,這是造成圖10b中試驗(yàn)數(shù)據(jù)有較大波動(dòng)的原因。橡膠懸置動(dòng)態(tài)特性的振幅相關(guān)性如圖11所示。由圖11a,在某一激振頻率下,動(dòng)剛度隨著激振振幅的增大而減小(即Payne效應(yīng)),動(dòng)剛度的計(jì)算值與試驗(yàn)值吻合較好,相對(duì)誤差也小于8%。由圖11b,滯后角隨著振幅的增大而增大,滯后角的計(jì)算值與試驗(yàn)值吻合程度一般,但基本能揭示滯后角隨激振振幅的變化趨勢(shì)。綜上所述,采用文中建立的超彈性-粘彈性-彈塑性本構(gòu)模型可以較好地表征橡膠材料動(dòng)態(tài)特性的振幅相關(guān)性和頻率相關(guān)性。雖然計(jì)算值與試驗(yàn)值存在一定程度的誤差,但對(duì)于橡膠材料,小于8%的計(jì)算誤差是可以被工程實(shí)際所接受的。6本構(gòu)模型的建立與分析針對(duì)橡膠隔振器動(dòng)態(tài)性能預(yù)測(cè)難度較大的問(wèn)題,將已有的橡膠材料本構(gòu)模型進(jìn)行簡(jiǎn)單地疊加,得到了填充橡膠材料的超彈性—粘彈性—彈塑性本構(gòu)模型。利用擬合得到的模型參數(shù)對(duì)某橡膠懸置的動(dòng)態(tài)特性進(jìn)行了計(jì)算,并與試驗(yàn)數(shù)據(jù)進(jìn)行了對(duì)比分析。(1)文中給出的超彈性—粘彈性—彈塑性本構(gòu)模型能較準(zhǔn)確地表征橡膠材料的動(dòng)態(tài)特性與激振振幅和激振頻率的相關(guān)性。該方法可以較大程度地提高對(duì)橡膠隔振器動(dòng)態(tài)特性的預(yù)測(cè)精度。(2)文中方法只須在已