abaqus扭轉屈服的案例
Abaqus扭轉仿真案例講解
Abaqus扭轉仿真案例講解
用ABAQUS軟件分析扭轉問題
1、引言
在處理扭轉問題時，常規的計算方法，往往會伴隨一些假設，這會降低了結果的準確程度。根據有限元理論，使用有限元軟件求解扭轉問題會大大提高求解的精確度，特別是對復雜的結構，效果更為明顯。本文以橡膠產品為例，討論的在ABAQUS軟件中，如何正確完成扭轉分析，并提取需要的分析結果。
2、問題描述
受扭轉件結構由鋼筒和橡膠筒組成，產品尺寸如圖1所示。
圖1產品結構簡圖
3、有限元建模
加載時內芯固定，在外圈施加扭轉位移。根據產品的CAD結構建立有限元模型如圖2所示：
圖2產品有限元模型圖
4、材料性質定義
鋼：彈性模量EX=2×105MPa，泊松比μ=0.3
橡膠：橡膠是一種超彈性材料，對于超彈性材料，不用楊氏模量和泊松比，而用應變勢能（U）來表達應力—應變關系。ABAQUS軟件中有兩種應變勢能可利用，分別是多項式模型和奧根（Ogden）模型，本例中使用多項式模型，表達式如下：
式中：U—應變勢能，Jel—彈性體積比；I1、I2—應變不變量；Di—定義材料的壓縮性；Cij—Rinvlin系數。本例中取N=1，以橡膠材料的單軸拉伸，單軸壓縮和平面剪切實驗數據為依據，并考慮到橡膠的不可壓縮性，輸入方程系數值：C01=0.36，C10=0.09，D1=0
注意事項: 橡膠的特性錯綜復雜，材料特性和幾何特性均呈非線性變化的。如果要準確預測模型中發生變形或應變部分的行為，那么提供的試驗數據的范圍要涵蓋計算模型中可能會出現的變形狀態和應變范圍。
5、加載求解
加載時，內鋼筒的內套固定，即UX=UY=UZ=0,將外鋼筒的最外層結點的坐標系定義為柱坐標系。在此柱坐標系中施加扭轉載荷。
展開 ABAQUS C30混凝土短柱扭轉仿真
最近做了混凝土短柱的扭轉案例，需要的可以下載。
Torsion column.rar
基于ANSA與ABAQUS的橡皮筋扭轉扇葉簡單實例
基于ANSA與ABAQUS的橡皮筋扭轉扇葉簡單實例
Abaqus 采用YLD2000-2D屈服準則的UMAT子程序
Barlat在2003年提出了專門針對平面應力問題的各向異性屈服準則，該屈服準則對于各向異性材料具有很高的精度，得到了廣泛的應用。
YLD2000-2D屈服面示意圖
Yld2000-2d屈服準則由下式給出
其中
矩陣X′和X″的元素分別由柯西應力的下列線性變換獲得
L′和L″的分量由下式求得
積分算法采用徑向返回算法，該方法是穩健和精確的。
當彈性試算超出屈服面時，則需要進行塑性修正
使其滿足
公式9可以通過牛頓法進行迭代求解。
計算的應力應變曲線如下圖所示
B, F. Barlat A , et al. "Plane stress yield function for aluminum alloy sheets—part 1: theory." International Journal of Plasticity 19. 9(2003):1297-1319.
王海波, 萬敏, 閻昱,等. 屈服準則在有限元軟件中實現的正確性驗證[J]. 固體力學學報, 2010, 031(002):173-180.
展開 針對平面應力問題的YLD2000-2D屈服準則及其在ABAQUS中UMAT子程序的實現
Barlat在2003年提出了專門針對平面應力問題的各向異性屈服準則，該屈服準則對于各向異性材料具有很高的精度，得到了廣泛的應用。
YLD2000-2D屈服面示意圖
Yld2000-2d屈服準則由下式給出
其中
矩陣X′和X″的元素分別由柯西應力的下列線性變換獲得
L′和L″的分量由下式求得
積分算法采用徑向返回算法，該方法是穩健和精確的。
當彈性試算超出屈服面時，則需要進行塑性修正
使其滿足
公式9可以通過牛頓法進行迭代求解。
計算的應力應變曲線如下圖所示
B, F. Barlat A , et al. "Plane stress yield function for aluminum alloy sheets—part 1: theory." International Journal of Plasticity 19. 9(2003):1297-1319.
王海波, 萬敏, 閻昱,等. 屈服準則在有限元軟件中實現的正確性驗證[J]. 固體力學學報, 2010, 031(002):173-180.
最后，有需要歡迎通過微信公眾號聯系我們。
微信公眾號：320科技工作室。
展開 ABAQUS中mises應力云圖顯示的最大值還不到屈服應力值為啥還有PEEQ值
ABAQUS中mises應力云圖顯示的最大值還不到屈服應力值為啥還有PEEQ值，PEEQ云圖有變形值
