復(fù)雜非金屬減震結(jié)構(gòu)的建模與等效研究

宋 勇，靖建全

(1. 中國船舶重工集團(tuán)公司第七一三研究所，河南 鄭州 450015；2. 北京理工大學(xué) 宇航學(xué)院，北京 100081)

復(fù)雜非金屬減震結(jié)構(gòu)的建模與等效研究

宋 勇1，靖建全2

(1. 中國船舶重工集團(tuán)公司第七一三研究所，河南 鄭州 450015；2. 北京理工大學(xué) 宇航學(xué)院，北京 100081)

針對(duì)減震系統(tǒng)中復(fù)雜非金屬減震結(jié)構(gòu)數(shù)值分析和等效處理的需求，利用有限元方法對(duì)預(yù)彎柱體非金屬減震器進(jìn)行建模和分析，獲得了與實(shí)驗(yàn)測試吻合良好的計(jì)算結(jié)果。在此基礎(chǔ)上，提出一種利用實(shí)物實(shí)驗(yàn)或數(shù)值實(shí)驗(yàn)數(shù)據(jù)等效復(fù)雜非金屬減震結(jié)構(gòu)的等效模型處理方法。該方法利用三維實(shí)體結(jié)構(gòu)等效復(fù)雜非金屬減震器外形和連接關(guān)系，以實(shí)際減震器的單軸壓縮實(shí)驗(yàn)和應(yīng)力松弛實(shí)驗(yàn)數(shù)據(jù)獲得等效模型的超彈和粘彈材料參數(shù)。實(shí)例分析表明，基于三維實(shí)體結(jié)構(gòu)的等效模型能夠有效表征實(shí)際減震器在復(fù)雜減震系統(tǒng)中的連接關(guān)系，并在主工作方向上有效模擬實(shí)際減震結(jié)構(gòu)的力學(xué)響應(yīng)特性。

減震器；超彈性；粘彈性；等效模型

0 引 言

非金屬材料如橡膠、聚氨酯發(fā)泡材料等制備的減震器由于具有獨(dú)特的材料屬性，以及易加工成型的工藝特點(diǎn)使得其在導(dǎo)彈發(fā)射裝置減震系統(tǒng)以及其他減震環(huán)節(jié)中得到了廣泛應(yīng)用。為滿足復(fù)雜振動(dòng)環(huán)境下的減震要求，這類減震器往往具有復(fù)雜的結(jié)構(gòu)特征，其力學(xué)性能的研究、模擬和等效一直是國內(nèi)外研究人員關(guān)注的重要內(nèi)容[1 – 4]。

近年來，張強(qiáng)等[5]為模擬某導(dǎo)彈發(fā)射裝置減震系統(tǒng)的減震性能，將減震結(jié)構(gòu)簡化為非線性彈性單元，研究了彈體遭受沖擊后的加速度響應(yīng)情況，結(jié)果表明由外筒、筒間減震塊、內(nèi)筒、適配器組成的兩級(jí)減震系統(tǒng)具有良好的減震效果，能大幅衰減傳遞到彈體上的沖擊。王平等[6]為提高機(jī)載光電吊艙的成像質(zhì)量和穩(wěn)定精度，設(shè)計(jì)了2種能限制一定自由度的小型金屬橡膠減振器并將其分別嵌入光電吊艙內(nèi)、外框架減振系統(tǒng)中，并對(duì)其減震性能進(jìn)行了實(shí)驗(yàn)和計(jì)算分析仿真結(jié)果表明光電吊艙兩級(jí)減振系統(tǒng)對(duì)25Hz以上的角振動(dòng)有很好的減振效果；Mattias[7]在考慮橡膠減震器的頻率與振幅相關(guān)特性的基礎(chǔ)上，提出了橡膠減震器的一維等效模型，并對(duì)其等效效果進(jìn)行了試驗(yàn)驗(yàn)證；和法家等[8]以研究某飛機(jī)發(fā)動(dòng)機(jī)橡膠減振器為研究對(duì)象，對(duì)其進(jìn)行了一系列實(shí)驗(yàn)，得到了振動(dòng)頻率、激勵(lì)振幅、預(yù)載荷與橡膠減震器動(dòng)剛度的主要關(guān)系。可以看出，在復(fù)雜減震系統(tǒng)中，有效表征和模擬減震器結(jié)構(gòu)的力學(xué)性能、并采用合適的方式進(jìn)行等效和簡化，是進(jìn)行復(fù)雜減震系統(tǒng)分析的關(guān)鍵內(nèi)容之一。

