摘要:文章对压力机用压力容器进行仿真分析,通过三维实体建模,得出了压力容器自重作用下的应力、位移分布,应力最大位移出现在压力容器入孔周围,支座与压力容器连接部应力最大。采用水作为施压介质,进行静水压力计算分析,因载荷的增加,所得到的分析结果更具真实性及可靠性。旨在通过对压力容器的有效分析,为其安全性提供合理保障。
加入收藏
静力学分析具有简明的操作流程。分为三个阶段:前处理、求解、后处理。前处理器建立二维、三维模型旭可与三维建模软件关联直接导人。结构中参数定义、网格划分也在此模块完成。求解模块施加载荷、约束,对求解形式进行设置。通用后处理器查看结果云图、输出数据。本文分析中,由于压力容器具有曲面结构,需要建立模型。软件在建模方面更加人性化,有更好的人机交换界面,适用于三维复杂零件的建模。本文中有限元分析的流程图如图1所示。
图1 有限元分析流程
本文以压力机使用的压力容器为例来计算分析其承载下的力学响应,压力容器筒体和封头均采用厚度为14ram的Q345R,支座设置为刚体。压力容器相关参数如表1所示。
1、压力容器有限元分析
1.1 建立压力容器的三维模型
本文中使用三维建模软件Pro—E对压力容器建立三维模型。因为有些如液位计所需的管口直径较小,对下面的分析结果影响很小,所以可将口径较小的管口去除,达到简化模型,提高有限元软件分析速度和效率的目的[1]。建立的三维模型如图2所示。
表1 压力容器结构尺寸;图2 压力容器三维模型
1.2 压力容器的网格划分
对于三维几何结构模型,ANSYS Workbench有下面的划分网络法:多域法(MuhiZone),这是Ansys Workbench网格划分的亮点之一,主要用来花划分六面体网格,(其特点就是)它本身具有几何体自动分解的功能,从而产生六面体网格圈;由于压力容器模型相对比较简单,而且Workbench软件的自动划分网格的质量和精度都能达到要求,故本文采用自动划分网格的功能对模型进行网格划分。当然,使用自动划分网格的功能时也可以对网格的数量和大小进行控制,主要功能选项在Sizing栏,包括:
Relevance:网格相关度,数值从一100至100,代表网格由疏到密;
Relevance Center:代表网格Coarse(稀疏)、Medium(中等)、Fine(细化);
SpanAngle Center:由于弯曲区域网格细分,包括Coarse(稀疏)、Medium(中等)、Fine(细化)。
本文划分网格时Relevance Center和SpanAngle Center选项均采用Fine级别进行划分,得出模型的网格划分如图3所示。
图3 压力容器网格划分
1.3 定义材料属性
在线性静力结构分析中,材料属性只要定义杨氏模量以及泊松比。同时必须明白:①假如有任何惯性载荷,必须要定义材料的密度;②热膨胀系数(可输入负热膨胀系数)和传热系数在有热载荷的时候才需要确定;③若要进行疲劳分析则需要FatigueModuleadd—onlicense。压力容器筒体与封头均采用Q345R为材料,其材料参数如表2所示。
表2 Q345R材料属性
1.4 施加边界条件和载荷
1.4.1 边界条件
因为容器为卧式压力容器,且安置在鞍座上,因此需定义鞍座的约束条件以及鞍座与压力容器间的接触类型。
(1) 约束条件。定义鞍座左支腿约束类型为FixedSupport(固定约束);右支腿约束类型为给定位移(Displacement)约束,规定在压力容器轴向即z轴方向位移为Free,其他方向即x轴和I,轴位移为0。
(2) 接触类型。Workbench软件中提供五种接触方式,其中包括:
Bonded(绑定):这是AWE中关于接触的默认设置。如果接触区域被设置为绑定,不允许面或线间有相对滑动或分离,可将此区域看做被连接在一起。因为接触长度/面积是保持不变的,所以这种接触可以用作线性求解。如果接触是从数学模型中设定的,程序将填充所有的间隙,忽略所有的初始渗透。
NoSeparation(不分离):这种接触方式和绑定类似。它只适用于面。不允许接触区域的面分离,但是沿着接触面可以有小的无摩擦滑动。
Frictionless(无摩擦):这种接触类型代表单边接触,即,如果出现分离则法向压力为零。只适用于面接触。因此,根据不同的载荷,模型间可以出现问隙。它是非线性求解,因为在载荷施加过程中接触面积可能会发生改变。假设摩擦系数为零,因此允许自由滑动。使用这种接触方式时,需注意模型约束的定义,防止出现欠约束。程序会给装配体加上弱弹簧,帮助固定模型,以得到合理的解。
Rough(粗糙的):这种接触方式和无摩擦类似。但表现为完全的摩擦接触,即没有相对滑动。只适用于面接触。默认情况下,不自动消除间隙。这种情况相当于接触体间的摩擦系数为无穷大。
Frictional(有摩擦):这种情况下,在发生相对滑动前,两接触面可以通过接触区域传递一定数量的剪应力。有点像胶水。模型在滑动发生前定义一个等效的剪应力,作为接触压力的一部分。一旦剪应力超过此值,两面将发生相对滑动。只适用于面接触。摩擦系数可以是任意非负值。
因为鞍座和压力容器大多数是通过焊接进行连接的,所以两者之问不会产生相对位移,因此综上可以设置接触方式为Bonded(绑定)。
1.4.2 载荷加载
设计工况包括压力容器空载自重,装满水时静水压力和内压载荷这几种载荷单独作用或组合作用时的情况,下文将分别就各种工况进行详细的分析计算。
1.5 设置需要的结果
在Workbench软件中的Mechanical模块中的后处理中可以得到各种不同的结果,如各个方向的变形和应力、应变、主应力应变或者应力应变不变量、接触输出、支反力等。用户只需选择插入所需的选项,软件分析结束后会以图解的形式将所需结果直观表现出来。
2、压力容器在自重下的力学响应
2.1 压力容器的有限元分析
由于压力容器筒体身长6000ram,内径1600mm,筒壁厚14ram,这样的压力容器白重下的力学响应分布是不容忽视的,通过在有限元软件ANSYSWorkbench中定义相关的量来计算说明,在ANSYSWorkbench中,有限元模型的体积会自动计算得出,同时根据设定的材料密度P可以计算出其质量,再定义重力加速度g就可以得到压力容器的重量(材料参数前面已给出,这里我们定义重力加速度g=9.8m/s2)。
在完成网格划分,材料参数输入,定义约束几个步骤后,通过工具条[intertial]选择[StandardEarthGravity]选项,对压力容器和鞍座施加重力加速度,并设置其方向垂直向下,即l,轴负方向。并插入总变形:[Deformation]一[Total],应力强度[Stress1-[Inten—sity]。然后运行求解选项[Solve],完成后得到应力应变和变形云图,分别如图4、图5所示。
图4 压力容器应力应变云图;图5 压力容器位移五图
从压力容器自重作用下的应力、位移分布看出:应力最大部位在支座与压力容器连接部位,为3.75MPa。这是因为支座与压力容器连接部位采用的是Bonded接触(相当于焊接在一起),在连接边缘处存在着应力集中现象;最大位移出现在压力容器中间部位,即人孔周围,最大位移为0.135mm。
2.2 压力容器自重作用下的结果分析
提取支座反力并和理论值比较:
(1) 根据理论计算压力容器自重:
G=mg=6610kgX9.8N/kg=64787N(1)
(2) Workbench中提取左右两个支座的支反力。左右支座支反力之和:
F=64869(N)
(3) 误差:
(pG)/G=(64869—64787)/64787=0.1265%(2)
Workbench计算的支反力和压力容器自身重力误差在5%以内,因此可以认为压力容器自重作用下网格有效,接下来的计算中采用同一套网格来计算。
3、压力容器在静水压力下的力学响应
因为液氨的密度会随着温度有较大的变化,所以为了更好的模拟压力容器装满液体时的工况,可以采用水来代替液氨进行静水压力分析计算,同时因为载荷的增加使得分析结果更具可靠性。
为了模拟该工况,在前处理完成后,施加静水压力载荷的同时应设置重力加速度。重力加速度设置如前一工况相同,施加静水压力载荷的过程为:首先,在[Definition]一[FluidDensity]中设置水的密度为le一6kg/mm3;然后在[HydrostaticAccelerationl-[Ycomponent]中设置y轴方向的加速度为9800mm/s2;最后在[FreeSurfaceLocation]一[YCoordinate]中设置液面高度为y轴正方向786mm。
液体由于重力会对容器不同深度处各个方向产生的压强不用,压强随着深度的增加而增大,如图6所示。液体压强的计算公式:
P=pgh(3)
图6 压强随深度变化示意图
设置完成后可以得到静水压力图如,由图可知,压力随着深度的增大而增大,完全符合上述公式,可知静水压力施加方法是正确的。在确定静水压力正确施加后,即可运行求解选项[Solve],得到应力应变和变形云图如图8、9所示。
图7 压力容器静水压力云图;图8 压力容器应力应变云图;图9 压力容器位移云图
由图8和图9可知,在装满水的工况下,压力容器所受应力最大的地方仍然与受自重载荷时相同,均是鞍座与筒体连接部位,与预计的结果相符。最大应力为15.55MPa。最大位移为0.544mm,位置仍然在人孔开口附近,由于压力容器的长度相较于筒径较长,故中间部位出现最大挠度是可以预期的。故总体分析结果仍然可靠。
4、小结
(1) 压力容器自重作用下支反力之和与压力容器自重相比,两者误差在1%以内,可以认为相等;压力容器内水压作用下支反力之和与压力容器中水的重力有效;
(2) 由不同节点的应力数值看出,内压对压力容器应力应变分布起着主要影响,压力容器自重和内部液体对压力容器应力应变分布影响很小;
(3) 压力容器开孔接管附近存在应力集中现象,由于应力集中产生的应力超过了材料的屈服强度,但这种屈服是局部的,压力容器整体仍处于弹性阶段,因此不会导致压力容器失效。
参考文献:
[1]潘家祯.压力容器材料实用手册[M].北京:化学工业出版社,2000.
[2]连耀东.ANSYS Workbench结构分析快速入门指南[M].北京:电子工、世出版社,2015.
陈启升,张红梅,卞光伟,等.压力机用压力容器的静力学研究[J].锻压装备与制造技术,2019,54(5):96-99.
分享:
气动伺服弹性是一门涉及空气动力学、结构动力学和控制理论的交叉学科,目前在航空领域与该学科有密切关联的研究为颤振主动抑制和阵风减缓[1,2]。而随着现代大型民机的柔性不断增加,导致阵风载荷的危害剧增,使得阵风载荷的主动抑制研究成为一种新的趋势[2]。
2020-08-10随着经济全球一体化的深入发展,高等教育国际化已成为一种必然趋势。为应对新时代对工科人才培养提出的高要求,国家适时提出了新工科建设,并在高校开展了工程教育专业认证工作。在新时代的背景下,作为传统的工科专业课程,理论力学在课程教学中也需要紧跟时代节奏,更新与改革迫切。
2020-08-10新工科是为适应人工智能、全智能时代,响应国家《中国制造2025》等战略需求,培养实践能力强、创新能力强、具备国际竞争力的高素质复合型人才,高等工程教育承担的职责与使命。“新工科”这一概念自2016年提出以来,在不到一年的时间里,教育部组织高校进行深入研讨,形成了“复旦共识”和“天大行动”[1]等,由此可以看出国家对新工科发展的全局考量。
2020-08-10“理论力学”课程是高等工科学校的技术基础课,是连接基础理论课与工程专业课的桥梁。课程的任务是使学生掌握质点系和刚体机械运动(包括平衡)的基本规律及其研究方法,为学习有关的后续课程打好必要的基础,初步学会应用理论力学的理论和方法分析解决一些简单的工程实际问题,同时结合课程的特点,培养学生的辩证唯物主义世界观,以及分析问题和解决问题的能力。
2020-08-10理论力学、材料力学课程是军队院校相关专业本科学员必修的专业技术基础课程,是连接基础知识和武器装备专业知识的桥梁。其具有双重属性,既有基础性的特点,理论经典、逻辑严谨、论证严密,又有应用性的广泛,能紧密联系实际,解决工程实际问题。正是由于课程的理论经典,且具有通用性,大部分教材列举的实例均是通用型例子。
2020-08-10理论力学是一门理论性较强的专业技术基础课。它是研究力学中最普遍、最基本的规律,是后续力学课程的重要基础[1]。受力分析是理论力学教学中极其重要的教学内容,是学好其它力学课程的根基。倘若这部分内容掌握不好,就不利于后续理论力学内容和其它力学课程的顺利学习。
2020-08-10理论力学是一门理论性较强的专业技术基础课。它是研究力学中最普遍、最基本的规律,是后续力学课程的重要基础[1]。受力分析是理论力学教学中极其重要的教学内容,是学好其它力学课程的根基。倘若这部分内容掌握不好,就不利于后续理论力学内容和其它力学课程的顺利学习。受力分析在中学和大学物理课程中已经学习过,学生有些基础。
2020-07-08大气动力学是气象学的一个分支,是应用物理学定律和数学方法研究天气和气候中与大气运动有关的现象,从理论上探讨大气环流、天气系统演变和其它大气运动过程的学科。该课程是本科阶段大气动力学(动力气象学)课程的后续课程,侧重于基础理论扩展和对更复杂的强迫耗散系统及非线性系统问题的讨论。
2020-06-04前人研究已揭示了封闭体系适用于Ⅲ型烃源岩,而Ⅰ、Ⅱ型烃源岩生烃特征研究时采用开放体系更为合理[2]。开放体系下生烃动力学研究方法是一种描述烃源岩生烃特征的有效方法,根据开放体系下烃源岩的生烃模拟试验结果,可获得动力学参数(活化能和指前因子),这两个参数可反映出烃源岩生烃的难易程度[3]。
2020-06-01土木工程不断地为人类社会创造崭新的物质环境,尤其是经历了近代工业技术革命的洗礼后,其便逐渐形成了一个比较完善的现代工程建设体系。随着土木工程项目的日益增多,生活中的建筑施工、道路工程、高铁项目等都可将其看作土木工程的内容。在这一过程中,土木工程施工的内容也日益丰富,项目工程的前期勘测、施工准备、技术选择、质量把控等都成为其实现的重要环节。
2020-05-29人气:6027
人气:3748
人气:3210
人气:2962
人气:2955
我要评论
期刊名称:应用数学和力学
期刊人气:2170
主管单位:重庆交通大学
主办单位:重庆交通大学
出版地方:重庆
专业分类:科学
国际刊号:1000-0887
国内刊号:50-1060/O3
邮发代号:78-21
创刊时间:1980年
发行周期:月刊
期刊开本:16开
见刊时间:一年半以上
影响因子:0.963
影响因子:1.730
影响因子:0.542
影响因子:0.744
影响因子:1.028
您的论文已提交,我们会尽快联系您,请耐心等待!
你的密码已发送到您的邮箱,请查看!