基于ABAQUS的内压厚壁圆筒的弹塑性分析论文.docVIP

基于 ABAQUS的内压厚壁圆筒的弹塑性分析 学院：航空宇航学院专业：工程力学 指导教师： 姓名： 学号： 1 1. 问题描述 一个受内压的厚壁圆筒（如图 1 ），内半径和外半径分别为 a 10mm 和 b 15mm（外径与内径的比值 b 15 1.5 1.2 ），受到均匀内压 p 。材料为理想 a 10 弹塑性钢材（如图 2），并遵守 Mises 屈服准则，屈服强度为 Y 380MPa ，弹性 模量 E 200GPa ，泊松比 0.3。 图 1 内压作用下的端部开口厚壁圆筒 图 2 钢材的应力 -应变行为 首先通过理论分析理想弹塑性材料的厚壁圆筒受内压作用的变形过程和各 阶段的应力分量， 确定弹性极限压力 pe 和塑性极限压力 p p ；其次利用 ABAQUS 分析该厚壁圆筒受内压的变形过程， 以及各个阶段厚壁筒内的应力分布， 与理论 分析的结果进行对比，验证有限元分析的准确性。 2. 理论分析 2.1 基本方程 由于受到内压 p 的作用，厚壁圆筒壁上受到径向压应力 r 、周向压应力 和 轴向应力 z 的作用，由开口的条件可推出 z 0 。因为这是一个轴对称问题， 所有的剪应力和剪应变均为零。平衡方程和几何方程用下式表示： d r r - 0 dr r r dur , ur dr r 弹性本构关系为： 1 * , 1 * r r E* r E*  （1） （2） （3） 2 由于此问题为平面应变问题，所以上式中 E* E 2 *1 1 相应的边界条件为： r r a , 0 （ ） p r r b 4 2.2 弹性阶段 根据弹性力学中的应力解法：取应力分量 r ， 为基本未知函数，利用平 衡方程和应力表示的协调方程联合求解，可得应力分量的通解 r C1 C2 r 2 C1 - C2 r 2 将边界条件带入可得应力分量为： a2 p b2 r - b2 a2 r 2 - 1 r a2 p b2 1 （5） b2 a 2 r 2 因为 a r b ，所以 r 0且 0 ，可以观察到： z 0 r ， 分析采用 Mises 屈服准则，表达为 2 2 2 6 2 2 2 2 2 （ 6） r rz z r rz z Y 该厚壁圆筒是轴对称平面应变问题，即 r rz z 0 ，由 Mises 屈服条 件其表达式可得到： 2 1.155 （7） r Y Y 3 当内压 p 较小时，厚壁圆筒处于弹性状态，在 r a 处， - r 有最大值， 筒体由内壁开始屈服，此时的内压为 pe ，由式（ 5）、（ 7）联立可求得弹性极限 压力为 1.155 b2 a2 Y （8） pe 2 2b 代入题目所给数据得到弹性极限压力为： 3 1.155 380 15 2 10 2 pe 2 152 121.92MPa 弹塑性阶段 p pe 时，圆筒处于弹性状态，当 p pe 的情况，在圆筒内壁附近出现塑 性区，产生塑性变形，随着内压的增大，塑性区逐渐向外扩展，而外壁附近仍为 弹性区。由于应力组合 - r 的轴对称性，塑性区与弹性区的分解面为圆柱面。 可用 r c （ a c b ）表示弹塑性边界，对于 a r c ，圆筒就处于塑性状态， 对于 c r b ，圆筒就处于弹性状态。 分别考虑两个变形区，同时在弹塑性边界 r c 上满足 r 的连续条件，这样 可以得到塑性区（ a r c ）内的应力分量为 r 1.155 Y ln r p a （9） 1.155 Y (1 ln r ) p a 弹性区（ c r b ）的应力分量为 r 1.155 Y c2 b2 1 2b2 r 2 1.155 Y c2 b2 （ 10） 1 2b2 r 2 2.4 塑性极限分析 随着压力的增加，塑性区不断扩大，当 c b 时，整个截面达到塑性状态， 即圆筒达到塑性极限状态，此时的压力不能继续增加，将 c b 带入式（ 9），可 得塑性极限压力为 pY 1.155 Y ln b （ 11） a 代入题目所给数据得到塑性极限压力为： pY 1.155 380 ln 15 177.96MPa 10 3. 有限元分析 4 3.1 有限元模型 此问题属于平面应变问题，采用二维有限元模型，选取厚壁圆筒的 1/4 部分 作为分析模型，其内径为 a 10mm、外径为 b 15mm，厚壁圆筒轴向无限延长， 则模型处于平面应变状态。 3.2 材料属性定义 圆筒材料为钢材，弹性模量 200Gpa，屈服强度 380Mpa，泊松比 0.3，截面 属性选用实体、匀质，采用理想弹塑性本构关系。 3.3 分析步的定义 由于是非线性分析， Step 中设置分析过程和输出要求选择静态分析， 最小分 析步取 0.05，最大分析步取 0.1，输出