理想塑性材料厚壁圓筒解析解與數(shù)值解對比研究

摘 要：文章針對理想塑性材料特性，選取厚圓筒壁進行了解析解與數(shù)值解對比分析。通過對比發(fā)現(xiàn)，有限元解和理論解相差很??；當筒體內(nèi)部處于塑性狀態(tài)，外層處于彈性狀態(tài)，當壓力卸除后，筒體內(nèi)層塑性區(qū)將有殘余變形存在，而外層彈性區(qū)受到內(nèi)層塑性區(qū)殘余變形的阻擋而不能完全恢復，結(jié)果使內(nèi)層塑性區(qū)受到外層彈性區(qū)的壓縮而產(chǎn)生殘余壓應力，而外層彈性區(qū)由于收縮受到阻擋而產(chǎn)生殘余拉應力。

關(guān)鍵詞：理想塑性材料；厚圓筒壁；解析解；數(shù)值解

1 計算工況

受均勻內(nèi)壓（p=12.5kg/cm）作用的理想塑性材料厚壁圓筒，其幾何參數(shù)為：內(nèi)徑Ri=10cm，外徑Re=20cm；材料參數(shù)為：E=86666.7kg/cm2，v=0.3，？滓s=17.32kg/cm2的理想塑性材料。厚壁筒計算模型長度取H=20cm， 在子午面上沿徑向劃分八個以上的八結(jié)點等參單元。從初始狀態(tài)開始，歷經(jīng)加載（內(nèi)壓到達p=12.5kg/cm）、然后完全卸載（p=0）。這一過程之后，求厚壁筒內(nèi)的殘余應力沿徑向r的應力（？滓r，？滓？茲）分布曲線。

2 數(shù)值解計算模型

建立有限元模型，在子午面上沿徑向劃分10個八結(jié)點等參單元。劃分單元以及結(jié)點如圖1所示。

3 計算結(jié)果及對比

將該厚壁圓筒的幾何參數(shù)代入理論解析解中可以得到，彈性極限載荷為7.4873kg/cm，塑性極限載荷理論解為Pp=13.8629kg/cm，塑性半徑為15.03cm。在題目中給出的均勻內(nèi)壓是12.5kg/cm，達不到塑性極限，但是超過彈性極限荷載，所以厚壁圓筒的一部分處于塑性狀態(tài)，一部分處于彈性狀態(tài)。通過計算結(jié)果可以發(fā)現(xiàn)第5個單元完全進入塑性，第6個單元都沒有進入塑性，所以，近似認為第5第6個單元交界處為塑性分界面，塑性半徑為15cm。通過有限元計算，加載、卸載后的結(jié)果如表1所示。

通過以上的比較可以看出，有限元解和理論解相差很小。當筒體內(nèi)部處于塑性狀態(tài)，外層處于彈性狀態(tài)，當壓力卸除后，筒體內(nèi)層塑性區(qū)將有殘余變形存在，而外層彈性區(qū)受到內(nèi)層塑性區(qū)殘余變形的阻擋而不能完全恢復，結(jié)果使內(nèi)層塑性區(qū)受到外層彈性區(qū)的壓縮而產(chǎn)生殘余壓應力，而外層彈性區(qū)由于收縮受到阻擋而產(chǎn)生殘余拉應力。

參考文獻

[1]朱菊芬，等.非線性有限元及其在飛機結(jié)構(gòu)設(shè)計中的應用[M].上海：上海交通大學出版社，2011.

[2]李元媛，淡勇，蔡睿賢.理想塑性軸對稱平面問題的解析解[J].機械工程學報，2009，7：270-273+277.