橡膠粘滯生熱的案例
abaqus橡膠熱仿真：減振橡膠疲勞黏滯生熱的仿真分析-源文件與子程序詳解
采用一種修正的Kraus模型定量描述了橡膠材料動態損耗模量隨溫度、載荷頻率和應變幅值的變化規律。得到了生熱率與溫度、載荷頻率和應變幅值的函數關系式。
利用依黏彈性理論得出的黏滯生熱率與溫度、載荷頻率和應變幅值的函數關系式，編制了相應的計算程序。建立了減振橡膠疲勞黏滯生熱的有限元分析方法。
通過將經典疲勞模型中用作疲勞壽命預測指標的最大主應變替換為穩態溫升，在冪律模型的基礎上開發了一種方法來快速評估橡膠結構的疲勞壽命。
08
—
源文件與操作步驟（沙漏試樣為例）
8.1分析流程
仿真分析主要包括三個環節：變形分析、熱源計算與熱分析。（1）在變形分析環節，對材料和減振元件施加設定的載荷歷史，采用超彈性本構描述橡膠材料的力學行為，求解每個加載時刻有限元模型中各積分點的應變狀態；（2）在熱源計算環節，對應每一加載時刻，將變形分析中對應的載荷頻率、應變狀態（動態應變幅值）以及熱分析中得到的溫度作為輸入變量，通過自編的Fortran語言子程序，計算得到各積分點的黏滯生熱率；（3）依已知的材料參數和問題的熱邊界條件進行Abaqus熱分析，得出溫度分布后再將溫度場數據返回到自編子程序，對黏滯生熱強度和溫度場進行迭代計算，從而得出橡膠材料和減振元件各位置的溫升歷程。
展開 減振橡膠疲勞黏滯生熱的仿真分析-源文件與子程序詳解
采用一種修正的Kraus模型定量描述了橡膠材料動態損耗模量隨溫度、載荷頻率和應變幅值的變化規律。得到了生熱率與溫度、載荷頻率和應變幅值的函數關系式。
利用依黏彈性理論得出的黏滯生熱率與溫度、載荷頻率和應變幅值的函數關系式，編制了相應的計算程序。建立了減振橡膠疲勞黏滯生熱的有限元分析方法。
通過將經典疲勞模型中用作疲勞壽命預測指標的最大主應變替換為穩態溫升，在冪律模型的基礎上開發了一種方法來快速評估橡膠結構的疲勞壽命。
08
—
源文件與操作步驟（沙漏試樣為例）
8.1分析流程
仿真分析主要包括三個環節：變形分析、熱源計算與熱分析。（1）在變形分析環節，對材料和減振元件施加設定的載荷歷史，采用超彈性本構描述橡膠材料的力學行為，求解每個加載時刻有限元模型中各積分點的應變狀態；（2）在熱源計算環節，對應每一加載時刻，將變形分析中對應的載荷頻率、應變狀態（動態應變幅值）以及熱分析中得到的溫度作為輸入變量，通過自編的Fortran語言子程序，計算得到各積分點的黏滯生熱率；（3）依已知的材料參數和問題的熱邊界條件進行Abaqus熱分析，得出溫度分布后再將溫度場數據返回到自編子程序，對黏滯生熱強度和溫度場進行迭代計算，從而得出橡膠材料和減振元件各位置的溫升歷程。
8.2建模
abaqus/cae操作
8.3賦材料屬性
鋼：
CAE操作：
inp文件：
*Material, name=ste*Conductivity43.
展開 巧用單元提高Abaqus計算效率：帶扭曲的軸對稱單元-橡膠阻尼器內摩擦生熱分析 ￥49.99
橡膠阻尼器的內摩擦生熱分析-節點溫度云圖
比如，假設上圖中的阻尼器不再是長方體，而是回轉體，且發生軸向扭曲變形，那么能不能用軸對稱單元來建模呢？
答案是可以的，在Abaqus的軸對稱單元系里還有一種可考慮Twist的單元，即帶字母G標識的那種類型，能夠在分析時充分考慮回轉體的整體扭轉變形。
首先，我們可以在part模塊使用Axisymmetric建立環形塊狀阻尼器的回轉截面；然后在mesh模塊劃分好四邊形網格；最后，定義單元類型為CGAX4T，即帶扭曲的4節點軸對稱位移-溫度耦合單元。
這里的橡膠阻尼器材料本構采用的是超彈性模型，應變能描述形式為Neo Hooke，再結合時域黏彈性Prony參數與非彈性變形能耗散比，來計算阻尼器周期性扭轉過程中的材料內摩擦生熱。
阻尼器上、下兩個端面的節點分別使用位于回轉軸上的兩個參考點來耦合，固定下端面參考點，并在上端面參考點施加軸向的周期性扭角位移。
阻尼器的回轉結構與網格-單元
雖然建模時只考慮了回轉截面，但是帶扭曲的軸對稱單元可以將回轉體發生扭轉時的整體結構響應考慮在內，這是因為這種單元多了一個扭轉自由度5，拿本例中的位移-溫度耦合單元CGAX4T來說，該單元的節點具有1、2、5和11四個自由度。
回轉截面上的節點溫度云圖
當然，我們還可以通過后處理設置，將阻尼器的三維實體以部分或全貌的形式顯示出來，這樣方便更直觀地看到整個結構的扭曲變形。
后處理定義回轉角度和網格堆疊
在周期性扭轉過程中，由于橡膠材料的黏彈性內摩擦耗散，導致阻尼器溫度逐漸升高，計算結果表明，30分鐘內阻尼器從開始時的室溫22℃升高到40℃。
阻尼器的升溫曲線
注：本案例的動圖為跳幀模式，直觀上不反映阻尼器真實的扭轉頻率。
展開 