matlab受力分析,基于Matlab的多支座蒸压釜的受力分析和强度计算

0前言蒸压釜是化工、建材行业中应用较为广泛的一种具有多支座的卧式压力容器。GB150-1989《钢制压力容器》只给出了双支座对称布置卧式容器的剪力、弯矩和应力计算方法,而对多支座卧式容器的计算方法仅在文献[2]中有一般性说明。HGJ16-1989《钢制化工容器强度计算规定》只给出了三支座卧式容器的设计和计算。本文从三弯矩方程出发,结合卧式容器的特点,基于材料力学理论基础,比照文献[2,3]的推导过程,导出了相应的弯矩和剪力计算公式,按齐克法给出应力校核式,并基于Matlab编制了相应的计算及校核程序。1公式推导多支座卧式蒸压釜可简化为受均布载荷的外伸梁,假设共有2n-1个支座,见图1。图中:L为圆筒两封头切线之间长度;h为封头内壁曲面深度;A为边支座中心线到近端封头切线的距离;q为单位长度上的载荷。图1n支座连续外伸梁受力图11支座截面弯矩的计算为求出卧式容器的各项应力,须首先求得各截面的剪力和弯矩。多支座卧式容器属静不定结构,需用三弯矩方程求解其剪力和弯矩。对于n支座,三弯矩方程一般性公式为[2]:Mi-1li+2Mi(li+li+1)+Mi+1li+1=-q(l3i+4li3+1)(1)在对多支座卧式容器进行结构设计时,为使其受力状况较好,通常将支座设计成等距布置。即有:l1=l2=l3=…=l2n-2=l其值为:q=L+(W4/3)h(2)故式(1)简化为:Mi-1+2Mi+Mi+1=-ql22(3)在支座n左右,由于对称性Mn-1=Mn+1,故2Mi-1+4Mi=-ql2/2(4)首先需要计算封头及其内装物料重量和作用于封头上的静载荷对封头切线与轴线交点的等效力矩:M=41q(R2m-h2)(5)根据文献[4]有:M1=M-q(32hA+21A2)(6)由式(4)、(5)、(6)组成联立方程如下:M1=M-q(23hA+21A2)M1+4M2+M3=-q2l2M2+4M3+M4=-q2l2……Mn-2+4Mn-1+Mn=-q2l22Mn-1+4Mn=-q2l2(7)令F1=M-q(23hA+21A2)F2=F3=…=Fn=-q2l2则上述方程组可以写成:LM=F(8)其中:L=1000……01410………0141…0……………410…………12M=[M1,M2,…,Mn]TF=[F1,F2,…,Fn]T由此方程组可解得各支座截面弯矩Mi(i=1,2,…,n),并由对称性得:Mn+i=Mn-i,(i=1,2,…,n-1)(9)12支座反力的计算当求得各支座截面弯矩之后,把该连续梁分解为2n-2个静定梁,如图2所示。左右端为均布载荷外伸梁、中间为均布载荷简支梁,从而求得支座反力。封头及其内装物料重量图2连续外伸梁分解图为:G=32hqR1=G+qAR2=21ql+M1-lM2……Rn=21ql+Mn-1l-Mn(10)R1=21ql+M2-lM1R2=21ql+M3-lM2……Rn=21ql+Mn+1l-Mn(11)R1=R1+R1R2=R2+R2……Rn=Rn+Rn(12)并由对称性得Rn+i=Rn-i,(i=1,2,…,n-1)13梁内剪力的计算根据简支梁的剪力计算公式可求出梁内的剪力如下:A段:F0~1=-G-qx,(x=0~A)l1段:F1~2=R1-qx,(x=0~l)l2段:F2~3=R2-qx,(x=0~l)……ln-1段:Fn-1~n=Rn-1-qx,(x=0~l)(13)求出A段以及1~n-1段的剪力,其它各段由对称性可得。14弯矩的计算首先根据简支梁弯矩的计算公式即可求出梁内弯矩,然后采用数学方法求出弯矩的最大值。141梁