基于abaquscae的压力容器有限元分析.docxVIP

基于abaquscae的压力容器有限元分析 压力容器广泛应用于岩石、能源、科学研究、军事、工业和其他经济部门。因此，压力容器的设计和运营可靠性尤为重要。对压力容器进行正确的应力分析是压力容器设计的关键环节,这对其设计、制造和安全使用有直接影响。但是,传统的压力容器设计方法只是通过简单的计算以及对制造工艺、结构安装、使用运行等采取一系列强制规定或根据经验采取相关措施来保障压力容器的安全,并没有对其进行系统的应力分析与强度计算。这种做法容易从某种程度上掩盖问题的本质,从而给设备的安全运行带来不安全因素。 压力容器的有限元应力分析是一个典型的非线性接触问题。接触问题经常涉及由大变形引起的材料非线性和几何非线性,还有接触界面的非线性。接触界面的非线性主要源于两个方面:① 接触边界的非线性。接触边界条件(包括接触界面大小、相互位置及相应的接触状态)不是在计算的开始就可以全部给出的,往往需要在求解过程中进行确定。② 接触条件的非线性。接触条件包括接触物体不能相互侵入、接触力的法向分量只能是压力以及发生切向接触的摩擦条件。 ABAQUS是应用广泛的有限元软件,其主要优点是具有功能强大的非线性分析功能。ABAQUS能够自动选择合适的载荷增量和收敛准则并在分析过程中不断调整这些参数值,以确保获得精确的结果。本文以ABAQUS软件为分析平台,对极限内压下某实例压力容器进行强度分析计算,并将软件分析结果与ASME规范的计算结果进行了比较。 1 有限分析 1.1 单元类型的建立 文中分析的某特定条件的压力容器主体分为顶盖、卡箍和筒体3部分,顶盖与筒体通过一对卡箍和螺栓连接,见图1。各部件结构示图分别见图2～图4。在3个部件中,卡箍与筒体的材料均为35CrMoVA,其常温下的屈服强度ReL=930 MPa,抗拉强度Rm=1 080 MPa,弹性模量E=210 GPa,泊松比ν=0.3;顶盖材料为40Cr,常温下的屈服强度为ReL=785 MPa,抗拉强度Rm=980 MPa,弹性模量E=210 GPa,泊松比ν=0.3。 忽略次要结构后,考虑整体结构的对称性,建立容器的1/4有限元模型。建模时保留所有应力集中点(内倒角、凹点),简化处理结构非应力集中点(外倒角、凸点)。3个部件建模后按定位组装状态进行单元划分,单元类型为八节点六面体线性完全积分单元。为保证计算精度,所有关键部位均进行了网格细化。模型共划分了92 084个单元,模型及网格划分见图5。 1.2 接触对的建立 压力容器从装配到正常使用,3个部件之间都是相互接触的,所以在分析过程中必须建立部件之间的接触对,即筒体上顶面与顶盖下底面、顶盖周边与卡箍内侧、筒体周边与卡箍内侧。接触对的建立遵守以下2个原则:① 从属表面应是网格划分得更精细的表面。② 若网格密度相近似,从属表面应由更为柔软的材料组成。基于这样的原则,为使在计算分析中结果便于收敛,分别选取筒体上顶面、卡箍内侧作为接触从面,从而形成3个接触对,见图6。接触面之间的摩擦因数设为0.15。 1.3 预紧密封比压计算 根据本例压力容器的对称性,在筒体的中截面和对称剖面、顶盖对称剖面加对称约束,另外在筒体中截面上一点上施加固定约束来限制结构整体的刚体位移。加载过程分为两步:① 施加螺栓预紧力。②在顶盖下侧和筒体内侧施加极限工作压力80 MPa(保证计算结果与实际工况一致)。 根据ASMEVIII-1,卡箍连接螺栓的预紧载荷必须使密封垫片达到足够的预紧密封比压。预紧载荷W按下式计算: W=2DGbytan(α+μ) (1) 式中,b为垫片有效密封宽度,DG为垫片压紧力作用中心圆直径,mm;y为预紧密封比压,MPa;α为卡箍倾斜面倾角,μ为摩擦角,(°)。本例垫片类型为石墨波形弹性垫片,b=9.8 mm,取比压力y=69.6 MPa。 1.4 顶盖、顶盖应力分析 对整体模型进行弹性分析, 得到整体结构的应力分布,其中卡箍的应力分布见图7,顶盖的应力分布见图8,筒体的应力分布见图9。 从图7中可以看出,危险点出现在卡箍与顶盖的接触面及卡箍内倒角附近,其Mises等效应力与最大主应力分别为964 MPa与974 MPa。卡箍最大主应力略大于Mises应力,原因在于应力三轴度较大,因而卡箍容易发生准脆性破坏。 从图8中可以看出, 顶盖的最大Mises应力出现在顶盖与筒体接触的底面内侧应力集中点,大小为798 MPa,而最大主应力点出现在顶盖与卡箍接触面根部,大小为587 MPa。 从图9可以看出,筒体的最大Mises应力点和最大主应力点分别出现在筒体与卡箍接触面及筒体反缘根部应力集中点,大小分别为642 MPa和680 MPa。 2 采用asme规范计算 2.1 密封垫片接枝公式 根据ASME Ⅷ-1,卡箍(图2)接触面的总载荷包括内压引起