本文針對(duì)復(fù)雜減震系統(tǒng)中的非金屬減震器結(jié)構(gòu)，利用有限元方法建立其計(jì)算分析模型，獲得了與實(shí)驗(yàn)測試吻合良好的計(jì)算結(jié)果。在此基礎(chǔ)上，提出利用三維實(shí)體結(jié)構(gòu)建立等效模型，并通過實(shí)物實(shí)驗(yàn)或數(shù)值實(shí)驗(yàn)方法獲得等效材料參數(shù)的復(fù)雜非金屬減震結(jié)構(gòu)簡化建模方法，為包含大量非金屬減震結(jié)構(gòu)的復(fù)雜減震系統(tǒng)三維建模和分析提供基礎(chǔ)。實(shí)例分析表明，文中提出的等效處理方法模型簡單、參數(shù)明確，在減震器的工作方向上與實(shí)際結(jié)構(gòu)的力學(xué)性能吻合良好。

1 復(fù)雜非金屬減震器有限元建模分析

1.1 有限元模型

復(fù)雜非金屬減震器的結(jié)構(gòu)形式多種多樣，本文以預(yù)彎柱體減震器為對(duì)象，對(duì)其有限元建模分析過程和等效處理方法進(jìn)行介紹。為便于對(duì)照，選取與文獻(xiàn)[9]中相同的預(yù)彎柱體減振器結(jié)構(gòu)，如圖1（a）所示。減震器的外形尺寸為L=105 mm，H=52 mm，W=105 mm，h=32 mm，b=7.5 mm，β=73°，材料為澆注型聚氨酯彈性體。

根據(jù)預(yù)彎柱體減震器的結(jié)構(gòu)尺寸建立如圖1（b）所示的三維有限元模型，在模型上表面施加0～20 mm的位移載荷，在下表面施加固定約束，在圖示B處施加對(duì)稱面約束，在減震器壓縮的過程中，1，2，3，4，5，6處會(huì)出現(xiàn)接觸，因此在這些位置施加對(duì)稱接觸進(jìn)行處理。接觸類型為點(diǎn)面接觸，接觸過程中不考慮摩擦。考慮到施加的變形尺寸相對(duì)原始尺寸較大，在有限元求解計(jì)算中采用了大變形模型。

針對(duì)聚氨酯材料的超彈特性，這里選用不可壓縮的九參數(shù)Mooney-Rivlin本構(gòu)模型對(duì)材料屬性進(jìn)行表征，其具體形式為：

式中：I1和I2分別為第一和第二Cauchy-Green應(yīng)變不變量；cij為與材料相關(guān)的參數(shù)，其中c10，c01，c20，c30分別為為2 200 000，300 000，–700 000，450 000，其余為0，單位為Pa。

1.2 有限元分析結(jié)果

利用商用有限元軟件Ansys對(duì)上述預(yù)彎柱體減震器有限元模型進(jìn)行求解計(jì)算，獲得隨加載位移變化的減震器響應(yīng)狀態(tài)。圖2給出了加載位移分別為5 mm，20 mm時(shí)的結(jié)構(gòu)變形和應(yīng)力分布。可以看出，在位移載荷為5 mm的情況下，減震器的響應(yīng)主要為中間柱體的彎曲變形，并在1，2，3，4位置處出現(xiàn)了局部接觸狀態(tài)，而5，6位置處未出現(xiàn)接觸，此時(shí)減震器的剛度特性主要與減震器中間柱體的抗彎能力有關(guān)，應(yīng)力最大值出現(xiàn)在接觸的位置。在位移載荷為20 mm的情況下，中間柱體產(chǎn)生了較大的彎曲變形，1，2，3，4，5，6處均出現(xiàn)了大面積接觸，且應(yīng)力最大值出現(xiàn)在中間柱體的彎曲拉伸部分，此時(shí)減震器的剛度特性不僅與聚氨酯材料的特性有關(guān)，還與減震器中間柱體的變形以及大面積接觸有關(guān)。

