用matlab求又向圈,基于MATLAB算法的外压容器环向加强圈优化设计

在石油和化工生产中,处于外压操作条件的容器是比较常见的,如石油分馏中的减压蒸馏塔、带蒸汽夹套的反应釜及各类真空贮槽等。由于这类容器的失稳往往是在强度破坏之前发生,所以稳定性就成为外压容器考虑的主要问题。由文献[1,2]的分析结果表明,外压圆筒在不改变直径尺寸和外压载荷的前提下,通过装置刚性构件来提高临界压力的方法,要比单独增加壁厚的方法更为有效、合理。因此,在外压圆筒上设置必要的加强圈,即成为保证外压容器安全运行的一个重要步骤。目前,虽然有不少关于加强圈设计计算的论述,但多是按照传统或常规方法[3]进行的。此方法虽然合理,但需要通过多次反复的试算才能得到一个较为满意的设计方案,而且结构尺寸的选取,主要是根据设计人员的经验来确定,缺乏严格科学的定量分析基础,故很难得到近乎最优的设计方案,导致结果偏于保守。优化设计方法则是建立在数学规划和现代计算机技术的基础上,它既能定量分析和解决产品的安全性,也能使产品的功能参数实现最优化,是解决复杂工程问题的一种有效工具。笔者运用优化设计理论和MAT-LAB算法语言,以一定体积的加强圈质量最小为设计目标,探讨环向加强圈尺寸优化设计的一般方法。1环向加强圈优化设计的数学模型[4,5]1.1确定设计变量与目标函数当材料一定时,环向加强圈的总质量主要由加强圈截面积As、截面中性轴直径Ds及加强圈间距Ls等参数确定。考虑Ds受限于筒体外径,故选择与加强圈面积有关的参数h、及Ls作为设计变量,即:X=[x1,x2,x3]T=[h,,Ls]T优化设计追求的目标为加强圈总质量最轻,于是目标函数的表达式为:F(X)=h(Do+h)(LLs-1)=x1x2(Do+x1)(Lx-31-1)(1)式中h、扁钢加强圈截面的高度和厚度,mm;Do外压圆筒外直径,mm;L计算长度,即筒体上两刚性构件(如封头、法兰、支座以及加强圈等)之间的最大距离,可根据文献[3]计算确定,mm;加强圈材料密度,kg/m3。1.2确定约束条件加强圈是一个刚性结构件,设计变量取值既要使其满足刚度要求,同时还要考虑尺寸规格和间距等多种因素的限制。1.2.1刚度条件为保证圆筒不致失稳,所设加强圈必须具有足够的刚性。如图1所示,焊接在外压圆筒上的加强圈可视为一个刚性较大的圆环,并与圆筒共同承担作用在其两侧各Ls/2范围内的载荷。由此计算得到加强圈与圆筒体当量截面所需的综合惯性矩[1,3]:I=LsDo2(e+As/Ls)10.9为保证加强圈安全工作,选用的加强圈与两侧圆筒有效段构成的组合截面实际惯性矩Is应大于并接近I,即IsI。图1环向加强圈结构图从图2可以看出,若加强圈与其两侧圆筒有效段组合截面的实际惯性矩[3,6]为:Is=I1+As(0.5h+0.5e-a)2+I2+Aca2则由Is=I1+As(0.5h+0.5e-a)2+I2+Aca2I得到约束条件:g1(X)=[D2ox3(e+x1x2/x3)/10.9]-[I1+x1x2(0.5x1+0.5e-a)2+I2+Aca22]0(2)式中Ac外压圆筒有效宽度b内的截面积,Ac=2be,mm2;a圆筒截面形心轴至组合截面形心轴的距离,mm;b外压圆筒的有效宽度,取b=0.55Doe,mm;I1加强圈截面对其形心轴的惯性矩,扁钢取I1=h3/12,mm4;I2外压圆筒截面对其形心轴的惯性矩,mm4;e外压圆筒有效厚度,e=n-C,mm;n、C外压圆筒的名义厚度和厚度附加量,mm;应变。图2扁钢加强圈及组合截面图1.2.2最大间距条件mpDoDoe由文献[2]可知,设置加强圈时,只有使计算长度L小于临界长度Lcr