基于matlab的多支座蒸壓釜的受力分析和強度計算

1、基于matlab的多支座蒸壓釜的受力分析和強度計算論文關(guān)鍵詞：蒸壓釜多支座臥式容器剪力彎矩atlab論文摘要：根據(jù)材料力學(xué)理論，推導(dǎo)出多支座臥式壓力容器支座彎矩、支座反力的計算公式，給出釜件應(yīng)力的校核方法，編制了基于atlab的計算程序，并附有五鞍座蒸壓釜的受力分析和強度計算實例。蒸壓釜是化工、建材行業(yè)中應(yīng)用較為廣泛的一種具有多支座的臥式壓力容器。gb15089?鋼制壓力容器?只給出了雙支座對稱布置臥式容器的剪力、彎矩和應(yīng)力計算方法，而對多支座臥式容器的計算方法僅在文2中有一般性說明。hgj1689?鋼制化工容器強度計算規(guī)定?只給出了三支座臥式容器的設(shè)計和計算。本文從三彎矩方程出發(fā)結(jié)合臥式容器

2、的特點，基于材料力學(xué)理論根底，比照文2，3的推導(dǎo)過程，導(dǎo)出了相應(yīng)的彎矩和剪力計算公式，按齊克法給出應(yīng)力校核式，并基于atlab編制了相應(yīng)的計算及校核程序。圖1n支座連續(xù)外伸梁受力圖1公式推導(dǎo)多支座臥式蒸壓釜可簡化為受均布載荷的外伸梁，假設(shè)共有2n-1個支座，見圖1。圖中l(wèi)為圓筒兩封頭切線之間長度;h為封頭內(nèi)壁曲面深度;a為邊支座中心線到近端封頭切線的間隔 ;q為單位長度上的載荷。1.1支座截面彎矩的計算為求出臥式容器的各項應(yīng)力，首先求得各截面的剪力和彎矩，多支座臥式容器屬靜不定構(gòu)造，需用三彎矩方程求解其剪力和彎矩。對于n支座，三彎矩方程一般性公式為2：(1)在對多支座臥式容器進(jìn)展構(gòu)造設(shè)計時，為

3、使其受力狀況較好，通常將支座設(shè)計成等距布置。即有l(wèi)1=l2=l3=l2n-2=l，其值為：(2)故公式(1)簡化為(3)在支座n左右，由于對稱，故(4)首先需要計算封頭及其內(nèi)裝物料重量和作用于封頭上的靜載荷對封頭切線與軸線交點的等效力矩：(5)根據(jù)文4有(6)由公式(4)、(5)、(6)組成連立方程如下：(7)令，那么上述方程組可以寫成：(8)其中：由此方程組可解得各支座截面彎矩i(i=1,2,，n)，并由對稱性得(i=1,2,，n-1)(9)1.2支座反力的計算圖2連續(xù)外伸梁分解圖當(dāng)求得各支座截面彎矩之后，把該連續(xù)梁分解為2n-2個靜定梁，如圖2所示。左右端為均布載荷外伸梁、中間為均布載荷簡

4、支梁，從而求得支座反力。封頭及其內(nèi)裝物料重量為：并由對稱性得(i=1,2,，n-1)1.3梁內(nèi)剪力的計算根據(jù)剪支梁的剪力計算公式可求出梁內(nèi)的剪力,如下：13)求出a段以及1n-1段的剪力，其它各段由對稱性可得。1.4彎矩的計算首先根據(jù)剪支梁彎矩的計算公式即可求出梁內(nèi)彎矩，然后采用數(shù)學(xué)方法求出彎矩的最大值。1.4.1梁內(nèi)彎矩的計算(14)求出a段以及1n-1段的彎矩，其它各段由對稱性可得。1.4.2最大彎矩的計算對方程組(14)求一階導(dǎo)數(shù)，以求最大彎矩所在點，如方程組15。把求得的x值代入原方程組(14)，即可求得各段的最大彎矩。在求得多支座臥式容器各支座處支反力和彎矩后，可作出其剪力圖和彎矩圖

