-
管板组合是换热设备中常见的一种结构形式,广泛应用于冷凝器和锅炉中.作为锅炉炉膛的焊制水冷壁,由管板交替排列组成,在工作中承受热流、炉膛烟气压力、自重等多种载荷作用.
对管板组合结构的静动力学分析,前人做了一些研究工作.文献[1-4]针对焊制水冷壁结构提出了一些具体的简化方法,并用于热应力分析.文献[5-6]利用二维有限差分和有限元法模拟了水冷壁中温度场的分布.文献[7]建立了锅炉水冷壁系统的有限元仿真模型.但从现有文献看,有限元法是分析管板组合结构的重要方法,但其复杂的建模,大的计算量,很难用于结构参数化设计.
由于焊制水冷壁由管板周期排列组合而成,使其整体表现出正交各向异性,文献[8-11]提出用各向异性板对这类结构进行简化,不仅可使计算量减小,还便于参数化设计.
本文将在文献[1]的基础上,将管板组合结构等效为正交各向异性板,通过变形协调方程和刚度分析得到其相关的等效参数,通过与3D有限元结果比较说明方法的正确性.
全文HTML
-
考虑图 1(a)所示周期性排列的管板组合结构,取图 1(b)所示截距为p的基本单元.其中,鳍片宽度l,管子外径r0,管子内径ri,鳍片厚度tf,管子厚度t,管子中面距管心距离rm.其等效正交各向异性板如图 2.
平面正交各向异性热力学问题,应力应变关系为
其中,σx,σy,σxy为应力分量,εx,εy,γxy为应变分量,Exx,Eyy,Exy,Gxy为弹性模量和剪切模量,αx=αy=α为热膨胀系数.
平面问题的应变协调方程为
将应力应变关系(1) 代入协调方程(2),并利用
${\sigma _x} = \frac{{{\partial ^2}\varphi }}{{\partial {y^2}}}, {\sigma _y} = \frac{{{\partial ^2}\varphi }}{{\partial {x^2}}}, {\tau _{xy}} =-\frac{{{\partial ^2}\varphi }}{{\partial x\partial y}}$ ,得到用应力函数表示的变形协调方程文献[1]给出了水冷壁结构在面内温度场作用下的变形协调方程为
其中,F为应力函数,T为温度场,α为材料的热膨胀系数,E为材料弹性模量,G为剪切模量,μ为泊松比. Fx,Fy,Fz为仅与结构几何有关的系数,表示为
比较式(3) 与式(4),并利用Exy=φExx以及F=φ·Fy可得
-
从管板结构几何构型可知,y方向的等效抗弯刚度为
其中,Dy1为管部分的抗弯刚度,Dy2为鳍片部分的抗弯刚度.由材料力学可得
-
管板结构沿x方向为变截面,整体抗弯能力主要由截面抗弯刚度最弱处决定.由薄板弯曲理论得其等效刚度为
式中,μxy和μyx为等效后正交各向异性板的泊松比,且μyxExx=Exy,μxyEyy=Exy,h为其等效厚度,将由式(8) 给出.
由文献[8]可得扭转刚度Dxy为
-
由薄板弯曲理论得y方向的抗弯刚度为
结合公式(6) 及上述Dy的表达式得等效厚度为
由等效前后质量相等,得等效后密度ρ1为
其中,ρ为管板材料密度.
1.1. 等效材料性质
1.2. 等效抗弯刚度
1.2.1. y方向的抗弯刚度
1.2.2. x方向的抗弯刚度及扭转刚度
1.3. 等效厚度与等效密度
-
为验证等效结果,分别用三统合有限元方法(3D-FEM)和本文方法进行计算.计算中取图 1(b)所示长为50倍节距长宽相等的管板组合结构为计算模型,基本单元参数如表 1.
-
考虑四边简支,受单位均布载荷的正交各向异性板.其挠度的经典Navier解[8]为
其中a,b为板的长和宽.计算中取上述级数前20项,结果如表 2.
-
对四边简支正交各向异性板,其频率公式为
其中,m表单位面积质量,等于其等效密度与等效厚度之积.计算管板结构前9阶固有频率,如表 3.
-
以上算例是针对规模为50p的结构进行的,实际工程中管板数量达到数百个,现讨论管板规模对等效模型精度的影响. 图 3给出了板中心受单位集中力该处挠度的相对误差随管板规模的变化情况. 图 4给出了四边简支管板前四阶固有频率的相对误差随管板规模的变化情况.
由图 3、图 4可见,随着管板规模的增加,等效板理论模型与3D有限元结果的相对误差逐渐减小,该简化模型可用于大规模管板结构计算.
2.1. 静态挠度
2.2. 动态频率
2.3. 管板结构尺寸规模对等效结果精度的影响
-
应用经典弹性理论将管板组合结构等效为正交各向异性板,得到其等效材料参量;使用3D有限元法分别从静力学和动态特性两方面对等效模型进行了验证;讨论了等效模型对管板规模的适用性.研究表明该等效模型可应用于大规模管板组合结构的静动力学计算.