式中 — 试验压力,MPa;
— 液柱静压力,MPa;
— 承压筒体的有效厚度,mm;
— 试验温度下材料的屈服极限,MPa。
3 优化分析前的准备-静态分析
3.1 筒体三维模型的建立
由于Solidworks Simulation若要进行优化分析,则必须先对分析对象作静态分析,因此此处先对承压容器的主筒体作静力有限元分析。Simulation在其结构分析中,与其他CAE分析过程一样,一般包括前处理、求解和后处理三个部分,即几何模型的建立、定义材料属性、网格的划分、加载、求解和结果分析。由于Solidworks Simulation分析模块与Solidworks是无缝集成的,因此在不出CAD环境下即可对所建模型进行相关分析,故本文的压力容器筒体即由三维制图软件Solidworks完成实体模型的建立,所得压力容器筒体的三维模型如图2所示。
图2 压力容器筒体三维模型
3.2 筒体约束、加载及网格的划分
在完成三维模型建立后,对压力容器模型进行约束、加载及划分网格,使其转为有限元模型。根据16MnDR材质,筒体材料属性定义为:弹性模量(EX)、泊松比值分别取为190GPa和0.3,屈服强度和抗拉强度分别为310MPa和490MPa。由于筒体内壁受均匀的圆周径向压力,故对筒体长端内壁加载设计压力(34Mp),顶端连接筒是起连接旋转盖和主筒体作用,不受水压,故顶端连接筒不加载,筒体俩端面部位为固定约束,由于主筒体内壁为受力部位,故需对筒体内壁进行局部细分,划分结果得节点数为15391个,单元数为7688个,网格划分结果及解算器信息如图3、图4所示。
图4 解算器信息
图3 筒体网格划分
3.3 结果分析
当模型、材料、网格、约束及加载均已完成后,即对筒体进行在此条件下的求解,本文采用的解法是FFEplus算法,并采用p-自适应方法调整。在求解后的的变形图显示(为米塞斯应力云图)中,为使变形结果较为清楚且又不夸张的显示,将变形比例定为放大100倍,结果如图5所示。由变形图可知,最大应力发生在筒体内壁处部位,并且变形量较大,最小应力发生在顶端连接筒处,与实际受压情形一致。
图5 筒体在加载后的应力云图 图6 筒体的安全系数图解
由图5可以看出,壁厚为25mm、材质为16MnDR的受压筒体在加34MPa静压力后,内壁最大应力为196.8MPa,远小于材料的屈服极限值310MPa,因此用材质为16MnDR、壁厚为25mm设计的高压容器筒体是安全的,此结论也可以从图6的安全系数图中得到进一步证实。根据GB 150-1998标准,对于钢制压力容器,当屈强比小于或等于0.7时,其安全系数取为1.6[3],由图6可以看出,其最小安全系数出现在筒体内壁处,值为1.69。考虑到原材料的尺寸偏差及高压容器筒体在使用过程中的腐蚀余量,依据此分析,可以将原壁厚增加3mm,以增加高压容器的安全性及可靠性。
4 优化分析
4.1 设计变量及约束条件数值的确定
综合以上所述及分析,根据Solidworks Simulation静态分析确定设计变量及约束条件的数值如下:
承压容器壁厚:12.5mm37.5mm
承压容器内径D1 :97.5mm292.5mm
承压容器受载筒长L :3000mm9000mm
4.2 结果分析
在对承压容器静态分析及各优化参量确定的基础上,对承压容器进行了优化分析,优化迭代次数为11次,优化后所得优化结果及设计变量变化过程和趋势如图7、图8所示。
DV1
DV2
图7 承压容器各参量优化结果
图8 设计变量在优化迭代过程中的变化趋势
由图7可以看出,在优化前后承压容器的质量m、容器壁厚及内径D1 均有所变化,结合图7、8得出优化设计所得结果与原设计结果相比:
A:承压容器质量变化率:
B:承压容器壁厚变化率:
C:承压容器内径变化率:
由以上各式可以很清楚的得出:在满足压力容器设计要求的前提下及确保高压容器长度基本不变的条件下,利用优化设计的方法对承压容器的外形结构尺寸进行优化,可使其承压容器壁厚和质量降低三分之一左右,结合实际使用情况;承压容器内径可以适当的扩大,从而增加容器的体积,由此可以一次对多个被测件进行压力试验。由此可以得出:在满足使用及设计要求的前提下,对高压容器的外形结构尺寸进行优化后,能够节约材料,降低了制造成本,这对于提高压力容器尤其是高压容器的经济性,使之符合当今市场经济条件下对容器实现低成本设计的要求具有非常重要的意义[5]。