用ABAQUS軟件分析扭轉(zhuǎn)問題

1、用ABAQUS軟件分析扭轉(zhuǎn)問題         1、引言    在處理扭轉(zhuǎn)問題時，常規(guī)的計算方法，往往會伴隨一些假設(shè)，這會降低了結(jié)果的準確程度。根據(jù)有限元理論，使用有限元軟件求解扭轉(zhuǎn)問題會大大提高求解的精確度，特別是對復(fù)雜的結(jié)構(gòu)，效果更為明顯。本文以橡膠產(chǎn)品為例，討論的在ABAQUS軟件中，如何正確完成扭轉(zhuǎn)分析，并提取需要的分析結(jié)果。    2、問題描述    受扭轉(zhuǎn)件結(jié)構(gòu)由鋼筒和橡膠筒組成，產(chǎn)品尺寸如圖1所示。圖1產(chǎn)品結(jié)構(gòu)簡圖

2、0;   3、有限元建模    加載時內(nèi)芯固定，在外圈施加扭轉(zhuǎn)位移。根據(jù)產(chǎn)品的CAD結(jié)構(gòu)建立有限元模型如圖2所示：圖2產(chǎn)品有限元模型圖    4、材料性質(zhì)定義    鋼：彈性模量EX=2×105MPa，泊松比=0.3    橡膠：橡膠是一種超彈性材料，對于超彈性材料，不用楊氏模量和泊松比，而用應(yīng)變勢能（U）來表達應(yīng)力應(yīng)變關(guān)系。ABAQUS軟件中有兩種應(yīng)變勢能可利用，分別是多項式模型和奧根（Ogden）模型，本例中使用多項式模型，表達式如下：

3、60;       式中：U應(yīng)變勢能，Jel彈性體積比；I1、I2應(yīng)變不變量；Di定義材料的壓縮性；CijRinvlin系數(shù)。本例中取N=1，以橡膠材料的單軸拉伸，單軸壓縮和平面剪切實驗數(shù)據(jù)為依據(jù)，并考慮到橡膠的不可壓縮性，輸入方程系數(shù)值：C01=0.36，C10=0.09，D1=0    注意事項: 橡膠的特性錯綜復(fù)雜，材料特性和幾何特性均呈非線性變化的。如果要準確預(yù)測模型中發(fā)生變形或應(yīng)變部分的行為，那么提供的試驗數(shù)據(jù)的范圍要涵蓋計算模型中可能會出現(xiàn)的變形狀態(tài)和應(yīng)變范圍。    5、

4、加載求解    加載時，內(nèi)鋼筒的內(nèi)套固定，即UX=UY=UZ=0,將外鋼筒的最外層結(jié)點的坐標系定義為柱坐標系。在此柱坐標系中施加扭轉(zhuǎn)載荷。    注意事項：在分析過程中最外層結(jié)點的柱坐標系的方向始終保持與結(jié)點原始柱坐標系的方向相同，并不隨著結(jié)點的移動而改變方向。因此，施加的位移量并不是結(jié)點始末位置間的弧線長度，而應(yīng)是結(jié)點始末位置間的直線距離。所以此例中所加的位移載荷應(yīng)為：UR=60*sin30°=30(mm)UF=60*cos30°-60=-8.038(mm)    6、結(jié)果分析

5、60;   將模型信息提交給求解器，經(jīng)迭代求解出結(jié)果：橡膠部分的VonMises 應(yīng)力云圖如圖3所示 ：圖3 橡膠的Von Mises應(yīng)力云圖    由上圖看出，橡膠上的應(yīng)力呈層狀分布，且隨著半徑的增加承遞減趨勢。其遞減情況如圖4所示。應(yīng)力的這種隨半徑遞減的分布趨勢，與實際使用中，產(chǎn)品率先在中間發(fā)生破壞的情況是相吻合的。    由分析求得的反作用力結(jié)果如下表所示        計算扭矩需要沿圓周切線方向的分量，而這里的F合是沿弦線的合力，其沿切線方向的分量應(yīng)為, F切=F合*cos15°。計算的扭矩結(jié)果與實驗結(jié)果的對照情況如下表所示        產(chǎn)生誤差的原因大致如下： 實測操作的本身；計算過程中模型的簡化，材料參數(shù)的選擇；在對連續(xù)的結(jié)構(gòu)做離散處理時。誤差在工程的許可范圍內(nèi)，分析計算結(jié)果是可信的。    7、結(jié)論    本文利用ABAQUS軟件分析了超彈性橡膠材料受扭時的性能，分析結(jié)果與實驗結(jié)果的對比表明，計算結(jié)果與實測結(jié)果吻合的很好，對確定產(chǎn)品的使用方案有指導(dǎo)意義。