圖3給出了計(jì)算獲得的減震器在位移載荷作用下的變形-受力曲線。圖中F為有限元模型下表面節(jié)點(diǎn)支反力之和，DISP為減震器在高度方向上的變形量，在很大程度下，該曲線反映了減震器在載荷作用下的剛度特性。從圖中可看出，在壓縮量較小時(shí)（小于8 mm），減震器承受的作用力隨壓縮量增加而快速增加；壓縮量在8～18 mm這一范圍內(nèi)，減震器承受的作用力隨壓縮量的變化較為緩慢。如前所述，這一階段減震器中間柱體發(fā)生彎曲變形，其整體剛度特性降低。在壓縮量超過18 mm后，減震器中間柱體彎曲后產(chǎn)生大面積接觸，使其剛度特性快速增加。從這一過程可以看出，對(duì)于復(fù)雜非金屬減震器，其整體剛度特性不但與材料屬性有關(guān)，也與其內(nèi)部結(jié)構(gòu)以及工作狀態(tài)有關(guān)，呈現(xiàn)出復(fù)雜的非線性狀態(tài)。在包含多個(gè)以至幾十個(gè)類似復(fù)雜減震器的減震系統(tǒng)中，采用完整的減震器結(jié)構(gòu)進(jìn)行減震性能分析在計(jì)算量上存在較大的困難，為此尋求合理的減震器等效模型，并準(zhǔn)確反映其在工作方向上的剛度特性，是進(jìn)行復(fù)雜系統(tǒng)分析的關(guān)鍵。

從圖3中還可以看出，在整個(gè)壓縮過程中，計(jì)算結(jié)果與文獻(xiàn)給出的實(shí)驗(yàn)數(shù)據(jù)吻合良好。在壓縮量小于18 mm的范圍內(nèi)，結(jié)果誤差小于5%；在壓縮量大于18 mm后，受減震器內(nèi)部結(jié)構(gòu)產(chǎn)生大面積接觸影響，剛度特性快速增加，計(jì)算與實(shí)驗(yàn)結(jié)果間的誤差也有所增加，但變化趨勢(shì)相同。由于數(shù)值計(jì)算與實(shí)物試驗(yàn)具有較好的一致性，因此在本文的研究中，主要采用數(shù)值試驗(yàn)的方法獲取結(jié)構(gòu)的響應(yīng)參數(shù)，并以此進(jìn)行等效過程的實(shí)例分析。

2 復(fù)雜非金屬減震器等效方法

2.1 等效原理

在多數(shù)情況下，尤其在本文涉及的導(dǎo)彈減震系統(tǒng)中，減震器工作狀態(tài)通常為壓縮變形和吸能，因此在減震器工作時(shí)，可將其靜態(tài)剛度特性考慮為具有超彈本構(gòu)特性的實(shí)心結(jié)構(gòu)在單軸壓縮時(shí)的響應(yīng)過程。對(duì)于不可壓縮的超彈性材料在單軸壓縮下的主應(yīng)力與應(yīng)變能函數(shù)的存在如下關(guān)系[10]：

式中：pe為靜水壓力；F為應(yīng)變梯度張量； B為左Cauchy-Green變形張量；W為應(yīng)變能函數(shù)；I1，I2為B的第一，第二基本不變量；σ11，σ22，σ33為真實(shí)應(yīng)力；λ為單軸壓縮方向的拉伸率。

將σ11-λ的關(guān)系轉(zhuǎn)化為工程應(yīng)力應(yīng)變的關(guān)系式為：