5、，本文不詳細(xì)表達(dá)。(15)1.5釜體應(yīng)力計算和校核多支座臥式容器的應(yīng)力計算可以按照齊克方法，依次求解1、2、38及，并按文獻(xiàn)1進(jìn)展應(yīng)力校合。1.5.1筒體軸向應(yīng)力計算(1)兩支座中間處的橫截面上:按各跨中點處的最大彎矩ax作用點處，計算橫截面的最高點和最低點的軸向應(yīng)力：最高點：(16)最低點：(17)式中，p為設(shè)計壓力，pa；r為釜體平均半徑，；=ax|i=1,2,，n。(2)支座處橫截面上筒體被加強的最高點或筒體不被加強的靠近中間程度平面處：最高點：(18)最低點：(19)式中k1、k2為計算應(yīng)力系數(shù)，根據(jù)ar12和支座包角按文獻(xiàn)2式85、式86計算，。(3)筒體軸向應(yīng)力的驗算(20)式中k

6、3為計算應(yīng)力系數(shù)，根據(jù)ar12和支座包角按文獻(xiàn)2式821進(jìn)展計算，=ax|i=1,2,，n。如0.8t,那么切應(yīng)力校核合格。1.5.3圓筒周向應(yīng)力計算按支座處無加強圈，先按鞍座墊板不起加強作用進(jìn)展計算。(1)支座處橫截面最低點：(21)式中k5為計算應(yīng)力系數(shù)，根據(jù)ar12和支座包角按文獻(xiàn)2式8-20計算，=ax|i=1,2,，n。(2)鞍座邊角處(22)式中k6為計算應(yīng)力系數(shù)，根據(jù)ar12和支座包角按文獻(xiàn)2式8-36計算，=ax|i=1,2,，n。(3)周向應(yīng)力驗算7t,81.252、atlab程序的編制根據(jù)上述推導(dǎo)公式，基于atlab語言編寫計算程序。本文利用atlab強大的矩陣計算功能進(jìn)展

7、方程組的求解，大大的簡化的計算的復(fù)雜性，且利用atlab的繪圖功能可以很方便的繪制出剪力圖與彎矩圖。程序的pad圖如圖3：shape*ergefrat圖3atlab程序pad圖3五鞍座蒸壓釜計算實例一臺五鞍座蒸壓釜構(gòu)造示意圖見圖4。：支座個數(shù)2n-1=5，設(shè)計壓力p=1.40pa，邊支座中心線到近端封頭切線的間隔 a=740，封頭內(nèi)壁曲面深度h=500，鞍座包角=150，圓筒平均半徑r=1006，圓筒長度(兩封頭切線之間)l=26680，設(shè)計溫度下容器材料的許用應(yīng)力t=170pa，筒體有效厚度0=11，容器殼體及充滿介質(zhì)時的總重量包括殼體、內(nèi)件、物料及保溫層=1.05106n。圖4五鞍座蒸壓釜

8、構(gòu)造示意圖將參數(shù)輸入atlab程序，得計算結(jié)果如下：支座彎矩支座反力跨間最大彎矩應(yīng)力校核n,(n)n,pa1=-1.32107r1=14463112ax=1.1611081=60.72=-1.663107r2=28439723ax=5.8471072=67.33=-1.153107r3=23581034ax=5.8471083=66.34=-1.663107r4=28439745ax=1.1611084=93.55=-1.32107r5=144631iax=1.1611085=-10.6ax=-1.32107rax=2843976=-151.21、2、3、4、5、6t，0.8t=10.86各項盈利校核均合格并且程序可以自動生成剪力圖和彎矩圖，如圖5、圖6所示。4結(jié)語本文所給出的受力分析和強度計算方法，雖然是針對詳細(xì)的蒸壓斧設(shè)計而得，但是具有非常普遍的意義，廣泛適用于其他多支座臥式容器。同時，在本文中，atlab強大的計算功能與繪圖功能得到充分表達(dá)，值得在設(shè)計與計算中廣泛推廣。參考文獻(xiàn)：1全國壓力容器標(biāo)準(zhǔn)化委員會，gb15089，鋼制壓力容器(一)北京：學(xué)院出版社，l9892全國壓力容器標(biāo)準(zhǔn)化委員會，gb15