变壁厚固顶型钢储罐在实测不均匀沉降下稳定性分析.docVIP

变壁厚固顶型钢储罐在实测不均匀沉降下的稳定性分析 一 分析模型 储罐几何参数与材料参数均与浮顶罐相同：高度与半径之比h/ tav =1，径厚比r/ tav =1000；弹性模量E=209GPa， 泊松比ν=0.3。 采用通用有限元程序ANSYS进行分析，采用壳体单元Shell63。有限元模型见图1。取柱面坐标系，x、y、z方向分别代表壳体的径向（即法向）、切向和纵向。由于固顶罐在其顶部具有固定的顶盖约束其变形，对罐顶、罐底均约束其径向和切向线位移，即。分析中考虑几何非线性，采用弧长法跟踪结构的全过程响应以完整了解储罐的屈曲和后屈曲行为。 图1 固顶罐有限元计算模型 二 固顶罐分析的沉降特点 由文献[1]分析可知，在固顶罐分析中除了考虑n≥2的不均匀沉降，还应考虑n＝1项。图2给出了由n≥1各项叠加得到的固顶罐沉降分布曲线。由于要考虑n＝1项，前述的整体不均匀沉降可进一步分为两种，一种是n＝1项明显大于其余各项，即储罐整体倾斜占绝对主导地位，这一类的储罐包括G1-G5罐；另一种是n＝1项与其他各项基本相当，甚至小于其它项，这一类的储罐包括G6、G7罐。[2] 图2 8个储罐的不均匀沉降（n≥1叠加） 三 整体不均匀沉降下的非线性分析 对于固顶罐，整体不均匀沉降可分为两类，下面就整体倾斜沉降，即n＝1占主导地位的储罐G1-G5进行分析。以G4罐为例，分析变壁厚固顶罐在不均匀沉降下的稳定性。图3给出了储罐G4的沉降位移曲线，所跟踪的节点为罐壁径向位移最大的节点，图4给出了变形发展过程中的几个典型变形图。从曲线上可以看到典型的跳跃屈曲现象，储罐从开始加载到A点（Smax / tav＝0.16）的响应基本是线性的，也没有出现明显变形（图4a）。此后，结构响应开始表现出一定的非线性，当沉降值增大至B点（Smax / tav＝0.25）时，出现一个荷载极值点以及典型的跳跃屈曲现象。沉降位移曲线在B点之后下降至C点（Smax / tav＝0.16）时，在罐壁上部出现明显的变形，结果是变形的幅值达到4.1tav 。之后，沉降位移曲线又开始上升，表明罐壁的局部屈曲并没有导致结构的整体失稳。总体上看结构的后屈曲性能是稳定的，但是结构剪切屈曲后，储罐出现了跳跃屈曲的现象，储罐所承受的沉降值下降后又有所上升，剪切变形不断发展，因此可以把出现剪切屈曲的临界值作为极限沉降值。与同等条件下的均匀壁厚储罐相比，二者的临界沉降值相差不大，但是当考虑后屈曲性能时，浮顶罐的极限沉降值大得多。和均匀壁厚的固顶罐相比，二者的后屈曲性能相差不大，但是变壁厚储罐的极限沉降值（Smax / tav＝0.25）小于浮定罐的（Smax / t＝0.50）。此外，从变形发展（图4）中可以看到，由于储罐壁厚从下到上递减，储罐过大变形主要集中在罐壁的上部区域。 图3 G4罐沉降－位移曲线 (a) A点 (b) B点 (c) C点 图4 G4罐的变形发展（变形放大30倍） 对G1、G2、G3、G5固顶型变壁厚储罐进行了同样的分析，发现各储罐在实测沉降下的结构响应均与G4罐类似，限于篇幅，这里不再给出各储罐的具体分析，仅在图5给出各罐的变形发展过程。 (a) G1 (b) G2 (c) G3 (d) G5 图5 G1、G2、G3、G5罐的变形发展（变形放大30倍） 再对整体不均匀沉降中n＝1项与其余各项相当的G6、G7罐进行分析。图6是G7罐的沉降－位移曲线，所跟踪的节点是罐壁径向位移最大点，图7给出了变形发展过程中的几个典型变形图。从曲线上可以看出，从储罐开始承受沉降到A点（Smax / tav＝0.19）时罐壁的响应基本上是线性的（图7a），A点之后开始表现出非线性，当沉降值达到B点（Smax/ t＝0.23）时罐壁开始发生屈曲变形（图7b）， B点对应的沉降值就是临界沉降值,此时结构变形并不明显。之后曲线一直呈现下降趋势，C点（Smax / tav＝0.22）为结构的极限沉降值。当沉降值达到D点（Smax / tav＝0.16）时，/ t已经从C点的0.89增加到3.41（图7c）。鉴于G7罐的后屈曲性能不是稳定的，所以它所能承受的极限沉降值就是它的临界沉降值，小于同等条件下浮顶罐的极限沉降值（Smax / tav＝0.61）相比要小很多。与G1~G5罐类似，由于壁厚从下到上的递减性变化，导致罐壁上部局部区域较下部更为