可以看出，σ11-λ的關(guān)系、ε的關(guān)系可由材料的應(yīng)變能函數(shù)W確定，通過合理設(shè)置應(yīng)變能函數(shù)的系數(shù)便可得到需要的應(yīng)力應(yīng)變關(guān)系。在等效過程中，只要選擇合適的應(yīng)變能函數(shù)W并通過由復(fù)雜非金屬減震結(jié)構(gòu)的單軸壓縮實(shí)物實(shí)驗(yàn)或數(shù)值實(shí)驗(yàn)獲得的工程應(yīng)力應(yīng)變關(guān)系確定其相應(yīng)的系數(shù)，就能得到與實(shí)際結(jié)構(gòu)相一致的力-位移關(guān)系，從而實(shí)現(xiàn)實(shí)心結(jié)構(gòu)對(duì)復(fù)雜非金屬減震器靜態(tài)剛度特性的等效。

復(fù)雜非金屬減震器除具有靜態(tài)剛度非線性外，通常還具有與應(yīng)變率相關(guān)的剛度特性，通常通過prony級(jí)數(shù)[11]定義材料的粘彈特性，無量綱后的應(yīng)力松弛模量可表示為：

通過復(fù)雜非金屬減震器的應(yīng)力松弛實(shí)驗(yàn)，可獲得時(shí)變工程應(yīng)力數(shù)據(jù)。在等效模型中，通過參數(shù)擬合，可確定其粘彈性材料的待定參數(shù)ei和li值。

2.2 等效方法

在包含多個(gè)復(fù)雜結(jié)構(gòu)減震器的減震系統(tǒng)中，在考慮減震性能時(shí)主要考察其整體在沖擊響應(yīng)過程中的力-位移曲線，而不關(guān)心復(fù)雜結(jié)構(gòu)減震器內(nèi)部所發(fā)生的變化。因此在對(duì)復(fù)雜結(jié)構(gòu)減震器進(jìn)行等效時(shí)，主要保證等效的減震器在工作方向上的響應(yīng)特性與復(fù)雜結(jié)構(gòu)減震器整體在沖擊響應(yīng)過程中的力-位移曲線相一致即可。在等效過程中，采用如下假設(shè)和方法：

1）等效模型為三維實(shí)心均勻材料結(jié)構(gòu)，材料各向同性；

2）利用超彈性材料參數(shù)表征由實(shí)際結(jié)構(gòu)材料以及結(jié)構(gòu)變形、接觸等因素引起的位移-載荷非線性變化過程；

3）利用prony級(jí)數(shù)形式的粘彈性材料參數(shù)表征實(shí)際結(jié)構(gòu)應(yīng)變率相關(guān)的剛度特性。

具體模型等效的思路和處理過程如圖4所示。

為了使等效模型最大程度地真實(shí)反映減震器在工作過程中所起的作用，等效模型的外形尺寸應(yīng)與減震器盡可能一致。在應(yīng)用過程中，可取等效模型在工作方向上的高度與減震器相同；在支撐平面上取面積與實(shí)際減震器支撐平面面積相同，并參考實(shí)際減震器外廓形式設(shè)置支撐平面形狀。

在通過單軸壓縮實(shí)驗(yàn)獲得實(shí)際減震器的力-位移響應(yīng)后，可通過如下公式直接計(jì)算其工程應(yīng)力和工程應(yīng)變：

式中：σc為復(fù)雜結(jié)構(gòu)減震器壓縮過程中的虛擬工程應(yīng)力；N為復(fù)雜結(jié)構(gòu)減震器的下表面支反力的和；A為等效實(shí)心結(jié)構(gòu)支撐面的面積；εc為復(fù)雜結(jié)構(gòu)減震器壓縮過程中的虛擬工程應(yīng)變；?L為復(fù)雜結(jié)構(gòu)減震器的壓縮位移；L為復(fù)雜結(jié)構(gòu)減震器的高度。

在實(shí)際應(yīng)用中，為通過實(shí)驗(yàn)數(shù)據(jù)獲得超彈性材料的本構(gòu)模型，通常需要將單軸壓縮數(shù)據(jù)轉(zhuǎn)變?yōu)榈入p軸拉伸數(shù)據(jù)，二者在理論上存在如下關(guān)系[12]：

式中：εb為等雙軸拉伸工程應(yīng)變；εc單軸壓縮工程應(yīng)變；σb為等雙軸拉伸工程應(yīng)力；σc為單軸壓縮工程應(yīng)力。

利用等雙軸拉伸實(shí)驗(yàn)數(shù)據(jù)，可通過參數(shù)擬合方式獲得常見超彈性材料本構(gòu)模型如多項(xiàng)式模型、Mooney-Rivlin模型的待定參數(shù)，即獲得等效模型的超彈性材料本構(gòu)參數(shù)。對(duì)于反應(yīng)率相關(guān)屬性的粘彈性模型，同樣可通過由實(shí)驗(yàn)獲得的時(shí)變應(yīng)力數(shù)據(jù)擬合出prony級(jí)數(shù)表征的粘彈性材料參數(shù)。

3 減震器的等效實(shí)例與分析

3.1 非線性彈性等效實(shí)例

在有限元建模分析中參考的減震器結(jié)構(gòu)在工作過程中主要表現(xiàn)為非線性彈性特征。為建立其等效模型，選擇其等效模型與實(shí)際減震器具有同樣的外廓尺寸，如圖5所示；設(shè)置反映非線性彈性屬性的本構(gòu)模型為5參數(shù)的Mooney-Rivlin本構(gòu)模型。

利用圖3所示的數(shù)值實(shí)驗(yàn)數(shù)據(jù)，通過工程應(yīng)力和工程應(yīng)變計(jì)算、單軸壓縮和雙軸拉伸應(yīng)力應(yīng)變轉(zhuǎn)換以及參數(shù)擬合，可獲得5參數(shù)Mooney-Rivlin超彈性本構(gòu)模型的系數(shù)如表1所示。

為校驗(yàn)等效模型，同樣利用有限元方法建立其計(jì)算模型，在其下表面施加固定約束，在其上表面施加0～20 mm位移載荷，獲得力-位移曲線如圖6所示。從圖中可看出，等效模型在壓縮過程中的力-位移曲線與實(shí)際減震器在壓縮過程中的力-位移曲線吻合良好。實(shí)際減震器在壓縮過程中的力-位移曲線的平臺(tái)效應(yīng)主要是由其內(nèi)部的復(fù)雜結(jié)構(gòu)引起；等效模型沒有復(fù)雜內(nèi)部結(jié)構(gòu)，但通過等效后的超彈性材料屬性，使其力-位移曲線仍然有平臺(tái)效應(yīng)，且與復(fù)雜結(jié)構(gòu)減震器的力-位移曲線相一致，說明等效實(shí)心結(jié)構(gòu)與復(fù)雜結(jié)構(gòu)減震器具有相同的靜態(tài)剛度特性。

3.2 率相關(guān)粘彈屬性的等效實(shí)例

以前述復(fù)雜非金屬減震器為例，為考慮相關(guān)粘彈特性及其等效過程，假定該減震器材料具有粘彈屬性，并可通過2參數(shù)prony模型進(jìn)行表示，模型參數(shù)[13]設(shè)置為a1=0.333 33，t1=0.4，a2=0.333 33，t2=0.2。同樣選擇其等效模型與實(shí)際減震器具有相同的外廓尺寸，并利用2參數(shù)prony模型表示其率相關(guān)的等效粘彈屬性。通過對(duì)復(fù)雜非金屬減震器施加20 mm的常位移載荷進(jìn)行數(shù)值實(shí)驗(yàn)并進(jìn)行轉(zhuǎn)換，獲得工程應(yīng)力隨時(shí)間變化曲線如圖7所示。利用參數(shù)擬合，獲得prony模型參數(shù)與假定的減震器粘彈性材料參數(shù)相同，表明該減震器在結(jié)構(gòu)上的阻尼特性并不顯著，其率相關(guān)特性主要體現(xiàn)為材料的粘彈特性。

為校驗(yàn)同時(shí)具有非線性彈性和率相關(guān)粘彈特性的減震器等效模型，利用有限元計(jì)算模型，在等效模型下表面施加位移約束，上表面施加的正弦位移載荷為：

表1 等效材料參數(shù)表Tab.1 Material parameters for the equivalent

經(jīng)過4個(gè)載荷周期計(jì)算，獲得等效模型計(jì)算穩(wěn)定后的力-位移曲線如圖8所示。可以看出在加、卸載過程中，力-位移曲線經(jīng)歷了不同的過程，呈現(xiàn)出典型的回滯環(huán)結(jié)構(gòu)，反映了結(jié)構(gòu)的率相關(guān)粘彈特征；在每次獨(dú)立的加載和卸載過程中，力-位移曲線均出現(xiàn)具有平臺(tái)效應(yīng)的非線性彈性狀態(tài)。值得注意的是，在最大位移20 mm處，由于粘彈性特性的應(yīng)力松弛效應(yīng)，同時(shí)具有超彈性與粘彈性特性的等效模型的支反力要比只有超彈性材料屬性等效結(jié)構(gòu)的支反力小；在位移為0時(shí)，由于粘彈性特性的阻尼效應(yīng)，同時(shí)具有超彈性與粘彈性性的等效模型支反力為負(fù)值。

4 結(jié) 語

本文針對(duì)復(fù)雜非金屬減震器結(jié)構(gòu)，利用有限元模型進(jìn)行了計(jì)算分析，獲得了與實(shí)驗(yàn)結(jié)果吻合良好的計(jì)算數(shù)據(jù)，表明基于有限元模型的數(shù)值實(shí)驗(yàn)?zāi)軌蜉^好的反映非金屬減震器的工作過程。

針對(duì)復(fù)雜大型減震系統(tǒng)簡化建模和分析的需求，提出一種利用材料屬性等效復(fù)雜非金屬減震器非線性響應(yīng)特性的等效建模方法。實(shí)例分析表明，等效模型在減震器的主要工作方向上能夠有效反映實(shí)際減震器結(jié)構(gòu)的非線性彈性和率相關(guān)粘彈特性，并能施加減震結(jié)構(gòu)的三維連接和約束關(guān)系，能夠?yàn)閺?fù)雜減震系統(tǒng)的三維建模和分析提供借鑒和參考。

參考文獻(xiàn)：

[ 1 ]BERG M. A model for rubber springs in the dynamic analysis of rail vehicles[J]. Proceedings of the Institution of Mechanical Engin, 1997, 211(2): 95–108.

[ 2 ]隋允康, 彭細(xì)榮, 杜家政. 橡膠緩振墊力學(xué)性能尺寸效應(yīng)分析[J]. 北京工業(yè)大學(xué)學(xué)報(bào), 2004, 30(4): 398–400.

SUI Yun-kang, PENG Xi-rong, DU Jia-zheng. Analysis of size effect on mechanics behavious of the rubber buffer for absorbing vibration[J]. Journal of Beijing University of Technology, 2004, 30(4): 398–400.

[ 3 ]王基, 朱石堅(jiān). V形橡膠減振器靜剛度研究[J]. 噪聲與振動(dòng)控制, 2003, 23(1): 37–40.

WANG Ji, ZHU Shi-jian. Study on the static stiffness of shape of V rubber Vibration - lsolator[J]. Noise and Vibration Control, 2003, 23(1): 37–40.

[ 4 ]周振凱, 徐兵, 胡文軍, 等. 橡膠隔振器大變形有限元分析[J].振動(dòng)與沖擊, 2013, 32(5): 171–175.

ZHOU Zhen-kai, XU Bing, HU Wen-jun, et al. Large deformation finite element analysis of rubber isolator[J]. Journal of Vibration and Shock, 2013, 32(5): 171–175.

[ 5 ]張強(qiáng), 吳明晨. 某型發(fā)射裝置水平減震系統(tǒng)仿真研究[J]. 艦船電子工程, 2012, 32(2): 73–75.

ZHANG Qiang, WU Ming-chen. Emulation research on lateral shock absorber system of certain launching system[J]. Ship Electronic Engineering, 2012, 32(2): 73–75.

[ 6 ]王平, 張國玉, 高玉軍, 等. 金屬橡膠減振器在機(jī)載光電吊艙復(fù)合減振系統(tǒng)中的應(yīng)用[J]. 振動(dòng)與沖擊, 2014, 33(5): 193–199.

WANG Ping, ZHANG Guo-yu, GAO Yu-jun, et al. Application of metal-rubber dampers in vibration reduction system of an airborne electro-optical pod[J]. Journal of Vibration and Shock, 2014, 33(5): 193–199.

[ 7 ]SJ?BERG M M, KARI L. Non-linear behavior of a rubber isolator system using fractional derivatives[J]. Vehicle System Dynamics, 2002, 37(3): 217–236.

[ 8 ]和法家, 盧曦. 某直升機(jī)橡膠減震器動(dòng)剛度特性試驗(yàn)[J]. 實(shí)驗(yàn)室研究與探索, 2014, 33(3): 41–44.

HE Fa-jia, LU Xi. Experimental research on dynamic stiffness characteristics of a Rubber-Shock-Absorber[J]. Research and Exploration in Laboratory, 2014, 33(3): 41–44.

[ 9 ]趙華, 王敏杰, 趙書德, 等. 預(yù)彎柱體減振器靜態(tài)特性的尺寸效應(yīng)[J]. 兵工學(xué)報(bào), 2009, 30(11): 1515–1521.

ZHAO Hua, WANG Min-jie, ZHAO Shu-de, et al. Size effect on static mechanics behaviors of isolator with Chevron-Shaped struts[J]. ACTA Armamentarii, 2009, 30(11): 1515–1521.

[10]YANG L M, SHIM V, LIM C T. A visco-hyperelastic approach to modelling the constitutive behaviour of rubber[J]. International Journal of Impact Engineering, 2000, 24(6): 545–560.

[11]傅志紅, 喻堅(jiān), 魏靈嬌. 聚碳酸酯的應(yīng)力松弛實(shí)驗(yàn)及數(shù)據(jù)處理分析[J]. 塑料工業(yè), 2014, 42(6): 89–92.

FU Zhi-hong, YU Jian, WEI Ling-jiao. The stress relaxation experiment and data processing analysis of polycarbonate[J]. China Plastics Industry, 2014, 42(6): 89–92.

[12]李樹虎, 賈華敏, 李茂東, 等. 超彈性體本構(gòu)模型的理論和特種試驗(yàn)方法[J]. 彈性體, 2011, 21(1): 58–64.

LI Shu-hu, JIA Hua-min, LI Mao-dong, et al. Theory and testing method of hyperelastic material constitutive model[J]. China Elastomerics, 2011, 21(1): 58–64.

[13]劉文武, 翁雪濤, 樓京俊, 等. 基于ANSYS對(duì)橡膠制品動(dòng)態(tài)分析的諧響應(yīng)法研究[J]. 武漢理工大學(xué)學(xué)報(bào)(交通科學(xué)與工程版), 2010, 34(5): 966–968.

LIU Wen-wu, WENG Xue-tao, LOU Jing-jun, et al. Study on the harmonic response method to analyse the dynamic characteristics of rubber products based on ansys[J]. Journal of Wuhan University of Technology (Transportation Science & Engineering), 2010, 34(5): 966–968.

Research on the equivalence and modelling of the nonmetallic isolator with complicate structure

SONG Yong1, JING Jian-quan2
(1. The 713 Research Institute of CSIC, Zhengzhou 450015, China; 2. School of Aerospace Eengineering, Beijing Institute of Technology, Beijing 100081, China)

In view of the requirement of numerical analysis and the equivalence for the nonmentallic isolator with complicate structure, The force-displacement relationship of compression was analyzed based on the 3dfinite element model with complicate structure, which show good agreement with the experiment result. A equivalent method based on the practicality experiment or numerical experiment was presented to reduce the modeling difficulty. The general shape and the connectivity of 3d equivalent structure is same with actual structure. The material parameter of the equivalent structure is determined according the stress relaxation test and the uniaxial compression experiment. Example analysis shows that the equivalent model can model the mechanical property and connectivity of the complicate structure in work direction.

isolator；hyperelastic；viscoelastic；equivalent model

TP535

A

1672 – 7619(2017)06 – 0088 – 06

10.3404/j.issn.1672 – 7619.2017.06.018

2017 – 04 – 10

宋勇(1985 – )，男，工程師，主要從事水下發(fā)射技術(shù)研究工作。